Модель UNIFAC [44-48]

4.8.1. Модель UNIFAC [44-48]

Метод UNIFAC (universal functional activity coefficient – универсальный функциональный коэффициент активности) был разработан в 1975 году Фреденслундом, Джонсом и Праузницем. Данный метод рассчитывает величины коэффициентов активности на основе концепции группового вклада. Предполагается, что взаимодействия между двумя молекулами являются функцией взаимодействий между группами. Данные межгруппового взаимодействия получаются путем обработки экспериментальных данных для пар компонентов. Число функциональных групп ограничено.

Метод UNIFAC основан на модели UNIQUAC, которая представляет избыточную энергию Гиббса (и логарифм коэффициента активности) как комбинацию двух эффектов. Таким образом, используется уравнение:

(4.8.1.)

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

(4.8.2.)

где:

(4.8.3.)

N_C – число компонентов;

N_G – число различных групп в смеси;

z – согласованное число для пространственной решетки, равное 10;

 – число функциональных групп типа k в молекуле i;

R_k – параметр объема для функциональной группы k;

Q_k – параметр площади для функциональной группы k;

x_i – мольная доля компонента i в жидкой фазе.

Параметры группового объема и площади определяются из атомной и молекулярной структуры:

(4.8.4.)

где : V_k, A_k = объем и площадь Ван дер Ваальса для группы k;

Остаточный член  определяется как:

(4.8.5.)

где: Г_k = остаточный коэффициент активности группы k в смеси;

 = остаточный коэффициент активности группы k в стандартном растворе, содержащем только молекулы, состоящие из групп типа i.

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

(4.8.6.)

где : m, n = 1, 2, ..., N_G.

Параметр τ_mn определяется как

(4.8.7.)

где : a_mn – параметр бинарного взаимодействия для групп m и n.

Принимается, что параметр бинарного взаимодействия a_mn постоянен и не зависит от температуры.

5. Постановка задачи

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

Для достижения поставленной цели необходимо:

•  выбрать математическую модель парожидкостного равновесия, адекватно описывающую экспериментальные данные;

•  синтезировать схемы разделения;

•  произвести расчет ректификации и выбрать решение, обеспечивающее минимальные энергозатраты.

6. Расчетная часть.

 

6.1. Методы и алгоритмы исследования

В настоящей работе для математического моделирования схем ректификации использован лицензионный программный комплекс PRO/II with PROVISION компании SIMSCI corp., обеспечивающий моделирование и расчет технологических схем ректификации.

PRO/II - компьютерный комплекс для инженерных расчетов процессов органического синтеза и нефтехимии, технологии полимеров и др. Он объединяет базы данных химических компонентов и расширенных методов расчета термодинамических свойств с гибкими методами расчета аппаратов. Программа обладает вычислительными средствами для выполнения расчетов всех материальных и энергетических балансов необходимых для моделирования большинства статических процессов. Экспертные системы, расширенная обработка входных данных и проверка ошибок обеспечивают его высокую эффективность и надежность.

6.2. Расчет ректификации

Все алгоритмы ректификации в программе PRO/П представляют собой строгие модели равновесных ступеней контакта. В каждой модели решаются тепловой и материальный балансы и уравнения равновесия жидкость - пар.

Программа PRO/II предлагает четыре различных алгоритма моделирования ректификационных колонн:

алгоритм Inside/Out (I/O),

алгоритм Sure,

алгоритм Chemdist и

алгоритм ELDIST.

Алгоритм I/O может быть использован для решения большинства задач нефтепереработки и обладает высоким быстродействием. В настоящей работе расчет колонн ректификации проводился по этому алгоритму.

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

Рис. 6.2.1.. Схема простой ступени контакта алгоритма I/O.

Тепловой баланс для ступени контакта j выглядит следующим образом:

(6.2.1)

где: V— расход пара, покидающего ступень контакта;

где: V— расход пара, покидающего ступень контакта;

L - расход жидкости, покидающей ступень контакта;

Ls — жидкостной боковой погон;

V_s - паровой боковой погон;

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

где К представляет собой соотношение равновесной фугитивности пар - жидкость, может быть записан как:

(6.2.2.)

где: - расход жидкого компонента;

v - расход парового компонента;

f- расход компонента сырья.

6.3. Объект исследования

 

В качестве побочных продуктов при синтезе хлорбензола получаются дихлорбензолы и трихлорбензолы. Общее содержание побочных продуктов в потоке приблизительно равно 1,5% масс. Так как в этой работе не ставилась задача получения всех побочных продуктов в качестве отдельных продуктовых потоков в рассматриваемой нами смеси все побочные продукты были заменены на самый легкокипящий из них – метадихлорбензол. Далее смесь продуктов синтеза хлорбензола будет рассматриваться как трёхкомпонентная зеотропная смесь: бензол - хлорбензол - метадихлорбензол. В качестве технических требований к готовым продуктам приняты требования по ГОСТ 646-73 для хлорбензола и по ОСТ 6-01-392-74 для дихлорбензола. Качество продуктов в мольных долях соответствующее стандартам приведено в таблице 6.3.1.

Составы потоков в мольных долях.

Таблица 6.3.1.

Компоненты	Поток питания	Бензол	Хлорбензол	Дихлорбензол
Бензол	0,7258	0,9965	0,0015	0,0
Хлорбензол	0,2652	0,0035	0,9985	0,026
Метадихлорбензол	0,009	0,0	0,0	0,974

Индивидуальные свойства компонентов представлены в таблице 6.3.2.

Физико-химические свойства чистых компонентов

Таблица 6.3.2.

	Ткип, ◦С	Плотность, кг/м3	Молекулярная масса, г/моль
Бензол	80,09	883,92	78,11
Хлорбензол	131,68	387,0	112.56
М-Дихлорбензол	172,81	1288,1	147,004

Коэффициенты уравнения Антуана.

Таблица 6.3.3.

	А	В	С
Бензол	15,9008	2788,51	-52,36
Хлорбензол	16,0676	3295,12	-55,6
М-Дихлорбензол	16,8173	4104,12	-43,15

Самый легколетучий компонент в исследуемой смеси - бензол, который является родоначальником углеводородов ароматического ряда. Химические свойства пределяются наличием в молекуле стабильной замкнутой системы -электронов. Для бензола характерна устойчивость к действию высоких температур и окислителей. Лишь выше 650 ⁰С он частично превращается в дифенил, выше 750 ⁰С разлагается на углерод и водород

Основное области применения бензола ( более 80 %) – производство этилбензола, кумола, циклогексана. Остальное количество бензола –для получения анилина, малеинового ангидрида, как компонент моторного топлива для повышения октанового числа , как растворители и экстрагент в производстве лаков, красок ,ПАВ и др.

Хлорбензол- бесцветная жидкость с характерным запахом. Хорошо растворим в органических растворителях, растворимость в воде 0,049% по массе (30 ⁰С), образует азеотропную смесь с водой.

Дихлорбензол растворим а этаноле и диэтиловом эфире. П-дихлорбензол так же растворим в бензоле, хлороформе и хлоре. При хлорировании в присутствии FeCl₃ дихлорбензол превращается в 1,2,4-трихлорбензол (из орто-дихлорбензола образуется также небольшое количество 1,2,3-трихлорбензола).Дальнейшее хлорирование дает Петро-, пента и гексахлорбензолы

В промышленности о-дихлорбензол и п-дихлорбензол выделяют из смеси полихлоридов бензолов, образующихся как побочные продукты в производстве хлорбензола. После отгонки основного количества хлорбензола темную жидкость, содержащую 3-4 % хлорбензола 55-60 % п-дихлорбензола, 35-38 % о-дихлорбензола, примеси трихлорбензолов и смолообразных веществ; осветляют дистилляцией (130 ⁰С,20 кПа),подвергают кристаллизации(0-5 ⁰С) и центрифугированием отделяют п-дихлорбензол.

6.4. Выбор адекватной модели.

 

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

Выбор модели описания парожидкостного взаимодействия для смесей бензол-хлорбензол производился на основании сопоставления экспериментальных и расчетных данных. Расчет парожидкостного равновесия проведен с использованием программного комплекса PRO/IL и её базы данных. На основании имеющегося опыта практического использования и общих рекомендаций по применению выбрано несколько моделей ПЖР, реализованных в программном комплексе:

I - NRTL, II - модель SRK, III – Wilson, IV – UNIQUAC, V – UNIFAC.

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

(6.4.1.)

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

Экспериментальные данные по парожидкостному равновесию пары бензол – хлорбензол.

Таблица 6.4.1

X мол.дол.	Y мол.дол.	t, °С	P, мм.рт.ст.
5.3	17.4	126.4	760
10.4	31.1	121.8
19.2	48.0	115.1
29.5	62.8	108.2
29.6	62.9	108.2
39.9	73.1	102.7
51.4	81.6	97.1
59.1	86.0	93.9
68.4	90.4	90.3
70.3	91.2	89.5
78.6	94.2	86.7
80.4	95.0	86.1
88.4	97.1	83.5

Результаты моделирования парожидкостного равновесия

Таблица 6.4.2.

Смесь	δ,%
	NRTL	UNIQUAC	Wilson	SRK	UNIFAC
Бензол - хлорбензол	2,21	0,665	0,649	1,48	0,627

Судя по среднему относительному отклонению от экспериментальных данных по парожидкостному равновесию, как видно в таблице 6.4.2., наиболее адекватно равновесие бинарной пары бензол-хлорбензол описывает модель UNIFAC. Данные для пар содержащих метадихлобензол в базе данных PROII есть только для модели UNIFAC, поэтому для дальнейших расчётов была принята модель UNIFAC.

  6.5. Синтез схем ректификации для разделения смеси бензол-хлорбензол-м-дихлордензол

Для разделения трехкомпонентной зеотропной смеси нами было предложено четыре технологические схемы (рис. 6.5.1. - 6.5.2.). Две из них являются последовательностью из простых двухсекционных колонн, две другие - являются комплексами из сложных колонн с боковыми секциями. В ряде случаев применение сложных колонн оправдано за счет приближения к термодинамической обратимости за счет структурных особенностей. Кроме того, такие схемы требуют меньшее число кипятильников и дефлегматоров, что может привести к значительному снижению энергозатрат на разделение. Но следует помнить, что такие технологические схемы целесообразно применять при невысокой четкости разделения и содержании среднекипящего компонента более 20%.

Рассмотрим предложенные схемы более подробно.

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

Технологическая схема 2 представлена на рис. 6.5.2. а, в. Первая колонна схемы 2 работает в режиме второго заданного разделения, что позволяет выделить в качестве кубового продукта тяжелокипящий метадихлорбензол. Дистиллат колонны 1 направляют на дальнейшее разделение на бензол и хлорбензол в колонну 2.

Используя широко развитые методы синтеза схем, основанные на теории графов [5, 6, 7, 8], нами синтезированы технологические схемы разделения трехкомпонентной смеси бензол - хлорбензол – метадихлорбензол, содержащие сложные колонны (рис.6.5.1. б, г, 6.5.2. б, г). Схемами-прообразами являются описанные схемы 1 и 2. Для трансформации схем 1 и 2 представим их в виде графов (рис. 6.5.1. в, г, 6.5.2. в, г). Структуры 3 и 4 получают путем стягивания по ориентированному ребру, эксплицирующему потоки между колоннами. Полученные схемы-образы являются структурами с частично (рис. 6.5.1.б,г, 6.5.2.б,г) связанными тепловыми и материальными потоками.

Синтезированная схема 3 (рис. 6.5.1. б) представляет собой сложную колонну с боковой укрепляющей секцией. Схема содержит два дефлегматора и один кипятильник. Схема 4 (рис. 6.5.2. б), напротив, отличается наличием двух кипятильников и одного дефлегматора, представляя тем самым сложную колонну со стриппинг-секцией.

Таким образом, нами предложено четыре схемы разделения, представленные на рис.6.5.3.(схема1, схема2) и рис. 6.5.4. ( схема3, схема4).