Автоматика. Информатика. Управление. Приборы
УДК 661.9
DOI: 10.17277^^.2017.04^.548-556
АВТОМАТИЗАЦИЯ ПРОЦЕССА АДСОРБЦИОННОГО РАЗДЕЛЕНИЯ ГАЗОВЫХ СМЕСЕЙ И ПОЛУЧЕНИЯ ВОДОРОДА
В. Г. Матвейкин1, А. А. Ишин2, С. А. Скворцов1, С. И. Дворецкий3
Кафедры: «Информационные процессы и управление» (1),
«Технологии и оборудование пищевых и химических производств» (3), ФГБОУВО «ТГТУ» г. Тамбов, Россия; sdvoretsky@tstu.ru;
ООО «Энерготехпроект» (2), г. Тамбов, Россия
Ключевые слова: водород; динамика; изотерма адсорбции; кинетика; короткоцикловая безнагревная адсорбция; математическая модель; многокомпонентная газовая смесь; оптимизация; управление; цеолитовый адсорбент.
Аннотация: Сформулирована задача оптимизации режимных переменных (давления на стадии адсорбции и времени цикла адсорбции) и разработана двухуровневая система адаптивной оптимизации процесса адсорбционного разделения газовых смесей и получения водорода максимальной чистоты в четырех-адсорберной установке короткоцикловой безнагревной адсорбции с гранулированным цеолитовым адсорбентом СаА. Проведено численное исследование эффективности функционирования системы адаптивного управления при ступенчатом изменении возмущающих воздействий (концентрации диоксида углерода, температуры и давления газовой смеси) в питании установки короткоцикловой безнагревной адсорбции.
Введение
Циклические адсорбционные процессы широко применяют в промышленности для разделения углеводородов, концентрирования водорода, кислорода, диоксида углерода и др. газов, безнагревной осушки газов, извлечения метана и т.п. Одной из актуальных задач в области адсорбционного разделения является получение водорода из водородсодержащих технологических потоков (газы конверсии и окисления углеводородов, нефтезаводские газы, синтез-газ и др.) [1]. Типичными веществами, сопутствующими водороду, являются азот, оксид и диоксид углерода, метан. Целью данной работы является анализ циклического процесса адсорбционного обогащения газовой смеси водородом как объекта управления и разработка системы автоматического управления данным процессом, обеспечивающей поддержание максимальной концентрации продукционного водорода при выполнении регламентных технологических условий (ограничений).
Анализ процесса получения водорода методом адсорбционного разделения газовых смесей как объекта управления
Анализ процесса короткоцикловой безнагревной адсорбции (КБА) обогащения газовой смеси водородом как объекта исследования позволил определить [2, 3]:
- входные переменные х (возмущающие воздействия) установки КБА (состав ут = (у1", у2" , Уз"), температуру Т"", расход О1" и давление Р1П газовой смеси в
питании установки, давления на продукционном P^ и на сбросном P£o2 выхо-
,,in rpin s^in r>in nQUt nQUt
H2 >PC02 } (рис. !);
дах установки КБА), то есть x = {yin, Tgn, Gin, Pin, P^
- режимные (управляющие) переменные u (давление Pads и время цикла адсорбции тads), то есть u = {Pads, тads} (рис. 1);
- выходные переменные y установки КБА (состав yoUt = (yfUt, yOUt, y^Ut), расходы GOUt, Gp и температуру T°Ut газовой смеси на выходе из установки
КБА, производительность установки Q = yOUt ■ Gp, где Gp = GOUt - Gdes, Gdes -объемный расход потока, обогащенного водородом, направляемого на регенерацию адсорбента), то есть y = {yOUt, GOUt, Gp, TgOUt, q} (рис. 1).
При математическом описании процесса обогащения газовой смеси водородом в установке КБА принимали следующие допущения [2, 4]: 1) исходная газовая смесь является трехкомпонентной (содержит H2 с концентрацией (65 ± 2) % об. (в формулах используется индекс «1»), СО2 с концентрацией (34 ± 2) % об. (индекс «2»), CO с концентрацией (1 ± 0,5) % об. (индекс «3»)) и рассматривается как идеальный газ, что вполне допустимо при давлении в адсорбере до 200 х 105 Па; 2) в качестве адсорбента используется гранулированный цеолит СаА с диаметром гранул 1,5 мм, коэффициентом пористости гранул ~ 0,394 и диаметром транспортных пор d < 0,5 х 10-3 м; 3) геометрические размеры адсорбционного слоя считаем постоянными в течение заданного срока эксплуатации ~105 ч; 4) слой адсорбента представляет собой сплошную среду с коэффициентом порозности е; 5) десорбционная ветвь изотерм сорбции H2, CO2 и CO на гранулированном цеолите СаА совпадает с адсорбционной и носит нелинейный характер. В соответст-
Рис. 1. Схема четырехадсорберной установки КБА как объекта управления при обогащении газовой смеси водородом
вии с принятыми допущениями математическое описание процесса концентрирования водорода в потоке газовой смеси включает уравнения покомпонентного материального и теплового балансов в потоке газовой смеси по высоте слоя адсорбента (параболические дифференциальные уравнения второго порядка); уравнение массопереноса адсорбтива (Н2, С02, СО) из газовой фазы в твердую фазу адсорбента (через границу раздела фаз); уравнение изотермы сорбции Ленгмюра -Фрейндлиха для расчета равновесной концентрации адсорбтива на внешней поверхности гранул адсорбента; уравнение Эргуна, связывающее изменение давления и скорости газовой смеси по высоте адсорбента; начальные и граничные условия для системы параболических дифференциальных уравнений второго порядка [2].
Постановка задачи оптимизации установившегося технологического режима
Задача оптимизации технологического процесса получения водорода по технологии КБА может быть сформулирована следующим образом: при заданных диапазонах изменения возмущающих воздействий
x = У
in (jin pin pout j}OMt
H,
CO,
требуется найти вектор оптимальных управляющих воздействий и = (Ра^, та^), при котором целевая функция (среднее значение концентрации упродукционного водорода на заданном периоде [0, р] функционирования установки КБА) достигает экстремального значения, то есть
I(u*)=
( 1 tpr ) — i yOut (u*)dt tpr i
pr 0
= max
u={Pаds , Tads j
1 tPr
— f yOut (u) tpr i
л
pr 0
(l)
выполняются связи в форме уравнений математической модели [2, 3] и ограничения:
- по концентрации продукционного водорода y
y°ut - y°ut < 0,
- производительности Qh2 установки КБА
Q < Q < Q,
- расходу газовой смеси Gm в питании установки КБА
Gin < Gin < Gin ,
- давлению на стадии адсорбции Pads
Р1П < Pads ,
(2)
(3)
(4)
(5)
где черта внизу (вверху) - нижние (верхние) предельно допустимые значения технологических переменных.
Сформулированная задача (1) - (5) относится к классу задач нелинейного программирования, для ее решения будем использовать методы штрафных функций и последовательного квадратичного программирования [5].
Автоматизация процесса получения водорода методом адсорбционного разделения газовых смесей
Для оперативного решения задачи оптимизации (1) - (5) и управления процессом получения водорода разработана двухуровневая система адаптивной оптимизации и управления.
Технологические переменные установки адсорбционного получения водорода подвержены случайным изменениям в ходе процесса адсорбции. Значения возмущений, которые представлены нерегулируемыми переменными исходной газовой смеси в питании установки КБА, в ходе процесса также меняются случайным образом. В этом случае необходимо поддерживать априори неизвестное максимальное значение заданной целевой функции оптимизации, что обусловливает применение системы адаптивной оптимизации и управления с изменяемыми заданиями автоматическим регуляторам, работающей по принципу управления по возмущению с эталонной моделью технологического процесса в контуре управления [6, 7].
В системе адаптивной оптимизации и управления осуществляется непрерывный контроль за текущими значениями возмущающих воздействий и при их отклонении от номинальных значений на верхнем уровне с помощью персонального компьютера оперативно решается задача оптимизации (1) - (5) и определяются текущие оптимальные задания й = (Рад, т а^) регуляторам автоматической системы регулирования технологического процесса получения водорода, функционирующей на нижнем уровне.
На основании полученного значения тпересчитывается циклограмма и работы клапанов и с помощью программного задатчика ПЗ2 и логического элемента ЛЭ осуществляется ее реализация в установке КБА. Текущее оптимальное значение Рас1я поступает в качестве задания ПИД-регулятору одноконтурной системы автоматического регулирования с обратной связью (рис. 2).
Таким образом, система адаптивной оптимизации и управления режимами функционирования установки КБА обеспечивает выполнение следующих функций:
- поиск и поддержание оптимального значения чистоты получаемого водорода у0*;
Рис. 2. Структурная схема системы управления:
у - степень открытия клапана; Д - датчик; Р - регулятор; ЛЭ - логический элемент; ПЗ - программный задатчик
- расчет текущего оптимального времени цикла адсорбции таЛ;
- расчет и реализация оптимальной циклограммы работы клапанов в установке КБА для получения водорода;
- расчет текущего оптимального давления на стадии адсорбции Рас^8;
- расчет и формирование управляющих воздействий на приводы клапанов.
Имитационное исследование системы управления
На рисунке 3 представлены результаты имитационных исследований функционирования системы управления при ступенчатом повышении температуры
исходной смеси от 30 до 50 °С.
Ор,нл/мин 800
1920 1940 1960 1980 2000 2020 2040 t, С
а)
У°", % Об. 100
Т-Т-Ч—!—г—I—f—I—1—I——I—Г—I—!
> I
I I 4 4 fHr f rrf
99 98 97 96 95 94 93
1920 2120 2320 2520 2720 6)
P , M Па
2
1.5
0.5
fw- 1 ч Pads
iV 11
t, с
Па
1,59
1,575
1,56
1,545
1,53
1,515
1,5
0
1920 1940 I960 1980 2000 2020 2040 Г,с в)
Рис. 3. Переходные процессы в системе управления при ступенчатом повышении температуры исходной смеси:
а - изменение расхода G продукционного водорода и оптимального времени цикла адсорбции т ads в зависимости от времени t; б - изменение концентрации yjut продукционного водорода в зависимости от времени при наличии 1 и отсутствии 2 адаптивной оптимизации; в - изменение давлений в адсорберах Aj - А4 установки
Анализ показывает, что система управления компенсирует ступенчатое возмущение уменьшением длительности та^ стадии адсорбции с 42 до 36 с (рис. 3, а). При этом система управления обеспечивает поддержание концентрации продукционного водорода на уровне ~ 99 % об. и расхода Ор продукционного водорода на минимально допустимом уровне (рис. 3, а), в то время как в установке КБА без адаптивной оптимизации концентрация продукционного водорода устанавливается на уровне ~ 94 % об. (рис. 3, б).
На рисунке 4 представлены результаты имитационных исследований функционирования системы управления при ступенчатом повышении содержания СО2 в исходной смеси у2" от 34 до 45 % об. В данном случае система управления
900 '
800 700 600 500 400
300
G р
А
1 п
^ads
7
43 41
37
1920 1940 1960 1980 2000 2020 Г, С
а)
уГ, % Об.
1920
2120 2320 2520 2720 б)
РМ1,МПа 1,59
1920 1940 1960 1980 2000 2020 /,С б)
Рис. 4. Переходные процессы в системе управления при ступенчатом повышении содержания С02 в исходной смеси:
а - изменение расхода ОР продукционного водорода и оптимального времени цикла адсорбции Т аа^ в зависимости от времени; б - изменение концентрации у'°и: продукционного водорода в зависимости от времени при наличии 1 и отсутствии 2 адаптивной оптимизации; в - изменение давлений в адсорберах А^А4 установки
Таблица 1
Результаты сравнительного анализа работы установки КБА
Возмущающее воздействие -ступенчатое увеличение * I , % об. I, % об. АI, % об.
Tgn от 30 до 50 °С 98,92 94,13 4,79
y2 от 34 до 45 % об. 95,76 89,46 6,3
PHut до 0,3 МПа H2 98,92 95,09 3,83
P£o до 0,1 МПа 96,51 90,22 6,29
Примечание: I , I - значения целевой функции при наличии и без адаптивной оптимизации в системе управления соответственно; А I - изменение целевой функции.
минимизирует снижение чистоты продукционного водорода, обеспечивая поддержание концентрации продукта на максимально возможном уровне для изменившихся характеристик исходной смеси. Система управления с адаптивной оптимизацией более эффективна по сравнению с системой без адаптивной оптимизации (рис. 4, б).
Результаты сравнительного анализа работы установки КБА при наличии и отсутствии адаптивной оптимизации в системе управления представлены в табл. 1.
Заключение
Математическое, информационное и алгоритмическое обеспечение двухуровневой системы адаптивной оптимизации и управления статическими режимами функционирования установки КБА могут быть использованы при проектировании новых автоматизированных процессов и адсорбционных технологических установок с циклически изменяющимся давлением для разделения и очистки многокомпонентных газовых смесей.
Работа выполнена в рамках проектной части государственного задания № 10.3533.2017/ПЧ.
Список литературы
1. Ruthven, D. M. Pressure Swing Adsorption / D. M. Ruthven, S. Farooq, K. S. Knaebel. - NewYork : VCH Publishers, 1994. - 352 p.
2. Mathematical Modeling of Hydrogen Production Process by Pressure Swing Adsorption Method / E.I. Akulinin [et al.] // Advanced Materials and Technologies. -2017. - No. 2. - P. 38 - 49. doi: 10.17277/amt.2017.02.pp.038-049
3. Ишин, А. А. Математическое моделирование и управление процессом получения водорода методом адсорбционного разделения газовой смеси : дис. ... канд. техн. наук : 05.13.06 / Ишин Андрей Анатольевич. - Тамбов, 2017. - 152 c.
4. Akulinin, E. Modeling and Automation of Adsorption Processes with Cyclically Changing Pressure During the Separation of Gas Mixture / E. Akulinin, A. Ishin, S. Skvortsov // 17 International Multidisciplinary Scientific Geoconference SGEM 2017 : Conference Proceedings. - 2017. - Vol. 17 : Energy and Clean Technologies, Issue 41. - P. 527 - 534.
5. Реклейтис, Г. Оптимизация в технике : в 2-х кн. : пер. в англ. / Г. Реклейтис, А. Рейвиндран, К. Рэгсдел. - М. : Мир, 1986. - 2 кн.
6. Симонов, М. С. Исследование динамических режимов пневматического концентратора кислорода / М. С. Симонов, С. А. Скворцов, М. М. Шлыков // Вестн. Тамб. гос. техн. ун-та. - 2012. - Т. 18, вып. 4. - С. 987 - 990.
7. Васильев, А. С. Исследование режимов концентраторов кислорода с учетом их конструктивных особенностей / А. С. Васильев, Т. А. Лапина, С. А. Скворцов // Вестн. Тамб. гос. техн. ун-та. - 2012. - Т. 18, вып. 4. - С. 1007 - 1011.
Automation of the Process of Gas Mixtures Adsorptive Separation and Hydrogen Production
V. G. Matveykin1, A. A. Ishin2, S. A. Skvortsov1, S. I Dvoretsky3
Department of Information Processes and Control (1), Department of Technology and Equipment for Food and Chemical Production (3), TSTU, Tambov, Russia; sdvoretsky@tstu.ru;
LLC "Energotehproekt" (2), Tambov, Russia
Keywords: adsorption isotherm; control; dynamics; hydrogen; kinetics; mathematical model; multicomponent gas mixture; optimization; short-cycle heatless adsorption; zeolite adsorbent.
Abstract: The problem of optimizing the regime variables (adsorption stage pressure and the adsorption cycle time) is formulated and the 2-level system for adaptive optimization of the process of gas mixtures adsorptive separation and hydrogen production with maximum purity in the 4-bed short-cycle heatless adsorption unit with a granular zeolite adsorbent CaA is developed. A numerical study of the efficiency of the adaptive control system with a step change of disturbances (carbon dioxide concentration, temperature and pressure of the gas mixture) on input of the short-cycle heatless adsorption unit was carried out.
References
1. Ruthven D.M., Farooq S., Knaebel K.S. Pressure swing adsorption, New York: VCH Publishers, 1994, 352 p.
2. Akulinin E.I., Ishin A.A., Skvortsov S.A., Dvoretsky D.S., Dvoretsky S.I. Mathematical Modeling of Hydrogen Production Process by Pressure Swing Adsorption Method, Advanced Materials and Technologies, 2017, no. 2, pp. 38-49, doi: 10.17277/ amt.2017.02.pp.038-049
3. Ishin A.A. Matematicheskoye modelirovaniye i upravleniye protcessom polucheniya vodoroda metodom adsorbtcionnogo razdeleniya gazovoy smesi [Mathematical modeling and control of the hydrogen production process by the adsorptive separation method for gas mixtures], PhD Dissertation (Engineering), Tambov, 2017, 152 p. (In Russ.)
4. Akulinin, E., Ishin, A., Skvortsov, S. Modeling and Automation of Adsorption Processes with Cyclically Changing Pressure During the Separation of Gas Mixture, 17th International Multiisciplinary Scientific Geoconference SGEM 2017, Conference Proceedings, 2017, vol. 17, Energy and Clean Technologies, issue 41, pp. 527-534.
5. Reklaitis G.V., Ravindran A., Ragsdell K.M. Engineering Optimization: Methods and Applications, New York: Wiley, 1983, 684 p.
6. Simonov M.S., Skvortsov S.A., Shlykov M.M. Issledovaniye dynamicheskih rejimov pnevmaticheskogo koncentratora kisloroda [Study of dynamic modes of the pneumatic oxygen concentrator], Transactions of the Tambov State Technical University, 2012, vol. 18, issue 4, p. 987-990. (In Russ., abstract in Eng.)
7. Vasilyev A.S., Lapina T.A., Skvortsov S.A. [The study of the modes of oxygen concentrators and their structural features], Transactions of the Tambov State Technical University, 2012, vol. 18, issue 4, p. 1007-1011. (In Russ., abstract in Eng.)
Automatisierung des Prozesses der Adsorptionstrennung von Gasmischungen und der Wasserstofferzeugung
Zusammenfassung: Es ist das Problem der Optimierung von Betriebsvariablen (Druck in der Adsorptionsstufe und Zeit des Adsorptionszyklus) formuliert und das 2-stufige System zur adaptiven Optimierung des Verfahrens der Adsorptionstrennung von Gasgemischen und Erzeugung von Wasserstoff höchster Reinheit in der 4-Adsorberanlage für kurzzyklische, nicht beheizte Adsorption mit körnigem Zeolith-Adsorptionsmittel CaA ist entwickelt.
Eine numerische Untersuchung der Effektivität der Funktionsweise des adaptiven Regelsystems mit einer schrittweisen Änderung der Störeffekte (Konzentration von Kohlendioxid, Temperatur und Druck des Gasgemisches) bei der Versorgung einer kurzzyklischen nicht erwärmten Adsorptionsanlage ist durchgeführt.
Automatisation du processus de la séparation d'adsorption des mélanges gazeux et de la production de l'hydrogène
Résumé: Est formulée la tâche de l'optimisation des variables de régime (pression de la phase d'adsorption et temps du cycle de l'adsorption); est conçu un système à 2 niveaux de l'optimisation d'adaptation du processus de la séparation d'adsorption des mélanges gazeux et de l'obtention de l'hydrogène de la pureté maximale dans une installation de 4 adsorptions à cycle court sans chauffage avec un adsorbant granulé zéolitique CaA. Est effectuée une étude numérique de l'efficacité du fonctionnement du système de la commande adaptative lors de la modification étagée des effets perturbateurs (concentration en dioxyde de carbone, température et pression du mélange gazeux) dans l'alimentation de l'installation d'adsorption à cycle court sans chauffage.
Авторы: Матвейкин Валерий Григорьевич - доктор технических наук, профессор, заведующий кафедрой «Информационные процессы и управление», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия; Ишин Андрей Анатольевич - кандидат технических наук, инженер-проектировщик, ООО «Энерготехпроект», г. Тамбов, Россия; Скворцов Сергей Александрович - кандидат технических наук, доцент кафедры «Информационные процессы и управление»; Дворецкий Станислав Иванович - доктор технических наук, профессор кафедры «Технологии и оборудование пищевых и химических производств», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.
Рецензент: Муромцев Дмитрий Юрьевич - доктор технических наук, профессор, проректор по научно-инновационной деятельности, ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.