The methods of calculating the efficiency of fracking column plates are considered. The efficiency of cap plates was calculated when separating the methanol-water mixture according to two models based on the use of the number of transfer units and mass transfer coefficients
rectification, poppet columns, efficiency of plates, transfer units, mass transfer coefficients
Ректификация является одним из важнейших способов разделения жидких смесей на отдельные фракции или чистые компоненты. С помощью её в химической, нефтехимической и смежных отраслях промышленности перерабатывают сотни миллионов тонн разнообразной продукции.
Процесс ректификации основан на взаимодействии пара и жидкости, движущихся противотоком. При этом пар обогащается белее летучими компонентами, а жидкость – менее летучими.
Наиболее распространёнными аппаратами для осуществления процессов ректификации являются тарельчатые колонны, которые допускают большие нагрузки по пару и жидкости, обеспечивают возможность отбора продуктов по высоте колонны и обладают рядом других достоинств. Особенность их устройства позволяет применять ступенчатый метод расчета, известный в литературе как расчет “от тарелки к тарелке” [1].
В качестве основного элемента выступает одиночная контактная тарелка, математическое описание которой включает систему балансовых, термодинамических и кинетических уравнений.
В процессе взаимодействия пара и жидкости происходит непрерывный переход компонентов смеси из одной фазы в другую. При продолжительном контакте фаз и постоянных условиях между паром и жидкостью устанавливается состояние равновесия, характеризующееся одинаковой скоростью перехода каждого компонента из фазы в фазу. При этом данной концентрации компоненты в жидкости х соответствует строго определенная равновесная концентрация его в паре у*. Известно, что системы самопроизвольно стремятся к равновесию и знание равновесных концентраций является необходимым условием расчета процесса ректификации.
В настоящее время накоплен большой объем экспериментальных данных по равновесию между жидкостью и паром и разработаны математические модели, позволяющие рассчитывать фазовое равновесие для бинарных и многокомпонентных смесей [2-4].
Сведения о равновесии между жидкостью и паром, о составе исходной смеси и продуктов позволяют выполнить термодинамический расчет колонн, то есть определить необходимое число тарелок, на которых достигается равновесие, так называемых теоретических тарелок nTT, необходимых для заданного разделения смеси.
На рисунке 1 показана тарелка с жидкостью, в которой концентрация летучего компонента равна х1. Снизу поступает пар с концентрацией этого компонента у1 и выходит из слоя обогащенным. В случае теоретической тарелки было бы достигнуто равновесие и концентрация пара на выходе составила у*.
а) б)
Рисунок 1 – Схема контактной тарелки с идеальным перемешиванием
жидкости (а) и фрагмент диаграммы у-х бинарной смеси
На практике в силу разных причин, одна из которых – недостаточное время контакта пара и жидкости, равновесное состояние системы не достигается и пар покидает тарелку с концентрацией у2 меньше чем у*. Отношение изменения концентрации пара, которое произошло на тарелке, к тому изменению, которое могло бы произойти при достижении равновесия, называют эффективностью тарелки Е.
Разумеется, что для достижения требуемого разделения смеси реальных тарелок потребуется больше чем теоретических. При этом существует простое соотношение теоретических nТТ и реальных nДТ тарелок, связанное с их эффективностью:
Следует отметить, что термодинамический расчет процесса ректификации, связанный с определением числа теоретических тарелок, при наличии данных по фазовому равновесию, не представляет больших затруднений. Основная трудность состоит в правильной оценке эффективности тарелки Е.
Эффективность тарелок, являясь интегральным кинетическим параметром процесса, зависит от множества факторов. На неё оказывают влияние физико-химические свойства пара и жидкости такие как плотность, вязкость, коэффициенты молекулярной диффузии. Очень большое значение имеет гидродинамическая обстановка, возникающая на тарелках, которая зависит от интенсивности потоков и конструктивных особенностей контактных устройств.
Попытки получить удовлетворительные корреляции эффективности от определяющих факторов [5, 6] не увенчались успехом. Одной из причин является стремление учесть единым уравнением две стадии процесса – перенос компонента из ядра жидкой фазы к границе раздела и перенос его от границы раздела в паровую фазу (рисунок 2).
Рисунок 2 – Схема переноса летучего компонента из жидкости в пар
Более плодотворным оказался подход к описанию процесса, связанный с раздельным анализом обстановки в каждой фазе. Количество компонента m, переносимое в каждой фазе описывается уравнениями массоотдачи:
, кмоль/с (1)
, кмоль/с (2)
где – коэффициенты массоотдачи в жидкой и паровой фазах, соответственно, м/c;
– плотности жидкой и паровой фаз, кмоль/м3; х, у, хГР, уГР – мольные доли компонента в жидкости, паре и на границе раздела фаз, соответственно; А – величина поверхности контакта фаз на 1 м2 активной площади тарелки, м2.
Выражения, стоящие в скобках, представляют собой движущие силы процессов переноса компонента в жидкой и паровой фазах.
Такое же количество компонента m примет паровая и потеряет жидкая фаза:
, кмоль/с (3)
, кмоль/с (4)
где V и – объемные расходы пара и жидкости, отнесенные к 1 м2 активной площади тарелки, м/с, соответственно.
Решая попарно уравнения массоотдачи и материального баланса получим соотношения для паровой и жидкой фаз:
(5)
(6)
Выражения, стоящие в левой части уравнений, показывают во сколько раз изменения концентрации потока больше движущих сил процесса и называются числами единиц переноса. Они включают в себя коэффициенты массоотдачи и, следовательно, характеризуют скорость переноса компонента в каждой фазе. Важно, что каждый из них зависит от свойств и режима движения лишь одной фазы и, следовательно, легче поддается описанию. Поэтому наибольшее число работ связано с опытным определением этих кинетических параметров и обобщением экспериментальных данных в виде расчетных зависимостей [7-12].
Рассчитав значения частных чисел NГ и NЖ, определяют общее число единиц переноса NОГ:
(7)
Фактор отклонения λ равный произведению тангенса угла наклона линии равновесия на соотношение потоков пара и жидкости
(8)
обеспечивает сопряжение чисел единиц переноса.
Зная общее число единиц переноса, можно определить эффективность тарелки*:
(9)
Использование чисел единиц переноса в расчетах способствовало получению более общих зависимостей для расчета эффективности тарелок. Однако точность и надежность их, вплоть до сегодняшних дней, остаются проблематич-
* Выражение (9) строго справедливо при полном перемешивании жидкости. В противном случае вводятся поправки, учитывающие степень перемешивания жидкости на тарелке.
ными. Сравнение результатов расчета эффективности ситчатых тарелок при ректификации смеси “изопропанол-вода” по шести имеющимся зависимостям для чисел единиц переноса показало большие расхождения [13]. Одной из причин этого является сложная структура чисел единиц переноса, которые включают в себя наряду с коэффициентами массоотдачи β величину межфазной поверхности.
Раздельный учет коэффициентов массоотдачи и поверхности контакта фаз представлялся логическим шагом в детализации описания процесса. Коэффициенты, отнесенные к действительной поверхности, определяют скорость переноса компонентов смеси. Именно для них устанавливается связь с молекулярными транспортными характеристиками и состоянием межфазной поверхности. К тому же, значение величины межфазной поверхности совместно с величиной газосодержания определяет средний поверхностно-объемный диаметр пузырьков пара, который можно считать наиболее характерным размером двухфазной системы. Его важно знать при обобщении результатов экспериментов методами теории подобия. По-видимому, такими соображениями руководствовались исследователи в шестидесятых годах прошлого века пытаясь измерить величину межфазной поверхности в газо- и парожидкостных слоях.
На тарелках ректификационных колонн поверхность контакта образована пузырьками, струями пара, каплями жидкости, которые находятся в непрерывном движении и трудно поддается определению. Тем не менее были разработаны методы её измерения.
Первым был использован метод фотографирования двухфазного слоя с последующей статистической обработкой фотографии [14]. В результате обработки фотографий определялся средний диаметр газовых образований dП и газосодержание слоя . Это позволяло определить удельную объемную поверхность контакта фаз а, м2/м3:
(10)
Умножив а на высоту двухфазного слоя hf, получим величину поверхности контакта А, м2/м2, отнесенную к 1 м2 активной площади тарелок.
В работе [15] величину межфазной поверхности определяли методом отражения светового потока от поверхности двухфазного слоя. Оба указанных метода позволяют оценить величину поверхности контакта пристенного слоя, в то время как внутренняя структура оказывается недоступной.
Химические методы измерения дают интегральное значение величины межфазной поверхности. Однако использование их ограничивается системами с необходимыми кинетическими и физико-химическими характеристиками. Для измерения поверхности наиболее широко применяют реакцию карбонизации щелочи в процессе поглощения углекислого газа из смеси его с воздухом и реакцию окисления сульфата натрия кислородом воздуха в присутствии катализатора [16-18].
Для измерения межфазной поверхности применяют также электроконтактные методы, фиксирующие замыкание и размыкание цепи при прохождении пузырьков газа через игольчатые электроды [19-20].
Нами [21] разработан метод измерения поверхности путем просвечивания двухфазного слоя монохроматическим поляризованным светом. Сравнение результатов измерения поверхности этим методом с результатами, полученными при поглощении СО2 растворами NaOH показали достаточно хорошую сходимость.
Знание величины поверхности контакта позволило определить коэффициенты массоотдачи βГ и βЖ и получить для них расчетные зависимости [22-25].
В производстве метиламинов реакционная смесь после стадии синтеза подвергается ректификации последовательно в пяти колоннах. Последняя из них предназначена для извлечения метилового спирта и примесей метиламинов из воды. Кубовая жидкость колонны с концентрацией метилового спирта не более 250 мг/дм3 поступает на биологическую очистку.
В настоящее время по рекомендации АнГТУ ведется проектная работа по замене желобчатых тарелок на колпачковые, и является актуальной оценка их эффективности.
В таблице 1 представлены основные параметры колпачковых тарелок. Для расчета нами были выбраны две характерные тарелки в нижней и верхней частях колонны. Гидравлические показатели их работы при номинальной нагрузке колонны по питанию представлены в таблице 2, а в таблице 3 приведены физико-химические свойства пара и жидкости.
Таблица 1
Основные параметры колпачковой тарелки диаметром D = 1,2 м
|
№ п/п |
Наименование параметра, его обозначение и размерность |
Значение |
|
1 |
Площадь сечения колонны SК, м 2 |
1,13 |
|
2 |
Площадь активной зоны тарелки, SАКТ, м2 |
0,67 |
|
3 |
Периметр слива, lW, м |
0,825 |
|
4 |
Доля свободного сечения тарелки, fС |
0,063 |
|
5 |
Высота сливной перегородки, hW, м |
0,06 |
В расчетах эффективности были использованы две модели. Одна из них базируется на числах единиц переноса в паровой и жидкой фазах [11], а другая учитывает величину поверхности контакта фаз и коэффициенты массоотдачи [22-25] (таблица 4). Результаты расчета эффективности тарелок по этим моделям представлены в таблице 5.
Таблица 2
Гидравлические показатели работы колпачковых тарелок при номинальной нагрузке
|
№ п/п |
Наименование параметра, его обозначение и размерность |
Номера тарелок (считая снизу) |
|
|
10 |
29 |
||
|
1 |
Скорость пара в рабочем сечении тарелки, wP, м/с |
0,58 |
0,51 |
|
2 |
Скорость пара в сечении колонны, wK, м/с |
0,5 |
0,44 |
|
3 |
Запас жидкости на тарелке, h0, м |
0,034 |
0,042 |
|
4 |
Газосодержание слоя φГ |
0,52 |
0,53 |
|
5 |
Высота двухфазного слоя, hf, м |
0,07 |
0,09 |
|
6 |
Гидравлическое сопротивление слоя жидкости, ΔРСТ, Па |
304 |
277 |
|
7 |
Объемные расходы жидкости, VЖ, м3/с |
2,5∙10-3 |
7,4∙10-4 |
|
8 |
Расход жидкости на 1 м2 активной площади, U, м3/(м2∙с) |
3,75∙10-3 |
1,12∙10-3 |
|
9 |
Расход пара, VП, м3/с |
0,44 |
0,38 |
|
10 |
Молярный расход жидкости, L, кмоль/ч |
457 |
72,3 |
|
11 |
Молярный расход пара, G, кмоль/ч |
71 |
59 |
Таблица 3
Физико-химические свойства системы и термодинамические параметры
|
№ п/п |
Наименование параметра, его обозначение и размерность |
Номера тарелок (считая снизу) |
|
|
10 |
29 |
||
|
1 |
Плотность жидкости, ρЖ, кг/м3 |
910 |
673 |
|
2 |
Плотность пара, ρП, кг/м3 |
0,8 |
1,38 |
|
3 |
Вязкость жидкости, μЖ, Па∙с |
0,236∙10-3 |
0,276∙10-3 |
|
4 |
Вязкость пара, μП, Па∙с |
9,4∙10-6 |
7∙10-6 |
|
5 |
Поверхностное натяжение, σ, Н/м |
0,057 |
0,025 |
|
6 |
Коэффициент молекулярной диффузии в жидкости, DЖ, м2/с |
6,46∙10-8 |
3,8∙10-8 |
|
7 |
Коэффициент молекулярной диффузии в паре, DП, м2/с |
1,64∙10-5 |
1,64∙10-5 |
|
8 |
Тангенс угла наклона линии равновесия, |
7 |
1,5 |
|
9 |
Температура жидкости на тарелке, t, °С |
109 |
70,6 |
|
10 |
Разность температур пара и жидкости, Δt, °С |
0,2 |
1 |
|
11 |
Теплоемкость жидкости, сЖ, Дж/(кг∙град) |
4100 |
3300 |
|
12 |
Удельная теплота испарения, r, Дж/кг |
2 200 000 |
1 200 000 |
Таблица 4
Модели расчёта эффективности колпачковых тарелок на основе чисел единиц переноса (1) и коэффициентов массоотдачи (2)
|
№ п/п |
Уравнения |
|
1 |
|
|
2 |
|
Таблица 5
Результаты расчета кинетических параметров процесса массообмена между паром и жидкостью на колпачковых тарелках (система “метанол-вода”)
|
Номера уравнений в таблице 4 |
Значения параметров |
|||||||
|
А |
паровая фаза |
жидкая фаза |
Е |
|||||
|
ShГ |
βГ |
NГ |
ShЖ |
βЖ |
NЖ |
|||
|
Тарелка № 10 |
||||||||
|
1 |
- |
- |
- |
1,57 |
- |
- |
12,7 |
0,74 |
|
2 |
26 |
19,1 |
0,037 |
1,92 |
162 |
1,25∙10-3 |
8,6 |
0,78 |
|
Тарелка № 29 |
||||||||
|
1 |
- |
- |
- |
1,86 |
- |
- |
25,2 |
0,83 |
|
2 |
26 |
19,8 |
0,038 |
1,94 |
122 |
5,4∙10-4 |
12,5 |
0,80 |
Из таблицы следует, что расчеты по обеим моделям дали близкие значения эффективности. Это указывает на надежность результатов. Однако окончательный ответ будет получен после пуска колонны в эксплуатацию и анализа её разделяющего действия.
1. Kasatkin, A.G. Osnovnye processy i apparaty himicheskoy tehnologii. – M., 1971. – 784 s.
2. Kogan, V.B. Ravnovesie mezhdu parom i zhidkost'yu / V.B. Kogan, V.M. Fridman, V.V. Kafarov. – M.-L.: Nauka, 1966. – 1426 s.
3. Kogan, V.B. Geterogennye ravnovesiya. – L.: Himiya, 1968. – 432 s.
4. Ueylis S. Fazovye ravnovesiya v himicheskoy tehnologii. – M.: Mir, 1989. – 662 s.
5. English, G.E. Efficiency of fractionating columns / G.E. English, M. Van Winkle // Jour. Chem. Eng. – 1963. – Vol. 70. – № 11. – P. 241.
6. Onda, K. Plate and columns efficiencies of continuous rectifying col-umns for binary mixtures / K. Onda, E. Soda, K. Takahashi, S.A. Mukhtar // Jour. A. I. Chem. Eng. – 1971. – Vol. 17. – № 3. – P. 1141.
7. Asano, O. Gas-liquid contacting area on bubbling trays / O. Asano, S. Fu-jita // Jour. Chem. Eng. Japan. – 1965. – Vol. 3. – № 1. – P. 85.
8. Ivanov, V.A. Issledovanie massoperedachi pri rektifikacii sme-sey etanol(metanol)-voda v kolonne s sitchatymi tarelkami / V.A. Ivanov, A.N. Planovskiy, A.A. Baraev // Him. prom. – 1967. – № 5. – S. 385.
9. Gerster, J.A. Tray efficiencies distillation columns: Final report from the university of Delaware // J.A. Gerster, A.B. Hill, N. Mochgraf // Am. Inst. Chem. Engrs. – 1958. – P. 118.
10. Weiss, S. Zum Stoffübertragung bei der Distillation in Bodenkolonnen // S. Weis, L. John // Chem. Techn. – 1973. – Bd. 25. – № 1. – S. 658.
11. Kasatkin, A.G. Massoperedacha na barbotazhnyh proval'nyh tarelkah / A.G. Kasatkin, D.M. Popov, Yu.I. Dytnerskiy // Him. prom. – 1962. – № 2. – C. 123.
12. Solomaha, G.P. Uravnenie massootdachi v gazovoy faze na kolpachkovyh tarelkah / G.P. Solomaha // TOHT. – 1970. – T. 4. – № 2. – C. 181.
13. Hübner, W. Experimentelle Untersuchungen zum Verstärkungsver-hälthis von Rektifizierböden / W. Hübner, E.U. Schlünder // Chem. Ing. Techn. – 1973. – Bd. 45. –№ 5. – S. 247.
14. Rodionov, A.I. O dispersnosti pennogo sloya i metode opredeleniya srednego razmera puzyr'ka / A.I. Rodionov, V.M. Radikovskiy. – V kn.: Massoobmennye processy himicheskoy tehnologii. – L, 1965. – C. 32.
15. Calderbank, P.H. Physical rate processes in industrial fermentation the interfacial area in gas-liquid contacting with agitation / P.H. Calderbank // Trans/Inst. Chem. Engrs. – 1958. – Vol. 36. – № 6. – P. 443.
16. Rodionov, A.I. Issledovanie himicheskim metodom poverhnosti kontakta faz na sitchatyh tarelkah / A.I. Rodionov, A.A. Vinter // Izv. vuzov. Ser. Himiya i him. tehnologiya. – 1966. – T. 9. – № 6. – s. 970; – 1967. – T. 10. – № 1. – s. 102.
17. Yoshida, F. Effective interfacial area in packed columns for absorption with chemical reaction / F. Yoshida, Y. Miura // Jour. A. I. Ch. E. – 1963. – Vol. 9. – № 3. – P. 331.
18. Porter, K.E. Interfacial areas and liquid-film mass-transfer coefficients for a 3 ft. diameter bubble-cap plate derived from absorption rates of carbon dioxide into water and caustic soda solutions / K.E. Porter, M.B. King, K.C. Varshney // Trans. Inst. Chem. Engrs. – 1966. – Vol. 44. – № 7. – P. 274.
19. Tarat, E.Ya. O vzaimodeystvii gazovyh puzyr'kov s datchikami v gazozhidkostnyh strukturah / E.Ya. Tarat, O.S. Kovalev, V.I. Schepin. – V kn.: Gidrodinamika i yavleniya perenosa v dvuhfaznyh dispersnyh sistemah. – Irkutsk, 1976. – C. 51.
20. Sabanin, V.A. Stereometricheskiy metod izmereniya udel'noy poverhnosti kontakta faz v dvuhfaznoy sisteme “zhidkost'-gaz” koleblyuschimsya datchikom / V.A. Sabanin. – V kn.: Gidrodinamika i yavleniya perenosa v dvuh-faznyh dispersnyh sistemah. – Irkutsk, 1976. – C. 74.
21. Avtorskoe svidetel'stvo № 213409 SSSR, MPK G01F23/292. Sposob opredeleniya poverhnosti kontakta faz : № 1126771/26-25 : zayavl. 17.01.1967 : opubl. 12.03.1968 : byulleten' № 10 / Rodionov A.I., Ul'yanov B.A., Kashnikov A.M., Stroganov E.F. ; zayavitel' MHTI im. D.I. Mendeleeva. – 2 s.
22. Vinter, A.A. Opredelenie poverhnosti kontakta faz na barbotazhnyh tarelkah: Dis. na soisk. uchen. step. kand. tehn. nauk / Vinter Al'bert Aleksandrovich ; MHTI im. D.I. Mendeleeva. – M, 1966. – 130 s.
23. Kochetov, N.M. Issledovanie gidrodinamiki i massoperedachi na kolpachkovyh i kolpachkovo-sitchatyh tarelkah: Dissertaciya na soiskanie uchenoy stepeni kandidata tehnicheskih nauk / Kochetov Nikolay Mihaylovich ; MHTI im. D.I. Mendeleeva. – M, 1972. – 116 s.
24. Ul'yanov, B.A. Poverhnost' kontakta faz i massoobmen v tarel'chatyh rektifikacionnyh apparatah / B.A. Ul'yanov. – Irkutsk : Izd. IGU, 1982. – 129 s.
25. Fereferov, M.Yu. Massoobmen v zhidkoy faze na tarelkah rektifikacionnyh kolonn : avtoreferat dissertacii na soiskanie uchenoy stepeni kandidata tehnicheskih nauk / Fereforov Mihail Yur'evich ; Angarskiy gosudarstvennyy tehnologicheskiy institut. – Angarsk, 1997. – 18 s.
