Просмотр кода страницы Ректификация

{{другие значения}}
'''Ректифика́ция''' (от {{lang-lat|rectus}} — прямой и ''facio'' — делаю) — процесс разделения двойных или многокомпонентных смесей за счёт противоточного массообмена между паром и жидкостью.
Ректификация — разделение жидких смесей на практически чистые компоненты, различающиеся температурами кипения, путём многократного частичного испарения жидкости и частичной конденсации паров.

== Применение ==
[[Файл:Тобольск-Нефтехим.jpg|thumb|300px|Ректификационные колонны разделения углеводородных фракций на [[Тобольск-Нефтехим]]]]
[[Файл:Fractional distillation lab apparatus-ru.svg|thumb|300px|Лабораторный аппарат для фракционной дистилляции]]
Ректификацию широко применяют в промышленности, например, для получения [[этанол]]а с отделением [[сивушные масла|сивушных масел]] и альдегидных фракций, для [[Дистилляты нефти|выделения]] [[бензин]]ов, [[керосин]]ов и других [[Фракционный состав нефти|фракций]] из [[нефть|нефти]], а также получения компонентов [[воздух]]а ([[кислород]]а, [[азот]]а, [[инертный газ|инертных газов]]).

Проводят ректификацию в башенных колонных аппаратах, снабжённых контактными устройствами (тарелками или насадкой) — [[ректификационная колонна|ректификационных колоннах]], в которых осуществляется многократный контакт между потоками паровой и жидкой фаз. Движущая сила ректификации — отличие фактических (рабочих) концентраций компонентов в паровой фазе от равновесных для данного состава жидкой фазы. Парожидкостная система стремится к достижению равновесного состояния. При контакте с жидкостью пар обогащается легколетучими (низкокипящими) компонентами — НК, а жидкость — труднолетучими (высококипящими) компонентами — ВК. Жидкость и пар движутся противотоком: пар — вверх, жидкость — вниз, поэтому при достаточно большой высоте колонны в её верхней части можно получить практически чистый целевой компонент. В зависимости от процесса, таковым компонентом может быть либо [[вещество]], либо [[азеотропная смесь]].

В зависимости от температур кипения разделяемых жидкостей ректификацию проводят под разным давлением: атмосферным — для кипящих при 30-150 °С, выше атмосферного — для жидкостей с низкими температурами кипения, например, сжиженных газов, в вакууме — для снижения температур кипения высококипящих. Ректификацию можно осуществлять непрерывно или периодически. Колонны для непрерывной ректификации состоят из двух ступеней: верхней — укрепляющей, где пар «укрепляется» — обогащается НК, и нижней — исчерпывающей, где жидкая смесь исчерпывается — из неё извлекаются НК и она обогащается ВК. При периодической ректификации в колонне производится только укрепление пара. Различают ректификацию бинарных (двухкомпонентных) и многокомпонентных смесей.

== Расчёт процесса ректификации ==

=== Ручной (графический) расчёт ===
[[Файл:Vapor-Liquid Equilibrium Diagram 2.jpg|frame|Равновесие пар — жидкость и рабочая линия процесса ректификации]]

Корректный ручной расчёт процесса ректификации возможен только для бинарных систем. Для этого по экспериментальным данным строят фазовую диаграмму пар — жидкость, на которой проводят рабочую линию процесса (на рисунке изображён периодический процесс) из точки на диагонали (линия жидкости), соответствующей составу дистиллята <math>x_d</math> до линии ординат в точке <math>B</math>,
: <math>B=\frac{x_p}{R+1}</math>, где <math>R</math> — флегмовое число, <math>R=\frac{F}{D}</math>, <math>F</math> — флегма, <math>D</math> — отбор дистиллята.
Из точки, соответствующей начальной концентрации перегоняемой смеси <math>x_p</math> (жидкость) строится вертикальный отрезок от рабочей линии до пересечения с линией пара и далее горизонтальный отрезок до пересечения с рабочей линией (это т. н. «ступень разделения» — такой метод расчёта и ввёл в оборот этот термин). Далее построение повторяется до точки <math>x_d</math>, флегмовое число и общее количество ступеней разделения являются минимально необходимыми для получения дистиллята заданного состава, ступени разделения предполагаются идеальными, с [[Коэффициент полезного действия|к.п.д.]] 100 % и бесконечно малой задержкой жидкости на тарелке. При этом состав жидкости в кубе постоянно меняется по мере отбора летучего компонента. Метод расчёта подробно изложен в<ref>''Касаткин Андрей Георгиевич''. Основные процессы и аппараты химической технологии. Учебник для вузов. 10-е издание, стереотипное, доработанное. Перепечатка с издания 1973 г. (Москва: Альянс, 2004), с.425</ref>, но в настоящее время считается устаревшим, так как не позволяет корректно рассчитать <math>x_p</math> и, соответственно, <math>x_d</math> через некоторое произвольное время <math>t</math>.

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

==== Материальный расчёт ====
Большинство методов основаны на модели [[теоретическая тарелка|теоретической тарелки]], из материального баланса которой следует: .
{{Нумерованная формула|:|<math>\frac{\mathrm{d}x_{i,\;k}}{\mathrm{d}t}=\frac{F(x_{i,\;k+1}-x_{i,\;k})+V(y_{i,\;k-1}-y_{i,\;k})}{L}</math>|1}}
Это уравнение тарелки ректифицирующей части колонны периодического действия для компонента <math>\mathbf{i}</math> на тарелке <math>\mathbf{k}</math><br>
Уравнения верхней тарелки (состав флегмы равен составу пара на тарелке)
{{Нумерованная формула|:|<math>\frac{\mathrm{d}x_{i,\;k}}{\mathrm{d}t}=\frac{F(y_{i,\;k}-x_{i,\;k})+V(y_{i,\;k-1}-y_{i,\;k})}{L}</math>|2}}
Уравнения куба — испарителя (<math>k=1</math>)
{{Нумерованная формула|:|<math>\frac{\mathrm{d}x_{i,\;k}}{\mathrm{d}t}=\frac{F(x_{i,\;k+1})-V(y_{i,\;k})}{LK-(V-F)t}</math>|3}}
Уравнения сборника:
{{Нумерованная формула|:|<math>\frac{\mathrm{d}x_{i,\;k}}{\mathrm{d}t}=\frac{(V-F)(y_{i,\;k-1}-x_{i,\;k})}{LD+(V-F)t}</math>|4}}
где

: <math>x_{i,\;k}</math> — концентрация компонента <math>\mathbf{i}</math> в жидкости на тарелке <math>\mathbf{k}</math> , [[Мольная доля|мольные доли]];
: <math>y_{i,\;k
}</math>, — концентрация компонента <math>\mathbf{i}</math> в паре на тарелке <math>\mathbf{k}</math>,м.д., считается;
: <math>F</math> — флегма, кол-во [[Моль (единица измерения)|молей]] в единицу времени, задаётся;
: <math>V</math> — пар, кол-во [[Моль (единица измерения)|молей]] в единицу времени, задаётся;
: <math>L</math> — задержка жидкости на тарелке, Моль, задаётся
: <math>LD</math> — начальное количество жидкости в сборнике, Моль, задаётся
: <math>LK</math> — начальное количество жидкости в кубе, Моль, задаётся
: <math>t</math> — общее время процесса (размерность произвольная, м.б. сек., мин., час, мес.). <math>F</math> и <math>V</math> задаются в тех же единицах.
{{начало скрытого блока |заголовок= Вывод уравнений '''(1 - 4)''' |стиль_тела= clear:both;width:65%; |выравнивание_заголовка = left |стиль_заголовка= background-color:#FFFAF0;  font-size:larger; }}
Материальный баланс тарелки предполагает соблюдение простого условия:<br>
Было компонента <math>i</math> на тарелке <math>k</math> + пришло за время <math>dt</math> = осталось + ушло. То есть для тарелки ректифицирующей части колонны
: <math> Lx_{i,\;k}+Vy_{i,\;k-1}dt+Fx_{i,\;k+1}dt = L(x_{i,\;k}+dx_{i,\;k})+V(y_{i,\;k}+dy_{i,\;k})dt+F(x_{i,\;k}+dx_{i,\;k})dt</math>
Уравнения верхней тарелки аналогичны с учётом того, что состав поступающей флегмы равен составу уходящего пара. <br>
Материальный баланс куба <math>(k=1)</math>: в момент времени <math>t</math> было компонента <math>i</math> + поступило за время <math>dt</math> :
: <math> (LK-(V-F)t)x_{i,\;k}+Fx_{i,\;k+1}dt</math>
Осталось в кубе через время <math>dt</math> + расход пара из куба за то же время
: <math> (LK-(V-F)(t+dt))(x_{i,\;k}+dx_{i,\;k})+(V(y_{i,\;k}+dy_{i,\;k})dt</math>
Материальный баланс сборника в момент времени <math>t</math>:
: <math> (LD+(V-F)t)x_{i,\;k}+(V-F)y_{i,\;k-1}dt = (LD+(V-F)(t+dt)(x_{i,\;k}+dx_{i,\;k})</math>

Из этих равенств (пренебрегая <math>\mathbf{dx\cdot dt}</math> и <math>\mathbf{dy\cdot dt}</math>) и выводятся уравнения ('''1''' — '''4''').<br>
Так же выводятся уравнения для колонны непрерывного действия, при этом уравнения укрепляющей части совпадают с ('''1'''), а уравнения исчерпывающей части включают ещё и питание колонны.
{{конец скрытого блока}}

Математическая модель ректификационной установки периодического действия представляет собой систему из <math>\mathbf{n\cdot m}</math> таких уравнений, где <math>\mathbf{n}</math> — число компонентов, <math>\mathbf{m}</math> — количество тарелок ректифицирующей части колонны. К ним добавляются <math>\mathbf{n}</math> уравнений куба-испарителя и столько же уравнений сборника, всего уравнений <math>\mathbf{n\cdot (m+2)}</math>.

Эта система интегрируется по времени с помощью программ стандартных библиотек, напр. [[Международная математическая библиотека подпрограмм|IMSL]]

В уравнениях '''(1 — 4)''' неизвестны концентрации компонентов в паре. Их корректный расчёт в общем случае невозможен, адекватные приближённые методы расчёта являются коммерческим секретом. Для смесей, незначительно отклоняющихся от закона Рауля расчёт равновесных концентраций компонентов в паре несложен. А для смесей, отклоняющихся от закона Рауля, в том числе азеотропных (напр. этанол — вода) такой расчёт может быть полезным для оценки динамики изменения концентраций примесей.<br>
Расчёт основан на [[Законы Рауля|законе Рауля]] и [[Закон Дальтона|законе Дальтона]].
: <math>y^*_{i,\;k}=\frac{P_{i,\;k}}{P_k}</math> где
: <math>y^*_{i,\;k}</math> равновесные концентрации компонента <math>\mathbf{i}</math> в парах на тарелке <math>\mathbf{k}</math>
: <math>P_{i,\;k}</math> — парциальное давление компонента <math>\mathbf{i}</math> на тарелке <math>\mathbf{k}</math>
: <math>P_k</math> — абсолютное давление на тарелке, <math>P_k=P_a + P_g + P_d</math>
: <math>P_a</math> — атмосферное давление;
: <math>P_g</math> — давление столба жидкости над тарелкой <math>\mathbf{k}</math>
: <math>P_d</math> — дополнительное давление, необходимое для преодоления [[Лобовое сопротивление|аэродинамического сопротивления]] массообменных устройств колонны.
В большинстве случаев последним слагаемым можно пренебречь без увеличения общей ошибки расчёта.<br>
Парциальное давление компонента <math>\mathbf{i}</math> на тарелке <math>\mathbf{k}</math>
: <math>P_{i,\;k}=P^0_{i,\;k}x_{i,\;k}</math> где
: <math>P^0_{i,\;k}</math> — давление чистого компонента <math>\mathbf{i}</math> при температуре на тарелке <math>\mathbf{k}</math>, считается по [[Уравнение Антуана|уравнению Антуана]]

: <math>\log_{10} P^0_{i,\;k} = A-\frac{B}{C+T_k}.</math>

Коэффициенты уравнения <math>A, B, C</math> для многих веществ приведены в<ref>{{Cite web |url=http://webbook.nist.gov/chemistry/ |title=NIST Chemistry WebBook |access-date=2019-04-03 |archive-date=2007-05-23 |archive-url=https://web.archive.org/web/20070523013553/http://webbook.nist.gov/chemistry/ |deadlink=no }}</ref>. Но в этом уравнении неизвестной является температура на тарелке <math>T_k</math>. Её находят из условия

: <math>\sum^{n}_{i=1} {P_{i,\;k}}={P_k}  </math> , то есть для тарелки <math>\mathbf{k}</math> считаются корни уравнения
: <math>f(T) = \sum^{n}_{i=1} {P_{i,\;k}} - {P_k} = 0 </math> ,то есть находят такую температуру, при которой сумма парциальных давлений компонентов на тарелке становится равной абсолютному давлению, после этого считают равновесные концентрации компонентов на тарелке и, с учётом к.п.д. тарелки — рабочие. Под к.п.д. тарелки здесь понимается отношение разности концентраций компонентов в жидкости и в паре к теоретически возможной:

: <math>\eta = \frac{y_{i,\;k}-x_{i,\;k}}{y^*_{i,\;k}-x_{i,\;k}} </math>

Таким образом становятся известны все правые части уравнений '''(1 — 4)''', для расчёта необходимо ещё задать начальные концентрации компонентов на тарелках.

Пример использования такого подхода для решения практической задачи приведён в статье<ref>{{Cite web |url=http://nasadkaspn.ru/teoreticheskie-osnovy/spiralno-prizmaticheskaja-nasadka/ |title=Влияние задержки жидкости на эффективность колонны |access-date=2019-12-10 |archive-date=2019-12-10 |archive-url=https://web.archive.org/web/20191210215313/http://nasadkaspn.ru/teoreticheskie-osnovy/spiralno-prizmaticheskaja-nasadka/ |deadlink=no }}</ref>.

==== Тепловой расчёт ====
Отсутствие учёта теплового баланса тарелки приводит к увеличению (часто значительному) ошибки расчёта. Преимуществом машинного расчёта (по сравнению с ручным) является наличие такой возможности. Тепловой баланс тарелки (как и материальный) основан на необходимом равенстве теплоты, поступившей на тарелку и расходуемой на ней.

Тепло поступает на тарелку (далее <math>Qp_k</math>)
* с более горячим паром с предыдущей тарелки;
* выделяется при конденсации паров, поступающих с предыдущей тарелки
<br>
Тепло расходуется на тарелке (далее <math>Qr_k</math>)
* на нагрев более холодной флегмы, поступающей сверху;
* на парообразование
* на тепловые потери (далее <math>Ql</math>)
Неравенство этих теплот приводит к изменению материального потока пара, так как других источников тепла на тарелке нет, при этом (без учёта тепловых потерь) тепловой баланс тарелки, за редким исключением, положителен, на тарелку поступает тепла больше, чем расходуется с неё по следующим причинам:
# при работе колонны с бесконечным флегмовым числом пар на тарелке обогащён более летучим компонентом по сравнению с паром, поступающем с тарелки ниже. Модель теоретической тарелки предполагает использование мольных концентраций, а мольная теплота парообразования летучих компонентов всегда меньше менее летучих (в гомологическом ряду, но это правило, за редким исключением, справедливо и для разнородных смесей).
# при работе колонны с конечным флегмовым числом (кроме п. 1) на тарелку поступает больше горячего пара, чем холодной флегмы, а разница температур на тарелках, в среднем по колонне, одинакова.

Единственная величина, которая может меняться на тарелке при ненулевом тепловом балансе — это поток пара. Если тепловой баланс тарелки положителен — на тарелке образуется дополнительное количество пара <math>\delta V_k</math>, если отрицателен — то <math>-\delta V_k</math>, при этом материальный баланс должен соблюдаться — изменение потока пара на тарелке изменяет и поток флегмы с неё (но это изменение на тепловой баланс тарелки не влияет).

Расчёт теплового баланса тарелок для всей колонны предполагает нахождение таких <math>\delta V_k</math>, при которых тепловые балансы всех тарелок будут нулевыми. Но тепловые балансы тарелок взаимосвязаны, их нельзя рассчитать по отдельности. Для расчёта составляется система уравнений вида <math>f_k\left(\delta V_k\right)=Qp_k - Qr_k = 0 </math>, <math>k=1, ..\mathbf{m}</math> и
: <math>Qp_k = (V+\delta V_{k-1})(Cp_{k-1}(t_{k-1}-t_k)+ H_{k-1})</math>
: <math>Qr_k = (F+\delta F_{k+1})(Cl_{k+1}(t_k-t_{k+1}))+(V+\delta V_k)H_k+Ql_k </math>
: <math> \delta F_k =  \delta F_{k+1} + \delta V_{k-1} - \delta V_k</math>
: <math>F=V-D</math>,
: <math>F</math> — флегма, поступающая в колонну из дефлегматора, моль в единицу времени;
: <math>V</math> — пар, поступающий в колонну из куба-испарителя, моль в единицу времени;
: <math>D</math> — отбор дистиллята, моль в единицу времени;
: <math>Cp_k, Cl_k</math> — теплоёмкость пара, жидкости на тарелке <math>k</math>, Дж*моль/град;
: <math>H_k</math> — теплота парообразования (конденсации) пара на тарелке <math>k</math>, Дж*моль;
: <math>t_k</math> — температура на тарелке <math>k</math>, град.
: <math>Ql_k</math> — тепловые потери на тарелке <math>k</math> в единицах мощности, Дж в единицу времени.
: <math> \delta F_k,\delta V_k</math> — изменение потока флегмы, пара на тарелке <math>k</math>, моль в единицу времени.

На верхнюю тарелку <math>(k=\mathbf{m})</math> поступает флегма такого же состава, что и уходящий пар и той же температуры, на тепловой баланс этот поток не влияет (разница температур равна нулю), величина потока флегмы из дефлегматора на верхнюю тарелку равна <math>F</math> и формально не зависит от <math>\delta V_m</math> (на верхнюю тарелку поступает жидкость в количестве, равном количеству уходящего пара за вычетом отбора <math>F+ \delta V_m = V + \delta V_m - D</math>), но всегда должно проверяться выполнение условия <math>V + \delta V_m > D</math> .<br>
На нижнюю тарелку <math>(k=1)</math> поступает пар в количестве <math>V</math> (<math>\delta V_0=0</math>). Тепловой баланс куба-испарителя не считается, так как предполагается, что куб получает тепло в количестве, необходимом для испарения <math>V</math> моль пара.

Тепловой расчёт производится после расчёта температуры на тарелках и концентраций компонентов в парах и до расчёта правых частей уравнений '''(1 — 4)''', которые требуют задания материальных потоков с учётом рассчитанных поправок (в уравнениях '''(1 — 4)''' не указаны).

Теплоёмкость и теплота парообразования смеси с достаточной для расчётов точностью [[Аддитивность (математика)|аддитивны]].

Система нелинейных уравнений <math>f_k</math> также решается с помощью стандартных программ.

Недостаток метода — большая [[вычислительная сложность]].

==== Начальные значения концентраций ====
При периодической ректификации перед началом отбора дистиллята рекомендуется работа колонны с бесконечным флегмовым числом для возможно большей концентрации летучих примесей. Так же поступают и в машинном расчёте: задаются произвольными начальными концентрациями жидкости на тарелках (обычно равными их концентрациям в кубе) при равных потоках флегмы и пара. Результаты вычислений (концентрации компонентов в колонне, рассчитанные для времени <math>t_1 <t_2 <t_3</math>) будут асимптотически приближаться к некоторому стационарному состоянию, но такой расчёт требует большого времени счёта и в нём нет необходимости — достаточно стабильных трёх значащих цифр.<br>
Так же нельзя задавать нулевое значение <math>LD</math>, в противном случае будет генерироваться {{iw|Ошибки времени исполнения|ошибка исполнения|en|Runtime error detection}}. Для исключения этой ошибки задаются некоторым ненулевым значением <math>LD</math> с произвольными концентрациями компонентов и по окончании расчёта вычитают первоначально заданное количество компонентов из рассчитанного.

==== Проверка корректности расчёта ====
Самый простой (но не самый быстрый) способ проверки корректности вычислений — упоминавшийся выше расчёт с бесконечным флегмовым числом (поток пара равен потоку флегмы). Если в программе есть ошибки, то результат расчёта концентраций компонентов в колонне для разных <math>t</math> при <math>t\to\infty</math> не будут постоянными.<br>
Кроме этого в любой момент времени <math>t</math> должны соблюдаться следующие условия:
* Материальный баланс колонны: общее количество каждого компонента, заданное изначально (в кубе, на тарелках и в сборнике) должно быть равным общему количеству каждого компонента в конце расчёта.
* Сумма концентраций компонентов в жидкости и в паре (в кубе и на каждой тарелке) должна быть равна единице.

Конечно, корректность модели не означает, что она будет точно соответствовать реальному процессу — допущения, принятые при расчёте равновесных концентраций компонентов в паре увеличивают ошибку расчёта. Так же на ошибку влияет (но в меньшей степени) и отсутствие учёта теплового баланса тарелки. Ещё одним источником ошибки являются завышенные ожидания эффективности массообменных устройств колонны (К.П.Д. тарелки — для тарельчатых колонн или насадки — для насадочных). Эти ошибки приводят к тому, что результаты расчёта будут завышены по сравнению с реальными (то есть модель ведёт себя как более эффективная реальная колонна с такими же параметрами процесса), что объяснимо — модель является идеальной абстракцией, даже с учётом К.П.Д тарелки или эффективности насадки (производители указывают эти характеристики как максимально возможные).

== См. также ==
* [[Дистилляция]]
* [[Фракционированная конденсация]]
* [[Зейгерование]]

== Примечания ==
{{примечания}}

== Литература ==
* Девятых Г. Г., Еллиев Ю. Е. Введение в теорию глубокой очистки веществ. — М.: Наука, 1981. — 320 с.

[[Категория:Методы разделения]]