WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


 


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

       

       

ТОПИЛЬНИКОВ ВЛАДИМИР ИВАНОВИЧ

Разработка модели процесса переработки продуктов синтеза Фишера-Тропша,

получаемых на кобальтовом катализаторе

02.00.13 – Нефтехимия

АВТОРЕФЕРАТ

диссертации на соискание учёной степени

кандидата технических наук

МОСКВА – 2012

Работа выполнена на кафедре газохимии Российского государственного университета нефти и газа имени И. М. Губкина.

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

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

Сосна Михаил Хаймович

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

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

Лысенко Сергей Васильевич

доктор технических наук

Шурупов Сергей Викторович

Ведущая организация: Московский государственный университет тонких химических технологий имени М. В. Ломоносова

Защита диссертации состоится  «___»    _ 2012 г. в « » часов на заседании диссертационного Совета 212.200.04 при Российском государственном университете нефти и газа им. И. М. Губкина по адресу: 119991, г. Москва, Ленинский просп., 65, аудитория ____.

С диссертацией можно ознакомиться в библиотеке РГУ нефти и газа им. И. М. Губкина

Автореферат разослан «___» _____________ 2012 г.

Учёный секретарь

диссертационного совета,

доктор технических наук, профессор                                                Сафиева Р. З. 

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ



Актуальность темы

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

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

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

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

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

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

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

Создана математическая модель процесса гидрокрекинга парафинов, учитывающая:

  • парожидкостное равновесие на входе и по высоте реактора;
  • свойства катализатора (количество кислотных центров на его поверхности), в результате чего появляется возможность прогнозировать пути совершенствования катализатора.

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

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

С помощью созданной математической модели можно рассчитать геометрические размеры реактора гидрокрекинга для катализатора Pt/ZSM-5.

Модель предоставляет возможность использовать свойства катализатора как параметры моделирования, выявлять и прогнозировать направления его совершенствования, уменьшив за счёт этого объём экспериментальной работы.

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

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

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

  1. Математическая модель гидрокрекинга парафинов для катализатора Pt/ZSM-5, учитывающая парожидкостное равновесие и Гауссовское распределение продуктов гидрокрекинга;
  2. Математическая модель, использующая свойства катализатора в качестве параметров моделирования.

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

Результаты работы докладывались на I Российском нефтяном конгрессе (14–16 марта 2011 г., Москва); III Международной конференции «АИСТ-2011», Альтернативные источники сырья и топлива (24–26 мая 2011 г., Минск, Белоруссия); II Международной конференции «Газохимия-2011» (5–6 октября 2011 г., Москва); Конференции молодых учёных по нефтехимии (13–16 сентября 2011 г., Звенигород); 9-ой всероссийской конференции молодых ученых, специалистов и студентов (4–7 октября 2011 г., Москва).

Публикации

По теме диссертации опубликовано 3 статьи.

Структура и объём работы

Диссертационная работа состоит из введения, 4 глав, основных результатов и выводов, библиографического списка использованной литературы из 101 наименования. Работа изложена на 96 страницах машинописного текста и содержит 35 рисунков и 19 таблиц.

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Во введении обоснована необходимость химического преобразования продуктов синтеза ФТ с помощью процесса гидрокрекинга для получения высококачественных моторных топлив и компонентов масел.

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

Процесс ФТ является сложным каталитическим процессом, включающим большое число последовательных и параллельных реакций.

Катализаторами являются переходные металлы VIII группы (в основном Fe, Co, Ru):

  • на железных катализаторах при давлении 2530 атм и температуре 230240 °C образуются олефины, парафины и кислородсодержащие соединения;
  • на кобальтовых системах при давлении 130 атм и температуре 170250 °C можно селективно синтезировать линейные алканы.

В синтезе ФТ получают спектр алифатических углеводородов – молекулярно-массовое распределение (Рис. 1), который описывается уравнением Андерсона-Шульца-Флори

,                                        (1)

где Wn – массовая доля углеводорода, содержащего n атомов углерода,  – коэффициент роста цепи.

Рис. 1. Молекулярно-массовое распределение продуктов синтеза ФТ

для различных значений

Для эффективного использования синтез-газа процесс ФТ необходимо проводить при больших значениях (регулируется температурой), это приводит к утяжелению фракционного состава смеси парафинов.

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

Из-за высокого содержания олефинов и оксигенатов продукты ФТ более активно вступают в химические превращения, чем компоненты сырой нефти.

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

В сводной табл. 1 приведено сравнение процессов нефтепереработки, применимых для переработки продуктов синтеза ФТ.

Таблица 1. Сравнение процессов нефтепереработки

Каталитический крекинг

Гидрокрекинг

преимущества

недостатки

преимущества

недостатки

Получение высокооктановых бензинов

Необходимость гидрооблагораживания продуктов крекинга

Улучшение низкотемпературных свойств топлив

Циркуляция большого количества водорода

Переработка тяжёлых парафинов

Потеря ~5 % сырья из-за образования кокса

Обеспечение необходимого цетанового числа и высокий выход дизельного топлива

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

Катализаторы гидрокрекинга состоят из трёх компонентов: кислотного, дегидро-гидрирующего и связующего, обеспечивающего механическую прочность и пористую структуру.

Можно выделить три основные группы катализаторов гидрокрекинга: сульфидные, оксидные и на основе благородных металлов.

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

Оксидные катализаторы имеют недостаточно высокую активность для разрыва связи углерод–углерод, поэтому к ним добавляют кислотные компоненты. Применение данного типа катализаторов обуславливается наличием в сырье гидрокрекинга значительного количества каталитических ядов. В сырье процесса ФТ количество каталитических ядов незначительно, поэтому применение данного типа катализаторов нецелесообразно.

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

При моделировании процессов переработки смесей углеводородов возникают трудности, характерные как для нефтепереработки, так и для переработки продуктов синтеза ФТ:

  • большое количество участников процесса;
  • большое количество однотипных реакций;
  • многостадийность процесса.

В работе описана математическая модель процесса гидрокрекинга нормальных парафинов (С4С70). Катализатором служит платина, нанесённая на аморфный алюмосиликат. Для упрощения моделирования исследователи объединили индивидуальные компоненты в псевдокомпоненты согласно существующим в нефтепереработке фракциям: C1-4, n-C5-9, iso-C5-9, n-C10-14, iso-C10-14, n-C15-22, iso-C15-22, n-C22+, iso-C22+; во фракцию C1-4 также включен образующийся изобутан. Процесс гидрокрекинга протекает в газовой фазе, несмотря на различия в температурах кипения фракций.





Схема процесса гидрокрекинга имеет следующий вид:

               (2)

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

Данную модель не следует использовать в качестве основы для дальнейшего моделирования, поскольку она не отражает реальный механизм и химизм гидрокрекинга парафинов.

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

[CnH2n] +[H2]  [CnH2n+2],                                                (3)

поэтому в модели олефины рассматриваются как парафины. 

Процесс гидрокрекинга парафинов протекает по следующему механизму:

[CnH2n+2] + [H2]  [CsH2s+2] + [CrH2r+2]                                (4)

где n = r + s.

Считается, что разрыв молекулы парафина происходит посередине. Если n – чётное, то s = r, если нечётное, то s = r + 1 = (n/2 + 0,5).

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

В другой модели процесс гидрокрекинга парафинов протекает по механизму (5): линейные парафины подвергаются изомеризации (равновесный процесс), а затем происходит крекинг изопарафинов:

               (5)

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

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

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

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

Механизм процесса

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


Реакции изомеризации

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

Рис. 2. Реакция изомеризации (pcp-изомеризация)

Реакции крекинга

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

В зависимости от структуры карбокатиона возможны пять вариантов реакции -разрыва, которые представлены на рис. 3.

Рис. 3. Типы реакции -разрыва. Точками представлены алкильные группы, символы p,s,t –первичный, вторичный и третичный карбокатион соответственно

Групповые объединения

В рассматриваемой модели учитываются изомеры парафинов, имеющие не больше трёх метильных заместителей. Учесть все структурные изомеры в модели сложно, поэтому структурные изомеры с одинаковым числом метильных заместителей объединяются в групповые объединения, и в дальнейшем уже они рассматриваются в качестве индивидуальных компонентов. Рассмотрим принцип групповых объединений на примере изомеров гептана (C7H16), объединённых в четыре типа: без метильных заместителей (NP), с одним (MBP), двумя (DBP) и тремя (TBP) метильными заместителями. Результаты объединения представлены в табл. 2.

Таблица 2. Групповые объединения для изомеров с числом атомов углерода, равным 7

Нормальный (NP)

н-гептан

Моноразветвлённый (MBP)

2-метилгексан

3-метилгексан

Диразветвлённый (DBP)

2,3-диметилпентан

2,4-диметилпентан

2,2-диметилпентан

3,3-диметилпентан


Триразветвлённый (TBP)

2,2,3-триметилбутан

Теория активированного комплекса

Согласно теории активированного комплекса (ТАК), значение константы скорости реакции зависит от изменения энтальпии и энтропии, как функции строения промежуточного активированного комплекса и компонентов, вступающих в реакцию.

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

Единичное событие

Рассмотрим принцип «единичного события» на примере реакции изомеризации н-гептана:

Рис. 4. Реакции изомеризации н-гептана

Реакции гидроизомеризации н-гептана могут протекать по двум параллельным направления с образованием 2-метилгексана и 3-метилгексана (Рис. 4). Каждая из этих реакций протекает со своей скоростью, исходя из этого, скорость конверсии н-гептана может быть записана в следующем виде:

–wн-С7 = k1[C н-С7] + k2[C н-С7],                                        (6)

где k – константы скорости реакции, [С] – концентрации н-гептана.

Скорости образования изомеров записываются в следующем виде:

w1 = k1[C н-С7]                                                        (7)

w2 = k2[C н-С7]                                                        (8)

Так как 2-метилгексан и 3-метилгексан входят в одно групповое объединение MB-C7, а н-гептан входит в объединение NP-C7, значение скоростей реакций k1 и k2 согласно ТАК равны (обозначим их как k0 – константа скорости единичного события), уравнение (6) можно записать как:

–wNP-С7 = 2k0[C NP-С7].                                        (9)

Число 2 в уравнении (9) – это так называемое число «единичных событий» – ne, которое показывает количество возможных направлений реакций данного типа и является функцией числа атомов углерода в молекуле парафина/карбокатиона. Таким образом, произведения числа «единичных событий» на константу скорости единичного события даёт нам значение константы скорости реакции для конкретного парафина/карбокатина. Исходя из этого, уравнение (9) можно записать в виде:

wMB-С7 = ne k0[C NP-С7]                                        (10)

В третьей главе представлены дополнения, внесённые в исходную модель для устранения её недостатков.

Конденсация из многокомпонентного газового потока

В принятой модели изначально используется двухплёночная модель массопереноса, использование которой сомнительно, поскольку она основывается на газо-жидкостном равновесии, т. е. предполагает, что все газообразные углеводороды находятся в закритической области и ведут себя как газы. Такое предположение неверно, так как только водород и лёгкие углеводороды C1–C4 находятся в закритической области и ведут себя как газы, в то время как остальные углеводороды ведут себя как пар. Для расчёта парожидкостного равновесия была выбрана методика конденсации из многокомпонентного газового потока. Перед началом расчёта предполагается, что все компоненты смеси находятся в газовой фазе, из которой они конденсируются. Конденсат считается идеальным раствором легких углеводородов и водорода в жидких углеводородах.

Часть входящего потока i-го компонента zi , которая конденсируется, определяется по уравнению:

,                                                (11)

где Pi – давление насыщенного пара i-го компонента, ат; Р – общее давление, ат; S – степень конденсации.

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

Расчёт парциальных давлений компонентов

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

,                                        (12)

где A, B, C, D – параметры уравнения Гарлахера, Т – температура (K), Pi – парциальное давление компонента (мм. рт. ст.).

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

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

Скорость реакции по компоненту i имеет вид:

,                                (13)

где ri – скорость реакции i, моль/м3/ч; – комплексная константа скорости реакции,

ч-1; m, w – тип исходного компонента и тип продукта реакции соответственно; i – число атомов углерода в молекуле; k – тип карбокатиона (вторичный или третичный); ne – число единичных событий; Ci – концентрация карбокатиона на поверхности катализатора, кмоль/м3.

В работах не приведена схема реакций процесса гидрокрекинга парафинов. С использованием последовательности реакций гидрокрекинга и правил составления уравнений гидрокрекинга была получена схема превращений (рис. 5) и система уравнений процесса гидрокрекинга парафинов (14–22), применимая для каждого члена гомологического ряда метана.

Рис. 5. Предложенный механизм реакций гидроизомеризации и гидрокрекинга, где NR(s) – вторичный карбокатион нормального строения,  MBR(s) – вторичный карбокатион c одним метильным заместителем, MBR(t) – третичный карбокатион c одним метильным заместителем, DBR(s) – вторичный карбокатион c двумя метильным заместителями, DBR(t) – третичный карбокатион c двумя метильным заместителями, TBR(t) – третичный карбокатион c тремя метильным заместителями, NO(s) – вторичный олефин нормального строения, MBO(t) – третичный олефин c одним метильным заместителем.

Запишем уравнения скоростей реакций PCP изомеризации (pcp):

                                       (14)

                                        (15)

                                        (16)

                                        (17)

                                        (18)

Уравнения скоростей реакций крекинга (Cr):

                                        (19)

                                        (20)

                                        (21)

                                        (22)


Распределение продуктов гидрокрекинга

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

Рис. 6 Распределение продуктов гидрокрекинга н-декана

Это предположение противоречит экспериментальным данным (рис. 7).

Рис. 7. Распределение продуктов гидрокрекинга н-C10H22 со степенью конверсии 85% на платиновом катализаторе при 3 МПа и соотношении H2:н-C10=4:1

Наложим на точки распределения продуктов гидрокрекинга н-C10H22 (рис. 7) графики непрерывных распределений: нормальное (Гауссовское), Коши, логистическое и лапласовское. Сравнение распределения продуктов гидрокрекинга н-C10 c непрерывными математическими распределениями приведено на рис. 8. Математическое ожидание принято равным 5 – это положение связи, по которой наиболее вероятен разрыв, а параметр масштаба распределения подобран таким образом, что максимумы распределений равны.

Рис. 8. Сравнение экспериментального распределения продуктов гидрокрекинга н-C10 с непрерывными распределениями

Таблица 3. Отклонение распределений от литературных данных

№ ат. C

Отклонение, %

Лаплас

Логистическое

Гаусса

Коши

2

96,42

26,27

1,83

67,54

3

-9,89

-9,32

-4,84

-15,97

4

-24,22

-4,88

-0,39

-15,41

5

0,00

0,00

0,00

0,00

6

-27,52

-9,02

-4,72

-19,08

7

-9,89

-9,32

-4,84

-15,97

8

96,42

26,27

1,83

67,54

Из табл. 3 видно, что минимальное отклонение даёт распределение Гаусса. В дальнейшем именно это распределение используется для описания распределения продуктов гидрокрекинга (данное предположение согласуется с выводами, сделанными в работе).

В четвёртой главе описаны результаты работы.

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

Рис. 9. Изменение количества н- и изопарафинов в зависимости от массы катализатора гидрокрекинга при T=304,4 °С, P=35,5 бар, H2/CnH2n+2=9 моль/моль, сырьё н-C16H34, V=0,43 м3,  конверсия н-C16H34 ~10-15% NP, Iso – мольные проценты нормальных и изопарафинов по литературным данным; NP*, Iso* – мольные проценты нормальных и изопарафинов, рассчитанные с помощью модели

Рис. 10. Изменение количества н- и изопарафинов в зависимости от массы катализатора гидрокрекинга при T=321,3 °С, P=35,5 бар, H2/CnH2n+2=9 моль/моль, сырьё н-C16H34, V=0,43 м3, конверсия н-C16H34 ~10-15%

Рис. 11. Изменение количества н- и изопарафинов в зависимости от массы катализатора гидрокрекинга при T=332,4 °С, P=35,5 бар, H2/CnH2n+2=9 моль/моль, сырьё н-C16H34, V=0,43 м3, конверсия н-C16H34 ~10-15%

Аналогичные расчёты, проведённые для смеси парафинов состава, представленного на рис. 12 условия T=304,5 °С, P=35,5 бар, H2/CH=9 моль/моль, объёмная скорость 1 ч-1.

Рис. 12. Мольное распределение парафинов, подвергаемых гидрокрекингу

В табл. 4 представлено сравнение литературных и расчётных данных:

Таблица 4. Сравнение литературных и расчётных данных гидрокрекинга смеси парафинов

 

Литературные

Расчётные

  Конверсия, % моль

83,28

82,58

  C1-C4, % масс

5,29

2,95

  C5-C9, % масс

19,65

15,63

  C10-C22, % масс

57,79

63,98

  C23+, % масс

17,28

17,44

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

Исследование работы модели процесса гидрокрекинга
в зависимости от кислотности катализатора

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

Чтобы использовать кислотность при моделировании, значение кислотности катализатора было вынесено из комплексной константы скорости реакции единичного события. В результате получено уравнение (23), в котором кислотность фигурирует как параметр моделирования:

  (23)

Из-за отсутствия экспериментальных данных для гидрокрекинга парафинов на цеолитном катализаторе, мы воспользовались данными, приведёнными в диссертационной работе. В ней приводятся зависимости селективности процесса гидрокрекинга (фракции твёрдых парафинов С21С34), проводимого в условиях термостатированного автоклава под давлением водорода с использованием каталитических систем Pd/SA-X (палладий на алюмосиликатах, где X-массовое содержание SiO2). Исследовалось влияние кислотности катализатора на селективность гидрокрекинга. Использование значений кислотности алюмосиликатного катализатора в нашей математической модели позволит выявить влияние кислотности как параметра моделирования процесса гидрокрекинга.

Вместо абсолютного значения кислотности катализатора был введён параметрический коэффициент кислотности KA, который находится как отношение кислотностей образца катализатора () и эталонного катализатора ():

.                                                        (24)

Расчёт процесса гидрокрекинга предлагается вести по следующему уравнению:

                               (25)

При условиях 320 °C, 40 бар, H2/CH=9, различных значений кислотности и состава сырья, представленного в работе, были проведены расчёты селективности по фракциям. Результаты расчётов представлены в табл. 5 и на рис. 13

Таблица 5. Сравнение литературных и расчётных данных селективности по фракциям (% масс)

Катализатор

Кислотность,
моль KOH/г кат

KA

Литературные данные

Результаты моделирования

С1-С4

С5-С9

С10-С20

С1-С4

С5-С9

С10-С20

Pd/SA-88

0,067

1

9,70

38,80

51,50

11,98

38,14

49,89

Pd/SA-31

0,095

1,42

11,40

41,00

47,60

13,58

42,14

44,28

Pd/SA-43

0,126

1,88

12,80

49,10

38,10

15,00

45,50

39,50

Pd/SA-69

0,140

2,09

13,90

48,50

37,60

15,56

46,76

37,69

Pd/SA-57

0,153

2,28

17,90

55,20

26,90

16,04

47,81

36,15

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

Рис. 13. Сравнение результатов моделирования и литературных данных
по селективности выхода топливных фракций

Из представленных результатов видно, что выбранный подход хорошо работает во всём интервале кислотности.

Рассмотрим, как параметрический коэффициент кислотности влияет на конверсию и выход целевых фракций. Для сырья, представленного на рис. 12 и условий: Pt/ZSM-5, 304,5 °C, 35,5 бар, [H2]/[CnHm]=35,5 моль/моль, проварьируем параметрический коэффициент кислотности катализатора и рассмотрим, как это повлияет на процесс гидрокрекинга. Результаты расчётов представлены в табл. 6:

Таблица 6. Результаты варьирования параметрического коэффициента кислотности

  Кислотность, абс

0,25KA

0,50 KA

0,75 KA

KA

1,25 KA

  Конверсия, % масс

22,69

47,16

68,20

82,58

91,14

  C1– C4, % масс

0,41

1,25

2,75

2,95

3,65

  C5– C9, % масс

2,04

6,19

10,89

15,63

20,23

  C10– C22, % масс

20,23

39,71

55,15

63,98

67,25

  C23+, % масс

73,32

52,85

31,21

17,44

8,87

Из табл. 6 видно, что изменение параметра кислотности в 4 раза снижает конверсию углеводородов на 60% и снижает выход целевых фракций – это свидетельствует о том, что кислотность является одним из ключевых свойств катализатора.

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

KA= KA0(1–t)        ,                                                (26)

где KA – текущее значение параметрического коэффициента кислотности, KA0 – исходное значение параметрического коэффициента кислотности,  – скорость отравления катализатора [KA/ч], t – время работы катализатора [ч].

Посчитаем влияние снижения параметрического коэффициента кислотности на процесс гидрокрекинга для сырья, представленного на рис. 12 и условий: Pt/ZSM-5, 304,5 °C, 35,5 бар, [H2]/[CnHm]=35,5 моль/моль, =0,00005 KA/ч (произвольное значение). Результаты моделирования представлены в  табл. 7

Таблица 7. Результаты моделирования процесса гидрокрекинга во времени

Время, ч

10

2000

4000

8000

KA

KA0

0,9 KA0

0,8 KA0

0,6 KA0

Конверсия, % масс

83

78

72

56

С1– С4, % масс

2,96

2,65

2,32

1,62

С5– С10, % масс

15,63

13,74

11,84

8,04

С11– С23, % масс

63,98

61,21

57,43

46,64

С23+, % масс

17,43

22,40

28,41

43,70

Из результатов, представленных в табл. 7, видно, что с течением времени конверсия и выход лёгких фракций снижается. Например, если требуется, чтобы реактор работал с конверсией не менее 70 %, то через 4000 часов работы необходимо провести регенерацию катализатора, либо поднять температуру процесса. С помощью созданной модели можно определить необходимую динамику роста температуры, позволяющую сохранять конверсию сырья на уровне 83 %.

В промышленном реакторе гидрокрекинга процесс протекает в слое катализатора. Необходимо знать, как будет протекать процесс гидрокрекинга в этом слое и на зерне катализатора.

Реакции гидроизомеризации и гидрокрекинга имеют первый порядок по карбокатиону, поэтому уравнение диффузии на зерне катализатора (27)

,                        (27)

где С – концентрация вещества на поверхности катализатора [кмоль/м3], r – радиус зерна [м], Di – коэффициент диффузии [м2/с], имеет аналитическое решение, которое находится с помощью модуля Тиле (28)

,                                                        (28)

где a – размер зерна [м].

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

                                               (29)

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

Коэффициент диффузии парафинов рассчитывался по уравнению (30):

,                                                (30)

где  – доля свободного объёма в пористой грануле; г – радиус пор; Mi – молекулярная масса парафина; – коэффициент извилистости (близкий к двум); T – температура процесса; R – газовая постоянная;  =3,14.

Рис. 14. Зависимость коэффициента использования зерна катализатора от числа атомов углерода в молекуле; для зерна катализатора радиусом 1, 3 и 6 мм; г=3 нм; =0,6; T=300 °С

Из рис. 14 видно, что коэффициент использования зерна катализатора снижается практически до 0 вне зависимости от диаметра зерна катализатора, это значит, что процесс гидрокрекинга протекает только на поверхности зерна катализатора, а внутренний объём катализатора недоступен для парафинов с числом атомов углерода в молекуле больше 4.

С помощью математической модели процесса гидрокрекинга парафинов (катализатор Pt/ZSM-5) рассчитаем габариты реактора. Для реакторов объёмом менее 100 м3 отношение общей высоты к диаметру равно ~3 (реже 4), а отношение цилиндрической части к диаметру от 1,7 до 2,3, такие значения характерны для реакторов с одним слоем катализатора. В качестве сырья примем смесь парафинов, представленную на рис. 12. В таблице 9 приведены условия протекания процесса и результаты расчётов для четырёх значений объёмной скорости.

Таблица 9. Варианты реактора гидрокрекинга парафинов

 

Варианты расчёта реактора гидрокрекинга

1

2

3

4

Расход смеси парафинов, м3

0,61

Объём катализатора, м3

0,75

0,63

0,5

0,37

Объёмная скорость, ч-1

0,8

0,97

1,2

1,63

Температура процесса, °C

321,3

Давление, атм

35,5

Отношение H2/C, нм3/м3

35,5

Диаметр реактора, м

0,78

0,73

0,69

0,62

Высота реактора, м

1,56

1,47

1,38

1,25

Конверсия, % масс

99,99

99,66

97,64

91,36

С1-С3, % масс

10,83

8,04

5,78

3,92

С4-С10, % масс

65,17

56,21

45,71

33,62

С10-С22, % масс

23,99

35,42

46,15

53,83

С23+, % масс

0,01

0,33

2,36

8,64

Из табл. 9 видно, что подбором объёмной скорости удалось получить продукт с различным соотношением бензиновой и дизельной фракции, сохранив при этом высокий уровень конверсии сырья. Таким образом, выбор объёмной скорости процесса необходимо делать на основе экономических расчётов.

Выводы

  1. Разработана математическая модель процесса гидрокрекинга парафинов, позволяющая рассчитать геометрические размеры реактора гидрокрекинга с катализатором типа Pt/ZSM-5;
  2. Предложено использовать в качестве параметра моделирования, характеризующего свойства катализатора, число кислотных центров (кислотность);
  3. Разработанная модель позволяет прогнозировать динамику процесса гидрокрекинга во времени;
  4. На основе выполненных расчётов показано, что процесс гидрокрекинга протекает в тонком приповерхностном слое зерна катализатора, поверхность во внутреннем объёме зерна не вовлечёна в процесс.

Публикации:

  1. Топильников В. И., Сосна М. Х. Об использовании энергии химических реакций // Химическая промышленность сегодня. № 8, 2010. – с. 26-29
  2. Топильников В. И., Сосна М. Х. Моделирование процесса гидрокрекинга парафиновых углеводородов. № 2, 2012. – с. 34-38
  3. Топильников В. И., Сосна М. Х., Лапидус А. Л. Разработка модели процесса гидрокрекинга нормальных парафинов  // Химия твердого топлива. № 2, 2012. – с. 25-31;
  4. I Российский нефтяной конгресс; 14–16 марта 2011, Россия, Москва // «СОЗДАНИЕ МОДЕЛИ ПРОЦЕССА ГИДРОКРЕКИНГА ПАРАФИНОВ»
  5. III Международная конференция «АИСТ-2011», Альтернативные источники сырья и топлива, 24–26 мая 2011, Белоруссия, Минск // «СОЗДАНИЕ МОДЕЛИ ПРОЦЕССА ГИДРОКРЕКИНГА ПАРАФИНОВ»
  6. II Международная конференция газохимия-2011, 5-6 октября 2011 г, Москва, Россия // «Моделирования процесса переработки продуктов процесса ФТ на Co-катализаторе» Топильников, Сосна
  7. Конференция молодых учёных по нефтехимии 13-16 сентября 2011, Звенигород, Россия  «Создание математической модели гидрокрекинга парафинов» Топильников.
 





© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.