Научный журнал

Фундаментальные исследования

ISSN 1812-7339

"Перечень" ВАК

ИФ РИНЦ = 1,674

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Коробочкин В.В. 1 Кравцов А.В. 1 Попок Е.В. 1

---

1 ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», Томск

Проведено исследование поверхностного механизма синтеза метанола на низкотемпературном Zn–Cu–Al-катализаторе. Составлена математическая модель процесса, учитывающая реакции, протекающие на поверхности катализатора. Разработанная математическая модель основана на допущении того, что процесс протекает в кинетической области и реактор работает в режиме идеального вытеснения. С использованием разработанной математической модели проведена оценка эффективности модернизации технологической схемы установки М-750, основанной на внедрении в технологическую схему дополнительного реактора предкатализа. Показано, что оптимальный объем катализатора в реакторе предкатализа составляет 30 м3, при этом возможно увеличение общей производительности установки синтеза на 9,8%. Показано, что в связи с изменением режима работы первой полки катализатора не будет происходить его преждевременной дезактивации.

Статья в формате PDF

508 KB

синтез метанола

математическое моделирование

1. Бесков В.С., Флокк В. Моделирование каталитических процессов и реакторов. – М.: Химия, 1991. – 256 с.

2. Брагинский О.Б. Мировая нефтехимическая промышленность. – М.: Наука, 2003. – С. 172–179.

3. Караваев М.М., Мастеров А.П. Производство метанола. – М.: Химия, 1973. – 160 с.

4. Кравцов А.В. О динамических особенностях механизма реакции гидрирования окиси углерода // Вопросы кинетики и катализа; межвузовский сборник. – Иваново: Наука, 1980. – С. 33–40.

5. Кравцов А.В., Попок Е.В., Юрьев Е.М. Изучение механизма превращения СО на Zn-Cu-Al-катализаторе с помощью квантово-химических расчетов // Известия Томского политехнического университета. – 2011. – Т. 3, №3. – С. 69–72.

6. Кузнецов В.Д., Шуб В.С., Темкин М.И. Кинетика синтеза метанола и гидролиза метанола на медьсодержащем катализаторе. (1. Экспериментальные результаты) // Кинетика и катализ. – 1984. – Т. 25, №3. – С. 606–613.

7. Патент USA №WO2010/146380 А1, 23.12.2010.

8. Попок Е.В., Юрьев Е.М., Кравцов А.В. Моделирование промышленного реактора низкотемпературного синтеза метанола // Фундаментальные исследования. – 2012. – №3(2). – С. 446–451.

9. Рынок метанола. – Электронный интернет-ресурс http:\\methaprocess.ru. – Дата обращения 26.12.2011.

10. Graaf G.H., ScholtensH., StamhuisE.J. Intra-particle diffusion limitations in low-pressure methanol synthesis // Chemical Engineering Science. – 1990. – Vol. 45, №4. –P. 773–783.

11. Graaf G.H., Stamhuis E.J., Beenackersz A.A. Kinetics of low-pressure methanol synthesis // Chemical Engineering Science. – 1988. – Vol. 43, № 12. – P. 3185–3195.

12. Lim H.W., Park M.J., Kang S.H. Modeling of the kinetics for methanol synthesis using Cu/ZnO/Al2O3/ZrO2 catalyst: influence of carbon dioxide during hydrogenation // Ind. Eng. Chem. Res. – 2009. – V. 23, № 48. – P. 10448–10455.

13. Weigel J., Koeppel R., Baiker A. Surface species in CO and CO2 hydrogenation over copper/zirconia: on the methanol synthesis mechanism // Langmuir. – 1996. – № 12. – P. 5319–5329.

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

По состоянию на 2011 год в мире производится более 45 млн т. метанола в год, на Россию приходится порядка 3 млн т. Производительность крупнейших предприятий синтеза метанола в России, работающих по низкотемпературному способу, составляет: ОАО «Метафракс» (г. Губаха, Пермская область) - 1 млн т метанола в год; ОАО «Тольяттиазот» (Самарская область) - более 900 тыс. т в год; ЗАО «Сибметахим» (г. Томск) - 750 тыс. т в год. [9].

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

Целью данной работы стало составление математической модели реакторасинтеза метанола на низкотемпературном Zn-Cu-Al-катализаторе, а также определение изменения производительности от возможных вариантов реконструкции установок синтеза метанола.

Объектом исследования является установка синтеза метанола М-750. Принципиальная технологическая схема приведена на рис. 1.

Свежий синтез-газ I, смешиваясь с циркуляционным газом, разделяется на два потока и поступает к основным реакторам синтеза. Первый поток, пройдя каскад рекуперационных теплообменников, поступает в основные реакторы синтеза 1. Второй поток направляется для создания промежуточного байпасного охлаждения после каждой полки катализатора. Газо-продуктовая смесь, пройдя каскад теплообменников, направляется дополнительно в аппарат воздушного охлаждения и дальнейшее разделение в сепараторе 2. Отделенный газ возвращается на смешение со свежим синтез-газом, смесь метанола и воды (метанол-сырец) направляется на дальнейшее разделение.

При составлении модели использовался поверхностный механизм синтеза метанола, основанный на образовании положительно заряженного хемосорбированного комплекса Ме-СО, приведенный в ра- боте [8].

Для подтверждения выбранного механизма синтеза и установления структуры активного центра авторами была проведена серия квантово-химических расчетов. Результаты расчетов приведены в работе [5]. Также были рассчитаны основные кинетические параметры реакций, протекающих на поверхности катализатора синтеза между адсорбированными веществами. Была подтверждена структура активного центра катализатора, основанная на внедрении в кристаллическую решетку оксида цинка атомов меди, выявлена лимитирующая стадия процесса - адсорбция водорода на поверхности катализатора.

Дальнейшим этапом стал выбор кинетической модели. В научной литературе предложено несколько видов кинетических моделей, с разной точностью описывающих основные закономерности синтеза метанола на лабораторных и промышленных установках. Одной из первых моделей была модель А.А. Караваева [3], в которой лимитирующей стадией была адсорбция водорода на поверхности катализатора. В дальнейшем Ф.С. Шуб и М.И. Темкин описали полученные результаты своими закономерностями [6]. В зарубежной литературе наиболее полное математическое описание процесса приведено авторами [10, 11]. В настоящее время также проводится разработка новых кинетических моделей и уточнение параметров для уже созданных.

В настоящей работе авторы рассматривают схему синтеза метанола с учетом реакций, протекающих на поверхности Zn-Cu-Al-катализатора, брутто-реакции которой имеют следующий вид:

1. СО + 2Н₂ = СН₃ОН;

2. СО₂ + H₂ = CO + H₂O;

3. CH₃OH + CH₃OH = CH₃OCH₃ + Н₂О

4. СО+3Н₂ = СН₄+Н₂О.

Для приведенного механизма синтеза метанола в дальнейших расчетах была выбрана кинетическая модель, приведенная в работе [12]. Значения скорости целевой реакции синтеза метанола (1), конверсии СО (2), синтеза диметилового эфира (3) и образования метана (4) для данной модели имеют следующий вид:

(1)

 (2)

 (3)

 (4)

где k - константа скорости соответствующей реакции; K - константа равновесия по соответствующему веществу; P - парциальное давление соответствующего компонента, Па; С - концентрация соответствующего компонента, моль/м³.

При составлении модели сделано допущение того, что процесс реализуется в кинетической области, что подтверждается литературными данными [3].

В пользу отсутствия внутридиффузионных осложнений говорят малые кинетические диаметры молекул исходных веществ и продуктов по сравнению с размерами пор типичных катализаторов синтеза, а также низкая температура процесса, ключевым образом влияющая на коэффициенты молекулярной и кнудсеновской диффузии в порах катализатора. Катализаторы синтеза метанола,судя по приведенным данным [7], имеют характеристический размер пор 50-70 нм, что значительно больше кинетических диаметров молекул H₂O, СО, СО₂, СH₃OН, которые колеблются от 2 до 5 Å. Таким образом, скорость процесса определяется скоростью химических реакций на поверхности катализатора, имеющих математическое выражение (1) - (4).

Модель реактора составлена с учетом того, что в слое катализатора реализуется режим идеального вытеснения. Выбор гидродинамической модели идеального вытеснения подтвержден высокими значениями диффузионных и тепловых критериев Пекле, принимающих в низкотемпературном синтезе метанола значения около 2000 [1]. Также выбор гидродинамического режима подтверждается промышленными данными с установки М-750, согласно которым не происходит значительного изменения температуры по поперечному сечению катализатора.

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

где C - мольная концентрация, моль/м³; u -линейная скорость, м/с; l - координата реакционной зоны, м; j - порядковый номер реакции; m - общее количество реакций в схеме превращения; ρ - плотность реакционной смеси, кг/м³; C_p - теплоемкость реакционной смеси, Дж/(кг К); T - температура реакционной смеси, К; Q - тепловой эффект реакции, Дж/моль.

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

С использованием промышленных данных был составлен алгоритм дезактивации низкотемпературного Zn-Cu-Al-ката­ли­затора. Одной из основных причин прекращения эксплуатации катализатора в реакторе синтеза метанола является уменьшение его производительности вследствие снижения активности катализатора. Основным фактором снижения активности является перепад температур по слою катализатора.

На созданной программе была проведена серия расчетов для оценки точности составления математической модели и разработки рекомендаций по увеличению производительности установки синтеза метанола. Для проверки математической модели на адекватность использовались проектные и технологические данные с установки типа М-750, результаты расчетов приведены в работе [8].

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

Далее в работе был рассмотрен вариантповышения производительности установки синтеза метанола без увеличения объема катализатора - внедрение в технологическую схему реактора предкатализа. Принципиальная схема установки с реактором предкатализа приведена на рис. 2.

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

Как видно из представленных данных, максимальный выход метанола-сырца наблюдается при объеме реактора предкатализа 30 м³. Дальнейшее увеличение объема катализатора приводит к уменьшению производительности установки, что связано с приближением состава газа на выходе из дополнительно реактора к равновесному и соответственно уменьшению производительности основных реакторов. Максимальное увеличение производительности установки 10522 кг/ч метанола-сырца, что составляет 9,8% от проектной производительности. Рост производительности достигается без дополнительной нагрузки по свежему синтез-газу и без увеличения общего объема катализатора. Использование для нагрева сырья и охлаждения продуктов рекуперационных теплообменных аппаратов позволяет снизить энергетические затраты и, в конечном итоге, увеличить эффективность производства.

Основные технологические параметры установки синтеза метанола при объеме катализатора в реакторе предкатализа 30 м³

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

Таким образом, предлагаемая схема обладает рядом преимуществ перед применяемой в настоящее время двухреакторной технологией М-750:

• процесс ведется дальше от равновесных условий, что позволяет увеличить производительность установки без увеличения объема катализатора;
• в связи с сохранением соотношения синтез-газ/объем катализатора температурный режим основных реакторов изменяется незначительно, что позволяет предотвратить преж­девременную дезактивацию катализатора;
• в технологическую схему добавляется только 2 аппарата, что не требует дополнительных технологических площадей, но вместе с тем существуют и затраты по монтажу и закупке дополнительного оборудования.

На основе современных представлений о механизме синтеза метанола была составлена математическая модель реактора синтеза. Для составления математической модели был проведен анализ существующих представлений о схемах превращения веществ на поверхности катализатора. Выбрана наиболее подходящая модель синтеза. Для выбранной модели проведена серия квантово-химических расчетов, направленных на вычисление кинетических параметров.

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

Рассчитан вариант реконструкции установки, с внедрением в технологическую схему реактора предкатализа, рассчитан оптимальный объем катализатора в данном аппарате. Показано, что в случае такой модернизации производительность увеличивается на 9,8% без увеличения объемного расхода перерабатываемого газа и объема катализатора. Также установлено, что изменение режима не приведет к преждевременной дезактивации катализатора в основных реакторах.

• Курина Л.Н., д.х.н., профессор кафедры Томского государственного университета, г. Томск;
• Сваровский А.Я., д.т.н., профессор Северского технологического института Национального исследовательского ядерного университета «МИФИ».

Работа поступила в редакцию 03.07.2012.

Библиографическая ссылка