Технологическая схема непрерывного алкилирования бензола, реализованная в среде комплекса программ CHEMCAD: модули расчета: 1 – модуль равновесного реактора; 6, 9, 15, 17, 19, 29, 31, 33 – модули теплообменников; 14, 16, 18, 28, 30, 32 – модули статических контроллеров; 3, 11, 13 – модули сепараторов; 4, 7, 8 – модули ректификационных колонн; 5 – модуль насоса; 2, 12, 22, 26, 27, 35 – модули смесителей, 20, 21, 23, 24, 25, 34 – модули делителей потоков

Авторы: А.Ю. Медведев, Ф.С. Советин, Т.Н. Гартман, А.С. Павлов, С.А. Сеннер (РХТУ им. Д. И. Менделеева).

Опубликовано в журнале Химическая техника №7/2015

Алкилароматические соединения широко используются в технологии основного органического и нефтехимического синтеза, а также в технологии тонкого органического синтеза для получения полимерных материалов, поверхностно-активных веществ, высокооктановых добавок к топливу и в качестве промежуточных продуктов [1–3]. Наибольшее значение имеют этилбензол и изопропилбензол как исходные вещества при получении таких мономеров, как стирол и β-метилстирол.

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

Большинство действующих крупнотоннажных производств алкилбензолов в России используют традиционные технологии алкилирования, где в качестве алкилирующих агентов используют олефины (этилен, пропилен и т. п.) и процесс проводится на катализаторе (хлористый алюминий или ортофосфорная кислота на носителе) [5–6]. Этилбензол в промышленности получают алкилированием бензола этиленом. Катализатором данного процесса выступают кислоты Льюиса, чаще всего AlCl3 [7]. При этом разделение смесей продуктов реакций (алкилаты) осуществляется ректификацией. Данные процессы являются одними из наиболее энергозатратных процессов в химической промышленности [8–9].

Модернизация указанных технологических схем с позиции энергои ресурсосбережения данных производств не представляется возможной без применения современных комплексов программ, позволяющих разрабатывать компьютерные модели технологических линий, создавать так называемые «виртуальные» производства. Кроме того, использование комплексов программ предоставляет возможность без значительных материальных и временных затрат проводить исследования указанных процессов на их моделях (так называемые вычислительные эксперименты) [10–11].

Работа [12] посвящена оценке эффективности нового смесительного устройства, а также оценке целесообразности реконструкции смесительного оборудования перед реакторным блоком путем установки алкилирования бензола этиленом. В указанной работе для решения поставленных задач использовали комплексы программ ABAQUS и Flow Vision, позволяющие в совокупности с высокой точностью решать задачи по моделированию гидродинамики смешения жидких и газовых потоков.

Интервальное оценивание параметров математического описания процесса жидкофазного алкилирования фенилацетонитрила в микроструктурном реакторе на основе вычислительного эксперимента произведено в работе [13]. В этой работе предлагаются методы определения интервальных оценок параметров математического описания химико-технологических процессов при моделировании, а также рассматриваются компьютерные технологии решения задач интервального анализа.

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

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

Описание технологической схемы процесса непрерывного алкилирования бензола этиленом

Этилбензол и диэтилбензол получают методом алкилирования бензола этиленом. Описание технологической схемы процесса непрерывного алкилирования бензола этиленом заимствовано из работ [3–4]. Технологическая линия данного процесса, реализованная в среде комплекса программ CHEMCAD, изображена на рисунке.

Процесс алкилирования бензола этиленом проводится в полом реакторе колонного типа в режиме барботажа [3–4]. В реактор вводятся потоки бензола и газообразного этилена. Процесс проводится при 100°С на катализаторе (хлористый алюминий или ортофосфорная кислота на носителе) в соответствии с уравнениями:

C6H6 + C2H4 → C6H5C2H5 (1)

C6H6 + 2C2H4 → C6H4(C2H5)2 (2)

C6H6 + 3C2H4 → C6H3(C2H5)3 (3)

C6H4(C2H5)2 + C6H6 → 2C6H5C2H5 (4)

На основе анализа экспериментальных данных, приведенных в работах [3–4], реакции образования полиалкилбензолов практически не протекают (содержание полиалкилбензолов в выходном потоке незначительно, поэтому реакциями их образования можно пренебречь).

Разделение смеси продуктов алкилирования (алкилаты) проводится на трехколонной ректификационной установке. Порядок отгонки компонентов проводится по порядку повышения их температур кипения: в первой колонне отгоняют непрореагировавший бензол, во второй – этилбензол, в третьей – диэтилбензол.

Для обогрева ректификационных колонн в кипятильники подаются пары высокотемпературных органических теплоносителей (ВОТ) – смесь дифенила и дифенилового эфира. Конденсаторы ректификационных колонн

Описание хода моделирования основных аппаратов и технологической схемы процесса непрерывного алкилирования бензола

Реактор алкилирования бензола моделируется комбинацией стандартных расчетных моделей комплекса программ CHEMCAD – модулем равновесного реактора и модулем сепаратора компонентов (см. рисунок, модули 1 и 3). Для модуля равновесного реактора (модуль 1) указаны следующие условия: стехиометрические данные реакций (1)–(4), конверсии базовых реагентов, изотермический тепловой режим (температура 100°С).

Значения заданных конверсий для базовых реагентов при моделировании процесса алкилирования бензола

Реагент Реакция Конверсия, %
C6H6 (1) 50
C6H6 (2) 21,5
C6H6 (3) 2
C6H4(C2H5)2 (4) 10

Указанные значения конверсий базовых реагентов, предлагаемые для моделирования процесса алкилирования бензола, определенные по экспериментальным данным, заимствованы из работ [3–4].

Для модуля сепаратора (модуль 3) заданы доли газообразных компонентов, равные 1. Подобная комбинация стандартных расчетных модулей комплекса программ CHEMCAD для моделирования сложных химических реакторов применялась в работе [9]. Результаты расчета процесса алкилирования бензола приведены в табл. 1. В этой же таблице проведено сравнение результатов рассчитанного потока алкилата с экспериментальными данными, заимствованными из работ [3–4].

Таблица 1

Сравнение результатов расчетов процесса алкилирования  бензола с экспериментальными данными

Параметр На входе – смесь этилена и бензола, подаваемая на алкилирование На выходе – алкилат
Расчетные данные Расчетные данные Экспериментальные данные
Давление, МПа 0,2 0,2
Температура, °С 25 100
Массовый расход, т/ч 106,17 105,703
Массовая доля компонента, % мас.
C6H6 73,24 19,64 19,6
C6H5C2H5 0,3 50,09 49,3
C6H4(C2H5)2 0,04 27,21 27,1
C6H3(C2H5)3 3,06 3,57
C2H4 26,42

Первая ректификационная колонна – колонна отгонки бензола (см. рисунок, модуль 4) моделируется модулем строгого расчета ректификации TOWR с использованием алгоритма Inside/Out [14]. Число тарелок – 20, № тарелки питания – 10, массовый расход дистиллята – 20 500 кг/ч, паровое число – 15. Поток теплоносителя (пары ВОТ), подаваемого в кипятильник колонны, моделируется модулем теплообменника (см. рис. 1, модуль 29).

Технологическая схема непрерывного алкилирования бензола, реализованная в среде комплекса программ CHEMCAD: модули расчета: 1 – модуль равновесного реактора; 6, 9, 15, 17, 19, 29, 31, 33 – модули теплообменников; 14, 16, 18, 28, 30, 32 – модули статических контроллеров; 3, 11, 13 – модули сепараторов; 4, 7, 8 – модули ректификационных колонн; 5 – модуль насоса; 2, 12, 22, 26, 27, 35 – модули смесителей, 20, 21, 23, 24, 25, 34 – модули делителей потоков
Технологическая схема непрерывного алкилирования бензола, реализованная в среде комплекса программ CHEMCAD:
модули расчета: 1 – модуль равновесного реактора;
6, 9, 15, 17, 19, 29, 31, 33 – модули теплообменников;
14, 16, 18, 28, 30, 32 – модули статических контроллеров;
3, 11, 13 – модули сепараторов; 4, 7, 8 – модули ректификационных колонн;
5 – модуль насоса; 2, 12, 22, 26, 27, 35 – модули смесителей,
20, 21, 23, 24, 25, 34 – модули делителей потоков

Процесс теплопередачи моделируется модулем статического контроллера (модуль 28), который передает расчетную тепловую нагрузку кипятильника колонны модулю теплообменника с учетом знака. Поток воды, подаваемой для охлаждения конденсатора колонны, моделируется модулем теплообменника (модуль 15). Процесс теплопередачи также моделируется модулем статического контроллера (модуль 14), который передает расчетную тепловую нагрузку конденсатора колонны модулю теплообменника. Применена процедура комбинирования стандартных модулей комплекса программ CHEMCAD для моделирования процессов ректификации, разработанная в работах [15–16]. Результаты расчета данной колонны приведены в табл. 2.

Таблица 2

Результаты расчета колонны отгонки бензола

 Параметр  На входе На выходе
Дистиллят Кубовый остаток
Давление, МПа 0,5 0,3 0,3
Температура, °С 100 121 193
Массовый расход, т/ч 105,703 20,5 85,203
Диаметр колонны, м 4 4 4
Высота колонны, м 10 10 10
Тепловая нагрузка кипятильника, МДж/ч 389 915 389 915 389 915
Тепловая нагрузка конденсатора, МДж/ч –367 182 –367 182 –367 182
Расход паров ВОТ для обогрева кипятильника, т/ч 1 000 1 000 1 000
Расход воды для охлаждения конденсатора, т/ч 2 500 2 500 2 500
Массовая доля компонента, % мас.
C6H6 19,64 98,25 0,73
C6H5C2H5 50,09 1,56 61,77
C6H4(C2H5)2 27,21 0,19 33,71
C6H3(C2H5)3 3,06 3,79

Вторая ректификационная колонна – колонна отгонки этилбензола (см. рисунок, модуль 7) моделируется модулем строгого расчета ректификации TOWR с использованием алгоритма Inside/Out. Число тарелок – 20, № тарелки питания – 10, массовый расход дистиллята – 53 500 кг/ч, паровое число – 9. Поток теплоносителя, подаваемого в кипятильник колонны, моделируется модулем теплообменника (модуль 31). Процесс теплопередачи моделируется модулем статического контроллера (модуль 30), который передает расчетную тепловую нагрузку кипятильника колонны модулю теплообменника с учетом знака. Поток воды, подаваемой для охлаждения конденсатора колонны, моделируется модулем теплообменника (модуль 17). Процесс теплопередачи также моделируется модулем статического контроллера (модуль 16), который передает расчетную тепловую нагрузку конденсатора колонны модулю теплообменника. Результаты расчета данной колонны приведены в табл. 3.

Таблица 3

Результаты расчета колонны отгонки этилбензола

 Параметр  На входе На выходе
Дистиллят Кубовый остаток
Давление, МПа 0,3 0,2 0,2
Температура, °С 193 162 217
Массовый расход, т/ч 85,203 53,5 31,703
Диаметр колонны, м 4 4 4
Высота колонны, м 12 12 12
Тепловая нагрузка кипятильника, МДж/ч 83 337 83 337 83 337
Тепловая нагрузка конденсатора, МДж/ч –85 498 –85 498 –85 498
Расход паров ВОТ для обогрева кипятильника, т/ч 150 150 150
Расход воды для охлаждения конденсатора, т/ч 1 000 1 000 1 000
Массовая доля компонента, % мас.
C6H6 0,73 1,15
C6H5C2H5 61,77 98,38
C6H4(C2H5)2 33,71 0,47 89,81
C6H3(C2H5)3 3,79 10,19

Третья ректификационная колонна – колонна отгонки диэтилбензола (см. рисунок, модуль 8) моделируется модулем строгого расчета ректификации TOWR с использованием алгоритма Inside/Out. Число тарелок – 20, номер тарелки питания – 10, массовый расход кубового остатка – 3 245 кг/ч, флегмовое число – 2. Поток теплоносителя, подаваемого в кипятильник колонны, моделируется модулем теплообменника (см. рис. 1, модуль 33).

Процесс теплопередачи моделируется модулем статического контроллера (модуль 32), который передает расчетную тепловую нагрузку кипятильника колонны модулю теплообменника с учетом знака. Поток воды, подаваемой для охлаждения конденсатора колонны, моделируется модулем теплообменника (модуль 19). Процесс теплопередачи также моделируется модулем статического контроллера (модуль 18), который передает расчетную тепловую нагрузку конденсатора колонны модулю теплообменника. Результаты расчета данной колонны приведены в табл. 4.

Таблица 4

Результаты расчета колонны отгонки диэтилбензола

 Параметр  На входе На выходе
Дистиллят Кубовый остаток
Давление, МПа 0,2 0,15 0,15
Температура, °С 217 199 286
Массовый расход, т/ч 31,703 28,46 3,24
Диаметр колонны, м 2,6 2,6 2,6
Высота колонны, м 12 12 12
Тепловая нагрузка кипятильника, МДж/ч 24 301 24 301 24 301
Тепловая нагрузка конденсатора, МДж/ч –25 388 –25 388 –25 388
Расход паров ВОТ для обогрева кипятильника, т/ч 350 350 350
Расход воды для охлаждения конденсатора, т/ч 1 000 1 000 1 000
Массовая доля компонента, % мас.
C6H4(C2H5)2 89,81 99,99 0,45
C6H3(C2H5)3 10,19 0,01 99,56

При расчете технологической схемы непрерывного алкилирования бензола для обеспечения сходимости рецикла используется метод простых итераций. Также использовано методическое обеспечение для решения задачи компьютерного моделирования энергои ресурсоемких химических производств, разработанное в работе [16].

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

Список литературы

  1. Топчиев А.В. Избранные труды. Кн. 2. Алкилирование. М.: Наука, 1965.
  2. Олбрайт Л.Ф., Голдсби А.Р. Алкилирование. Исследование и промышленное оформление процесса. Пер. с англ. Под ред. Липовича В.Г. М.: Химия, 1982.
  3. Орехов В.С., Дьячкова Т.П., Субочева М.Ю., Колмакова М.А. Технология органических полупродуктов. Ч. 1. Тамбов: ТГТУ, 2007.
  4. Бесков В.С. Общая химическая технология. М.: Академкнига, 2005.
  5. Новые процессы органического синтеза/Под. Ред. Черных С.П. М.: Химия, 1989.
  6. Байрамов М.Р., Гаджиев М.М., Зейналов Н.Ю., Алиева Ф.Д., Гусейнова Р.А. Разработка технологии процесса алкилирования бензола олигомерами бутенов-1//Химия и химическая технология. 2006. Т. 49. Вып. 7.
  7. Серафимов Л.А., Тимофеев В.С., Тимошенко А.В. Принципы технологии основного органического и нефтехимического синтеза. М.: «Высшая школа», 2010.
  8. Гартман Т.Н., Советин Ф.С. Применение пакетов программ CHEMCAD для моделирования процессов многокомпонентной ректификации в тарельчатых колоннах при получении синтетического жидкого топлива//Химическая техника. 2010.№2.
  9. Советин Ф.С., Гартман Т.Н., Сафонова В.Д., Шумакова О.П., Тамбовцев И.И. Применение комплекса программ CHEMCAD для разработки компьютерной модели технологической схемы хлорирования бензола по Беркману//Химическая техника. 2014. №10.
  10. Гартман Т.Н., Клушин Д.В. Основы компьютерного моделирования химико-технологических процессов. М: Академкнига, 2008.
  11. Гартман Т.Н., Советин Ф.С. Аналитический обзор современных пакетов моделирующих программ для компьютерного моделирования химико-технологических систем//Успехи в химии и химической технологии. 2012. Т. 26. №11 (140).
  12. Ивашкина Е.Н., Хлебникова Е.С., Беккер А.В., Белинская Н.С., Ткачев В.В. Исследование процесса смешения реагентов в технологии алкилирования бензола этиленом с использованием методов вычислительной гидродинамики//Химическая промышленность сегодня. 2014. №8.
  13. Боровинская Е.С., Решетиловский В.П., Холоднов В.А. Интервальное оценивание параметров математического описания химико-технологических процессов на основе вычислительного эксперимента//Известия ОрелГТУ. Фундаментальные и прикладные проблемы техники и технологии: информационные системы и технологии. 2009. №2 (52)/563.
  14. Гартман Т.Н. Анализ и синтез непрерывной ректификации на основе автоматизированного накопления и классификации информации. Дис … д-ра техн. наук. М.: 2000.
  15. Советин Ф.С., Гартман Т.Н. Логико-вычислительные процедуры разработки блочных компьютерных моделей реакторных и ректификационных процессов//Изв. ТулГУ. Технические науки. 2011. Вып. 5. Ч. 3.
  16. Советин Ф.С. Разработка и применение методического обеспечения блочного компьютерного моделирования энергоресурсоемких химико-технологических систем c применением инструментальных комплексов программ. Дис. … канд. техн. наук. М.: РХТУ им. Д. И. Менделеева. 2011.