Математические модели управляемых массоИ теплообменных процессов в технологическом комплексе систем «Абсорбция-десорбция» Текст научной статьи по специальности «Химические технологии»

Radionov Ivan Alekseevich

Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”.

E-mail: radionov87@gmail.com.

2, Chehov Street, Ap. 403, Taganrog, 347900, Russia.

Phone: +78634360707.

The Department of Synergetics and Control; Postgraduate Student.

УДК 681.5.017:66.08.012-52

C.E. Абрамкин, C.E. Душин, H.H. Кузьмин

МАТЕМАТИЧЕСКИЕ МОДЕЛИ УПРАВЛЯЕМЫХ МАССО-И ТЕПЛООБМЕННЫХ ПРОЦЕССОВ В ТЕХНОЛОГИЧЕСКОМ КОМПЛЕКСЕ СИСТЕМ АБСОРБЦИЯ-ДЕСОРБЦИЯ»

Приведены основные ограничения и допущения, положенные в основу получения математических моделей систем «Абсорбция» и «Десорбция». Разработаны нелинейные динамические модели взаимосвязанных управляемых массообменных и теплообменных процессов замкнутого технологического комплекса систем «АБСОРБЦИЯ-ДЕСОРБЦИЯ». Компьютерное моделирование проводилось с использованием программного средства MATLAB/Simulink. Особое внимание уделено получению математических моделей подсис-« » « ».

Математическая модель; структурная схема; природный газ; абсорбция; десорбция; ; ; .

S.E. Abramkin, S.E. Dushin, N.N. Kuzmin

MATHEMATICAL MODELS OF CONTROLLED MASS- AND HEAT EXCHANGE PROCESSES IN TECHNOLOGIC COMPLEX OF SYSTEMS «ABSORPTION-DESORPTION»

Paper shows the results of analysis of the gas dewatering process. The basic restrictions and the assumptions taken as a principle of reception of mathematical models of systems absorption and desorption are presented. Nonlinear dynamic models interconnected controlled mass and heat exchange processes of the closed technological complex of the systems «ABSORPTION-DESORPTION» are developed. Computer modeling was carried out with use of software MATLAB/Simulink. The special attention is given reception of mathematical models of subsystems «Evaporation» and «Air cooling».

Mathematical model; structure scheme; natural gas; absorption; desorption; mass exchange; heat exchange; controlled process.

В настоящее время газовая промышленность составляет значительную часть российской экономики. Необходимость удовлетворения как внутреннего, так и внешнего спроса на газ определяет актуальность проблем управления процессами в добыче, подготовке и транспортировке газа. В программах ОАО «Г^пром» декларируется переход к малолюдным технологиям. Такой переход может быть осуществлен в результате разработки современных автоматических систем управления процессами добычи, подготовки и транспортировки природного газа. В рамках данного направления развития газовой отрасли проводится реконструкция технологических аппаратов и внедрение автоматизированных систем управления технологическими процессами на газовых промыслах. Проведение данных работ определено высокими требованиями к качеству подготовки природного газа для его

транспортировки, повышением надежности технологического оборудования в связи с непрерывным характером производства и стремлением к оптимизации режимов работы скважин и установок комплексной подготовки газа (УКПГ).

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

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

Технологический процесс подготовки природного газа методом абсорбционной осушки на УКПГ, осуществляемый по замкнутому циклу (по жидкой поглощающей фазе), можно представить как системный комплекс «Абсорбция - Десорбция», в который входят две взаимосвязанные системы «Абсорбция газа» и «Десорбция абсорбента» [1]. Процесс физической абсорбции газа протекает при низкой температуре и высоком давлении, а процесс десорбции, напротив, происходит при высокой температуре и низком давлении.

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

Система «Десорбция абсорбента» также состоит из ряда подсистем, основными из которых являются: подсистема «Ректификация» (колонна регенерации), подсистема «Выпаривание» (испаритель) и подсистема «Воздушное охлаждение», предназначенная для охлаждения флегмы (рефлюкса) в аппарате воздушного ох.

комплекса на процесс осушки.

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

Всестороннее изучение вопроса абсорбционной осушки газа [2], [3] показало, что при построении адекватных математических моделей (ММ) процессов абсорбции и десорбции важная роль должна отводиться нелинейным динамическим мо, ( ) « ». Только нелинейные ММ позволяют выявить и исследовать взаимосвязь массо- и

.

Работа замкнутого комплекса систем «Абсорбция-Десорбция» может сопровождаться тремя видами нестационарных процессов, обусловленных:

1) установлением стационарного рабочего режима;

2) изменением внешних условий (например, состоянием газового пласта, климатическими и погодными возмущениями и т. п.);

3)

.

Нестационарные процессы первого вида возникают во время подачи абсорбента в абсорберы при настройке систем «Десорбция» и «Абсорбция» на рабочие , . процессы первого вида возникают при переключениях УКПГ с одной технологической линии на другую.

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

« ».

В рассматриваемых условиях функционирования комплекса при установленных ограничениях и правильной его эксплуатации, нестационарные процессы , , .

При разработке ММ технологических процессов важной задачей является правильный выбор физических принципов и закономерностей, допущений и ограничений, накладываемых на модель. В результате анализа поведения систем «Аб» « » -:

♦ жидкость (абсорбент) в насадках колонн (абсорбционной и регенерирую-

) ( );

♦ в основе получения ММ массообменных процессов лежит гидродинами-

( );

♦ режим продольного пере мешивания фаз отсутствует;

♦ скорости газа ь?, пара ьп и жидкости ьа одинаковы по сечениям аппа-

;

♦ концентрации газа С,, пара Сп и абсорбента Са одинаковы по сечениям

;

♦ массопередача на насадках эквимолярная, т. е. через поверхность раздела фаз в противоположных направлениях переносится одинаковое количество целевого компонента;

♦ абсорбент, проходя через насадку, не задерживается, что исключает образование зон простоя в насадках, а также исключается зависание абсорбента в насадках;

♦

( ) ( закона Г енри), для процесса десорбции - как линейной, так и нелинейной;

♦ температурные профили газа и абсорбента по всей высоте абсорбера при.

С учетом принятых допущений нелинейная динамическая ММ управляемого массообменного процесса в абсорбере может быть представлена системой дифференциальных уравнений в частных производных (ДУЧП):

ЭС, ЭС ,

^ ~д^ - R [С - С$ (Сца)];

ЭС ЭС

= f (и)-£ + R„[С,,. - С'.(С„)], (1)

0< г <1а , t >0,

где Сиг, Сиа - концентрации целевого компонента в газе и жидкости (абсорбенте); Сир„, С1а - равновесные концентрации целевого компонента; \)иг - скорость целевого компонента в газе, которая не меняется по высоте абсорбера; \а - высота абсорбера. Постоянные коэффициенты Л, и Ла определяются физическими

свойствами фаз и геометрическими особенностями аппарата. Управление процессом абсорбции осуществляется изменением подачи жидкого компонента в абсорбер, что отражено в модели введением сомножителя / (и) . Соответствующая

ММ структурная схема приведена на рис. 1.

При выполнении закона Генри система уравнений (1) приобретает вид:

дС дС

= ІГ-Л (С- -Е„С„);

д Ъг (2)

дС , ЪС

д = Г + К ~ ЕСа ),

где Е - коэффициент фазового равновесия, пропорциональный коэффициенту

Генри, для массовых концентраций. В силу присутствия управляющего воздействия во втором уравнении модель (2) сохраняет нелинейный характер.

Рис. 1. Структурная схема ММ управляемого массообменного процесса в

абсорбере

(2)

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

Получение ММ подсистемы «Ректификация» системы «Десорбция» основывается на ранее принятых допущениях, что позволяет представить процессы мас-сообмена нелинейной системой ДУЧП:

ді “ ді

ЪС , ,ЪС

- кЬ,»-с(сЛ'

ип

д! + К (3)

0 < г < I.., ! > 0,

к 7 7

где Сш , Сип - концентрации целевого компонента в потоках абсорбента и пара; С1а - равновесная концентрация целевого компонента в абсорбенте; ьа , ьп -скорости потоков абсорбента и пара; вп - температура пара; /к - высота регенерационной колонны. Коэффициенты К и К определяются физическим состоянием фаз и геометрическими размерами аппарата. При этом скорость пара ьп и коэффициент К зависят от температуры пара.

При выполнении закона Генри система уравнений (3) принимает вид:

дС„„ дС

'"--Ч,^-К (С„ - Б(в„)С„,);

д! '' дг (4)

дС дС

-др = Ч (в„ )-£- + К (в К, - Е (в, )С,„,

причем коэффициент фазового равновесия Е(ва) зависит от температуры абсор-. , (4) , -

пературы абсорбента ва и пара вп постоянны. Модель дополняется граничными и .

Для подсистемы «Ректификация» в переходных режимах температурные профили пара и абсорбента по высоте колонны не постоянны. В связи с этим было проведено моделирование тепловых процессов, протекающих в колонне:

= ~Ьп (в )~д~ ~ (в )(0« -®а ),

д! дг (5)

1Г = +Кв^-9*)

где вп , ва - температуры пара и абсорбента; ьп, ьа - скорости тепловых потоков пара и абсорбента; , Ква - теплофизические коэффициенты. В данной мо-

дели учитывается влияние температуры пара на его физические характеристики.

(5) -

филями распределения температур вдоль колонны. Структурная схема подсистемы « » . 2.

,

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

Управление взаимосвязанных массообменных и теплообменных процессов в подсистеме «Ректификация» осуществляется за счет изменения температуры в

подсистеме «Выпаривание» и поддержания температуры флегмы в подсистеме «Воздушное охлаждение».

Рис. 2. Структурная схема ММ подсистемы «Ректификация»

Объектом моделирования в подсистеме «Выпаривание» является испаритель. К топочной трубе испарителя подсоединяется газовая горелка, при помощи которой происходит управление тепловыми процессами в аппарате. В топочную трубу испарителя подается пламя от газовой горелки, которое нагревает трубное пространство и в результате нагрева образуются дымовые газы, поступающие затем в . . В межтрубное пространство испарителя из насадочной секции ректификационной колонны стекает частично регенерированный абсорбент. Дымовые газы нагревают стенки жаровых труб. В свою очередь, за счет отдачи тепла от стенок жаровых труб осуществляется нагревание абсорбента до заданной температуры. В результате нагревания абсорбента происходит выпаривание из него паров воды. В испарителе теплообмен между фазами происходит при движении по противоточной схеме, так как абсорбент течет в одном направлении, а дымовые газы - в противопо-.

При построении динамической ММ теплообменных процессов в подсистеме «Выпаривание» учитывались следующие допущения:

♦ продольное перемешивание в каждом из потоков теплоносителей отсут-, . . вытеснения;

♦ поперечное перемешивание в потоках теплоносителей считается идеальным. Это предположение позволяет разрабатывать одномерную ММ, в которой учитывается изменение температуры только по длине аппарата;

♦ термическое сопротивление стенок теплообменника сравнительно мало.

, .

Принимая во внимание эти допущения, была получена модель с распределенными параметрами в виде системы дифференциальных уравнений первого по,

, :

где $дг,6а,9ст - температуры дымовых газов, абсорбента и стен к и; > >>>>, >> а -скорости дымовых газов и абсорбента; ^, Ка, Ядгст, Каст - физико-

, -, .

Для ММ (6) были определены начальные распределения температур теплоносителей и граничные условия, характеризующие исследуемый процесс. Начальные распределения температур теплоносителей:

Однако при исследовании ММ (6) средствами Ма1;ЬаЪ/81шиИпк было выявлено, что при определенных условиях температура абсорбента недопустимо уменьшалась, так как температура дымовых газов не достигала значения, необходимого для нагрева абсорбента до заданной температуры. Для решения данной проблемы была разработана нелинейная ММ управляемого теплообменного процесса в испарителе:

где Стг - управляемый расход топливного газа.

Исследование полученной ММ (7) показало, что управляемый теплообменный процесс в испарителе позволяет поддерживать температурный режим и концентрацию абсорбента на требуемом уровне. А поддержание концентрации абсорбента в установленных техническим заданием пределах и является конечной це-

« ».

Выполнение этого условия означает, что температура стенки

(6)

0*0 (х ) = Мх * і ^ Є«0 (х )= * I=0.

:

и*. в„(х, * ),=о=вг().

(7)

ст ________

( - ®ст ) - Кет (( - ).

/ (и )= О,,

Подсистема «Воздушное охлаждение» также представляет собой важнейший элемент в системе «ДЕСОРБЦИЯ абсорбента». Объектом исследования в данной подсистеме являются технологические процессы в аппарате воздушного охлаждения (АВО), который предназначен для конденсации паров, выходящих с верхней наса-дочной секции колонны регенерации. На выходе АВО, после конденсации, образу-( ), -ной части колонны. Одновременно с процессом конденсации в рабочем пространстве АВО происходит накопление воздуха и других неконденсирующихся газов, которые выделяются из парожидкостной смеси, а также проникают через неплотности в соединениях аппаратуры из окружающего воздуха. Накопление неконденсирующихся газов и рост их парциального давления приводят к уменьшению разрежения в .

заданном уровне производится непрерывный отвод из АВО неконденсирующихся газов. Для этого в системе «ДЕСОРБЦИЯ» применяется вакуум-насос, который также служит для предотвращения колебаний давления, обусловленных изменением

( ).

При составлении математической модели (ММ) теплообменных процессов в АВО были приняты следующие допущения:

♦ в конденсационных секциях среда движется в режиме идеального вытес-

♦ поперечное перемешивание считается идеальным. Данное допущение приводит к составлению одномерной модели, в которой учитывается динамическое распределение температуры только по длине конденсационных секций АВО;

♦ принимается, что теплообменный процесс сопровождается конденсацией

;

♦ в ММ учитывается накопление теплоты в трубах конденсационной секции;

♦ трубы в конденсационных се кциях выполнены с оребрением.

Принимая во внимание вышеперечисленные допущения, была получена динамическая ММ теплообменного процесса в конденсационных секциях АВО:

где внп, вса, вср - температуры насыщенного пара, стенки теплообменных секций

АВО и окружающей среды; ьип - скорость насыщенного пара; К1т, Яс1, Яс2 -

физико-технологические коэффициенты, зависящие от физических свойств насыщенного пара и материала стенки теплообменных секций АВО.

Второе уравнение в системе (8) является обыкновенным дифференциальным уравнением первого порядка. Для его однозначного решения необходимо знать

начальное условие вса в момент времени ? —0.

в :

где 0 0 - начальная средняя температура окружающей среды при неработающем вентиляторе; 0 - установившееся значение температуры окружающей среды;

;

(В)

в = f q),

(9)

О„.

максимальный расход воздуха, нагнетаемый вентилятором; Ц - коэффи-

циент, принадлежащий интервалу значений, близких к единице.

Граничные условия для температурного профиля пара ММ (8) определяются :

'I,—0 — «:(') 0.»(х.гI,—— 67 (г),

где 1ат - длина трубы теплообменника.

Начальные распределения температур вдоль трубы АВО задаются выраже-:

г 1—0 —К,(,\ »„(,,г),.0 — «“(,)

Для обеспечения требуемой температуры флегмы на выходе из аппарата, что является основной задачи регулирования на основе ММ (8), была разработана ММ управляемого теплообменного процесса в конденсационных секциях АВО:

Э0.,„ Э0.

д*

О»

йї

■ = -и.,

дх

■- Кп -9«Л

(10)

= Кс1/(и ) + К2^- К^,

где К = Яс1 + Яс2, /(и) - управляющее воздействие.

Структурная схема ММ управляемого теплообменного процесса в конденсационных секциях АВО представлена на рис. 3.

Рис. 3. Структурная схема ММ управляемого теплообменного процесса

Исследование ММ (10) проводилось с применением программного средства Ма1ЬаЪ/81шиИпк. Было исследовано поведение управляемой величины (температуры флегмы в выходном коллекторе АВО) с применением линейных типовых . -

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

, -

вые математические модели взаимосвязанных управляемых массо- и теплообменных процессов для комплекса систем «АБСОРБЦИЯ-ДЕСОРБЦИЯ», проблемноориентированные на задачи управления профилями концентраций и температур

. , -

ские модели являются замкнутыми и учитывают нелинейный характер взаимодей-.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Абрамкин С.Е. Анализ технологического процесса абсорбционной осушки газа как объекта управления // Изв. СПбГЭТУ <^ЭТИ». - 2007. - Вып. 3. - С. 24-31.

2. Абрамкин С.Е., Душин С.Е. Математическая модель массотеплообменных процессов технологического комплекса абсорбционной осушки газа // Материалы VI научной конференции «Управление и информационные технологии». - СПб., 2010. - С. 220-224.

3. Абрамкин СЕ.,Душин С.Е. Разработка математической модели технологического комплекса Абсорбция - Десорбция» // Изв. СПбГЭТУ <^ЭТИ». - 2011. - В ып. 1. - С. 29-33.

Статью рекомендовал к опубликованию д.т.н., профессор Г.Е. Веселов.

Абрамкин Сергей Евгеньевич

- « ».

E-mail: abrsergey@yandex.ru.

. - , , . 4, . 11.

Тел.: +79217767124.

Кафедра автоматики и процессов управления; аспирант.

Душин Сергей Евгеньевич

E-mail: dushins@yandex.ru.

Тел.: +79219704631; 88122343798.

; . . .; .

Кузьмин Николай Николаевич

E-mail: nnkuzmin@mail.eltech.ru.

Тел.: +79591482; 88122343798.

Заведующий кафедрой автоматики и процессов управления; к.т.н.; профессор. Abramkin Sergey Evgen’evich

Saint-Petersburg State Electrotechnical University «LETI».

E-mail: abrsergey@yandex.ru.

4-11, Kotelnikova Avenue, Saint-Petersburg.

Phone: +79217767124.

The Department of Automation and Control Processes; Post-graduate Student.

Dushin Sergey Evgen’evich

Phone: +79219704631; +28122343798.

E-mail: dushins@yandex.ru.

The Professor of the Department of Automation and Control Processes; Dr. of Eng. Sc.; Associate Professor.

Kuzmin Nikolai Nikolaevich

E-mail: nnkuzmin@mail.eltech.ru.

Phone: +79591482; work phone: +78122343798.

Head of Automation and Control Processes Department; Dr. of Phil. Sc.; Professor.