Definition parameters of mathematical model of the distillation process ethanolamine

Full Text

Этаноламины являются базовым сырьем и входят в состав рецептур готовых продуктов во многих отраслях народного хозяйства: от нефте- и газоперерабатывающей, цементной и металлургической до косметической, фармацевтической и сельскохозяйственной. К качеству этаноламинов предъявляются высокие требования. В промышленности этаноламины получают в результате синтеза оксида этилена и аммиака. Заключительной стадией производства является разделение трехкомпонентной смеси моноэтаноламина (МЭА), диэтаноламина (ДЭА) и триэтаноламина (ТЭА) путем ректификации в двух колоннах насадочного типа. Управление работой последней из этих колонн по разделению ДЭА и ТЭА вследствие того, что оба ее выходных потока являются целевыми, требует регулирования состава как дистиллята, так и кубовой жидкости. Применение двух прямых или косвенных регуляторов состава дистиллята и кубовой жидкости из-за наличия перекрестных связей по управляющим воздействиям этих регуляторов приводит к их неустойчивой работе. Для решения данной проблемы предлагается использование математической модели в системе управления ректификационной колонной, позволяющей предсказывать состав дистиллята (содержание ДЭА в дистилляте) и состав кубовой жидкости (содержание ДЭА в кубовой жидкости) в зависимости от величины управляющих воздействий L (расход флегмы) и G (паропроизводительность кипятильника) соответственно. Таким образом, задача управления процессом ректификации может быть поставлена как задача оптимального управления, заключающаяся в поиске значений L и G с использованием математической модели, при которых величины и соответствуют технологическому регламенту. Найденные таким образом величины L и G выдаются в качестве заданий регуляторам, стабилизирующим эти технологические параметры [1]. Технологическая схема процесса бинарной ректификации представлена на рис. 1. Рис. 1. Технологическая схема ректификации смеси диэтаноламин - триэтаноламин: 1 - выносной кипятильник; 2 - колонна; 3 - дефлегматор Составим требуемую для решения задач управления процессом ректификации математическую модель, а затем определим ее параметры. Расчет состава выходных потоков ректификационной колонны обычно проводится итерационным методом [2] с уточнением показателя состава кубовой жидкости . По технологическому регламенту эта величина меняется от 0,2 до 0,3 кмоль/кмоль; в качестве начального приближения принималось среднее регламентное значение мольной доли низкотемпературного компонента (ДЭА) в кубовой жидкости кмоль/кмоль. . (1) В соответствии с общепринятым подходом гидродинамика процесса в кубовой части колонны описывается с помощью модели идеального смешения, поскольку в ней происходит кипение кубовой жидкости. Состав пара, поступающего из куба колонны на насадку, определяется уравнением равновесия пар - жидкость , (2) где - содержание ДЭА в паровой фазе, кмоль/кмоль; - содержание ДЭА в кубовой жидкости, кмоль/кмоль; - коэффициенты формализованной зависимости. Гидродинамика процесса в насадке колонны описывается как модель идеального вытеснения, поскольку площадь поперечного сечения материальных потоков в насадке сравнительно мала. Модель идеального вытеснения представляется двумя дифференциальными уравнениями для жидкости и пара, определяющими изменения состава этих потоков по высоте насадки. Граничные условия для интегрирования дифференциальных уравнений определяются из уравнения материального баланса для кубовой части колонны. При значении координаты высоты сечения насадки h = 0 мольная доля низкотемпературного компонента в паровом потоке , а в потоке жидкости , (3) где G - расход пара, кмоль/ч; W - расход кубовой жидкости, кмоль/ч; L - расход флегмы, кмоль/ч; F - расход потока питания колонны, кмоль/ч. Дифференциальное уравнение материального баланса для потока пара: , (4) где - объемный коэффициент массопередачи, кмоль/(); - равновесная мольная доля низкотемпературного компонента в паровой фазе, кмоль/кмоль; h - высота бесконечно малого элемента насадки, м. Уравнение (4) справедливо для всей насадки. Дифференциальное уравнение материального баланса потока жидкости для части насадки, расположенной ниже подвода питания, принимает вид . (5) Аналогичное уравнение справедливо и для другой части насадки: . (6) Уравнения (4) и (5) интегрируются до сечения по высоте насадки, соответствующего подаче на насадку питания (h = H1), а затем интегрируются уравнения (4) и (6) до верхнего сечения насадки (h = Н). При интегрировании дифференциальных уравнений (4), (5) и (6) величина равновесной мольной доли низкотемпературного компонента в паровой фазе определяется при переменных значениях Х по уравнению равновесия пар - жидкость, аналогичному уравнению (2): . (7) Значение величины , рассчитанной при h = Н, соответствует составу пара, поступающего с насадки на дефлегматор. В приближении полной конденсации пара в дефлегматоре состав конденсата (флегма и дистиллят) имеет показатель состава . Уравнение материального баланса для всей колонны позволяет уточнить состав кубовой жидкости (найти расчетное значение) : , (8) где - содержание ДЭА в питании, кмоль/кмоль; D - расход дистиллята, кмоль/ч. Для уточнения значения величины , принятой в качестве начального приближения, используется улучшенная итерационная процедура с инерционным множителем μ = 0,01 - 0,05: . (9) Если погрешность итерационного расчета не соответствует допустимой погрешности , , то итерационный расчет повторяется начиная с расчета по уравнению (2). После достижения допустимой погрешности итерационный расчет завершается. Уравнения (1) - (9) и порядок их применения (моделирующий алгоритм) представляют собой математическую модель процесса ректификации, реализованного в насадочной колонне. Эта модель при наличии данных о значениях неуправляемых технологических параметров F и , а также данных о значениях управляющих воздействий L и G позволяет определить значения выходных технологических параметров процесса D, W, и. Блок-схема модели процесса ректификации представлена на рис. 2. Порядок применения уравнений на схеме при расчете математической модели представлен снизу вверх. Для практического использования математической модели процесса ректификации с целью управления им необходимо определить ее параметры. К таким параметрам относится объемный коэффициент массопередачи , используемый в уравнениях (4), (5) и (6). Кроме того, ввиду отсутствия данных по равновесию системы пар - жидкость компонентов ДЭА и ТЭА коэффициенты уравнений (2) и (7) (одинаковые по величине для обоих уравнений) также являются неизвестными параметрами математической модели процесса ректификационного разделения бинарной смеси ДЭА - ТЭА. Общее количество параметров рассматриваемой модели, подлежащих определению, составляет пять. Рис. 2. Блок-схема модели процесса ректификации этаноламинов Наличие математической модели, включающей пять неизвестных параметров, позволяет определить эти параметры с использованием опытных данных по промышленной эксплуатации процесса ректификации смеси ДЭА и ТЭА при наличии данных о работе колонны в различных режимах (не менее пяти), отличающихся значениями входных технологических параметров. Комплект опытных данных должен включать значения всех входных и выходных технологических параметров работы ректификационной колонны. В данной работе для определения параметров модели собрана информация о двенадцати режимах функционирования ректификационной колонны (обучающая информация), разделяющей смесь ДЭА и ТЭА (табл. 1). Кроме обучающей информации здесь приводится информация еще о пяти режимах работы этой колонны для тестирования адекватности модели при найденных значениях ее параметров. Данная информация представлена в табл. 2, там же приводятся данные расчета показателей состава выходных потоков и , выполненного с использованием модели колонны. Для снижения погрешности опытных данных, представленных в табл. 1 и 2, выполнялся предварительный отбор таких данных, по которым соблюдался материальный баланс по расходам компонентов с погрешностью не более 1 % от величины расхода. Экспериментальные данные, приведенные в табл. 1 и 2, после получения на производстве этаноламинов были предварительно преобразованы путем изменения размерности величин так, чтобы расходы потоков были выражены в размерности Кмоль/час и все показатели состава потоков - в мольных долях низкотемпературного компонента диэтаноламина. Мольная доля высокотемпературного компонента в соответствующем потоке определяется вычитанием из единицы мольной доли низкотемпературного компонента. Обработка опытных данных с целью определения параметров модели процесса ректификации проводилась следующим образом. Для каждого из двенадцати режимов работы ректификационной колонны по значениям входных технологических параметров рассчитывались с использованием математической модели процесса величины выходных технологических параметров и (расчетные значения), которые сравнивались с экспериментальными значениями этих параметров для данного режима работы колонны. Сумма квадратов расхождений между расчетными и экспериментальными значениями выходных параметров по всем опытам рассматривалась как критерий оптимальности при выборе значений параметров математической модели: . Оптимальный выбор пяти параметров математической модели выполнялся путем минимизации критерия оптимальности Q методом Розенброка [3] с учетом ограничений на расчетные величины выходных технологических параметров. Эти ограничения относятся к расчетным величинам мольных долей низкотемпературного компонента в выходных потоках колонны и состоят в том, что они должны находиться в диапазоне от нуля до единицы: Результаты определения параметров модели для процесса ректификации в насадочной колонне представлены в табл. 3, там же даны стартовые значения (начало поиска) этих параметров и величины суммы квадратов расхождения между расчетными и опытными значениями выходных технологических параметров Q по всем опытам. Уточнение значений параметров математической модели для процесса ректификации в результате применения метода наименьших квадратов оказалось существенным. Таким образом, поиск оптимальных значений параметров математической модели для процесса ректификации бинарной смеси ДЭА и ТЭА прошел успешно. Таблица 3 Результаты определения параметров модели для процесса ректификации в насадочной колонне Приближение Парамeтры модели для процесса ректификации Q Стартовое 0,50781 0,1003 4,3613 -7,4548 4,1831 0,9479 Финальное 9,4597 0,6432 2,2979 -4,00032 2,104253 0,000582 Однако для установления адекватности составленной модели ректификации этаноламинов с найденными здесь значениями ее параметров необходимо сравнить дисперсии расхождения опытных и расчетных данных, полученных при обучении (см. табл. 3) и тестировании адекватности модели (см. табл. 2). Дисперсия, достигнутая при обучении: , где - сумма квадратов расхождений, достигнутая при обучении; - число степеней свободы дисперсии обучения. Согласно данным табл. 3: . Значение определяется количеством опытных данных в табл. 1, количеством выходных параметров ( и ), количеством параметров модели, найденных по данным табл. 1: , где n - количество опытных данных (n = 12); m - количество выходных параметров модели (m = 2); k - количество неизвестных параметров модели (k = 5). Таким образом, дисперсия, достигнутая при обучении, равна . Сумма квадратов расхождений между опытными и расчетными данными по величинам и , полученная при тестировании адекватности модели (согласно данным табл. 2): . Дисперсия, достигнутая при тестировании: . Дисперсионное отношение: . Критическое значение дисперсионного отношения при числе степеней свободы числителя , числе степеней свободы знаменателя и уровне значимости равно [4]. Поскольку , можно считать, что тестирование модели показало ее адекватность с вероятностью 0,95. Таким образом, представленная модель может быть использована для управления технологическим процессом ректификации этаноламинов в насадочной колонне. Результаты исследования переданы в ООО «Синтез Ока» (г. Дзержинск Нижегородской области) для дальнейшего внедрения на предприятии. Полученные алгоритмы могут быть применены и на других производственных площадках для проведения ректификации бинарных смесей в колоннах насадочного типа с целью получения двух целевых продуктов.

About the authors

Nikolay P Peskov

Dzerzhinsky Polytechnic Institute (branch) of the federal state budgetary educational institution of higher professional education "Nizhny Novgorod State Technical University name R.E. Alexeyev"

Email: peskov87@bk.ru
Senior Lecture 49, Gaidar st., Dzerzhinsk, Nizhny Novgorod region, 606026, Russian Federation

Ernst M Moncharzh

Dzerzhinsky Polytechnic Institute (branch) of the federal state budgetary educational institution of higher professional education "Nizhny Novgorod State Technical University name R.E. Alexeyev"

(Ph.D (Chem.)), Associate Professor 49, Gaidar st., Dzerzhinsk, Nizhny Novgorod region, 606026, Russian Federation

Anatoly A Popov

Dzerzhinsky Polytechnic Institute (branch) of the federal state budgetary educational institution of higher professional education "Nizhny Novgorod State Technical University name R.E. Alexeyev"

(Ph.D. (Techn.)), Associate Professor 49, Gaidar st., Dzerzhinsk, Nizhny Novgorod region, 606026, Russian Federation

Nataliya O Kuligina

Dzerzhinsky Polytechnic Institute (branch) of the federal state budgetary educational institution of higher professional education "Nizhny Novgorod State Technical University name R.E. Alexeyev"

Senior Lecture 49, Gaidar st., Dzerzhinsk, Nizhny Novgorod region, 606026, Russian Federation

References

Supplementary files