Динамический материальный баланс линеаризуется путем замены нелинейного соотношения парожидкостного равновесия его локальной линейной формой с последующим приведением членов состава. Для ступени (n) баланс компонента в жидкой фазе имеет вид

[ H_n\frac{dx_n}{dt}

Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]

Если соотношение равновесия аппроксимируется локально как (y_i=m_i x_i) с постоянным наклоном (m_i), уравнение ступени превращается в линейное дифференциальное уравнение, связанное с соседними ступенями. Эта форма является основой для вывода передаточных функций, изучения распространения возмущений и проектирования систем управления составом для пилотных установок абсорбции и дистилляции.

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

Начиная с динамического баланса ступени

Накопление компонента в жидкости

При постоянной удерживающей способности жидкости (H_n) накопление абсорбируемого или дистиллируемого компонента на ступени (n) составляет

[ H_n\frac{dx_n}{dt}. ]

Здесь (x_n) — мольная доля компонента в жидкой фазе на ступени (n). Единицы измерения должны быть согласованы: если (H_n) — мольная удерживающая способность, то (L) и (V) должны быть мольными расходами.

Вклад потоков

Поток жидкости, входящий на ступень (n), вносит (Lx_{n-1}), а выходящий поток жидкости уносит (Lx_n).

Поток пара, входящий с соседней ступени, вносит (Vy_{n+1}), а выходящий со ступени (n) пар уносит (Vy_n). Таким образом,

[ H_n\frac{dx_n}{dt}

Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]

Этот баланс предполагает постоянство (L), постоянство (V), постоянную удерживающую способность, отсутствие реакций и пренебрежимо малое накопление в паровой фазе.

Применение аппроксимации локального равновесия

Линейное соотношение равновесия

Составы пара и жидкости связаны функцией равновесия

[ y=f(x). ]

Вокруг выбранной рабочей точки это соотношение можно аппроксимировать как

[ y_i \approx y_{i,0}+m_i(x_i-x_{i,0}), ]

где

[ m_i=\left.\frac{dy}{dx}\right|{x{i,0}}. ]

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

[ y_i=m x_i. ]

Наклон (m) может представлять локальную летучесть или соотношение Генри, в зависимости от системы.

Подстановка в баланс ступени

Используя (y_n=mx_n) и (y_{n+1}=mx_{n+1}), получаем

[ H_n\frac{dx_n}{dt}

Lx_{n-1}-Lx_n+mVx_{n+1}-mVx_n. ]

Сгруппировав члены, содержащие (x_n), получаем

[ \boxed{ H_n\frac{dx_n}{dt}

Lx_{n-1}

(L+mV)x_n + mVx_{n+1} } ]

или, после деления на удерживающую способность,

[ \boxed{ \frac{dx_n}{dt}

\frac{L}{H_n}x_{n-1}

\frac{L+mV}{H_n}x_n + \frac{mV}{H_n}x_{n+1} }. ]

Это стандартное линеаризованное динамическое уравнение для идеальной равновесной ступени при допущении о постоянстве расходов.

Что означает «линеаризованный» в этой модели

Линеаризация соотношения равновесия

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

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

Линеаризация материального баланса

Сам материальный баланс уже является линейным по (x), когда (H), (L) и (V) постоянны. Дополнительные нелинейности возникают, когда эти величины меняются в зависимости от состава, температуры, давления, уровня или общего расхода.

Для более общей модели записывают

[ H_n(x_n)\frac{dx_n}{dt}

L(x_{n-1}-x_n) + V(y_{n+1}-y_n), ]

затем разлагают модель вокруг установившегося состояния:

[ x_i=x_{i,0}+\delta x_i, \qquad y_i=y_{i,0}+\delta y_i. ]

Возмущение равновесия аппроксимируется как

[ \delta y_i=m_i\delta x_i. ]

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

Важное уточнение по индексации

Необходимо сохранить состав соседнего пара

Частая ошибка индексации — подстановка (y_{n+1}=mx_{n-1}). При обычной противоточной нумерации ступеней пар, входящий на ступень (n), — это (y_{n+1}), который связан с составом жидкости на ступени (n+1):

[ y_{n+1}=mx_{n+1}. ]

Следовательно, правильно проиндексированное уравнение имеет вид

[ H_n\frac{dx_n}{dt}

Lx_{n-1}

(L+mV)x_n + mVx_{n+1}. ]

Выражение

[ H\frac{dx_n}{dt}

Lx_{n-1}-(L+mV)x_n+mVx_{n-1} ]

верно только в том случае, если переменная, обозначенная как (x_{n-1}), была намеренно определена как эквивалентный состав жидкости выше по потоку для входящего пара. Это не является обычной подстановкой равновесия с индексацией по ступеням.

Почему это различие важно

Сохранение (x_{n+1}) создает связанную модель колонны. Состав на одной ступени реагирует на обе соседние ступени, что необходимо для представления распространения возмущений через тарельчатую колонну или дискретизированную насадочную колонну.

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

От уравнений ступеней к передаточным функциям

Форма в области Лапласа

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

[ H_n sX_n(s)

LX_{n-1}(s)

(L+mV)X_n(s) + mVX_{n+1}(s). ]

Перегруппировав, получаем

[ \boxed{ \left[H_n s+(L+mV)\right]X_n(s)

LX_{n-1}(s)+mVX_{n+1}(s) }. ]

Это уравнение представляет собой одну строку динамической передаточной функции колонны или модели в пространстве состояний.

Одноступенчатая аппроксимация

Если состав входящего пара рассматривается как известное входное воздействие, а не как состояние, связанное с (x_{n+1}), определим

[ u_n=Lx_{n-1}+Vy_{n+1}. ]

Баланс принимает вид

[ H_n\frac{dx_n}{dt}

u_n-(L+mV)x_n. ]

Соответствующая передаточная функция

[ \boxed{ \frac{X_n(s)}{U_n(s)}

\frac{1}{H_n s+L+mV} } ]

с постоянной времени

[ \boxed{ \tau_n=\frac{H_n}{L+mV} }. ]

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

[ \frac{X_n(s)}{X_{n-1}(s)}

\frac{L}{H_n s+L+mV}. ]

Применение к абсорбции и дистилляции

Газожидкостная абсорбция

При абсорбции растворенное вещество обычно переходит из газовой фазы в жидкую, когда фактический состав газа превышает его равновесное значение:

[ y>y^*(x). ]

Наклон равновесия (m) описывает, как равновесный состав газа изменяется в зависимости от состава жидкости. Большее значение (mV) увеличивает силу связи паровой фазы и изменяет реакцию ступени.

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

Дистилляция

При дистилляции соотношение равновесия обычно представляется как

[ y_i=y^*(x_i), ]

где (y^*(x)) может быть нелинейным. Вокруг рабочей точки состава

[ \delta y_i \approx \left.\frac{dy^*}{dx}\right|{x{i,0}}\delta x_i. ]

Локальный наклон может меняться от ступени к ступени, поэтому дистилляционная колонна обычно имеет коэффициенты (m_n), (L_n), (V_n) и (H_n), а не одно универсальное (m), (L) и (V).

Граничные условия

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

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

Расширение модели на переменную удерживающую способность и расход

Общий материальный баланс

Когда уровень жидкости или удерживающая способность меняются, необходимо включить общий баланс. Общий баланс удерживающей способности жидкости имеет вид

[ \frac{dH_n}{dt}

L_{n-1}-L_n+V_{n+1}-V_n, ]

со знаками и определениями потоков, скорректированными под выбранную конвенцию ступеней.

Баланс компонента должен быть записан через удерживающую способность компонента:

[ \frac{d(H_nx_n)}{dt}

\text{компонент входящий}

\text{компонент выходящий}. ]

Раскрытие производной дает

[ H_n\frac{dx_n}{dt} + x_n\frac{dH_n}{dt}

\text{компонент входящий}

\text{компонент выходящий}. ]

Дополнительный член (x_n,dH_n/dt) важен, когда уровни жидкости на тарелках или удерживающая способность кипятильника и конденсатора заметно меняются.

Энергетические и барические эффекты

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

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

Понимание компромиссов

Преимущества линеаризованной модели

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

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

Ограничения постоянного наклона равновесия

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

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

Ограничения постоянного расхода и удерживающей способности

Допущение о постоянстве (L), (V) и (H) исключает динамику уровней и гидравлические эффекты. Это может привести к оптимистичным прогнозам скорости реакции и скрыть важные взаимодействия между управлением расходом, давлением и составом.

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

Путаница рабочих линий с динамическими балансами

Рабочая линия описывает соотношение установившегося материального баланса между составами фаз. Она сама по себе не описывает, как составы развиваются во времени.

Динамическая модель требует членов удерживающей способности, таких как (H,dx/dt). Удаление этих членов превращает модель в установившееся соотношение, а не в переходное.

Правильный выбор для вашей цели

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

Если ваша основная цель — проектирование контроллера: Линеаризуйте (y^*(x)) вокруг намеченной рабочей точки, сохраните связь соседних ступеней и соберите уравнения в форме пространства состояний.
Если ваша основная цель — простая реакция одной ступени: Рассматривайте состав входящего пара как известное входное воздействие и используйте передаточную функцию первого порядка (1/[Hs+(L+mV)]).
Если ваша основная цель — распространение возмущений по всей колонне: Сохраняйте (x_{n-1}) и (x_{n+1}) в каждом уравнении ступени и применяйте правильные граничные условия для конденсатора, кипятильника, питания, флегмы и продукта.
Если ваша основная цель — крупные переходные процессы на пилотной установке: Используйте нелинейное равновесие, переменную удерживающую способность, меняющиеся расходы потоков и связанные балансы энергии и давления.

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

Сводная таблица:

Ключевой элемент	Инсайт линеаризованной модели
Динамический баланс	`H_n dx_n/dt = Lx_{n-1} - Lx_n + Vy_{n+1} - Vy_n`
Аппроксимация равновесия	Замена `y=f(x)` на `delta y=m delta x` вблизи установившейся рабочей точки
Связанное уравнение ступени	`H_n dx_n/dt = Lx_{n-1} - (L+mV)x_n + mVx_{n+1}`
Основные допущения	Постоянная удерживающая способность, расходы жидкости и пара, пренебрежимо малое накопление в паровой фазе
Форма передаточной функции	`[H_n s + (L+mV)]X_n = LX_{n-1} + mVX_{n+1}`
Постоянная времени ступени	`tau_n = H_n/(L+mV)`
Критическая точка индексации	Член входящего пара связан с `x_{n+1}`, сохраняя связь соседних ступеней
Лучшие области применения	Настройка контроллеров, анализ возмущений, передаточные функции и моделирование в пространстве состояний
Когда нужны нелинейные модели	Большие возмущения, переменная удерживающая способность или расход, сильная неидеальность, изменения давления или энергетические эффекты

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