Разработка энергосберегающих технологий процесса ректификации продуктов синтеза хлорбензола

CELLPADDING=8 CELLSPACING=0>
США ФРГ Япония

Производство о – и n-дихлорбензолов,

тыс. тон

56 34 25

Мощность по о-дихлорбензолу,

тыс. тон

36,24 - - Потребность 20,4 20,3 20,5 Области потребления ортодихлорбензола, % :

- органический синтез 90 - - - другие цели 5 - -

Таблица 3.2.3.

Мощность по n-дихлорбензолу , тыс. тон	55.89
Потребность п- дихлорбензола на ,тыс. т
-1998	33.5
-1999	34.9
-2000	38.8
. Структура потребления, %:
-специальные дезодоранты	34
-1,2,4-трихлорбензол	9
-полифениленсульфидиые смолы	27

4. Моделирование парожидкостного равновесия

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

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

Модели локальных составов (Wilson, NRTL, UNIQUAC);

Уравнения состояния (SRK, Peng Robinson);

Групповая модель (UNIFAC).

4.1. Модели локальных составов

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

Согласно этой теории раствор рассматривается как упорядоченная структура.

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

4.2. Модель Вильсона [33-35]

Уравнение Вильсона было первым уравнением, в котором была применена концепция локального состава. Основная идея ее состоит в том, что из-за разницы в межмолекулярных взаимодействиях локальный состав вблизи конкретной молекулы в растворе будет отличаться от состава жидкости. Для бинарной пары два параметра связаны со степенью, в которой каждая молекула влияет на состав своего локального окружения. Выражение для коэффициента активности представлено ниже:

(4.2.1.)

где:, (aij, °K); (4.2.2.)

– мольный объем жидкости компонента i.

Параметры aij представляет энергию взаимодействия между молекулами i и j.

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

К недостаткам модели можно отнести то, что уравнение Вильсона не может описывать локальные максимумы или минимумы коэффициента активности, а также не подходит для описания равновесия с частично смешивающимися жидкими фазами.

4.3. Модель NRTL [36-37]

Уравнение NRTL (non-random two-liquid – неслучайное двужидкостное) было разработано Реноном и Праузницем с целью использовать концепцию локального состава в тех случаях, когда уравнение Вильсона неспособно к предсказанию разделения фаз жидкость – жидкость. Модель NRTL дает хорошие результаты для широкого круга систем, в частности для смесей в высокой степени неидеальных и для частично несмешивающихся систем.

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

(4.3.1.)

, (когда единицей измерения является °K); (4.3.2.)

(4.3.3.)

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

4.4. Модель UNIQUAC [38-41]

Уравнение UNIQUAC (universal quasi-chemical – универсальное квазихимическое) было развито Абрамсом и Праусницем на основании статистически-механических положений и решеточной квазихимической модели Гуггенхайма. Каждую молекулу характеризует два параметра: объем и площадь (поверхность) .

Избыточная энергия Гиббса (и соответственно логарифм коэффициента активности) делится на комбинаторную и остаточную части. Комбинаторная часть зависит только от размеров и форм отдельных молекул, она не содержит бинарных параметров. Остаточная часть, которая учитывает энергетические взаимодействия, имеет два регулируемых бинарных параметра. Выражение для коэффициента активности имеет следующий вид:

(4.4.1.)

(4.4.2.)

(4.4.3.)

(4.4.4.) (4.4.5.)

(4.4.6.) (4.4.7.)

(4.4.8.) (4.4.9.) (4.4.10.)

где , – объем и площадь Ван-дер-Ваальса молекулы i

С и R – комбинаторный и остаточный вклады в коэффициент активности жидкости.

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

Уравнение UNIQUAC адекватно описывает широкий ряд систем, оно применяется для неэлектролитических смесей, содержащих полярные или неполярные компоненты, также подходит для частично смешивающихся систем.

4.5. Уравнения состояния [42-43]

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

4.6. Модель SRK

В 1972 г. для улучшения предсказания парового давления чистых компонентов и парожидкостного равновесия многокомпонентных смесей Соав предложил следующую температурную зависимость:

(4.6.1.)

(4.6.2.)

где – приведенная температура, Т/Тс

– ацентрический фактор

Константы в уравнении (4.6.2.) были получены из преобразования данных по паровому давлению для ограниченного числа обычных углеводородов. Эти пределы использования уравнения состояния СРК ограничены неполярными компонентами.

Уравнение состояния Соав-Редлих-Квонга – это модификация уравнения состояния Редлиха-Квонга (которое основано на уравнении Ван-дер-Ваальса). Соав заменил член более общей температурной зависимостью а(Т). Выражение получило следующий вид:

(4.6.3.)

где (4.6.4.)

, – критические температура и давление для i-го компонента

(4.6.5.)

(4.6.6.) (4.6.7.)

(4.6.8.)

(4.6.9.)

– ацентрический фактор для компонента i

– константа бинарного взаимодействия для компонентов i и j

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

Использование формулировки Соава для предсказания свойств смеси включает в себя два этапа. Во-первых, для каждого из компонентов подбирается ацентрический фактор компонента (wi) таким образом, чтобы точно предсказать давление паров компонента. Во-вторых, из экспериментальных данных для бинарных систем с компонентами i и j, для которых достигается фазовое равновесие, определяется параметр kij.

4.7. Модель Peng-Robinson

Уравнение состояния Пенга-Робинсона было опубликовано в 1976 году и является модификацией уравнения Редлиха-Квонга. Во многих отношениях оно похоже на уравнение SRK, но все же он разработано для улучшения предсказания плотности жидкости. Как и в случае уравнения Соава-Редлиха-Квонга, в уравнении Пенга-Робинсона присутствует температурная зависимость а(Т) и выражение имеет следующий вид:

(4.7.1.)

где: (4.7.2.)

(4.7.3.)

, – критические температура и давление компонента i

(4.7.4.)

(4.7.5.) (4.7.6.)

(4.7.7.)

(4.7.8.)

– ацентрический фактор для компонента i

– константа бинарного взаимодействия для компонентов i и j

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

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

4.8. Групповые модели

Групповые модели основаны на предположении об аддитивности вкладов различных химических групп в термодинамические свойства компонентов и смесей. Наибольшее распространение получила модель UNIFAC. Именно её мы использовали для создания псевдоэкспериментальных данных парожидкостного равновесия. Рассмотрим модель более подробно.

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

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

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

(4.8.1.)

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