Определение радиуса робастной устойчивости автоматической системы регулирования давления в стабилизационной колонне установки каталитического риформинга
- Авторы: Джамбеков А.М.1, Дмитриевский Б.С.2, Терехова А.А.2
-
Учреждения:
- Астраханский колледж вычислительной техники
- Тамбовский государственный технический университет
- Выпуск: Том 30, № 3 (2022)
- Страницы: 31-44
- Раздел: Информатика, вычислительная техника и управление
- URL: https://journals.eco-vector.com/1991-8542/article/view/108975
- DOI: https://doi.org/10.14498/tech.2022.3.3
- ID: 108975
Цитировать
Полный текст
Аннотация
Предложен учет влияния неконтролируемых возмущений на регулирование давления в стабилизационной колонне установки каталитического риформинга изменением коэффициентов полинома знаменателя передаточной функции автоматической системы регулирования давления на основе методов робастной устойчивости. Для описания задачи в общем виде получен общий вид полинома знаменателя передаточной функции автоматической системы регулирования давления в стабилизационной колонне установки каталитического риформинга. Цель исследования – определение максимального размаха изменения коэффициентов полинома знаменателя передаточной функции автоматической системы регулирования давления в стабилизационной колонне установки каталитического риформинга, при котором сохраняется робастная устойчивость. Методы исследования – методы робастной устойчивости, теории автоматического управления. Рассмотрено семейство полиномов знаменателя передаточной функции автоматической системы регулирования давления в стабилизационной колонне установки каталитического риформинга с параметрами неопределенности, изменяющимися в квадрате. Задача определения области робастной устойчивости семейства полиномов сведена к определению областей устойчивости четырех реберных полиномов и пересечения всех областей для определения области робастной устойчивости семейства. В общем виде описана задача определения робастной устойчивости каждого реберного полинома. Для определения области робастной устойчивости семейства полиномов знаменателя передаточной функции автоматической системы регулирования давления в стабилизационной колонне установки каталитического риформинга заданы три полинома, включая номинальный. На основе общего описания задачи определен радиус устойчивости каждого реберного полинома с помощью графического критерия и построения годографа Цыпкина – Поляка. Результаты – на основе пересечения областей устойчивости четырех реберных полиномов определена область робастной устойчивости семейства полиномов знаменателя передаточной функции автоматической системы регулирования давления в стабилизационной колонне установки каталитического риформинга с радиусом устойчивости, равным 1,5. Заключение – данный результат может быть использован при разработке автоматической системы регулирования давления в стабилизационной колонне установки каталитического риформинга и выборе настроек регуляторов.
Полный текст
Введение
Октановое число продуктовой смеси (стабильного катализата) на выходе из стабилизационной колонны (СК) блока стабилизации катализата (БСК) определяет качество всей цепочки процесса каталитического риформинга (КР), протекающего в установке КР [1, с. 8; 2, с. 9; 3, с. 90]. Для эффективного управления БСК помимо разработки эффективных алгоритмов управления, обеспечивающих достижение экономического или иного эффекта, необходимо обеспечение устойчивого режима работы автоматических систем регулирования (АСР) технологических параметров (температуры, давления, расхода и пр.)[4, с. 55; 5, с. 11]. Одним из основных технологических параметров БСК является давление в СК [6, с. 91; 7, с. 67]. На регулирование давления в СК оказывают влияние неконтролируемые возмущения: давление поступающего сырья (нестабильного катализата), перепад температуры в СК [8, с. 7; 9, с. 71].
Данные возмущения необходимо учитывать при разработке АСР давления в СК [10, с.107]. Коэффициенты полинома знаменателя передаточной функции АСР давления в СК постоянны [11, с.127]. В настоящей работе предложен учет влияния возмущений путем изменения коэффициентов полинома знаменателя передаточной функции АСР давления в СК на основе методов робастной устойчивости [12, с. 76; 13, с. 5]. Предполагается, что данные коэффициенты изменяются под влиянием возмущений с течением времени [14, с. 757; 15, с. 172]. Важной является задача исследования робастной устойчивости АСР давления в СК [16, с. 18].
Цель исследования – определение максимального размаха изменения коэффициентов полинома знаменателя передаточной функции АСР давления в СК установки КР, при котором сохраняется робастная устойчивость.
Методы исследования – методы робастной устойчивости, теории автоматического управления.
Постановка задачи
Для исследования робастной устойчивости АСР давления в СК рассмотрим типовую схему, представленную на рис. 1 [17, с. 40].
Рис. 1. Структурная схема АСР давления в СК
На рис. 1 используются обозначения: r(t) – задающее воздействие; e(t) – ошибка управления; u(t) – управление; y(t) – выход; R – регулятор (ПИД-регулятор); OU – объект управления (процесс в СК). В качестве регулируемого параметра y(t) рассмотрено давление в СК. Задающим воздействием r(t) является ступенчатое изменение положения (хода) регулирующего органа (задвижки) на линии (трубопроводе) сброса газов из емкости орошения СК [18, с. 878].
Полином знаменателя передаточной функции АСР давления в СК, определенной на основе экспериментальной переходной кривой для установки КР Л-35-11/1000, имеет вид
. (1)
Для описания задачи в общем виде получим общий вид полинома (1):
, (2)
где ai – постоянные коэффициенты.
Предполагается, что при воздействии возмущений коэффициенты полинома (2) изменяются с течением времени. Для робастной устойчивости АСР давления в СК при воздействии возмущений необходимо определить максимальный размах изменения коэффициентов полинома (2), при котором сохраняется робастная устойчивость, т. е. радиус устойчивости [19, c. 56].
Для рассмотрения наиболее общего случая изменения коэффициентов полинома имеет место следующее.
При изменении коэффициентов полином (2) преобразуется в семейство полиномов (3) знаменателя передаточной функции АСР давления в СК
(3)
с параметрами, изменяющимися в квадрате
, (4)
где ξϵR2 – вектор неизвестных параметров, который принадлежит заданному множеству допустимых значений ξϵΞ (множеству неопределенности); P0(s) – номинальный полином, равный (2); |ξ|2 – евклидова норма вектора ξ
(5)
Одномерное семейство вида (6) носит название реберного полинома, количество полиномов для семейства (3) равно 4:
(6)
Вершинными полиномами называют полиномы вида (7), их количество для семейства (3) равно 4:
(7)
Геометрически вершинные (6) и реберные (7) полиномы соответствуют вершинам и ребрам квадрата (4), т. е. реберный полином «соединяет» два «соседних» вершинных полинома (соответствующих соседним вершинам квадрата) (рис. 2) [20, c. 267].
На рис. 2 изображены четыре вершины: A) ; B) , ; С) ; D) , и четыре ребра: , (AB), (CD); , (AD), (BC).
Пусть имеет место (8) и (9) и полином P0(s) устойчив:
(8)
(9)
где aki (k=0,3; i=1,2) – коэффициенты при sk полиномов семейства (3).
Рис. 2. Вершины и ребра множества неопределенности
Исходя из реберной теоремы для робастной устойчивости семейства (3) необходима и достаточна устойчивость всех реберных полиномов. При этом задача определения границ робастной устойчивости семейства (3) сводится к определению устойчивости четырех интервальных семейств полиномов (реберных полиномов). Пересечение всех областей устойчивости реберных полиномов дает область робастной устойчивости (4) семейства (3).
Задача заключается в определении области неопределенности (4) семейства полиномов знаменателя передаточной функции АСР давления в СК (3) при заданных полиномах P0(s) (устойчивого), P1(s), P2(s), робастной устойчивости четырех реберных полиномов и подтверждении (неподтверждении) выполнения (8), (9), на основании которого делаем вывод о робастной устойчивости (неустойчивости) семейства (3).
Для этого получим реберные полиномы для семейства (3).
Пусть полиномы P1(s), P2(s) в (3) заданы выражениями
, (10)
, (11)
где , – постоянные коэффициенты.
Исходя из (6) получим четыре реберных полинома. При получаем первый и второй реберные полиномы:
, (12)
. (13)
При получаем третий и четвертый реберные полиномы:
, (14)
(15)
Опишем задачу определения робастной устойчивости каждого из четырех реберных полиномов.
Робастная устойчивость первого реберного полинома
В общем виде опишем задачу определения робастной устойчивости первого реберного полинома (12). С учетом (2), (10), (11) интервальное семейство полиномов (12) принимает вид
(16)
Семейство (16) перепишем в виде
(17)
Семейство (17) представим в виде
, (18)
где – коэффициенты семейства (18), зависящие от ξ2 и γ.
Для семейства (18) получим номинальный полином
, (19)
где – постоянные коэффициенты.
Исходя из (18), (19) получим неравенство
. (20)
С учетом (20) запишем интервальное семейство (18) в форме
, (21)
где – масштабы изменения коэффициентов ; γ – размах неопределенности.
Определим радиус устойчивости семейства (21).
Для этого необходимо ввести функции
, , ,
и построить годограф Цыпкина – Поляка
, .
Исходя из графического критерия для робастной устойчивости семейства (21) необходимо и достаточно, чтобы выполнялось условие
, (22)
и годограф при изменении частоты от 0 до ∞ проходил последовательно через 3 квадранта против часовой стрелки и не пересекал квадрата с вершинами .
Построив годограф , определяем радиус устойчивости интервального семейства (12) по формуле
,
где – размер наибольшего квадрата , вписанного в годограф ; , .
Последние равенства вытекают из условия (22).
Аналогично опишем задачу определения робастной устойчивости для остальных реберных полиномов (13)–(15).
Робастная устойчивость второго реберного полинома
В общем виде опишем задачу определения робастной устойчивости второго реберного полинома (13). С учетом (2), (10), (11) интервальное семейство (13) принимает вид
(23)
Семейство (23) перепишем в виде
(24)
Семейство (24) представим в виде
(25)
где – коэффициенты семейства (25), зависящие от и .
Для семейства (25) получим номинальный полином
(26)
где – постоянные коэффициенты.
Исходя из (25), (26) получим неравенство
(27)
С учетом (27) запишем интервальное семейство (25) в форме
(28)
где – масштабы изменения коэффициентов di; γ – размах неопределенности.
Определение радиуса устойчивости семейства (28) выполняется также, как и для семейства (21).
Робастная устойчивость третьего реберного полинома
В общем виде опишем задачу определения робастной устойчивости третьего реберного полинома (14). С учетом (2), (10), (11) интервальное семейство (14) принимает вид
(29)
Семейство (29) перепишем в виде
(30)
Семейство (30) представим в виде
(31)
где di(ξ1)=ai+γci+ξ1bi – коэффициенты семейства (31), зависящие от ξ1 и γ.
Для семейства (31) получим номинальный полином
(32)
где di0=ai – постоянные коэффициенты.
Исходя из (31), (32) получим неравенство
(33)
С учетом (33) запишем интервальное семейство (31) в форме
(34)
где – масштабы изменения коэффициентов di; γ – размах неопределенности.
Определение радиуса устойчивости семейства (34) выполняется так же, как и для семейства (21).
Робастная устойчивость четвертого реберного полинома
В общем виде опишем задачу определения робастной устойчивости четвертого реберного полинома (15). С учетом (2), (10), (11) интервальное семейство (15) принимает вид
(35)
Семейство (35) перепишем в виде
(36)
Семейство (36) можно представить в виде
(37)
где di(ξ1)=ai-γci+ξ1bi – коэффициенты семейства (37), зависящие от ξ1 и γ.
Для семейства (37) получим номинальный полином
(38)
где di0=ai – постоянные коэффициенты.
Исходя из (37), (38) получим неравенство
, . (39)
С учетом (39) запишем интервальное семейство (37) в форме
(40)
где – масштабы изменения коэффициентов di; γ>0 – размах неопределенности.
Определение радиуса устойчивости семейства (40) выполняется так же, как и для семейства (21).
Определение робастной устойчивости первого (12) и третьего (14) реберных полиномов совпадает: для обоих полиномов масштабы изменения коэффициентов равны . Также совпадает определение робастной устойчивости второго (13) и четвертого (15) реберных полиномов: для обоих полиномов масштабы изменения коэффициентов равны .
Определим робастную устойчивость реберных полиномов знаменателя передаточной функции АСР давления в СК.
Робастная устойчивость реберных полиномов знаменателя передаточной функции автоматической системы регулирования давления в стабилизационной колонне
С учетом номинального полинома знаменателя передаточной функции АСР давления в СК (1) зададим полиномы (2), (10), (11):
,
,
.
Не существует конкретных правил по выбору коэффициентов полиномов P1(s), P2(s). В данной работе они задаются в виде малых положительных чисел, примерно равных 10 % от значений коэффициентов полинома P0(s).
Необходимо определить область неопределенности семейства (3) в виде квадрата .
Определяем коэффициенты полиномов (2), (10), (11): a0=1, a1=1,73, a2=1,27, a3=0,61; b0=0,12, b1=0,19, b2=0,15, b3=0,08; c0=0,08, c1=0,15, c2=0,11, c3=0,04.
Исследуем робастную устойчивость первого (12) и третьего (14) реберных полиномов, которым соответствуют интервальные семейства (21) и (34) соответственно.
Получим коэффициенты номинальных полиномов (19), (32) для интервальных семейств (21), (34): d00=a0=1, d10=a1=1,73, d20=a2=1,27, d30=a3=0,61.
Определяем размах изменения коэффициентов семейств (21), (34):
Необходимо определить, возможна ли робастная устойчивость семейств полиномов (21), (34) с номинальными полиномами (19), (32).
Проверим, существует ли γ>0, обеспечивающее выполнение неравенства (22):
.
Таким образом, существует 0<γ<5.
Введем величины
и построим годограф Цыпкина – Поляка:
Найдем точки пересечения с осями:
,
,
,
.
Результаты моделирования представлены на рис. 3.
Рис. 3. Годограф Цыпкина – Поляка для первого и третьего реберных полиномов
Из рис. 3 видно, что существует квадрат со сторонами, равными γϵ(0;5), вписанный в годограф z(jω); значит, семейства полиномов (21), (34) обладают робастной устойчивостью.
Определяем радиус устойчивости интервальных семейств (21), (34):
,
где γ*≈1,5; γ0=d00/α0=5, γ3=d30/α3=5,08.
Квадрат со сторонами, равными радиусу устойчивости 1,5, показан на рис. 3.
Исследуем робастную устойчивость второго (13) и четвертого (15) реберных полиномов, которым соответствуют интервальные семейства (28) и (40) соответственно.
Получим коэффициенты номинальных полиномов (26), (38) для интервальных семейств (28), (40):
d00=a0=1, d10=a1=1,73, d20=a2=1,27, d30=a3=0,61.
Определяем размах изменения коэффициентов семейств (28), (40):
Необходимо определить, возможна ли робастная устойчивость семейств полиномов (28), (40) с номинальными полиномами (26), (38).
Проверим, существует ли , обеспечивающее выполнение неравенства (22):
.
Таким образом, существует .
Введем величины
и построим годограф Цыпкина – Поляка
Найдем точки пересечения с осями:
,
,
,
.
Результаты моделирования представлены на рис. 4.
Рис. 4. Годограф Цыпкина – Поляка для второго и четвертого реберных полиномов
Из рис. 4 видно, что существует квадрат со сторонами, равными γϵ(0;15,25), вписанный в годограф z(jω); значит, семейства полиномов (28), (40) обладают робастной устойчивостью.
Определяем радиус устойчивости интервальных семейств (28), (40):
,
где γ*≈8,5; γ0=d00/α0=25, γ3=d30/α3=15,25.
Квадрат со сторонами, равными радиусу устойчивости 8,5, показан на рис. 4.
Пересечение областей устойчивости четырех реберных полиномов дает радиус устойчивости семейства (3) γ=1,5.
Поскольку выполняются условия (8), (9), полином P0(s) устойчив и четыре реберных полинома обладают робастной устойчивостью, делаем вывод о робастной устойчивости семейства полиномов (3) с параметрами неопределенности, изменяющимися в квадрате .
Заключение
Таким образом, на основе определения робастной устойчивости четырех реберных полиномов рассмотренное семейство полиномов знаменателя передаточной функции АСР давления в СК обладает робастной устойчивостью с радиусом устойчивости, равным 1,5. Данный результат может быть использован при разработке АСР давления в СК и выборе настроек регуляторов.
Об авторах
Азамат Матифулаевич Джамбеков
Астраханский колледж вычислительной техники
Автор, ответственный за переписку.
Email: azamat-121@mail.ru
кандидат технических наук, преподаватель
Россия, 414056, г. Астрахань, пер. Смоляной, 2Борис Сергеевич Дмитриевский
Тамбовский государственный технический университет
Email: azamat-121@mail.ru
профессор кафедры «Информационные процессы и управление», доктор технических наук
Россия, 392000, г. Тамбов, ул. Советская, 106/5, помещение 2Анастасия Андреевна Терехова
Тамбовский государственный технический университет
Email: azamat-121@mail.ru
аспирант
Россия, 392000, г. Тамбов, ул. Советская, 106/5, помещение 2Список литературы
- Ancheyta J. Modeling and simulation of catalytic reactors for petroleum refining. Hoboken: Wiley, 2011. 528 p.
- Smith J.M. Chemical engineering kinetics. St. Louis: McGraw-Hill, 1981. 676 p.
- Дерюгина О.П., Мечик С.В., Трапезников Е.А. Процессы каталитического риформинга и компаундирования как способы повышения октанового числа в бензинах, применяемые в промышленных масштабах // Известия высших учебных заведений. Нефть и газ. 2020. № 3. С. 89–99.
- Gumen M.I. Increasing of the efficiency of the reforming LG-35-11/300 // Petroleum Processing and Petrochemistry. 2001. № 11. P. 54–57.
- Pinheiro C.I.C., Fernandes J.L., Domingues L. Catalytic cracking (FCC) process modeling, simulation and control // Industrial I Engineering Chemistry Research. 2012. № 51 (1). P. 1–29.
- Weekman V.A. Model of catalytic cracking conversion in fixed, moving and fluidbed reactors // Industrial and Engineering Chemistry Process Design and Development. 1968. № 7 (1).P. 90–95.
- Mircea C., Agachi S., Marimoiu V. Simulation and model predictive control of a UOP fluid catalytic cracking // Chemical Engineering and Processing. 2003. V. 42. P. 67.
- Шариков Ю.В., Ткачев И.В. Анализ литературного обзора по моделированию процессов реакторного узла каталитического риформинга // Современные научные исследования и инновации. 2018. № 5 (85). С. 7.
- Дюсембаева А.А., Вершинин В.И. Моделирование каталитического риформинга прямогонного бензина: влияние технологических параметров на состав продуктов риформинга // Катализ в промышленности. 2018. № 5. С. 70–75.
- Кузичкин А.А. Статическая оптимизация процесса каталитического риформинга // Естественные и технические науки. 2017. № 10 (112). С. 106–112.
- Альмохсен М.А.К. Модификация процесса каталитического риформинга бензиновой фракции // Вестник магистратуры. 2018. № 12–3 (87). С. 127–128.
- Ефремов В.В., Шелудько А.Г. Расчет робастно-устойчивости автоматической системы регулирования температуры водообогреваемой спецодежды водолаза // Дизайн и технологии. 2013. № 34 (76). С. 75–79.
- Поляк Б.Т., Щербаков П.С. Робастная устойчивость и управление. М.: Наука, 2002. 303 с.
- Nguang S.K. Robust stabilization of a class of time-delay nonlinear systems // IEEE Transactions on Automatic Control. 2000. V. 45. № 4. P. 756–762.
- Wu L., Zhou W. Delay-dependent robust stabilization for uncertain singular systems with discrete and distributed delays // Journal of Control Theory and Applications. 2008. № 6. P. 171–176.
- Хлебников М.В., Квинто Я.И. Условия робастной устойчивости для семейства линейных дискретных систем с неопределенностями // Проблемы управления. 2020. № 5. С. 17–21.
- Цирлин А.М., Борисов Г.Б. Робастная степень устойчивости и настройки регуляторов для технологических объектов с запаздыванием // Приборы и системы. Управление, контроль, диагностика. 2021. № 10. С. 39–45.
- Ахроров А.А.У., Комилов М.З., Тиллоев Л.И., Кобилов А.Б.У. Эффективность протекания процесса каталитического риформинга // Научный аспект. 2018. Т. 7. № 4. С. 877–879.
- Щеглова А.А. К вопросу о сверхустойчивости интервального семейства дифференциально-алгебраических уравнений // Автоматика и телемеханика. 2021. № 2. С. 55–70.
- Гайворонский С.А., Езангина Т.А., Хожаев И.В., Несенчук А.А. Определение вершинных полиномов для анализа степени робастной устойчивости интервальной системы // Мехатроника, автоматизация, управление. 2019. Т. 20. № 5. С. 266–273.
Дополнительные файлы
![](/img/style/loading.gif)