Управление температурным режимом химического реактора

УДК: 66.011:681.51

DOI: 10.30987/artide_5c387d62f2a733.72337613

В.Ю. Невиницын, А.Н. Лабутин, Г.В. Волкова

УПРАВЛЕНИЕ ТЕМПЕРАТУРНЫМ РЕЖИМОМ ХИМИЧЕСКОГО РЕАКТОРА

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

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

V.Yu. Nevinitsyn, A.N. Labutin, G.V. Volkova TEMPERATURE MODE MANAGEMENT OF CHEMICAL REACTOR

This aticler solves the problem of synthesizing a nonlinear algorithm for stabilizing the temperature mode in a liquid-phase chemical reactor, ensuring invariance to perturbations, covariance with mastering effects and asymptotic stability of a closed system. The control laws are synthesized using the method of analytical design of aggregated regulators.

Keywords: analytical design of aggregated regulators, synergetic control theory, chemical reactor, control system, computer simulation.

Введение

Центральным звеном производственного процесса превращения исходных веществ в конечные продукты является реакторный узел, который должен быть кибернетически организованным. Это требование означает, что на стадии проектирования решается задача оптимального синтеза реакторного узла, а на стадии эксплуатации подзадача организации оптимального функционирования объекта в условиях действия параметрических и сигнальных возмущений [1]. Вторая подзадача может быть решена только путем разработки системы управления объектом.

Основной особенностью химических реакторов, как объектов управления, является их многомерность, нелинейность и многосвязность. Задача синтеза оптимальных систем управления такого рода динамическими объектами в полной мере не решена [2]. Вопросы разработки систем автоматической стабилизации и управления реакторами рассматривались в основном в линейной постановке [3] и лишь для ряда простых случаев - в нелинейной постановке [4].

Одним из перспективных способов решения рассматриваемой задачи, является метод аналитического конструирования агрегированных регуляторов (АКАР), базирующийся на синергетических принципах целевой самоорганизации нелинейных динамических систем [5].

Использование идей синергетики в задачах управления требует перехода от непредсказуемого поведения системы по алгоритму диссипативной структуры к направленному движению к целевому аттрактору (желаемому инвариантному многообразию) и дальнейшему движению вдоль многообразия в конечное состояние. Данный способ самоорганизации называется направленным или целевым [5]. Эффективность алгоритмов, синтезированных с применением синергетической теории управления, показана в ряде работ [6 - 9].

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

1. Описание технологического процесса

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

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

Л+В —^^Р1, А+Р1—Р2, А+Р2 —^Р3,

где А и В - исходные реагенты; Р1, Р2, Р3 - продукты реакции; к1, к2, к3 - константы скоростей стадий. Целевым компонентом является вещество Р2. В общем случае при изменении спроса на продукты реакции целевым может быть и другой компонент или смесь ряда веществ. Аппарат функционирует в политропическом режиме. Исходные реагенты А и В подаются в аппарат раздельными потоками. Смесь из реактора забирается насосом. Поскольку в аппарате протекает экзотермическая реакция, в рубашку реактора подается хладоагент для охлаждения реакционной массы.

Рис. 1. Принципиальная схема химического реактора

На рис. 1 введены следующие обозначения: х^, х^ - концентрации исходных

в х1 в х2

реагентов; и1, и2 - расход исходных реагентов; х6 , х6 - температуры потоков исходных

реагентов; х", ху - температуры хладоагента на входе и выходе из аппарата; охл - расход

хладоагента на входе и выходе из аппарата; х6 - температура реакционной смеси в аппарате;

V - расход реакционной смеси на выходе из аппарата; х1, х2, х3, х4, х5 - концентрации

компонентов А, В, Р1, Р2, Р3 в реакторе; V - объем реакционной смеси в аппарате; Уш -

объем хладоагента в рубашке.

Задача управления химическим реактором заключается в стабилизации температуры смеси в аппарате на заданном уровне хб в условиях действия возмущений. Управляющим

воздействием является расход хладоагента, подаваемый в рубашку.

Система дифференциальных уравнений, характеризующая модель динамики объекта, имеет вид:

dZ = Rl + MA-h2 Xl -Ьз Xl, dz

dx-2 = R2 + MB h2 X2 h3 X2 ,

dz

dx3 = R

= R3 h2 x3 h3 x3,

dz

4 = R4 -b2 x4 -h3x4,

¿X,

йт

¿Г = Я5 Ъ2 Х5 Ъ3 Х5,

йт

йX

— ОС^к^Х1Х2 ^ ОС2^2Х^Хз ^ ^3^3Х\Х4 ^ ЬХ6 ^ Д1Х7 (Д1 ^ ^2)Х6 ^ (Х6

йт

— Д(Х6 -х7) + Ъ^х76х-х7)м,

йт

где — к х1 Х2 к2 х1 Х3 кз х1 Х4, ^2 к1 Х1Х2, Я — кх X Х2 Х^ Х3 , ^4 — к2 Х^ Х3 кз Х1Х4 , Я — к3ххх4 - скорость реакции по компонентам; МА — и1х¡х / V; Мв — и2хв2х / V; Ъ1 — ИУл;

Ъ^—и^/У; Ъъ —и2/V; а,. — АН,. /(рС), I — 1.....3; Д — ВД/(рСУ); Д =КТГТ /(ршСшУш);

АН,., . —1,...,3 - тепловой эффект соответствующей стадии реакции; р, С - плотность и теплоемкость реакционной смеси; рл, Сл - плотность и теплоемкость хладоагента; Кт -коэффициент теплопередачи; ¥Т - поверхность теплообмена; к1 — кю • ехр(— Е1 / Я(х6 + 273)), . —1,...,3 - константы скоростей стадий; к,0, . —1,...,3 - постоянный множитель (предэкспонента) констант скоростей стадий; Е1, . —1,...,3 - энергия активации соответствующей стадии реакции; Я - универсальная газовая постоянная; и — ил -регулирующее воздействие.

2. Аналитическое конструирование регулятора температуры

Поскольку математическая модель объекта (2) содержит одно внешнее управляющее воздействие и — и, необходимо использовать метод АКАР на основе последовательной совокупности инвариантных многообразий [5]. Анализ уравнений (2) показывает, что управляющее воздействие и — ил воздействует на переменную х6 через переменную х7 . Таким образом, канал управления температурой смеси в аппарате в развернутом виде запишется: и ^ х7 ^ х6 .

Согласно методу АКАР [5], в фазовом пространстве динамических систем можно построить ряд многообразий, к которым притягиваются фазовые траектории. Отсюда следует возможность конструирования такой совокупности притягивающих инвариантных

многообразий (X) — 0, ^ — 1, т, когда изображающая точка системы, начав двигаться из произвольного начального положения в фазовом пространстве, последовательно перемещается от одного многообразия к другому, пока не попадет на последнее ц/т (X) — 0, приводящее в заданное конечное состояние. Таким образом, изображающая точка сначала сближается с многообразием ц/1(X) = 0, затем с у/2(X) — 0 и т.д. При использовании т притягивающих многообразий размерность каждого .-го многообразия будет на единицу меньше предыдущего, вследствие чего происходит сжатие фазового объема и динамическая декомпозиция задачи.

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

у = х7 + у^) = 0,

где у,(х6) - неизвестная функция от х6. Закон управления синтезируется таким образом, чтобы осуществлялся перевод изображающей точки системы в фазовом пространстве из произвольного начального положения в окрестность многообразия у = 0. Изменение агрегированной макропеременной, играющей роль параметра порядка, должно подчиняться функциональному уравнению

Ту + у = 0.

Уравнение (4) с учетом (3) примет вид

t-

dxy ôv, dx6

--i----

d— ôx6 d—

В силу уравнений объекта (2) это выражение запишется

8 у,

T

ß2

ôx6

+ x7 + v, = О,

где /6 = а1к1х1х2 + а2к2х,х3 + а3к3х,х4 + Ь2хв6л - (Д + Ь2)х6 + (х^2 - х6)Ь3. Из (5) получаем выражение для закона управления

и =

Т1Ь1(х7х - х7) Ь,(х7вх - х7) 8х6 Ь,(х7вх - х7).

Управление и переводит изображающую точку системы в окрестность многообразия у = 0, на котором реализуется связь х7 = -у, и наблюдается эффект «сжатия фазового пространства», т.е. снижение размерности системы уравнений (2). Уравнения декомпозированной системы с учетом соотношения х7 = -у, примут вид:

= Я + МА -Ь2х -Ь3х,,

а—

^ = Я + Мв - Ь2 х2 - Ь3 х2, ат

—— = Я3 - Ь2 х3 - Ь3 х3 ,

аахт

-— = Я4 - Ь2 х4 - Ь3 х4 ,

а—

-х5 ри и

— = Я5 - Ь2 х5 - Ь3 х5 ,

а—

= /6 -Ду.

d—

Функцию у(х6) в декомпозированной системе (7) можно рассматривать как «внутреннее» управление, под воздействием которого происходит движение объекта (7) вдоль многообразия ц1 — 0. На втором шаге процедуры синтеза закона управления осуществляется поиск выражения для у(х6) . Для этого вводится в рассмотрение цель движения системы (7) в форме инвариантного многообразия, отражающего технологическое требование к системе

ц2 — Х6 —Х6 —

Макропеременная ц2 удовлетворяет решению функционального уравнения Т2ц2 + ц2 — 0, которое в развернутом виде с учетом выражения (8) в силу модели декомпозированной системы (7) примет вид:

Т2

«Внутреннее» управление в соответствии с выражением (9) запишется

1 Т2Д1 Д

Закон (10) обеспечивает асимптотическое приближение изображающей точки ко второму притягивающему многообразию ц2 — х6 — Х6 — 0.

Окончательное выражение для закона внешнего управления и( Х) можно получить путем подстановки в (6) функции У1 и ее частной производной ду/ дх6. Частную производную ду / дх6 можно получить как аналитическим способом, так и с помощью численных методов. Параметрами настройки закона управления, влияющими на качество динамики процессов в замкнутой системе «реактор - управляющее устройство», являются постоянные времени Т1, Т2.

Из структуры закона управления температурой (6), (10) видно, что параметр Т1 оказывает непосредственное влияние на регулирующее воздействие и — и,, а параметр Т2 оказывает влияние на «внутреннее» управление У1, которое, в свою очередь, влияет на регулирующее воздействие и — и,. В работе [5] показано, что время попадания изображающей точки системы из произвольного начального положения в окрестность заданного инвариантного многообразия приближенно определяется следующим выражением:

ц*

где Т - параметр основного функционального уравнения (4).

Таким образом, выбор настроечных параметров Т1 , Т2 осуществляется на основании требуемого времени перевода изображающей точки системы в окрестность инвариантных многообразий ц1 — 0 и ц2 — 0, т.е. на основании заданного времени переходного процесса в соответствии с оценкой (11). В исследуемом варианте системы управления время переходного процесса (время регулирования) будет определяется в соответствии с выражением:

тр *

3. Компьютерное моделирование и результаты

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

Моделирование проводилось при технологических и конструктивных параметрах, обеспечивающих оптимальный режим работы химического реактора [10]: V = 500 л; Ухл = 290 л; х1х = 19.74 моль/л; х" = 10.93 моль/л; ц = 1.5 л/мин, и2 = 3.5 л/мин, и = 5.0 л/мин, ихл = 3.84 л/мин; ле6хХ = 20 °С; хвх2 = 30 °С; х7вх = 20 °С; Кт = 12 кДж/(м2 мин К); ^ = 2.9 м2; р = 0.9 кг/л; С = 2 кДж/(кг К); рш = 1 кг/л; Схл = 4.18 кДж/(кг К); АН1 = АН2 = АН3 = 80 кДж/моль; Е1 = 48635 Дж/моль; к10 = 109860 л/(моль мин); соотношения констант скоростей последовательных стадий к 2 /к1 = 2.0, к3 /к1 = 2.5. Параметры закона управления (6), (10): значение постоянных времени Т1 = Т2 = 15 мин; заданное значение температуры смеси в аппарате х6 = 140 °С.

На рис. 2, 3 приведены примеры переходных процессов в замкнутой системе «химический реактор - регулятор температуры» при начальном отклонении переменных состояния объекта от статики (Ал1 = —0.2х0) и ступенчатом возмущении по расходу и равном Аи1 = —0.1ц. Допустимая погрешность регулирования температуры в статике (точность в установившемся режиме) составляет ± 2 °С. Для наглядности движение системы (переходные процессы) до момента приложения входного воздействия (т = 50 мин) приводится в статическом режиме.

150

140 О 130

о

120

110 100

50

/

X

s

Ц

а

¡0

50 100 150 200 т, мин

250

300

50 100 150 200 т, мин

250

300

Рис. 2. Изменение регулируемой переменной и управляющего воздействия при начальном отклонении

переменных состояния от статики

143

142

141

и

о 140

о

* 139

13 8

13-7

'0

s s

.s

s

50 100 150 200 т, мин

250

300

5 5 0

4

'3 Г

2 *

50 100 150 200 т, мин

250

300

Рис. 3. Изменение регулируемой переменной и управляющего воздействия при ступенчатом возмущении

по расходу иг

8

2

0

0

0

0

Заключение

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

Список литературы:

1. Лабутин А.Н. Оптимизация гибких многопродуктовых реакторных систем непрерывного типа // Изв. вузов. Химия и хим. технология. 1999. Т. 42. № 1. С. 117-122.

2. Красовский А.А. Развитие и становление современной теории управления // Синергетика и проблемы теории управления / Под ред. А.А. Колесникова. М.: ФИЗМАТЛИТ, 2004. С. 13-34.

3. Рей У. Методы управления технологическими процессами: Пер. с англ. М.: Мир, 1983. 368 с.

4. Кван Н.В., Семичевская Н.П. Нелинейное робастное управление двухэтапным химическим реактором // Информатика и системы управления. 2011. № 4(30). С. 133-141.

5. Колесников А.А. Синергетическая теория управления. - М.: Энергоатомиздат, 1994. 344 с.

6. Labutin A.N., Nevinitsyn V.Y. Analytical Synthesis of Chemical Reactor Control System // Theor. Found. Chem. Eng. - 2014. Vol. 48. N. 3. P. 296-300. DOI: 10.1134/S0040579514030105.

7. Лабутин А.Н., Невиницын В.Ю. Синтез нелинейного алгоритма управления химическим реактором с исполь-зованием синергетического подхода // Изв. вузов. Химия и хим. технология. 2017. Т. 60. № 2. С. 38-44.

8. Labutin A.N., Nevinitsyn V.Y. Analytical Synthesis of Chemical Reactor Control System // International Journal of Advanced Studies. 2016. Vol. 6. N. 1. P. 27-37. DOI: 10.12731/2227-930X-2016-1-27-37.

9. Лабутин А.Н., Невиницын В.Ю., Деветьяров А.Н., Волкова Г.В. Синтез эффективного комплекса "реактор - управляющая система" с использованием синергетического подхода // Химическая промышленность. 2014. Т. XCI. № 2. С. 63-67.

10. Невиницын В.Ю., Лабутин А.Н., Волкова Г.В., Деветьяров А.Н. Системный анализ химического реактора как объекта управления // Изв. вузов. Химия и хим. технология. 2017. Т. 60. № 9. С. 92-99.

References:

1. Labutin A.N. (1999). Optimization of flexible multi-product reactor systems of continuous type, Izvestiya vuzov. Chemistry and Chemical technologies, Vol. 42. (1), pp. 117-122. [in Russian language]

2. Kolesnikova A.A. (Ed.), Krasovsky A.A. (2004). Development and formation of modern management theory. Synergetics and problems of control theory. Moscow: FIZMATLIT, pp. 13-34. [in Russian language]

3. Rey U. (1983). Advanced process control. Moscow: Mir. [in Russian language]

4. Kwan N.V., Semichevskaya N.P. (2011). Nonlinear robust control of a two-stage chemical reactor. Informatics and control systems. 30 (4), pp. 133-141. [in Russian language]

5. Kolesnikov A.A. (1994). Synergetic theory of management. Moscow: Energoatomizdat. [in Russian language]

6. Labutin A.N., Nevinitsyn V.Y. (2014). Analytical Synthesis of the Chemical Reactor Control System. Theor. Found. Chem. Eng, Vol. 48, (3), pp. 296-300. DOI: 10.1134 / S0040579514030105. [in Russian language]

7. Labutin AN, Nevinitsyn V.Yu. (2017). Synthesis of a nonlinear algorithm for controlling a chemical reactor using a synergetic approach. Izvestiya vuzov. Chemistry and Chemical technologies, (2), pp. 38-44. [in Russian language]

8. Labutin A.N., Nevinitsyn V.Y. (2016). Analytical Synthesis of the Chemical Reactor Control System. International Journal of Advanced Studies, vol. 6.p, (1), pp. 27-37. DOI: 10.12731 / 2227-930X-2016-1-27-37. [in Russian language]

9. Labutin AN, Nevinitsyn V.Yu., Devetyarov A.N., Volkova G.V. (2014). Synthesis of the effective «reactor-control system» complex using a synergetic approach. Chemical Industry, T. XCI, (2), pp. 63-67. [in Russian language]

10. Nevinitsyn V.Yu., Labutin A.N., Volkova G.V., Devetyarov A.N. (2017). A System Analysis of a Chemical Reactor as a Control Object. Izvestiya vuzov. Chemistry and Chemical technologies, Vol. 60, (9), pp. 92-99. [in Russian language]

Сведения об авторах:

Невиницын Владимир Юрьевич

кандидат технических наук, доцент кафедры «Техническая кибернетика и автоматика» Ивановского государственного химико-технологического университета. Служебный адрес: 153000, г. Иваново, Шереметевский проспект, 7. Тел. рабочий: 8-4932-32-72-26. Тел. сот: 8-915-837-94-53. E-mail: nevinitsyn@gmail.com

Лабутин Александр Николаевич

доктор технических наук, профессор кафедры «Техническая кибернетика и автоматика» Ивановского государственного химико-технологического университета. Служебный адрес: 153000, г. Иваново, Шереметевский проспект, 7. Тел. рабочий: 8-4932-32-72-26. Тел. сот: 8-910-985-43-05. E-mail: lan@isuct.ru.

Волкова Галина Витальевна

кандидат технических наук, доцент кафедры «Техническая кибернетика и автоматика» Ивановского государственного химико-технологического университета. Служебный адрес: 153000, г. Иваново, Шереметевский проспект, 7. Тел. рабочий: 8-4932-32-72-26. E-mail: konf gv@mail.ru

Статья поступила в редколлегию 10.09.18. Рецензент: д.т.н., доцент Брянского государственного технического университета

Петрешин Д.И. Статья принята к публикации 12.12.18.

Information about authors: Nevinitsyn Vladimir Yurievich

Candidate of Engineering Sciences, Associate Professor of the Department of «Technical Cybernetics and Automation» of the Ivanovo State University of Chemistry and Technology.Office address: 153000, Ivanovo, Sheremetevsky Avenue, 7. Tel. working: 8-4932-32-72-26. Tel. mob.: 8-915-837-94-53. E-mail: nevinitsyn @ gmail .com

Labutin Alexander Nikolaevich

doctor of Engineering Sciences, Professor of the Department of «Technical Cybernetics and Automation» of the Ivanovo State University of Chemistry and Technology.

Office address: 153000, Ivanovo, Sheremetevsky Avenue, 7.

Tel. working: 8-4932-32-72-26. Tel. mob.: 8-910-985-43-05. E-mail: lan@isuct.ru

Volkova Galina Vitalievna

Candidate of Engineering Sciences, Associate Professor of the Department of «Technical Cybernetics and Automation» of the Ivanovo State University of Chemistry and Technology.

Office address: 153000, Ivanovo, Sheremetevsky Avenue, 7.

Tel. working: 8-4932-32-72-26. E-mail: konf gv@mail.ru

Учредитель: Федеральное государственное бюджетное образовательное учреждение высшего образования "Брянский государственный технический университет" Адрес редакции и издателя: 241035, Брянская область, г. Брянск, бульвар 50 лет Октября, 7 ФГБОУ ВО «Брянский государственный технический университет» Телефон редакции журнала: (4832) 56-49-90. E-mail: aim-pu@mail.ru Вёрстка A.A. Алисов. Корректор А.Ю. Малюкина. Сдано в набор 10.12.2018. Выход в свет 27.12.2018. Формат 60 х 88 1/8.Бумага офсетная. Усл. печ. л. 5,88.

Тираж 500 экз. Свободная цена. Отпечатано в лаборатории оперативной полиграфии Федерального государственного бюджетного образовательного учреждения высшего образования "Брянский государственный технический университет" 241035, Брянская область, г. Брянск, ул. Институтская, 16