Синтез системы управления процессом каталитического крекинга нефти с использованием прогнозирующей модели
Б.А. Пашаева Каталитический крекинг нефти
Каталитический крекинг является важнейшим процессом в переработке нефти с целью получения бензина. Установка каталитического крекинга предназначена для переработки вакуумного дистиллята с целью получения компонента высокооктанового бензина. Результаты каталитического крекинга определяются в целом такими показателями, как глубина превращения (конверсии) сырья, выход целевых продуктов и их качество. Целевыми продуктами процесса являются бензин и сжиженный газ. Кокс, хотя и фигурирует в материальном балансе процесса (вместе с потерями), но не выводится из установки и полностью сгорает в регенераторе, обеспечивая тепловой баланс реакторного блока.
Управление процессом каталитического крекинга является сложной задачей. Сложность задачи обусловлена следующими проблемами:
-процесс имеет нелинейный характер;
-объект управления является многосвязным с существенными перекрестными связями;
-разница в постоянных времени для различных подсистем значительна
-существует большое количество неизмеряемых возмущений [1,2].
Одним из подходов к анализу и синтезу системы управления каталитическим крекингом является управление с использованием прогнозирующих моделей [4,5].
В работе рассмотрены вопросы синтеза математической модели процесса каталитического крекинга с использованием прогнозирующих моделей.
Установка каталитического крекинга состоит из двух основных частей: из реактора и регенератора (рис.1). Центральной частью установки каталитического крекинга является реактор, функционирующий следующим образом: сырье проходит через нагреватель, смешивается с катализатором и поступает в вертикальную трубу (райзер), ведущую в нижнюю часть большого сосуда (отстойная часть реактора). Время пребывания сырья в реакторе — несколько секунд, реакция протекает мгновенно. Основной задачей реактора является отделение углеводородов от катализатора. Это происходит в отстойной зоне реактора. Паровая фаза (прореагировавшее сырье) поднимается вверх и проходя через циклоны направляется в ректификационную колонну для дальнейшей очистки и обработки. Твердая фаза (закоксованный катализатор) за счет разности статических напоров катализатора в реакторе и регенераторе самотеком по наклонной транспортной линии поступает в регенератор.
Часть углеводородов, которая во время крекинга превращается в кокс, оседает в виде отложений на катализаторе. Когда поверхность катализатора покрывается отложениями, катализатор становится неактивным (отработанным). Чтобы удалить эти углеродные отложения, отработанный катализатор подают в регенератор, где его смешивают с горячим воздухом. В результате происходит окисление кокса. Процесс, протекающий в регенераторе, называют выжигом кокса.
Восстановленный катализатор выходит из нижней части регенератора. Его можно снова смешать с сырьем и направить в реактор. Таким образом, катализатор находится в непрерывном движении, проходя по циклу крекинг—регенерация.
Рис. 1. -Установка каталитического крекинга
Математическая модель каталитического крекинга нефти
Химические процессы, протекающие в реакторе каталитического крекинга, являются сложными и сопровождаются рядом одновременно протекающих химических реакций. Согласно редуцированной модели реактора [3,6,7], эти реакции можно разделить на три вида:
Реакция 1: Г —-—> О Реакция 2: О —Ь
к 3
->ь
Реакция 3: Г -
где Б - это сырье (гидроочищенный вакуумный дистиллят, иначе газойль), О -целевой продукт (бензин), Ь -легкие газовые фракции (в частности кокс). Первая реакция желаемая, так как бензин является целевым продуктом. Реакции 2 и 3 являются побочными, следовательно, нежелательными реакциями.
Лифт-реактор считается реактором идеального вытеснения, так как соотношение длины реактора к его диаметру велико, а реакция протекает в течение нескольких секунд.
Ниже приведена математическая модель реактора на основе материальных и тепловых балансов.
Материальный баланс по сырью:
dy f 0 2 - Ef
-У- = -R * tc * [COR] * Ф 0 * yf * exp f
dz
R *T0*(1 + 0)
: exp(-a * tc * z)
Материальный баланс по выходному продукту:
dyg , 2 - Ef
= -R1 * tc * [COR] * Ф0 * y f * exp ^ ф f—— * exp(-a * tc * z) -
dz
R *T0*(1 + 0)
- R20 * tc * [COR] * Ф 0 * yf * exp
- Eg
R *T0*(1 + 0)
: exp(-a * tc * z)
Где yf - массовая доля газойля в парах в реакторе; Z - безразмерная переменная длинны; R0 - постоянная скорости реакции образования газойля; Rj - постоянная скорости реакции образования газолина; R2 - постоянная скорости реакции получения углерода; COR - массовое соотношение катализатора к нефти; Ф0 - активность катализатора во входе в реактор; Ef - энергия активации для крекинг газойля; Т0 -
температура сырья при поступлении в реактор; 0 - безразмерная переменная
температуры; tc - время нахождения катализатора в реакторе; а - коэффициент старения
катализатора; у8 - массовая доля газолина в парах в реакторе; Е& - энергия активации для
крекинг газолина.
Тепловой баланс:
М =_________________А* АИ1 * Г_________________* ^у£.
Г0*(Г; * Є,. + А* Г0 * Сро + (1 -А)* * СрВ) й2
где А - массовая доля газойля в сырье ;Л#/ - теплота реакции крекинга газойля; Г0 -
подача сырья (нефть + пар) в реактор; Г. - скорость циркуляции катализатора; Ср, -
теплоемкость катализатора; Сро - теплоемкость газойля; СрВ - теплоемкость пара, ¥в -
скорость подачи пара.
Регенератор можно представить в виде реактора идеального смешения.
Ниже представлена модель регенератора [3,8]:
Материальный баланс по коксу:
4 ( Ш * С„с ) = Г * (Сс - С„с ) - к * у. * С„с * ш аї
Материальный баланс по кислороду: а і+15* а
- Ш * у.) = *Л *( У і - у.)^-г * к * у. * СКС * Ш
а Мс *(1 + а)
Тепловой баланс:
- ЧТГ1 *(Ш * Ср.)] = Ті * Г, * С. + Т *Га*Сра - гГ! *(Г, * Ср. + Г, * Ср,) -
ж и * *
ст *Atr ч*к**Q, *W
со
-(АИШ + ---*АИсо)*
^яс
1 + ст 2 Mc
где W - масса катализатора в регенераторе; Wa - масса воздуха в регенераторе; CR массовая доля кокса в восстановленном катализаторе; Ra - подача воздуха в регенератор; к - коэффициент выжига кокса; CSc - массовая доля кокса в использованном катализаторе; ух - мольная доля кислорода в регенераторе; ст - соотношение CO2/CO в дымовых газах; Mc - молярная масса кокса; T - температура катализатора на выходе из регенератора; Tri - температура на выходе из реактора; Ta - температура воздуха, подаваемого в регенератор; Fa - массовый расход воздуха в регенератор; Са -теплоемкость воздуха; АИсо - тепло, выделяемое при образовании CO; АИсо - тепло,
выделяемое при образовании CO2.
Управление процессом каталитического крекинга с использованием прогнозирующих моделей
Управление с использованием прогнозирующих моделей является одним из современных формализованных подходов к анализу и синтезу управления. Прогнозирующий регулятор состоит из двух компонентов: модель процесса и модуль оптимизации (рис. 2) [10].
Рис. 2. -Блок управления с использованием прогнозирующих моделей С точки зрения алгоритмов регулирования, предпочтительно представление модели в виде вход-выход. Для определения упрощенной модели для процесса
каталитического крекинга в работе идентификационная схема представленная на рис. 3.
Рис. 3. -Схема идентификации
Метод идентификации модели изложен в [9]. Фактически, предложенный авторами метод идентификации состоит из трех основных этапов.
Первый шаг заключается в составлении структурной схемы процесса, которая приведена на рисунке 4, где Тсырья - это температура сырья, тр, Т - температуры в
реакторе и регенераторе соответственно, Сошс - массовая доля кокса на катализаторе, дсырья,
дкат, (д - расходы сырья, катализатора и воздуха соответственно.
г
Рис. 4. Структурная схема процесса каталитического крекинга нефти
Второй шаг заключается в определении вида передаточной функции для каждого канала связи. В этой работе, принимаем, что каждый канал связи может быть охарактеризован апериодическим звеном первого порядка без чистого запаздывания: к
Ж (л) =----т—
Тл * У +1
Где кт -коэффициент усиления; Тл -временная постоянная.
Передаточные функции, связанные с реактором обозначены через букву Н, связанные регенератором обозначены через О, а общие передаточные функции связанные процессом в целом обозначены буквой Б. Передаточные функции реактора и регенератора были определены для статического режима, а передаточная функция всего процесса была найдена путем арифметических вычислений.
Таблица 1
Передаточные функции упрощенной модели процесса каталитического крекинга нефти
^""'-^Выход Вход Температура в реакторе Температура в регенераторе
Темпе ратура сырья ж (л) = 0 23 ттг , . 0.0045л2 + 0.105л + 0.165
0.064л +1 '' 2 1 (л ) 4 3 2 0.00064л4 + 0.029л3 + 0.386л 2 + 1.163л +1
Темпе ратура использован ного катализатора ж (л) = 0 79 ттг , . 0.075л2 + 0.39л + 0.61
) 0.064л +1 '' 22 (у ) 4 3 2 0.00084л4 + 0.045л3 + 0.4л2 + 1.177л +1
Расхо д сырья Ж (л) =- 811*10-4 /ч -2.059*10-4л2 -0.00045л-0.00072
) 0.062л +1 ' 23 \ / 4 3 2 0.00072л 4 + 0.031л 3 + 0.39л 2 + 1.165л +1
Расхо д использован ного катализатора Ж л 1751*10-4 ттг -2.16*10-7л4 -4.12*10-6л3 -3.62*10-5л2
) 0.062л +1 ” 2^ " 0.00037л5 + 0.056л4 + 0.43л3 + 0.89л2 + 1.75л +1 - 8.13*10-5 л - 6.04*10-5 0.00037л5 + 0.056л4 + 0.43л3 + 0.89л2 + 1.75л +1
Расхо д воздуха Ж15(л) = 0 ж„(.) = 00028 2 0.48л +1
Последний шаг заключается в определении параметров передаточной функции для каждого канала, которые корректно будут отражать поведение процессов. Метод определения параметров передаточных функций описан в [9].
В таблице 1 приведены связи, которые описывают упрощенную модель процесса после идентификации.
Схема системы управления с использованием прогнозирующих моделей процесса каталитического крекинга нефти приведена на рис. 5.
Рис. 5. -Схема управления процессом каталитического крекинга нефти с
использованием прогнозирующих моделей
Основная задача каталитического крекинга нефти заключается в максимизации удельного выхода целевого продукта. Эта цель достигается, когда в реакторе протекает реакция с хорошей степенью превращения, а в регенераторе происходит хороший выжиг кокса. На практике, для достижения наилучшего результата протекания реакций в реакторе регулируется температура реактора, а для достижения наилучшего выжига кокса, соответственно температура в регенераторе.
Моделирования процесса каталитического крекинга нефти Моделирование системы управления с использованием прогнозирующих моделей производилось для различных значений управляющих параметров (температуры в реакторе и регенераторе) и возмущений. Для регулятора по оценке с прогнозом с несколькими входами и выходами использовались следующие параметры:
-горизонт прогнозирования р=100 шагов;
-период дискретизации Т=4сек;
-горизонт управления М=19.
Проводились три типа экспериментов:
В первом эксперименте поведение системы исследуется при скачкообразном изменении входных параметров регулятора (температура реактора и температура регенератора). На рисунках 6 и 7 представлены переходные процессы для температур реактора и регенератора, вместе с графиками изменения регулирующих параметров (подача катализатора в реактор и подача воздуха в регенератор). Из графиков видно, что многосвязный регулятор по оценке с прогнозом обеспечивает достаточную точность и быстродействие. Система отслеживает уставки без перерегулирования и с нулевой установившейся ошибкой.
Рис. 6. -Динамические характеристики изменения температуры и расхода
катализатора в реактор
Рис. 7. -Динамические характеристики изменения температуры и расхода воздуха в регенераторе
Во втором эксперименте изменяются возмущающие параметры. На рисунках 8 и 9 приведены графики изменения входных и выходных параметров во времени, при изменении температуры сырья. Из полученных результатов можно увидеть, что система управления не чувствительна к возмущениям (внешние возмущения компенсируются системой управления).
Рис. 8. -Динамические характеристики изменения температуры и расхода
использованного катализатора в реакторе
хІО1
Рис. 9 -Характеристики изменения температуры и расхода воздуха в регенераторе
Рис. 10. -Динамические характеристики изменения температуры и расхода
использованного катализатора в реакторе
0.6
і
хІО5
Рис. 11. -Динамические характеристики изменения температуры и расхода воздуха в регенераторе
Последний эксперимент заключается в изменении параметров регулятора. При изменении периода дискретизации с 4сек до 8сек увеличивается время переходного процесса.
Заключение
Рассмотрены особенности функционирования системы каталитического крекинга
нефти.
Получена математическая модель процесса каталитического крекинга нефти.
Синтезирована система управления для каталитического крекинга нефти.
Разработанная система управления с использованием прогнозирующих моделей демонстрирует удовлетворительные показатели качества процесса.
Система управления не чувствительна к возмущениям, возникающим во время протекания процесса, однако увеличение периода дискретизации может привести к увеличению времени переходного процесса.
Список литературы:
1. Ахметов С. А, Ишмияров М.Х., Веревкин А.П., Докучаев Е.С., Малышев Ю.М. Технология, экономика и автоматизация процессов переработки нефти. - М.: Химия, 2005. - 670с.
2. Гаврилов А.И., Пашаева Б.А. Интеллектуальная система управления
каталитическим крекингом нефти// Интеллектуальные системы: Труды девятого
международного симпозиума/ Под ред. К.А.Пупков 2010. - С. 637-641.
3. Pashayeva B. Mathematical model of the fluid catalytic cracking for work in testing control systems for the cracking plant/ PCI, Baku, Azerbaijan - Vol.1 - 2010. pp. 328331.
4. Mircea C., Agachi S., Marimoiu V. Simulation and Model Predictive Control of a UOP Fluid Catalytic Cracking/ Chemical Engineering and Processing - 2003. - Vol. 42 - 67p.
5. Loeblein C. and Perkins J.D. Structural Design for On-line Process Optimization: Application to a Simulated F^/ AIChe Journal - 1999. - Vol 45 - 1015p.
6. Weekman V. A Model of Catalytic Cracking Conversion in Fixed, Moving and Fluid-Bed Reactors/ Industrial and Engineering Chemistry Process Desing and Development -1968. - 90p.
7. Weekman V. and Nace D.M. Kinetics of Catalytic Cracking Selectivity in Fixed Moving and Fluid bed reactors/ AIChE Journal - 1970. - 397p.
8. Errazu A.F., DeLasa H.I. and Sarti F. A Fluized Bed Catalytic cracking Regenerator model, Grid Effects/ Canadian Journal of Chemical Engineering - 1978. - 191p.
9. Coleman B., Babu J. Techniques of Model-Based Control - Pretice Hall PTTr, 2002. - 576p.
10. Bequette W. Process Control Modeling Desing and Simulation - Pretice Hall PTTr, 2003. - 564p.