НАУКА и ОБРАЗОВАНИЕ

УДК 004.942

Институт нефтехимии и катализа РАН, г. Уфа

OOO «БашНИПИнефть», г. Уфа

IrekMars@mail.ru

SabitovaKamila@rambler.ru

LindUB@bashneft.ru

Введение

Термическое гидроалюминирование олефинов до сих пор является одним из наиболее известных и широко применяемых способов синтеза высших алюминийорганических соединений. Однако этот метод отличается рядом специфических свойств, ограничивающих его применение – реакцию необходимо проводить при повышенной температуре, и она применима не ко всем олефинам. В лабораториях Института нефтехимии и катализа РАН под руководством директора института, чл.-корр. РАН У.М. Джемилева проводится экспериментальное изучение альтернативного способа – каталитического гидроалюминирования олефинов изобутилаланами.

Исследование механизма сложных реакций, протекающих в многокомпонентных системах, требует решения двух фундаментальных задач: 1) экспериментальное исследование химического объекта; 2) математическая обработка экспериментальных данных с целью построения и идентификации кинетической модели. При этом большое значение приобретает скорость расчета. Член-корреспондент РАН М.Г. Слинько в одной из своих последних публикаций пишет «Чтобы кинетическая модель могла служить опорой для решения практических вопросов, предварительная кинетическая модель должна быть получена не более, чем за месяц, а полная модель не более, чем за 4 месяца». Для достижения этого требования необходимо использование при расчете современных вычислительных технологий, таких как параллельные вычисления на высокопроизводительных многопроцессорных вычислительных системах (суперкомпьютерах).

1.      Химическая постановка задачи

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

HAlBuⁱ₂ – диизобутилалюминийгидрид (ДИБАГ),

ClAlBuⁱ₂ - диизобутилалюминийхлорид (ДИБАХ),

AlBuⁱ₃ - триизобутиалюминий, ТИБА

и олефинами

гептен-1,

октен-1,

нонен-1,

децен-1.

На этой основе был детализирован общий механизм реакции. В работе рассмотрен общий механизм гидроалюминирования олефинов с ДИБАХ.

Экспериментаторами Института нефтехимии и катализа РАН было предложено несколько механизмов протекания реакции. На начальном этапе изучения реакции удалось выделить частные реакции с АОС и олефинами в виде итоговых уравнений. Далее велись работы по детализации этого механизма до элементарных стадий. В табл. 1 представлен первый этап детализации частных реакций на основании квантово-химических расчетов (схема 1). Например, был сделан вывод, что димер (X₁) до вступления в реакцию с АОС разлагается до мономера (X₂), который уже реагирует с АОС. Дальнейшие исследования показали, что не только мономер реагирует с АОС, но и сам димер вступает в реакцию с алкеном. Так была получена вторая (более подробная) детализация (табл. 1, схема 2). На основании второй детализации был предложен обобщенный механизм гидроалюминирования олефинов алкилаланами, катализируемого  в присутствии ДИБАХ, на основе элементарных стадий (табл. 1, схема 3).

Таблица 1

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

Здесь приняты следующие обозначения:

X₁=[Cp₂ZrH₂·ClAlBu₂]₂;   X₂=[Cp₂ZrH₂·ClAlBu₂];

X₃=CH₂CHR;    X₄=Cp₂ZrCl(CH₂CH₂R);

X₅=HAlBu₂ (ДИБАГ);   X₆=Bu₂Al(CH₂CH₂R);

X₇=Cp₂ZrHCl;    X₈=[Cp₂ZrH₂·HAlBu₂·ClAlBu₂];

X₉=ClAlBu₂ (ДИБАХ); X₁₀=[Cp₂ZrHCl·ClAlBu₂];

X₁₁=l²AlBu;    X₁₂=[Cp₂ZrHBu·ClAlBu₂];

X₁₃=C₄H_8;    X₁₄=AlBu₃ (ТИБА);

X₁₅=Сp₂ZrCl_2;   X₁₆=[Cp₂ZrH₂·HAlBu₂·2(ClAlBu₂)];

X₁₇=[Cp₂ZrH₂·HAlBu₂·ClAlBu₂]; X₁₈=Cp₂ZrClBu;

X₁₉=ClBuAl(CH₂CH₂R); X₂₀=Cp₂ZrHBu∙ClAlBu₂, Bu=C₄H₉, Cp=C₅H₅.

2.      Математическое описание химической реакции

При математическом моделировании химических реакций необходимо учитывать, что для одинаковых стадий в разных механизмах реакций кинетические параметры должны быть одинаковыми (допускается отклонение 10% [1]).

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

                                                    (1)

где  – стехиометрические коэффициенты;  – скорость j-ой стадии, 1/с;  – количество катализатора, моль;  – объем реакционной смеси, л; J – количество стадий.

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

Таблица 2

Кинетические уравнения реакции гидроалюминирования олефинов олефинов

Квантово-химические исследования [6] показали существование ограничений на константы скоростей стадий реакции, т.е. для первой детализации

                                                               (2)

для второй детализации

                                                    (3)

Эти соотношения необходимо учитывать при проведении расчетов.

3.      Организация вычислительного процесса

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

Более подробное изучение химической реакции, ее детализация до элементарных стадий приводят к усложнению задачи. Поэтому необходимо разрабатывать новые способы обработки информации в зависимости от сложности решаемых задач. Эти методы включают в себя: 1) принципы математического моделирования; 2) численные методы и алгоритмы; 3) системы управления базами данных; 4) программирование; 5) проведение многочисленных вычислительных экспериментов. Усложнение задач ведет и к усложнению расчетов, к увеличению времени работы программ. Это вызывает необходимость обновления технических средств обработки информации, увеличения скорости вычислений и объемов хранения данных, использования технологий параллельных вычислений. Все эти процессы являются компонентами информационно-аналитической системы обратных задач химической кинетики (ИАС ОЗХК) (рис. 1) [5].

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

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

Из рис. 1 видно, что входные и выходные информационные потоки образуют единую базу данных. Так как информация о химических реакциях чаще всего предоставляется химиками в виде таблиц [5], используются реляционные базы данных. Для организации базы данных необходимо составить инфологическую модель предметной области. Для этого используется язык ER-диаграмм  (рис. 2).

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

На рис. 2 построена диаграмма, отражающая связи между основными сущностями в сложной химической реакции. Так, для экземпляра сущности «Реакция» имеется много или один экземпляр сущности «Вещество». Для экземпляра сущности «Вещество» может быть много или ноль экземпляров «Реакция», так как вещество может участвовать во многих реакциях или ни в одной. Для каждой реакции может быть поставлено много или ноль экспериментов, но каждый эксперимент поставлен только для одной реакции. По каждому эксперименту проводится много или ноль замеров, но каждый замер соответствует одному эксперименту. Для каждой реакции можно составить много или ноль механизмов ее протекания, но каждый механизм соответствует только одной реакции. Каждый механизм имеет много или хотя бы одну стадию, а каждая стадия может участвовать во многих или ноль механизмах. Каждая стадия имеет один или несколько (для обратимой стадии) кинетических параметров (констант скоростей, энергий активации), который может участвовать в много или ноль стадиях.

Для упрощения доступа к данным при распараллеливании вычислительного процесса произведена детализация базы данных. Т.е., для реакции каталитического гидроалюминирования олефинов при изучении макростадий база данных была общей. После выделения итоговых уравнений по АОС и олефинам единая база данных была разбита на две базы: для частных схем для АОС и олефинов, и для обобщенной схемы реакции. Дальнейшее детализирование схем позволило разделить базу данных на семь баз – по трем АОС и олефинам для частных реакций и по каждому АОС для обобщенной реакции. Такое разделение базы данных позволяет вести параллельный расчет по этим реакциям с использованием ранее полученных результатов (рис. 3).

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

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

– существует несколько наборов экспериментальных данных, полученных при разных условиях (обычно более пяти); для учета погрешности эксперимента обрабатываются все и выбираются два-три лучших, то есть наиболее корректно описывающих экспериментальные данные;

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

Для последовательного решения всех этих задач иногда необходимы существенные затраты времени (от нескольких месяцев до года), что привело авторов к идее последовательно-параллельного ведения расчета [7]. Предлагается сгруппировать задачи, по которым можно вести параллельный расчет, а расчет по этим группам вести последовательно (рис. 4).

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

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

Уровень 2. Для каждого механизма реакций можно параллельно обрабатывать все предложенные эксперименты. Определив константы скоростей при разных температурах, из уравнения Аррениуса [1] методом наименьших квадратов можно получить значения энергий активации для каждой стадии, на основе которых определяются константы скоростей для любой температуры.

Уровень 3. Реализуется принцип геометрического параллелизма, который предполагает декомпозицию расчетной области изменения кинетических констант на подобласти соответственно числу процессоров. Так, для двух кинетических констант строится двумерная плоскость (соответственно, для n констант – n-мерная). Область делится на заданное число подобластей, которые распределяются между процессорами. Каждый процессор решает обратную задачу (минимизирует невязку в своей области методом параболического спуска [2]) только в своих подобластях. После нахождения в каждой подобласти значений констант, минимизирующих функционал невязки [2], исключаются те из них, в который этот функционал имеет наибольшие значения. Из новой области формируются подобласти меньшего размера и процедура повторяется до тех пор, пока не определятся константы, наилучшим образом описывающие эксперимент.

4.      Кинетическая модель реакции гидроалюминирования олефинов

После обработки экспериментальных данных для детализации  частных реакций были найдены кинетические параметры стадий. Выписаны приведенные значения кинетических параметров относительно начальной мольной плотности реакционной смеси (табл. 3, 4). Результаты первой детализации удовлетворяют условиям (2). Здесь и далее величины k_i, k_–i измеряются в мин^–1, а величина E_j  - в  ккал/моль.

Таблица 3

Кинетические параметры стадий частных реакций с АОС по первой детализации

Таблица 4

Кинетические параметры стадий частных реакций с олефинами (значения констант k₇ (k₈) и энергий активаций E₇ (E₈)) по первой детализации

На рис. 5 приведены результаты сравнения экспериментально заданных и рассчитанных по модели (1) с учетом (2) концентраций наблюдаемого вещества для частных реакций с АОС и олефинами по первому этапу детализации. Среднее отклонение расчета от эксперимента для предложенной кинетической модели не превышает 10%. Сравнивая скорости прямой и обратной реакции образования мономерного комплекса X₂ из X₁, можно сделать вывод, что в интервале температур -60÷50⁰С равновесие смещено в сторону димерной формы, однако, с увеличением температуры доля X₂ в смеси возрастает. Самой быстрой оказалась стадия перехода комплекса X₂ в тригидридный комплекс X₈, а равновесие X₂ ↔ X₈ сильно смещено в сторону тригидридного комплекса, особенно при низких температурах. В реакции гидрометаллирования олефинов комплексом X₁ лимитирующими оказались равновероятные стадии взаимодействия мономерного комплекса X₂ с алкенами. По-видимому, высокая скорость гидрометаллирования обеспечивается существованием обратной реакции перехода комплекса X₈ в активный по отношению к олефинам мономер.

а) a)   a)	б)
в)	г)
Рис. 5. Экспериментальные и расчетные результаты (детализация I): а) реакция с НAlBui2, T= -600C; б) реакция с ClAlBui2, T = -200C; в) реакция с AlBui3, T = -200C; г) реакция с деценом-1, T = 70C

Соотношение (2) для детализации I выполняется в рассматриваемом экспериментатором интервале температур. Построив на основе уравнения Аррениуса зависимость ln(k) от 1/T для k₁ и k_-1, можно увидеть, для каких температур это соотношение будет нарушено (рис. 6).

Из рис. 6 видно, что графики функций ln(k₁), ln(k_-1) пересекаются в точке , левее которой соотношение (2) нарушается. При этом T = 400 K = 127°C, т.е. для температур выше 127°C соотношение (2) нарушается. Поэтому был сделан вывод, что первая детализация частных реакций не в полной мере характеризует экспериментальные и теоретические данные.

Для второй детализации частных реакций выписаны приведенные значения кинетических параметров относительно начальной мольной плотности реакционной смеси (табл. 5, 6).

Таблица 5

Кинетические параметры стадий частных реакций с АОС по второй детализации

Таблица 6

Кинетические параметры стадий частных реакций с олефинами по второй детализации

Результаты второй детализации удовлетворяют условиям (3): константа  обратной реакции значительно больше константы прямой реакции, k₄ больше чем k₂, k₁₁ много больше k₁₄, что корректно описывает эксперимент.

Результаты сравнения эксперимента и расчета для всех АОС и олефинов представлены на рис. 7.

Отклонение расчетных результатов от экспериментальных не превышает 8,5 %.

налогично предыдущим вычислениям, для выяснения пределов соблюдения соотношения (3) построим зависимость ln(k) от 1/T для k₁ и k_-1 (рис. 8).

Из рис. 8 видно, что графики функций ln(k₁) и ln(k_-1) пересекаются в точке , правее которой соотношение (3) нарушается.

При этом T=173 K=-100°C, т.е. для температур ниже (-100)°C соотношение (3) нарушается. Такие температуры встречаются крайне редко при постановке экспериментов, поэтому можно считать, что вторая детализация частных реакций достаточно хорошо характеризует экспериментальные и теоретические данные.

Для общей реакции гидроалюминирования олефинов в присутствии ДИБАХ построены кинетические модели без учета количества катализатора и с учетом количества катализатора. Сделан вывод о том, что модель (1), учитывающая количество катализатора, лучше описывает реакцию, и практически все кинетические параметры стадий, определенные по этой модели, совпадают с соответствующими кинетическими параметрами частных реакций (исключение составляет k₉(E₉)). В табл. 7 представлены кинетические параметры стадий общей реакции, полученные с использованием модели (1). Представлены приведенные константы относительно начальной мольной плотности реакционной смеси. Результаты расчета и эксперимента для общей реакции гидроалюминирования олефинов в присутствии диизобутилалюминийхлорида  по модели c учетом катализатора (x_кт – катализатор) представлены на рис. 9, 10. Отклонение расчетных результатов от экспериментальных данных не превышает 10%.

Таблица 7

Кинетические параметры стадий реакции с учетом количества катализатора

На рис. 11 представлена кинетическая модель обобщенной реакции каталитического гидроалюминирования олефинов с помощью ДИБАХ. Значения кинетических параметров обобщенной реакции, полученные посредством математического моделирования, подтверждают выводы, полученные квантово-химическим подходом [6]. Таким образом, в рассматриваемом случае макрокинетика (экспериментально-расчетные исследования на уровне молекул) обосновывает выводы микрокинетики (теоретические исследования на уровне атомов). Построенная математическая модель реакции позволяет определить кинетические параметры реакции в условиях, которые сложно создать экспериментально.

Из рис. 11 видно, что первоначально Сp₂ZrCl₂ взаимодействует с молекулой ClAlBu₂ⁱ с образованием Cp₂ZrClBu (E=31 ккал/моль). Эта стадия обладает наибольшим активационным барьером. В начальный момент времени наблюдается низкая скорость межлигандного обмена между исходными веществами, благодаря чему в реакции гидроалюминирования олефинов с ДИБАХ наблюдается индукционный период – необходимо время на образование комплекса (15), который является ключевым и отвечает за гидрометаллирование олефинов. Далее включается каталитический цикл. При увеличении начального количества катализатора увеличивается скорость межлигандного обмена, что ускоряет образование (15) для которого включается каталитический цикл, что приводит к уменьшению индукционного периода реакции гидроалюминирования олефинов с ДИБАХ.

Реакция димеризации комплекса (15) (E=7,02 ккал/моль) и взаимодействие его же с очередной молекулой ДИБАГ (E=7 ккал/моль) практически равновероятны. Но для гидроалюминирования олефинов с ДИБАХ ДИБАГ не является исходным веществом, а образуется в ходе реакции, в малых количествах. Поэтому более вероятна димеризация комплекса (15), в результате чего образуется комплекс (6). Соединение (15) находится в мономер-димерном равновесии с бициркониевым комплексом (6). Комплекс (15) и (6) могут взаимодействовать с молекулой ДИБАГ с образованием тригидридного комплекса (7).

Принимая во внимание активационные барьеры реакций взаимодействия комплекса (12) с ДИБАГ (E=6,62 ккал/моль) и ДИБАХ (E=26,7 ккал/моль), можно сделать вывод, что более вероятно протекание реакции Cp₂ZrHCl с HAlBu₂ⁱ, что приводит к образованию мономера Cp₂ZrH₂·ClAlBu₂ⁱ.

Гидрометаллирующая активность комплексов уменьшается в ряду (12)(E=6,6 ккал/моль)>(15)(E=7 ккал/моль)>(7)(E=26 ккал/моль). Тогда в каталитическом процессе скорость будут определять эффективная концентрация и взаимное отношение компонентов (12), (15), (7).

Учитывая активационные параметры реакций, первоначально происходят превращения (12)→(15)→(7). Однако комплекс (7) обладает меньшей активностью, так как активационный барьер его реакции с олефином значительно выше (E=26 ккал/моль) относительно (12) (E=6,6 ккал/моль) и (15) (E=7 ккал/моль). Ранее, в первой детализации, комплекс (7) считался неактивным по отношению к олефинам.

5.      Математическое моделирование индукционного периода

На графиках, представленных на рис. 9, 10, наблюдается индукционный период. Для определения характера индукционного периода был проведен вычислительный эксперимент. Взяв за эталонные значения констант скоростей при X_kt =  0,3 ммоль (табл. 7), был определен индукционный период для различного начального количества катализатора. Исследования проводились для следующих значений количества катализатора, ммоль: 0,01<=X_kt<=0,09 с шагом 0,01; 0,1<=X_kt<=0,9 с шагом 0,1; 1<=X_kt<=5 с шагом 0,5.

Для любого начального количества катализатора (X_kt) константы скоростей определяются следующим образом:

Зная константы скоростей, можно определить концентрацию наблюдаемого вещества X₃  в любой момент времени (t_l), решив прямую задачу, т.е. задачу Коши для системы обыкновенных нелинейных дифференциальных уравнений с заданными начальными условиями с фиксированными и неизвестными . Обозначив зависимость концентрации вещества X₃ от времени через x(t), получим

где t_l – время, с.

Исходя из определения индукционного периода [8] и производной функции, можно считать, что индукционный период функции равен t, если

 или

Как видно из рис. 12, имеет место высокая параметрическая чувствительность продолжительности индукционного периода от начального количества катализатора. С уменьшением количества катализатора закономерно увеличивается индукционный период. Количество катализатора, большее 2 ммоль, несущественно влияет на индукционный период.

На рис. 13 показана скорость расходования наблюдаемого вещества в зависимости от количества катализатора. Здесь

Из рис. 13 следует, что с увеличением количества катализатора скорость расходования наблюдаемого вещества (X₃) увеличивается. Следуя определению индукционного периода, амплитуда полученных графиков отображает Δl. Чем больше количество катализатора, тем меньше Δl, то есть тем меньше продолжительность ускорения расходования наблюдаемого вещества.

На рис. 14 представлена зависимость максимальной скорости расходования наблюдаемого вещества Z_max от T_ind. С увеличением длительности индукционного периода величина Z_max уменьшается по гиперболическому закону.

Вычислив коэффициенты этого уравнения, прологарифмировав его и проведя расчеты для двух значений количества катализатора, можно определить значение Z_max для любого T_ind и наоборот (рис. 15).

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

6.      Заключение

В ходе исследования выявлен внутренний параллелизм задачи построения механизма реакции каталитического гидроалюминирования олефинов, на основе которого предложен трехуровневый метод распараллеливания вычислительного процесса. Разработана информационно-аналитическая система обратных задач, включающая базы данных детализированных и обобщенных стадий. Построена кинетическая модель частных реакций гидроалюминирования олефинов, по которой для первого и второго этапа детализации частных реакций определены константы скоростей и энергии активации. На основе второй детализации частных реакций определены кинетические константы общего механизма гидроалюминирования олефинов с ДИБАХ, в том числе для экспериментов с индукционным периодом. Определен характер индукционного периода для общей реакции гидроалюминирования олефинов с Д ИБАХ. Установлено и количественно оценено увеличение индукционного периода с уменьшением концентрации катализатора. Определены максимальные значения скорости расходования наблюдаемого вещества в зависимости от количества катализатора.

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

Литература

1. Слинько М.Г. Основы и принципы математического моделирования каталитических процессов. Новосибирск: Ин-т катализа им. Борескова СО РАН, 2004. 488 с.

2. Спивак С.И., Губайдуллин И.М., Вайман Е.В. Обратные задачи химической кинетики. Уфа: РИО БашГУ, 2003. 110 с.

3. Парфенова Л.В., Печаткина С.В., Халилов Л.М., Джемилев У.М. Исследование механизма гидроалюминирования олефинов алкилаланами, катализируемого Cp₂ZrCl₂ // Изв. РАН, Сер.хим. 2005. Т. 2. С. 311-322.

4. Крылов О.В. Гетерогенный катализ. Москва: ИКЦ «Академкнига», 2004. 679 с.

5. Губайдуллин И.М., Спивак С.И. Информационно-аналитическая система обратных задач химической кинетики // Системы управления и информационные технологии. 2008. ╧1.1(31). С. 150-153.

6. Панкратьев Е. Ю., Тюмкина Т. В., Халилов Л. М., Парфёнова Л. В., Хурсан С. Л., Джемилев У. М. DFT моделирование образования Al,Zr-комплекса как ключевого интермедиата в реакции каталитического гидроалюминирования олефинов // Сборник тезисов V всероссийской конференции «Молекулярное моделирование». Москва, 2007. С. 84.

7. И.М. Губайдуллин, Ю.Б. Линд. Информационно-аналитическая система решения обратных задач химической кинетики на основе современных высокопроизводительных вычислений // Вестник Омского университета. 2010. ╧ 4(58). С. 137-146.

8. Хилько А.В., Спивак С.И., Губайдуллин И.М., Парфенова Л.В. О математическом моделировании индукционного периода химических реакций // Системы управления и информационные технологии. 2008. ╧1.2(31). С. 264-267.