|
|
11 |
|
дг t дг дг x Rдг дг cт ; |
|
||
ж |
t ж ж |
x Rж ст ж ; |
(3) |
d ст |
dt Rдгст дг ст Rжст ст ж , |
|
|
где дг , ж , ст – температуры дымовых газов, жидкой фазы и стенки; |
дг , ж – |
||
скорости дымовых |
газов и |
абсорбента; Rдг , Rж , Rдгст , Rжст – |
физико- |
технологические коэффициенты, зависящие от физических свойств дымовых газов, жидкости и материала стенки жаровой трубы.
Граничные условия: |
|
x |
|
вх t , |
|
|
x |
|
|
вх t . |
|
|
|
|
|
||||||
|
дг |
|
x l |
дг |
|
ж |
|
|
x 0 |
ж |
|
|
|
|
|
Начальные условия: дг0 x дг x,t t 0 , ж0 x ж x,t t 0.
Для подсистемы «Воздушное охлаждение» получена динамическая ММ теплообменного процесса в конденсационных секциях АВО:
нп x, t t нп нп нп x, t x Rнп нп x, t ст x, t , |
(4) |
|
d ст dt Rс2 нп ст nRс1 ст ср , |
||
|
где нп , ст , ср – температуры насыщенного пара, стенки теплообменных секций АВО и окружающей среды; Rнп , Rс1, Rс2 – физико-технологические коэффициен-
ты, зависящие от физических свойств насыщенного пара и материала стенки теплообменных секций АВО; n – количество трубопроводов в АВО.
Граничные условия: |
нп |
x |
|
|
вх |
t ; |
|
нп |
x |
|
|
вых t , |
|
|
|
|
|||||||||
|
|
|
x0 |
нп |
|
|
|
|
x lаво |
нп |
||
|
|
|
|
|
|
|
|
|
|
|
где lаво – длина трубы теплообменника.
Начальные условия: |
|
x |
|
|
вх |
x ; |
|
|
x |
|
|
вх x . |
|
|
|
|
|
|
|||||||
|
нп |
|
|
t 0 |
нп |
|
|
ст |
|
|
t 0 |
ст |
|
|
|
|
|
|
|
Результаты моделирования подсистем «Выпаривание» и «Воздушное охлаждение» представлены на рисунке 3.
вых |
|
|
|
|
|
θж , °С |
|
|
|
|
|
170 |
|
|
|
|
|
165 |
|
|
|
|
|
160 0 |
|
|
|
|
t, ×10^3, |
1 |
2 |
3 а 4 |
5 |
6 |
Рисунок 3. Графики переходных
вых
θнп, °С
45
40
35
с |
|
|
|
|
300 |
5 |
10 |
б |
15 |
|
|
|
|
процессов в испарителе (а) и АВО (б)
t, с 20
На рисунке 3, а приведен график переходного процесса в испарителе при повышении температуры жидкой фазы на входе в испаритель, на рисунке 3, б − график переходного процесса в АВО при повышении температуры охлаждающего воздуха. Время установления нового равновесного состояния составляет около 1,5 ч и 20 с соответственно.
12
Исследование моделей подсистем «Выпаривание» и «Воздушное охлаждение» позволило сделать заключение, что установившиеся значения, полученные при проведении моделировании, либо не достигают заданного уровня, либо превышают его. Эти результаты подтверждают необходимость применения систем автоматического регулирования тепловых процессов для подсистем «Выпаривание» и «Воздушное охлаждение», которые в свою очередь оказывают влияние на течение тепломассообменных процессов в подсистеме «Ректификация».
В четвертой главе сформулированы и решены задачи управления массо- и теплообменными процессами при абсорбционной осушке природного газа. Исходя из сформулированных задач были разработаны ММ управляемых массо- и теплообменных процессов в подсистемах «Абсорбция газа», «Выпаривание» и «Воздушное охлаждение».
Для получения управляемой модели подсистемы «Абсорбция газа» (2) во второе уравнение введена управляющая функция f (u) zж . Управление процес-
сом абсорбции осуществляется изменением подачи жидкой фазы в абсорбер. Для этих целей предложено использовать многорежимное регулирование соотношения расходов «газ-абсорбент». Многорежимный нелинейный П-регулятор реализован в программной среде MATLAB/Simulink.
Эффективная работа абсорбера осуществляется при изменении расхода газа от 250000 м3/ч до 500000 м3/ч. При этом скорость газа может варьироваться в пределах от 0,29 м/с до 0,57 м/с. Результаты моделирования в системе с многорежимным регулятором при подаче возмущающих воздействий различного типа по скорости газа на входе абсорбера представлены на рисунке 4. Графики переходных процессов соответствуют минимальному уровню возмущающих воздействий.
По результатам вычислительных экспериментов можно сделать вывод, что система автоматического регулирования, реализованная в виде многорежимного нелинейного П-регулятора, приводит объект в заданное состояние.
вых |
3 |
|
Cцг, г/м |
|
|
0,027 |
вых |
3 |
|
Cцг, г/м |
|
0,025 |
0,022 |
|
|
|
|
0,023 |
|
|
|
|
|
|
|
|
|
|
t, с ×10^4 |
|
|
|
|
|
|
|
t, с ×10^4 |
|
0,0210 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
0,020 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
|
|
|
|
а |
|
|
|
|
|
|
|
|
б |
|
|
|
|
Рисунок 4. Графики переходных процессов по концентрации целевого компонента в газе при прямоугольном импульсном (а) и ступенчатом (б) возмущениях
Анализ концептуальной модели КТС показывает, что влиять на массотеплообменные процессы в ректификационной колонне возможно посредством управления тепловыми процессами в испарителе и АВО.
Регулирование теплообменных процессов в подсистеме «Выпаривание» осуществляется за счет изменения расхода топливного газа или газовоздушной сме-
13
си, подаваемых в горелку с целью обеспечения требуемой температуры абсорбента на выходе испарителя. Регулирование реализовано введением в первое уравнение ММ подсистемы «Выпаривание» (3) функции управления f u дг Gтг .
Система регулирования предусматривает коррекцию по температуре дымовых газов. В качестве исполнительного устройства выбирался клапан-регулятор с мембранным исполнительным механизмом и позиционером. Для регулирования процессов применялся ПИ-закон. Требуемая величина температуры абсорбента на выходе испарителя составляет 160 С. Вычислительные эксперименты показали, что система регулирования температуры абсорбента в испарителе позволяет поддерживать заданную температуру.
Основной задачей системы регулирования теплообменных процессов в подсистеме «Воздушное охлаждение» является поддержание заданной температуры жидкой фазы на выходе из аппарата. Регулирование обеспечивается введением во второе уравнение ММ подсистемы «Воздушное охлаждение» (4) функции управления f (u) ср Gв . Процесс регулирования температуры флегмы на выходе
АВО производится изменением расхода охлаждающего воздуха, обдувающего конденсационные секции аппарата. В качестве исполнительного устройства выбирался частотно-регулируемый электропривод. Регулирование производилось по ПИ-закону. Требуемая величина температуры флегмы на выходе АВО составляет 34 С. Результаты вычислительных экспериментов показали, что система регулирования обеспечивает требуемую величину температуры.
Как следует из анализа концептуальной модели КТС «АБСОРБЦИЯ– ДЕСОРБЦИЯ», его главной особенностью является замкнутость технологической схемы циркуляции раствора абсорбента (замкнутый цикл по жидкой фазе). В связи с этим возникает задача исследование процессов в замкнутом контуре. «Замыкание» ММ двух подсистем «Абсорбция газа» и «Ректификация» производилось по величине концентрации ЦК в жидкой фазе.
Дискретно-непрерывная управляемая модель КТС «АБСОРБЦИЯ–
ДЕСОРБЦИЯ» записывается в виде следующей системой уравнений: |
||
dСцгi 1 dt г h Rг г Сцгi 1 Rг г EрСцжn-i г h Cцгi ; |
||
dСцжi 1 dt f u h Rж Eр Сцжi 1 |
RжCцгn-i f u h Cцжi ; |
|
|
|
(5) |
dСцпi 1 dt п h Rп Сцпi 1 г h Cцпi Rп EСцжn-i ; |
||
dСцжi 1 dt ж h Rж ж E Сцжi 1 ж h Сцжi Rж ж Сцпn-i , |
||
|
|
|
где i 0,...,n 1, |
г , п , ж – скорости |
газовой, паровой и жидкой фаз; |
Cцг ,Cцп ,Cцж – концентрации ЦК в газовой, паровой и жидкой фазах; Eр , E – коэффициенты фазового равновесия; f u – управляющее воздействие по скорости
абсорбента в абсорбере; R , R , R , R – физико-технологические коэф-
г г ж ж п ж
фициенты по газообразной и жидкой фазам, h – шаг дискретизации.
Для регулирования взаимосвязанных массообменных процессов КТС «АБ- СОРБЦИЯ–ДЕСОРБЦИЯ» так же применяется многорежимное регулирование.
14
На рисунке 5 представлена КМ управляемых массообменных процессов КТС «АБСОРБЦИЯ–ДЕСОРБЦИЯ». Вычислительные эксперименты с КМ проводились в условиях приложения ступенчатых и последовательности прямоугольных импульсных возмущающих воздействий по скорости газа на входе абсорбера.
|
|
|
|
|
|
подсистема |
|
|
|
|
|
|
|
«Ректификация» |
|
|
n=υа/υг |
|
|
|
kтр |
Out1 |
|
|
|
|
|
|
|
In1 |
|
|
1-D |
|
|
|
|
|
|
|
|
|
|
|
Out2 |
|
|
|
T(u) |
|
|
|
|
|
|
|
подсистема |
|
|
|
|
||
|
|
|
|
|
|
||
|
|
«Абсорбция газа» |
|
РА |
Out3 |
|
|
|
|
|
|
Газ |
|
|
|
|
|
|
|
|
|
|
|
|
υг |
In1 |
Out1 |
|
|
Out4 |
In2 |
|
|
In2 |
Out2 |
|
|
Out5 |
|
|
|
|
|
|
|
||
|
|
|
Out3 |
|
|
|
|
|
Cцг0 |
In3 |
|
|
|
|
|
|
Out4 |
|
|
Out6 |
|
||
|
|
|
НА |
Пар |
|
||
|
|
In4 |
|
|
In3 Cцп0 |
||
|
|
Out5 |
|
|
|
||
|
|
|
|
|
Out7 |
|
|
|
× |
In5 |
Out6 |
|
|
|
|
υна |
|
|
|
|
|||
|
|
|
|
|
|
|
|
Рисунок 5. Компьютерная модель управляемых массообменных процессов КТС «АБСОРБЦИЯ–ДЕСОРБЦИЯ»
На рисунке 6 представлены результаты компьютерного моделирования управляемых массообменных процессов КТС «АБСОРБЦИЯ–ДЕСОРБЦИЯ» в случае подачи ступенчатого возмущения по скорости газа на входе абсорбера. При изменении скорости газа от номинальной до минимальной наблюдается плавное отклонение выходных величин от первоначальных значений и постепенный переход к заданным значениям по ЦК в газовой фазе на выходе аппарата.
вых |
3 |
|
|
|
вых |
|
|
|
|
|
Cцг, г/м |
|
|
|
|
Cцж, %вес. |
|
|
|
|
|
0,025 |
|
|
|
|
8 |
|
|
|
|
|
0,023 |
|
|
|
|
7,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
0,021 |
|
|
|
|
6,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,019 |
|
|
|
|
6 |
|
|
|
|
|
вых |
|
а |
|
|
вых |
|
|
б |
|
|
Cцп, %вес. |
|
|
|
Cцж, %вес. |
|
|
|
|
||
99 |
|
|
|
|
4 |
|
|
|
|
|
98 |
|
|
|
|
3,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
97 |
|
|
|
|
3 |
|
|
|
|
|
96 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
95 |
|
|
|
t, с ×10^4 |
2,5 |
|
|
|
|
t, с ×10^4 |
|
|
|
|
|
|
|
|
|||
940 |
|
|
|
2 |
|
|
|
|
||
5 |
в |
10 |
15 |
0 |
5 |
г |
10 |
15 |
||
|
|
|
|
|
|
|
|
|
||
Рисунок 6. Графики переходных процессов в абсорбере и РК по газовой/паровой (а,в) и жидкой (б,г) фазам при ступенчатом снижении скорости газа
15
Следует отметить, что значения, полученные при вычислительном эксперименте, для содержания ЦК в жидкой фазе лежат в области значений, выявленных при наблюдении за ТП на реальном производственном объекте. Для регенерированного диэтиленгликоля эти значения лежат в области [0,71÷4,1] %вес, для насыщенного диэтиленгликоля – [2,45÷9,9] %вес. В случае ступенчатого снижения скорости газа содержание ЦК в регенерированном диэтиленгликоле установилось в новом значении 3,3 %вес., а в насыщенном – 7,1 %вес. В случае подачи возмущения в виде прямоугольных импульсов содержание ЦК в регенерированном диэтиленгликоле устанавливалось в значении 1,2 %вес., а в насыщенном – 4,9 %вес.
Сравнительный анализ результатов вычислительного и натурного экспериментов позволяет сделать выводы что, применение многорежимного регулирования обеспечивает поддержание параметров массобменных процессов КТС «АБ- СОРБЦИЯ–ДЕСОРБЦИЯ» в требуемом интервале. Разработанная математическая модель КТС «АБСОРБЦИЯ–ДЕСОРБЦИЯ» адекватна реальному технологическому процессу.
ЗАКЛЮЧЕНИЕ
Результаты диссертационного исследования состоят в следующем.
1. На основании анализа технологической, структурно-потоковой, укрупненной функциональной схем и физики процессов разработана концептуальная модель взаимосвязанных массотеплообменных процессов комплекса технологических систем «АБСОРБЦИЯ – ДЕСОРБЦИЯ». Установлены процессы в аппаратах комплекса, которые оказывают доминирующее влияние на абсорбционную осушку природного газа, что определило направления разработки и исследования математических моделей комплекса технологических систем «АБСОРБЦИЯ– ДЕСОРБЦИЯ».
2.В результате анализа состояния проблемы управления объектами с распределенными параметрами были определены нерешенные задачи управления пространственными массо- и теплообменными процессами абсорбционной осушки газа.
3.На основе анализа технологических процессов, происходящих в абсорбере, ректификационной колонне, испарителе и аппарате воздушного охлаждения, определены допущения и ограничения, составляющие базу для разрабатываемых математических моделей.
4.С учетом допущений и ограничений разработаны динамические математические модели подсистем «Абсорбция газа», «Ректификация», «Выпаривание», «Воздушное охлаждение», представленные системами дифференциальных уравнений в частных производных с начальными и граничными условиями.
5.На основе непрерывных моделей для исследования поведения подсистем численными методами были получены дискретно-непрерывные модели.
6.Для проведения компьютерного моделирования в программной среде
MATLAB/Simulink были разработаны:
− методика расчета физико-технологических коэффициентов абсорбции для номинального установившегося режима, используемых в математической модели подсистемы «Абсорбция газа»;