CC BY
174
паундирования товарных бензинов, позволяющая проводить оценку октановых чисел товарного бензина на выходе после узла смешения. Поскольку полученные расчетные кинетические коэффициенты индивидуальны для каждого углеводорода, программа не требует постоянного их пересчета и остается чувствительной к изменениям состава сырья. Использование программы позволяет опти-
мизировать компаундируемые потоки и выбирать необходимое количество присадок и добавок для получения товарного бензина заданной марки. Погрешность расчетов компьютерной системы не превышает одного пункта, что соответствует требованиям ГОСТ 511-82 о воспроизводимости результатов определения октанового числа по моторному методу.
СПИСОК ЛИТЕРАТУРЫ
1. Лисицын Н.В., Гошкин В.П., Поздяев В.В., Кузичкин Н.В. Методология построения системы оптимального компаундирования товарных нефтепродуктов // Химическая промышленность. - 2003. - № 8. - С. 15-20.
2. Поздяев В.В., Сомов В.Е., Лисицын Н.В., Кузичкин Н.В. Оптимальное компаундирование бензинов // Нефтепереработка и нефтехимия. - 2002. - № 10. - С. 53-57.
3. Киселев Б.Д., Горелова Н.Л. Расчет октановых чисел прямогонных бензинов // Нефтепереработка и нефтехимия. - 1982. - № 3. - С. 10-12.
4. Жоров Ю.М., Гуреев А.А., Смидович Е.В. Производство высокооктановых бензинов. - М.: Химия, 1981. - 219 с.
5. Левин И.А., Попов А.А., Энглин Б.А. Определение октановых чисел бензинов прямой перегонки по их физико-химическим показателям // Нефтепереработка и нефтехимия. - 1985. -№5. - С. 10-12.
6. Кравцов А.В., Иванчина Э.Д. Компьютерное прогнозирование и оптимизация производства бензинов. Физико-химические и технологические основы. - Томск: STT, 2000. - 192 с.
Поступила 27.11.2008 г.
УДК 66.011
ФОРМАЛИЗАЦИЯ СХЕМЫ ПРЕВРАЩЕНИЙ УГЛЕВОДОРОДОВ В ПРОЦЕССЕ ДЕГИДРИРОВАНИЯ ВЫСШИХ АЛКАНОВ С9-С14 НА ПОВЕРХНОСТИ ПЛАТИНОВЫХ КАТАЛИЗАТОРОВ
Е.В. Францина, Ю.И. Афанасьева, А.А. Функ, Е.Н. Ивашкина, А.В. Кравцов
Томский политехнический университет E-mail: [email protected]
На основе квантово-химических расчетов показана термодинамическая вероятность протекания химических реакций в процессе дегидрирования высших алканов С9-Си и предложен уровень формализации механизма превращений на платиновом катализаторе. Составленная схема превращений стала основой кинетической модели процесса дегидрирования высших алканов С9-Сн на платиновом катализаторе, позволяющей учесть влияние химического состава сырья на эффективность его работыы.
Ключевые слова:
Дегидрирование, углеводородыы, квантово-химические расчеты!, термодинамика, моделирование, линейный алкилбензол.
Начиная с 60-х гг. XX в. в мире резко возросли потребности в различных моющих веществах, которые используются для бытовых нужд и в промышленности. В связи с этим увеличилось производство основного компонента синтетических моющих средств. В настоящее время таковыми являются линейные алкилбензолы и линейные алкилбен-золсульфонаты [1]. Для их производства используют до 75 % от общего потребления алканов [2].
Основным методом производства линейных ал-килбензолов, используемым в промышленности, является дегидрирование н-алканов в алкены в присутствии нанесенных платиновых катализаторов с последующим алкилированием ими бензола. Полученные этим способом алкилбензолы и ал-килбензолсульфонаты обладают удовлетворительной скоростью биохимического разложения в
сточных водах, вследствие чего не загрязняются водоемы и не нарушается работа очистных сооружений.
Для повышения эффективности процесса получения н-алкенов разрабатываются новые катализаторы дегидрирования н-алканов, увеличивающие длительность рабочего цикла, а также проводиться поиск оптимальных технологических параметров проведения процесса (температура, давление, мольное соотношение водород/сырье и т. д.). Перечисленные выше приёмы относятся к типу экспериментальных и требуют крупных экономических вложений и длительного времени. В настоящее время для решения подобных задач используют методы математического моделирования и разработанные на их основе технологические моделирующие системы [3]. Успех применения этих методов
определяется, в первую очередь, правильностью и глубиной знаний о механизме протекания реакций на поверхности катализаторов, условиями термодинамического равновесия и количественными кинетическими закономерностями с оценкой мак-рокинетических осложнений. Ранее была разработана математическая модель процесса дегидрирования высших алканов на платиновом катализаторе. Предложенная авторами [3] формализованная схема учитывала превращение алканов на платиновом катализаторе до алкенов и алкадиенов, а также их изомеризацию, дегидроциклизацию и крекинг. В то же время в схеме была учтена побочная реакция коксообразования.
Вместе с тем, выбранный уровень формализации не позволил учесть образование алкенов и ал-кадиенов с двойной связью в различных положениях. Кроме того сделано предположение о термодинамической вероятности протекания реакции изомеризации алканов и не учтено их присутствие в качестве сырьевых компонентов. Это внесло погрешность в расчеты, проведенные с использованием разработанной на основе выбранной схемы математической модели. В связи с этим, целью данной работы стало уточнение схемы превращений веществ в процессе дегидрирования, а также механизма протекающих реакций для увеличения чувствительности математической модели к составу перерабатываемого сырья и повышения ее адекватности. Первым этапом исследований стала оценка термодинамических параметров реакций, протекающих в ходе рассматриваемого процесса.
Для предсказания молекулярных свойств химических систем широко используют методы молекулярно-механического и квантово-химического (полуэмпирического и неэмпирического) моделирования. Они не имеют жесткой привязки к узким классам органических соединений и способны учесть множество эффектов, таких как колебательные и вращательные движения атомов, конфигурацию электронных орбиталей, эффекты сопряжения двойных связей и др. [4]. Данные методы позволяют на основе расчетов электронной структуры молекул предсказать различные молекулярные свойства химических систем: энергию, структуру, спектральные характеристики, термохимические параметры и др. [5].
Оценка термодинамических характеристик реакций и составление формализованной схемы превращений в процессе дегидрирования высших алканов С9-С14 на платиновых катализаторах проводилась с использованием программного пакета Gaussian, содержащего процедуру PM3 (полуэмпи-рический метод расчета) метода NDDO [6]. Отличительной особенностью используемого метода является то, что он с высокой точностью воспроизводит строение и энергетику гипервалентных соединений и обеспечивает достаточную точность для качественного воспроизведения многих физико-
химических свойств молекул.
Анализ литературных данных [7] показал, что существующие до настоящего времени представления о механизме каталитической дегидроциклизации алканов можно разделить на две группы: первая рассматривает ароматизацию алканов как кон-секутивную реакцию с образованием алкенов в качестве промежуточного продукта, а вторая считает образование алкенов и аренов параллельными процессами. По мнению некоторых авторов, ароматизация на оксидных катализаторах протекает через алкены, в то время как на металлах образование алкенов и аренов происходит параллельно [8].
По мере накопления экспериментального материала по дегидроциклизации индивидуальных н-алканов был установлен ряд общих закономерностей в образовании алкенов и аренов, анализ которых свидетельствует о протекании ароматизации алканов через образование алкенов как на оксидных, так и на металлических катализаторах.
Реакции дегидрирования, сопровождающиеся удалением четырех и более атомов водорода, являются ступенчатыми. Промежуточный продукт, как правило, выделяется в объем прежде, чем претерпевает дальнейшее превращение. В некоторых случаях (например, дегидрирование циклогексана в бензол) преобладающим может стать и другое направление реакции - непосредственное превращение в конечный продукт. Таким образом, дегидроциклизация протекает через последовательные стадии десорбции продуктов дегидрирования и адсорбции этих же продуктов в другой форме [8].
Ароматизация н-алканов С6+ происходит через образование алкенов, алкадиенов и (в меньшей степени) триолефинов. Циклизация происходит на более поздних стадиях реакции. Такое направление реакции является основным как на оксидных, так и на металлических катализаторах. Образование промежуточных пяти- или шестичленных алканов составляет значительную долю ароматического продукта, по-видимому, лишь при дегидрировании некоторых разветвленных углеводородов.
Если природу объемных промежуточных продуктов дегидрирования можно считать в большинстве случаев установленной, то для суждения о поверхностных промежуточных соединениях имеется еще мало данных. Схема дегидрирования может быть представлена как последовательность объемных и поверхностных промежуточных соединений. Например, для ароматизации алкена или ал-кана:
спи2п ^ (Спи2п.,) ^ епи2п.2 ^ (Спи2п.2.г) ^
п>6 0<г<2 0<Г]<2
^ СпИ2п.4 ^ (СпН^) ^ СпИ2п.б
С<2<2
(в скобках - поверхностные соединения). Каждая из стадий приведенной схемы может в свою оче-
редь состоять из двух или большего числа ступеней. И наоборот, не все указанные стадии являются обязательными: при определенных условиях, возможно одновременное отщепление более двух атомов водорода от дегидрируемой молекулы, или превращение на некоторых стадиях может протекать ступенчато, но без выделения в газовую фазу промежуточных соединений.
Вместе с тем при составлении математической модели промышленного процесса невозможно учесть все стадии поверхностного механизма. Для сокращения списка параметров модели необходима формализация механизма. Однако построенная на основе выбранной формализованной схемы превращений математическая модель должна оставаться чувствительной к составу перерабатываемого сырья. Лишь в этом случае она будет обладать прогнозирующей способностью.
На основе имеющихся представлений о механизме дегидрирования была составлена предположительная схема превращений веществ в данном процессе. Используя имеющуюся в Gaussian процедуру PM3, были рассчитаны электронные структуры молекул и основные термодинамические параметры этих процессов, таблица.
Таблица. Средние значения термодинамических характеристик реакций (AG, - изменение энергии Гиббса; АН, - изменение энтальпии; AS - изменение энтропии) в процессе дегидрирования (при 753 К, Р=0,20 МПа)
Реакция AGr, кДж/моль АН,, кДж/моль AS, кДж/(моль’К)
1. Алкан ^ алкен-1+Н2 -33,52 49,89 0,11
2. Алкан ^ алкен+Н2 -55,92 33,04 0,12
3. Алкен ^ алкадиен-1,2-2,3+Н2 оо оо 90,37 0,11
4. Алкен ^ алкадиен-2,4-2,(п-2)+Н2 -54,84 32,49 0,12
5. Алкен ^ алкадиен-2,(п-1)+Н2 -33,27 49,89 0,11
6. Изоалкан ^ изоалкен+Н2 -66,72 24,50 0,12
7. Изоалкен ^ изоалкадиен+Н2 -53,49 33,38 0,12
8. Алкен ^ изоалкен 12,15 0,94 -0,02
9. Алкан ^ изоалкан 17,48 2,51 -0,02
10. Алкан ^ циклоалкан+Н2 -49,79 -33,12 0,02
11. Алкан ^ арен+4Н2 -331,21 -58,68 0,36
12. Алкен ^ циклоалкан 6,56 -64,51 -0,09
13. Алкен арен+3Н2 -274,87 -90,07 0,25
14. Циклоалкан арен+3Н2 -318,70 -25,52 0,34
Результаты расчетов показали, что реакции ароматизации являются наиболее термодинамически вероятными (Дфа—300 кДж/моль). Реакции циклизации алкенов, а также изомеризация алканов и ал-кенов не протекают в ходе данного процесса. Также термодинамически маловероятны реакции образования диенов с кумуллированными двойными связями (Д0^9 кДж/моль). Наличие изомерных алке-
нов в продукте обусловлено протеканием процессов дегидрирования изоалканов, присутствующих в сырье. Все остальные возможные реакции термодинамически вероятны и имеют примерно одинаковое значение Дф»-50 кДж/моль.
На основании проведенных термодинамических расчетов была составлена следующая схема превращений (рис. 1). При формализации схемы превращений вещества были объединены в группы псевдокомпонентов по реакционной способности, которая оценивалась по ДОг
Рис. 1. Формализованная схема механизма процесса дегидрирования
Принятая формализованная схема превращений углеводородов на платиновых катализаторах позволяет описать процесс дегидрирования н-де-кана и изодекана следующей стадийной схемой:
С,„И
к
k-1
±С,оИ2о+И2
н-декан 1 н-децен
к
СюН,
■^10 20 ^ н-децен
С10И22 <“ н-декан
С И
10 20 бутилциклогексан
С И ___________:
10 22
*ад8+И2
± С10И20 +И2
3 бутилциклогексан
к
> ВД4 +ЗИ2
бутилбензол
-> СюИ,4 +4И2
н-декан бутилбензол
К
С И
10 20
С,0И
> С10И14 +ЗИ2
бутилбензол
■^10 22 ^
изодекан
± С10И20+И2
С.0И2
■^10 20 ^
изодецен
± СюИ18 +И2
изодецидиен
Тогда кинетическая модель процесса дегидрирования по закону действующих масс запишется следующим образом:
-2
-8
dC
dCcnH2n+2________________k C — kC — k C + k C C + k C C
dt 1 CnH2n+2 (алкан) 3 CnH2n+2 (алкан) 5 CnH2n+2(aлкaн) -1 С^^^кен) H —3 C1H2n (циклoaлкaн) H2
CnH2n(aлкен) _— k C C — k C — k C + k C + k C C
dt -1 CnH2n (алкен) H2 'Ь^С^^кен) ^б^С^^кен) ^ *1 ^„Н^^к ен) ^л—2^CnH2n-2
dCCnH2n-2 (aлкaДиен) _ kC — k C
dt 2 ^Н^^кен) —2 CnH2n-2(aлкaдиен)
dCCnH2n (циклoaлкaн) _ j c 1C 1C C
_ k3CCnH2n+2 (aлкaн) — k4 CCnH2n (циклoaлкaн) — k-3 CCnH2n(циклoaлкaн) C H2
dCcnH2n-6 (apен) _ kC + kC + kC
dt 4 CnH2n (циклoaлкaн) 5 C nH2n+2(алкан) 6 CnH2n(aлкен)
dC
dC
CnH2n+2 (изoaлкaн) _ — k C + k C C
df 7CnH2n+2 (изоалкан) /C—7 ^CnH2n (изоалкен) 2
CnH2n(изoaлкен) _ kC — kC — k C C + k C C
dt 7 CnH2n+2 (изoaлкaн) 8 CnH2n(изoaлкен) -7 С „H2n (изоалкен) 2 —8 Cn H2n-2 (изoaлкaдиен) H2
dC
CnH2n-
dt
CnH2n-2 _ kC — k C C
~ Ка^CnH2n (изoaлкен) К-8 ^CnH2n-2(изoaлкaдиен) ^H 2
dCH2 _
dt _ k1CCnH2n+2(aлкaн) + k2 CCnH2n ^кен) + k3 CCnH2n+2(aлкaн) + 3 ’ k4 CCnH2n(циклoa лкaн) +
+4 • k C + 3 • k C +
^5^^+^^)^J 'Ьз^С^^кен) ^
+ k7CCnH2n+2 (изoaлкaн) +k8 CCnH2n (изoaлкен) k-1 CC „H 2n(aлкен) ^H 2 k-2 CC „H 2n-2 (aлкaциен)CH2 k— 3CCn H2n (циклoaлкaн) CH2
k—7CCnH2n (изoaлкен) CH2 k— 8 CC nH2n-2 (изoaлкaдиен) ^H 2
Начальные условия /=0, Сг-=С0/, где i - соответствующий углеводород (алкан, алкен, алкадиен и т. д.).
Полученная кинетическая модель является формализованной и квазигомогенной, следовательно, константы ^-^ и ^-^ будут эффективными, т. е. представляют собой комбинацию констант всех промежуточных стадий с учетом стадий адсорбции реагентов и десорбции продуктов.
Модель процесса была проверена на адекватность (рис. 2, 3). Показана удовлетворительная сходимость с экспериментальными данными, взятыми с реальной установки. Погрешность расчета концентрации целевого продукта (алкена) составила 5...7 %, побочного (диена) - до 10 % (чем ниже концентрация компонента в смеси, тем больше ошибка при его определении). Погрешность расчетов сопоставима с погрешностью метода хроматографического анализа, который используется для определения концентраций компонентов в продуктовой смеси.
Таким образом, проведенные расчеты показали, что при дегидрировании фракции н-алканов С9-С14
не протекают реакции изомеризации компонентов сырья и образующихся продуктов, а также реакции циклизации алкенов. Также в ходе процесса образуются алкадиены преимущественно с сопряженными и изолированными двойными связями, тогда как образование алкадиенов с кумуллированными связями термодинамически запрещено при данных условиях.
Рис. 2. Расчетные и экспериментальные значения выхода алкенов
UM
0 0.85--------------------------------------------------------Ч-
1 0.8----------------------Ш----------------------------------Лг
10.55 ------------------------—-------- Дата
U 0,5 1--------------1-------1-------1-------1--------1---------1-----1
л& # A' <S>- -d if-' qv (?■ *?■
л#
Ö
<5К *>■ &
----Расчет а Эксперимент
Рис. 3. Расчетные и экспериментальные значения выхода диенов
Предложенный уровень формализации схемы превращений заключается в объединении углеводородов в реакционные серии по показателю ДО. При этом реакционная способность веществ оценивается квантово-химическим методом N000, содержащим процедуру РМ3. Выбранный уровень формализации схемы превращений с одной стороны не загромождает математическое описание и в то же время позволяет учесть образование алкенов и алкадиенов с двойной связью в различных положениях, а также превращение изоалканов в качестве сырьевых компонентов.
Проверка на адекватность составленной формализованной схемы превращений показала удовлетворительную сходимость расчетных и экспериментальных данных.
СПИСОК ЛИТЕРАТУРЫ
1. Баннов П.Г. Процессы переработки нефти. - М.: ЦНИИТЭ-нефтехим, 2000. - 224 с.
2. Пантелеева Е. Линейный алкилбензол: технологии, рынок, перспективы [Электронный ресурс]. - режим доступа: http://www.newchemistry.ru/letter.php. - 12.12.2008.
3. Кравцов А.В., Иванчина Э.Д., Сизов С.В., Ивашкина Е.Н. Разработка формализованного механизма дегидрирования высших парафинов С10-С13 на Pt-катализаторах // Известия Томского политехнического университета. - 2006. - Т. 309. - № 5. - С. 81-85.
4. Баскин И.И., Полюлин В.А., Зефиров Н.С. Прогнозирование энтальпий образования алифатических полинитросоединений // Вестник Московского университета. Серия 2. Химия. -2001. - Т. 42. - № 6. - С. 387-389.
5. Полещук О.Х., Кижнер Д.М. Химические исследования методами расчета электронной структуры молекул. - Томск: Изд-во ТПУ, 2006. - 146 с.
6. Хурсан С.Л. Квантовая механика и квантовая химия. - Уфа: чП Раянов, 2005. - 164 с.
7. Кравцов А.В., Иванчина Э.Д., Юрьев Е.М., Ивашкина Е.Н. 11-технологии в решении проблем промышленного процесса дегидрирования высших парафинов. - Томск: 8ТТ, 2008. - 230 с.
8. Скарченко В.К. Дегидрирование углеводородов. - Киев: Нау-кова думка, 1981. - 328 с.
Поступила 20.01.2009 г.
УДК 66.011
РАЗРАБОТКА КИНЕТИЧЕСКОЙ МОДЕЛИ ПРОЦЕССА АЛКИЛИРОВАНИЯ БЕНЗОЛА ОЛЕФИНАМИ
И.О. Шнидорова, В.А. Фетисова, Е.Н. Ивашкина, Э.Д. Иванчина, А.А. Функ
Томский политехнический университет E-mail: [email protected]
С использованием квантово-химических методов расчета термодинамических функций оценена термодинамическая возможность протекания целевых и побочных реакций в процессе жидкофазного каталитического алкилирования бензола высшими олефинами. Составленная на основе результатов расчета формализованная схема превращений углеводородов положена в основу кинетической модели процесса.
Ключевые слова:
Алкилирование, высшие моноолефины, квантово-химический расчет, кинетическая модель, схема превращений, линейный алкилбензол.
Производство синтетических моющих средств является одним из важных направлений нефтехимической промышленности. Наибольшее распространение среди анионных синтетических моющих веществ получили алкиларилсульфонаты, в первую очередь алкилбензосульфонаты натрия. Эти вещества получили название сульфонолов [1].
Линейные (мягкие) алкилбензосульфонаты получают из линейных алкилбензолов (ЛАБ) (с алкильным радикалом С10-С14) путем их сульфирования оксидом трехвалентной серы или олеумом с последующей нейтрализацией. Такие поверхностно-активные вещества обладают высокими моющими свойствами и хорошей биологической разлагаемостью.
