CC BY
58
Автоматика. Информатика. Управление. Приборы
УДК 661.9
Б01: 10.17277/уе81тк.2018.04.рр.556-568
ОПТИМИЗАЦИЯ И УПРАВЛЕНИЕ ЦИКЛИЧЕСКИМ ПРОЦЕССОМ АДСОРБЦИОННОГО ОБОГАЩЕНИЯ ВОЗДУХА КИСЛОРОДОМ
В. Г. Матвейкин1, С. А. Скворцов1, Е. И. Акулинин2, С. И. Дворецкий2
Кафедры: «Информационные процессы и управление» (1); «Технологии и оборудование пищевых и химических производств» (2), ФГБОУ ВО «ТГТУ», г. Тамбов, Россия; ipu@ahp.tstu.ru
Ключевые слова: воздух; динамика; изотерма адсорбции; кинетика; кислород; короткоцикловая безнагревная адсорбция; математическая модель; оптимизация; управление; цеолитовый адсорбент.
Аннотация: Поставлена и исследована задача оптимизации (в смысле определения максимума степени извлечения кислорода) процесса обогащения воздуха кислородом в двухадсорберной установке короткоцикловой безнагревной адсорбции (КБА) с учетом выполнения требований по чистоте кислорода, производительности установки и ресурсосбережению гранулированного адсорбента. Разработаны алгоритм и двухуровневая система адаптивного управления процессом обогащения воздуха кислородом, реализация которых обеспечит автоматический контроль возмущающих воздействий, оперативное решение задачи оптимизации (на верхнем уровне) и коррекцию заданий программному задатчику и регуляторам автоматической системы управления (на нижнем уровне). Использование адаптивной оптимизации в системе управления позволит повысить значение степени извлечения кислорода в установке КБА в среднем на 10 % (по сравнению с использованием традиционной системы управления без адаптивной оптимизации) при выполнении требуемых ограничений по чистоте кислорода, ресурсосбережению гранулированного адсорбента и производительности установки.
Обозначения
ak - величина сорбции ^го компонента в адсорбенте, моль/м3; а\ - концентрация адсорбата на поверхности раздела фаз или равновесная текущей величине концентрации Ck адсорбтива в газовой фазе, моль/м3; Ь - вектор параметров изотермы сорбции;
Ck - мольная концентрация ^го компонента газовоздушной смеси, моль/м3;
сp , сp - удельные теплоемкости адсорбента и газовоздушной смеси соответственно, Дж/(кг К);
Dа, dgГ - диаметры слоя и гранул адсорбента соответственно, м;
G1П - расход атмосферного воздуха в питании установки КБА, л/мин; О - расход на выходе продукционной линии установки КБА, л/мин;
QOUt
расход газовоздушнои смеси,
обедненной кислородом, на отводноИ линии установки КБА;
h/tds - теплота адсорбции k-го компонента газовой смеси, Дж/моль; Mg - молярная масса газовоздушной
смеси, кг/моль; nin nin
"ads, "des - давление на стадиях адсорбции и десорбции соответственно; R - универсальная газовая постоянная, Дж/ (моль-К);
Q - производительность установки КБА, л/мин;
5"ud - коэффициент удельной поверхности гранул адсорбента, м2/м3; /ц, /щ - продолжительность цикла и
полуцикла «адсорбция - десорбция» периодических стадий адсорбционного обогащения воздуха кислородом соответственно, с;
/уст - время выхода установки на установившийся режим, с; /ф - время функционирования, с;
Pout - давление на продукционном выходе установки КБА, Па; „out
P - давление на линии отвода газовоздушной смеси, обедненной кислородом, Па;
Ta, Tg, ТОС - температуры адсорбента,
воздуха и окружающей среды соответственно, K;
Tgou/ - температура на выходе продукционной линии установки КБА, К; У ос - состав воздуха окружающей среды;
.out
y - состав на выходе продукционной линии установки КБА; yout - состав газовоздушной смеси, обедненной кислородом, на отводной линии установки КБА; а - коэффициент теплоотдачи от поверхности гранул адсорбента к потоку газовоздушной смеси, Вт/(К-м2); L - высота (длина) слоя адсорбента, м; вк - эффективный коэффициент массопе-реноса k-го компонента, отнесенный к концентрации адсорбата (кинетический коэффициент);
е - коэффициент порозности адсорбента с учетом пористости гранул, м3/м3; Y - тип адсорбента;
1 a, 1g - коэффициенты теплопроводности адсорбента и газовоздушной смеси соответственно, Вт/(м-К); X - коэффициент обратного потока обогащенного кислородом воздуха; е0 - порозность слоя адсорбента без учета пористости частиц, м3/м3; vg - скорость воздуха, м/с;
pa - плотность адсорбента, кг/м3; р g - мольная плотность газовоздушной смеси, моль/м3;
^ g - динамическая вязкость газовоздушной смеси, Па-с;
q - коэффициент сферичности гранул
адсорбента; с - форма адсорбента; ук. - закон открытия i-го регулирующего клапана;
П - степень извлечения кислорода, %.
Введение
В последние десятилетия все более распространенным становится использование циклических адсорбционных процессов для разделения газовых смесей и концентрирования целевых продуктов (короткоцикловая безнагревная адсорбция (КБА)), которые не предполагают наличия внешнего источника тепла [1, 2].
Короткоцикловые процессы адсорбционного разделения смесей газов широко применяют в промышленности для извлечения водорода, обогащения воздуха кислородом, безнагревной осушки газов, разделения углеводородов, концентрирования диоксида углерода, извлечения метана и т.п. Одной из актуальных задач в области адсорбционного разделения является обогащение воздуха кислородом. Типичными
веществами, сопутствующими кислороду, являются азот, аргон и диоксид углерода. Особенностью адсорбционного концентрирования кислорода является тот факт, что азот имеет более высокие значения адсорбционной селективности и является более сорбируемым компонентом.
Анализ многочисленных работ зарубежных и российских ученых в области адсорбционного разделения многокомпонентных газовых смесей и концентрирования целевого продукта (водорода, кислорода, азота, углекислого газа и др.), оптимизации режимных переменных технологического процесса, управления и автоматизации позволил определить место данной статьи среди других работ, ее актуальность и перспективность.
Цель работы - разработка алгоритма оптимизации и системы автоматического управления циклическим адсорбционным процессом обогащения воздуха кислородом, обеспечивающих поддержание максимальной степени извлечения кислорода при выполнении требований технологического регламента по чистоте кислорода, производительности установки и ресурсосбережению гранулированного адсорбента.
Анализ процесса концентрирования кислорода методом адсорбционного разделения воздуха как объекта управления
Анализ процесса адсорбционного разделения воздуха в двухадсорберной установке позволил определить (рис. 1):
- в качестве управляющих переменных и(/): /ц; /^(О, Р&з(/); или законы
Ук1 (t),Ук4 (t), t е [0;0,5/ц], ¥к2 (0, ¥к3 (t), t е [0,5/ц;/ц] открытая регулирующих
ц
V/' ' ^ Г' ТК^^ ^ ;гц
«впускных» К], К2 и «сбросных» К3, К4 клапанов; коэффициент х обратного потока обогащенного кислородом воздуха, подаваемого на регенерацию адсорбента параллельного адсорбера через клапан К5;
- в качестве возмущающих воздействий г(/), / е [0, /ф]: Уос, Тос, Вос ;
на отрезке времени [0, /ф]; i|out = Вос, то есть z = {уос, Тос, Вос, P°ut, P°ut};
Gin
;
- выходные переменные установки КБА (y°ut, G°ut, Tg°ut, Q, n, y°ut, G°ut).
Будем считать, что существенные изменения возмущающих воздействий z(t), t е [0, /ф] происходят в моменты времени ti, i = 1,2,... и далее на интервалах
времени [ti, ti+1], i = 1,2, ... значения возмущений не изменяются. Причем, на данных интервалах времени намного превышают отрезок времени [0, /уст], на котором, собственно, и осуществляется решение задачи оптимизации неустановившихся режимов работы установки КБА, вызванных возмущением z(t) в момент
времени ti. В этом случае возмущения z(t), t е [0, /ф] на интервале времени [0, /ф] удобно представлять в форме кусочно-постоянных функций z(tj) = z1 = const, i = 1,2,3... Аналогично, непрерывные функции управлений Ук1 ((),УK4 (t), t е [0;0,5/ц], ук2 (t), ¥к3 (t), t е [0,5/ц;/ц] на отрезке времени [0, /ц] будем аппроксимировать кусочно-постоянными функциями:
¥1 (t) = уK1 = ¥к1 (tj, ), ¥4(t) = ¥КК:4 = ¥k4 (tj, ), tj е [0; 0,5/ц], j = 1, m; V2(t) = ViK'2 = ¥k2 (tj, z), ¥3'(t) = ¥K3 = ¥к3 (tj, Z), tj е [0,5/ц;гц], j = 1m
Окружающая среда
е
f G" t
t* г ц к х <u ч
а
а с
PП
1 ad
P"
de
X
Л X
ч о х m
У = У
ОС
V
«
s X <u
в
B P0*
ОС
rrnn _ ГТ1
l ГГ Г»
V pout = в
л
r out
G01
П
T out
>У1
GO
<u
3
И
и
<u
<u ft <u С
У *
f Л X
4 о х 3
pq
J
Конструктивные параметры
Рис. 1. Схема двухадсорберной установки КБА для концентрирования кислорода:
А - адсорбер; Вос - барометрическое давление воздуха окружающей среды;
К,- - г'-й клапан; Пр - побудитель расхода; Р - ресивер
При математическом описании процесса обогащения воздуха кислородом в установке КБА принимали следующие допущения: 1) атмосферный воздух является трехкомпонентной газовоздушной смесью (содержит 1 - О2 с концентрацией 20,3 - 20,8 об. %; 2 - N с концентрацией 78,2 об. %; 3 - Аг и другие примеси с концентрацией 1,0 - 1,5 об. %) и рассматривается как идеальный газ, что вполне допустимо при давлении в адсорбере до 200-105 Па [3]; причем концентрации аргона и других примесей в газовоздушной смеси не изменяются при обогащении воздуха кислородом; 2) диаметр и высота насыпного слоя цеолитового адсорбента МаХ составляли 0,05 и 0,5 м соответственно; диаметр гранул -2 мм, коэффициент пористости гранул ~ 0,394 и диаметр транспортных пор ~ 0,5^10 м; 3) геометрические размеры адсорбционного слоя считали постоянными в течение заданного срока эксплуатации ~ 105 ч [4, 5]; 4) слой адсорбента представляет собой сплошную среду с коэффициентом порозности е, учитывающим пористость гранул; 5) адсорбционное равновесие (изотерма адсорбции) описывается уравнением Ленгмюра-Фрейндлиха; 6) десорбционные ветви изотерм адсорбции 02, N на цеолите №Х совпадают с адсорбционными [6, 7]; 7) температура воздуха в ресивере равна температуре воздуха на выходе из адсорбера, тепловые потери в окружающую среду пренебрежимо малы.
В соответствии с принятыми допущениями математическое описание циклического процесса адсорбционного разделения воздуха и концентрирования кислорода в установке КБА включает следующие уравнения [8].
1. Уравнения покомпонентного (для О2 - к = 1; N2 - к = 2) материального баланса в потоке газовоздушной смеси по высоте слоя адсорбента
дек (x, t) + (1 - s) дак + d(vgCk (х, t)) = q (1)
dt s dt дх
В уравнении (1) первое слагаемое описывает скорость накопления k-го компонента смеси в газовой фазе; второе - скорость накопления k-го компонента в адсорбенте; третье - конвективный перенос вещества в слое адсорбента.
2. Уравнение кинетики адсорбции:
^Г = вк (ак - ак ), к = 1,2. (2)
dt
Уравнение (2) представляет собой описание кинетики адсорбции для смешанно-диффузионной области массопереноса адсорбтива (O2, N2) через границу раздела фаз.
Кинетические коэффициенты вк определялись путем решения регуляризо-
ванной экстремальной задачи идентификации с использованием экспериментальных данных, полученных методом физического моделирования установки
КБА в лабораторных и опытно-промышленных условиях: Ро2 = 5,776 с-1; р n2 = 1,925 с-1.
Величина адсорбции ак, равновесная текущей концентрации адсорбтива ск в потоке на внешней поверхности гранул, рассчитывается по уравнению изотермы Ленгмюра-Фрейндлиха [9, 10]
к b1,кек exP (ь2,к/Та ) b5,кек exP (b6,к /Ta ) ,
ак =-2-+-2-, к =1,2.
1 + Z b3, jej exP(4,к /Ta ) 1 + X b3, jej exP(b4,к/Ta ) j=1 j=1
3. Уравнение, описывающее распространение тепла в потоке газовоздушной смеси по длине х адсорбента:
dTg(x,t) dTg(x,t) a г i
e$ P^-^— + ср Pg v^-Jdx---^dPk(x,t) -Tg(x,t)] = 0, 0 < x < L , (3)
где Sud = (1 - s)
"gr
В уравнении (3) первое слагаемое описывает накопление тепла в газовой фазе; второе - конвективную составляющую переноса тепла; третье - теплоотдачу от газовой фазы к твердой (адсорбенту).
4. Уравнение, описывающее изменение температуры в адсорбенте
с>а + [(х,,)-(х,,)]-X^ ^ = Ха . (4)
01 к 01 дх
В уравнении (4) первое слагаемое описывает энтальпию твердой фазы (адсорбента); второе - теплоотдачу от твердой фазы (адсорбента) к газовой; третье -выделение теплоты сорбции компонентов газовой смеси; четвертое - теплопроводность в адсорбенте по координате х .
5. Уравнение Эргуна, связывающее изменение давления и скорости газовой смеси по длине адсорбента [11]:
д P дх
(150 (1 - е 0 )2 , (1 - е о ) 2 А v 0/ рgVg +1,75MgPg^-v 2
(d gr tf
b2
d gr 0
(5)
где уравнение состояния идеального газа имеет вид: Р(х, г) = ЯТё (х, г ак (х, г).
Начальные и граничные условия для уравнений (1) - (5), формулы для расчета коэффициентов математической модели приведены в работе [8], поэтому в данной статье не приводятся.
Объемный расход потока, направляемого в ресивер Р (см. рис. 1), определялся как Оои1 = Оа°Ц - О^ , где О^ = хООи Ре1 - объемный расход потока, обогащенного кислородом, направляемого на десорбцию преимущественно 02, N2 и отвод в атмосферу.
Таким образом, уравнения (1) - (5) с начальными и граничными условиями представляют собой математическое описание циклических процессов адсорбции -десорбции, осуществляемых в адсорберах Аь А2 установки КБА (см. рис. 1). Для решения системы дифференциальных уравнений в частных производных с соответствующими начальными и граничными условиями использовали метод конечных элементов в программной среде МЛТЬЛБ.
Далее концентрации кислорода и азота будем обозначать через У = ((1 У2) об. %. Проверка адекватности математической модели проводилась по результатам независимых экспериментов, то есть отличных от экспериментальных данных, использованных для идентификации модели. Рассогласование расчетных по модели и экспериментальных данных оценивали по формуле
out,e
5max = max i=1,d
Уъ
out ■ y1 V'ads.
(fads,i)
,out,e
Уъ
где уои1 (гай8>г-) - экспериментальное и расчетное значения концентрации ки-
слорода в продукционном потоке на выходе установки КБА в /-м опыте, об. %, 1 = 1, й, й - число экспериментальных
У\ .об.%
40
20 4-
д А
L □
/ 1 J
I ■. ■. д J
точек; гаЙ8,/ - продолжительность стадии адсорбции в /-м опыте, с.
Максимальная относительная погрешность 8шах рассогласования расчетных у°и1 и экспериментальных
У°и1,е данных на стадии адсорбции не
превышает 13,2 % (рис. 2), что позволяет с приемлемой для практики точностью использовать математическую Рис. 2. Экспериментальные и расчетные модель для технологического расчета, по модели зависимости концентрации оптимизации и управления процессом пр°дукционног° кисл°р°да при давлении обогащения воздуха кислородом мето- на выходе компрессора: дом адсорбционного разделения газо- 1 - 2,710 Па ; 2 - Ъ,110 Па ; гщ = 0,5гц; воздушной смеси. А, □ - эксперимент; — - расчет по модели
k
Постановка задачи оптимизации режимных параметров процесса обогащения воздуха кислородом
Сформулируем задачу оптимизации циклических режимов работы установки КБА при разделении атмосферного воздуха и концентрировании кислорода: при фиксированных значениях возмущающих воздействий и нагрузки по сырью
, in - in ~7n /">in ^ /^in , Т-"1" г *
yi < yi < yi , G < G < G требуется определить управления u - продолжительность цикла /Ц «адсорбция - десорбция», степени открытия
!* IJ* !* U* Г/Л Г\ Г+ 1 !* !J* !* !J*
¥1 = , ¥4 = VK4 на отрезке времени [0;0,5/ц] и у2 = ¥K2, ¥3 = ¥K3 на отрезке времени [0,5/ц; /ц ] впускных K1, K2 и сбросных клапанов K3, K4 такие, что среднее значение степени извлечения (концентрирования) кислорода п на отрезке времени [0, /уст] достигает максимума, то есть
1 ^уст(и*) 1 ^уст(и)
I (u*, Z) =- f п(и*, z' )dt = max - f n(u,z' )dt (6)
t ^ J,i J,i J,i J,i ■ 1 \ t ^
уст 0 u=(tц, ¥1 , ¥3 , ¥2 ' ¥4 , J=1'™) уст 0
при связях (l) - (5) с начальными и граничными условиями и ограничениях:
- на «установившееся» значение концентрации продукционного кислорода
у°"ад -Уоуст(u,z)<0,i = 1,2,3, ... ; (7)
- производительность установки
бзад-6(u,zi, yfu^) < 0; (8)
- скорость газовоздушной смеси в «лобовом» слое адсорбента
ma
fe[ö,
ax I vg lu, z
(,,,<)
- v+< 0; (9)
перепад давления в слое адсорбента
max Pa°dSt(u,zi, L, t) - pdSO^ ,0, t) < Д/^; /е[0;0,5/ц ]
max PdSW, L, t) - Ou^' ,0, t) <APd+eS; (10)
/Е[0,5/ц;/Ц]
расход атмосферного воздуха в питании установки КБА
max Gin(t) < Gin+; (11)
/Е[0;0,5/ц ]
- управления
tц <tn <tц; (12)
0 < yj < 1, X = 1,4; j = 1m; (13)
Pands < PdS < P'ands, (14)
где у°зад - заданное значение концентрации кислорода в продукционном потоке воздуха; Qзад - заданное значение производительности (расхода продукционного
потока воздуха при концентрации У^д кислорода) установки КБА; - максимально допустимое значение скорости газовоздушной смеси в слое адсорбента; ДРа+ъ, ДРаез - максимально допустимые значения перепада давления в адсорбционном слое на стадиях адсорбции и десорбции; гц, гц минимальное и максималь-
ное значения продолжительности цикла «адсорбция - десорбция»; Р^з, Раеэ -
минимальное и максимальное допустимые давления на стадии адсорбции соответственно.
Критерий оптимальности (6) и ограничения (1) - (10) вычисляются через решение уравнений динамики (1) - (5) до наступления нового стационарного периодического режима работы установки КБА. Ограничение максимального значения скорости потока воздуха в «лобовом» слое адсорбента (9) необходимо учитывать при решении задачи оптимизации для предохранения гранулированного адсорбента от разрушения вследствие аэродинамического удара при резком (скачкообразном) открытии регулирующих клапанов К - К4.
Сформулированная аргументная задача оптимизации (6) - (14) относится к классу задач нелинейного программирования, для решения которой будем использовать метод последовательного квадратичного программирования [12].
Разработка алгоритма и системы управления циклическим адсорбционным процессом обогащения воздуха кислородом
Для управления циклическим адсорбционным процессом концентрирования кислорода в установке КБА разработана двухуровневая система управления, в которой на нижнем уровне находятся: контроллер, измерительные преобразователи, исполнительные механизмы; на верхнем уровне: 8СЛЭЛ-система на персональном компьютере.
При разработке алгоритма управления короткоцикловым процессом концентрирования кислорода будем считать, что возмущающие воздействия Т являются измеримыми в моменты времени г/ и действуют в течение интервала времени [гг-, г/+1 ], ' = 0,1,2... Алгоритмом управления предусматривается непрерывный
контроль текущих значений возмущений Т = {уОс , ТО с , Б'0 с , Рои1,/} и при их отклонении от номинальных значений осуществляется решение задачи оптимизации (6) - (14) методом последовательного квадратичного программирования. Результатом решения задачи оптимизации являются: оптимальное значение продолжительности цикла гц и законы оптимального изменения степени открытия клапанов
К!, К4 - = ^К'*, ^4* = ^К'4*, 1 = 1,т на отрезке времени [0;0,5гц] и К2, К3 -
^2* = ^К2*, < = < на отрезке времени [0,5гц; гц]. Далее на основании рассчитанного оптимального управления производится пересчет циклограммы переключений клапанов К1 - К4. Программно-логическое управление клапанами реализуется посредством программного задатчика.
Пусть при tj = 1600 с в окружающей среде скачкообразно изменяется температура 7Ос от номинального значения 293 до 323 K, а также состав воздуха: концентрация кислорода >}" скачкообразно уменьшается от 20,8 до 20,3 об. %, концентрация примесей >3" увеличивается от 1,0 до 1,5 об. %. Для этих условий оптимальными управлениями u* являются tk = 31 с, у^* - у4*, j = 1,10 (рис. 3).
Из анализа графиков на рис. 4 следует, что ограничение по ресурсосбережению адсорбента - скорости газовоздушной смеси в «лобовом слое» ( vg < 0,08 м/с
при х = 0) выполняется.
Анализ других результатов показывает, что при адаптивной оптимизации режимов функционирования установки КБА степень извлечения кислорода увеличивается в среднем на 9,6 % и выполняются ограничения по чистоте кислорода
>out > 39 об. % и производительности установки Q > 2 нл/мин (табл. 1).
100 80 60 40 20
w , V2 ,%
1600
1650
1700
1750
1800 t, С
а)
100 80 60 40 20
V4 ,Vз ,%
Кз
1800 t, С
б)
j* j* j* j*
Рис. 3. Динамика степени у[ , у^ , ¥4 , ¥3 открытия впускных Kj, K2 (а)
и сбросных K3, K4 (б) клапанов
vg, м/с
_о,о1 1600 1610 1620 1630 1640 1650 г, с
Рис. 4. Динамика скорости воздуха в слое адсорбента при скачкообразном возмущении по температуре окружающей среды (293 ... 323 К):
1 - х = 0; 2 - х = Ь/2; 3 - х = Ь
0
0
1650
1750
Таблица 1
Сравнительный анализ работы установки с адаптивной оптимизацией (+) в системе управления и без адаптивной оптимизации (-)
Показатели установки КБА Оптимальные значения для номинального режима Случайные значения возмущений
Тос = 233 K Тос = 323 K Тос = 323 K y'{ = 20,3 %
- + - + - +
П* , % 41,5 30,9 40,3 33,1 42,2 31,7 41,9
, об. % 39,0 58,7 39,0 25,5 39,0 23,2 39,0
Gout*, нл/мин 2,16 1,71 2,12 4,35 2,20 4,41 2,18
Заключение
С использованием методов теории автоматического управления, системного анализа и математического моделирования получены новые научные результаты для проектирования автоматизированных процессов короткоцикловой адсорбции и систем управления ими с учетом выполнения требований технологического регламента по степени извлечения и чистоте кислорода, производительности установки КБА и сбережению гранулированного цеолитового адсорбента от разрушения.
Работа выполнена в рамках проектной части государственного задания № 10.3533.2017/ПЧ.
Список литературы
1. Шумяцкий, Ю. И. Промышленные адсорбционные процессы: учебное пособие для студентов химико-технологических специальностей вузов / Ю. И. Шумяцкий. - М. : КолосС, 2009. - 184 с.
2. Ruthven, D. M. Pressure Swing Adsorption / D. M. Ruthven, S. Farooq, K. S. Knaebel. - University of New Brunswick Frederiction, NB, Canada. - VCH Publishers, Inc, 1994. - 352 p.
3. Кириллин, В. А. Техническая термодинамика / В. А. Кириллин, В. В. Сычев, А. Е. Шейндлин. - 5-е изд., перераб. и доп. - М. : Издательский дом МЭИ, 2008. - 496 с.
4. Акулов, А. К. Моделирование разделения бинарных газовых смесей методом адсорбции с колеблющимся давлением : дис. ... д-ра техн. наук : 05.17.08 : защищена 22.02.96 : утв. 08.11.96 / Акулов Аркадий Клавдиевич. - СПб., 1996. - 304 с.
5. Акулинин, Е. И. Исследование процессов тепло- и массообмена при обогащении воздуха кислородом методом короткоцикловой адсорбции / Е. И. Акулинин, Д. С. Дворецкий, С. И. Дворецкий // Вестн. Тамб. гос. техн. ун-та. -2016. - Т. 22, № 3. - С. 411 - 419.
6. Кельцев, Н. В. Основы адсорбционной техники / Н. В. Кельцев. - 2-е изд., перераб. и доп. - М. : Химия, 1984. - 592 с.
7. Modelling the Process of Air Oxygenation by Nitrogen Adsorption in a Pressure Swing Adsorption Unit / E. I. Akulinin [et al.] // Вестн. Тамб. гос. техн. ун-та. - 2012. -Т. 18, № 1. - С. 182 - 196.
8. Васильев, А. С. Оптимизация и управление циклическим процессом адсорбционного концентрирования кислорода : дис. ... канд. техн. наук : 05.13.06 / Васильев Александр Сергеевич. - Тамбов, 2018. - 160 с.
9. Farooq, S. Numerical Simulation of a Pressure Swing Adsorption Oxygen Unit / S. Farooq, D. M. Ruthven, H. A. Boniface // Chemical Engineering Science. - 1989. -Vol. 44, No. 12. - P. 2809 - 2816.
10. Versatile Process Simulator for Adsorptive Separations / R. Kumar [et al.] // Chemical Engineering Science. - 1994. - Vol. 49, No. 18. - P. 3115 - 3125.
11. Белоглазов, И. Н. Основы расчета фильтрационных процессов / И. Н. Бе-логлазов, В. О. Голубев. - М. : Издательский дом «Руда и Металлы», 2002. - 210 с.
12. Реклейтис, Г. Оптимизация в технике. В 2-х кн. Кн. 2 / Г. Реклейтис, А. Рейвиндран, К. Рэгсдел : пер. с англ. : В. Я. Алтаева, В. И. Моторина. - М. : Мир, 1986. - 320 с.
Optimization and Control over the Cyclic Process of Adsorption Enrichment of Air with Oxygen
V. G. Matveykin1, S. A. Skvortsov1, E. I. Akulinin2, S. I Dvoretsky2
Departments of Information Processes and Control (1), Technologies and Equipment for Food and Chemical Industries(2), TSTU, Tambov, Russia; ipu@ahp.tstu.ru
Keywords: air; dynamics; adsorption isotherm; kinetics; oxygen; short cycle adsorption; mathematical model; optimization; control; zeolite adsorbent.
Abstract: The paper investigates optimization of the process of air enrichment with oxygen (i.e. determination of the maximum degree of oxygen extraction) in a dualadsorption plant for pressure swing adsorption (PSA), taking into account the requirements for oxygen purity, installation efficiency and resource-saving granular adsorbent. An algorithm and a two-level system of adaptive control of air enrichment with oxygen have been developed; their implementation will provide automatic control of disturbing influences, promptly solving the optimization problem (at the upper level) and correcting tasks for the program generator and regulators of the automatic control system (at the lower level). The use of adaptive optimization in the control system will increase the value of the degree of oxygen extraction in the PSA plant by an average of 10 % (compared to the traditional control system without adaptive optimization) while fulfilling the required limitations on oxygen purity, resource saving of the granular adsorbent and plant capacity.
References
1. Shumyatskiy Yu.I. Promyshlennyye adsorbtsionnyye protsessy [Industrial Adsorption Processes], Moscow: KolosS, 2009, 184 p. (In Russ.)
2. Ruthven D.M., Farooq S., Knaebel K.S. Pressure Swing Adsorption, VCH Publishers, Inc, 1994, 352 p.
3. Kirillin V.A., Sychev V.V., Sheyndlin A. Ye. Tekhnicheskaya termodinamika [Technical Thermodynamics], Moscow: Izdatel'skiy dom MEI, 2008, 496 p. (In Russ.)
4. Akulov A.K. PhD Dissertation (Technical), St. Petersburg, 1996, 304 p. (In Russ.)
5. Akulinin Ye.I., Dvoretskiy D.S., Dvoretskiy S.I. [Study of the Processes of Heat and Mass Transfer During Air Enrichment with Oxygen by the Method of Short-Cycle Adsorption], Transactions of the Tambov State Technical University, 2016, vol. 22, no. 3, pp. 411-419. (In Russ., abstract in Eng.)
6. Kel'tsev N.V. Osnovy adsorbtsionnoy tekhniki [Basics of Adsorption Technology], Moscow: Khimiya, 1984, 592 p. (In Russ.)
7. Akulinin E.I., Dvoretsky D.S., Dvoretsky S.I., Tugolukov E.N. Modelling the Process of Air Oxygenation by Nitrogen Adsorption in a Pressure Swing Adsorption Unit, Transactions of the Tambov State Technical University, 2012, vol. 18, no. 1, pp. 182-196. (In Eng., abstract in Russ.)
8. Vasil'yev A.S. PhD Dissertation (Technical), Tambov, 2018, 160 p. (In Russ.)
9. Farooq S., Ruthven D.M., Boniface H.A. Numerical Simulation of a Pressure Swing Adsorption Oxygen Unit, Chemical Engineering Science, 1989, vol. 44, no. 12, pp. 2809-2816.
10. Kumar R., Fox V.G., Hartzog D.G., Larson R.E., Chen Y.C., Houghton P.A., Naheiri T. Versatile Process Simulator for Adsorptive Separations, Chemical Engineering Science, 1994, vol. 49, no. 18, pp. 3115-3125.
11. Beloglazov I.N., Golubev V.O. Osnovy rascheta fil'tratsionnykh protsessov [Basics of Calculation of Filtration Processes], Moscow: Izdatel'skiy dom «Ruda i Metally», 2002, 210 p. (In Russ.)
12. Rekleytis G., Reyvindran A., Regsdel K. Optimizatsiya v tekhnike [Optimization in Engineering], Moscow: Mir, 1986, 320 p. (In Russ.)
Optimierung und Steuerung des zyklischen Prozesses der Adsorptionsbereicherung der Luft mit Sauerstoff
Zusammenfassung: Gestellt und erforscht ist das Optimierungsproblem (im Sinne der Bestimmung des maximalen Grads der Sauerstoffextraktion) des Luftanreicherungsprozesses mit Sauerstoff in einer zwei-Wege-Anlage der Druckwechseladsorption (CBA) im Hinblick auf die Anforderungen an die Sauerstoffreinheit, Arbeitsleistung der Anlage und die Ressourcenschonung des granulierten Adsorptionsmittels. Es ist der Algorithmus und das zweistufige System der adaptiven Steuerung des Prozesses der Anreicherung der Luft mit Sauerstoff entwickelt. Ihre Realisierung stellt die automatische Steuerung der störenden Effekte, die operative Lösung des Problems der Optimierung (auf der obersten Ebene) und die Korrektur der Aufgaben an den Software-Sensor und die Regler des automatischen Steuerungssystems (auf der unteren Ebene) zur Verfügung. Die Nutzung der adaptiven Optimierung im Steuersystem wird ermöglichen, den Wert für den Grad der Sauerstoffextraktion bei der CBA-Anlage durchschnittlich um 10 % zu erhöhen (im Vergleich mit der Verwendung des herkömmlichen Managementsystems ohne adaptive Optimierung), wenn die erforderlichen Einschränkungen für die Reinheit von Sauerstoff, die Ressourceneinsparung des granulierten Absorbers und die Leistung der Anlage erfüllt sind.
Optimisation et contrôle du processus cyclique de l'adsorption de l'oxygène de l'air
Résumé: Est fixée et étudiée la tâche de l'optimisation (en ce qui concerne la détermination du degré maximum d'extraction de l'oxygène) du processus d'enrichissement de l'air en oxygène dans une installation à deux étages d'adsorption à
cycle court sans chauffage (ACS), en tenant compte des exigences de la pureté de l'oxygène, de la productivité de l'installation et de la conservation des ressources de l'adsorbant granulaire. Sont mis au point un algorithme et un système de contrôle adaptatif à deux niveaux pour l'enrichissement de l'air en oxygène ce qui permettra un contrôle automatique des effets perturbateurs, une solution d'optimisation opérationnelle (au niveau supérieur) et une correction des tâches aux réglages programmatiques et aux régulateurs du système de contrôle automatique (au niveau inférieur). Grâce à l'optimisation adaptative du système de contrôle, la valeur de l'extraction de l'oxygène dans l'installation ACS est augmentée en moyenne de 10 % (par rapport au système de contrôle traditionnel sans optimisation adaptative), tout en respectant les limites de pureté requises, en économisant les ressources l'adsorbant granulaire et la productivité de l'installation.
Авторы: Матвейкин Валерий Григорьевич - доктор технических наук, профессор, заведующий кафедрой «Информационные процессы и управление»; Скворцов Сергей Александрович - кандидат технических наук, доцент кафедры «Информационные процессы и управление»; Акулинин Евгений Игоревич - кандидат технических наук, доцент кафедры «Технологии и оборудование пищевых и химических производств»; Дворецкий Станислав Иванович - доктор технических наук, профессор кафедры «Технологии и оборудование пищевых и химических производств», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.
Рецензент: Гатапова Наталья Цибиковна - доктор технических наук, профессор, заведующий кафедрой «Технологические процессы, аппараты и тех-носферная безопасность», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.
