Семенов С.А. Планирование эксперимента в химии и химической технологии. Том 3, страница 3

Уравнение (6) адекватно эксперименту: Воспользовавшись свойствами ортогональности использованного плана эксперимента, определим среднеквадратичные ошибки полученных констант: 135 10 0.703 ~ь, Лк, 135 10 0.238 ~ь, Лу,> (7) ~ь, К 135 10 0.00198 23 4.186 Лк, 1де 0.00008 ' ~Око =135 Из выражений (7) следует, что ошибки констант зависят от интервалов варьирования факторов. За счет расширения интервалов варьирования можно уменьшить ошибки определения констант. Рассмотренный метод линеаризации кинетических уравнений приводит к нарушению условий максимума правдоподобия и получению смещенных оценок для констант, так как константы определятся из условия минимума квадратичной формы: Ф, =,') (у., — у.,) =~ (18%, — 184,), анена Кроме того, зависимость (1) достаточно формальна. Опираясь на стадийный механизм реакций, было предложено определять константы скоростей каждой стадии из условия минимума квадратов откпонений [7): 20 аачи.

пйп1.ги~е-1 ~Ьгаг т и Ф, =,'~ ,'à —,(с„, .— Е,,) 3=1 1=1 Кс. (8) ~г ' С кз В+А->0 †>Е система кинетических уравнений имеет вид: Ж, = = (Е, + 1с, ИА1 + 1~,ГАПВ1, 4А1 от где и — число опытов; гп — число компонентов; я, - ошибка воспроиз- 3 водимости определения концентрации 1-го компонента. Установлено, что при определении концентрации веществ без систематической ошибки оценки констант, минимизирующие квадратичную форму Фз, будут несмещенными. Вычисление концентраций ~, производится или на основе интегральной формы кинетического уравнения, или численным интегрированием системы кинетических уравнений. Минимум квадратичной формы Фз определяют каким-либо методом нелинейного программирования.

Сходимость решения зависит от начального приближения. Для реакций, кинетика которых описывается системой обыкновенных дифференциальных уравнений с правыми частями в виде полиномов по концентрациям реагирующих веществ, предлагается для получения оценок констант пользоваться дифференциальными уравнениями, а не их интегральными формами (8~. Например, для последовательно-параллельной реакции вида 4с! А,2 Й %4 ~з! Апв! 1с4ГО! 4Р! Й ж = =к 1Щ, с1[е! й Из ~9) следует, что уравнения кинетики для рассматриваемого класса реакций в общем виде соответствуют квадратичному полиному, обычно применяющемуся для описания почти стационарной области: 1 Ь + ХЬР3 + Хьщ~ ~ + Хьз~ в котором независимыми факторами х являются концентрации реагирующих веществ, а откликом у — скорость образования того или иного компонента.

Константы скоростей реакций первого порядка интерпретируются при таком подходе как линейные эффекты Ь;, константы скоростей смешанного второго порядка — как эффекты взаимодействия Ь„; и константы скоростей реакции второго порядка — как квадратичные эффекты Ь„-. Если схема реакции заранее неизвестна, но известны все те вещества, которые могут участвовать в реакции, то план эксперимента составляется таким образом, чтобы по результатам опытов можно было найти независимые оценки всех коэффициентов полиномов второго порядка. Далее проводится проверка значимости коэффициентов и определяется схема реакции.

При протекании реакции концентрации действующих веществ непрерывно меняются, поэтому при планировании эксперимента в качестве факторов берут начальные концентрации реагентов, а в качестве отклика — начальные скорости реакции. Для определения по- ииилтий1. го~в-! ~ Ьгаг следних обычно проводят графическое или численное дифференцирование начальных участков кинетических кривых. Пример 2 ~81. Изучалась кинетика реакций каталитического окислительного дегидрирования бутенов в дивинил.

На основании предварительных исследований был принят механизм реакций, в котором исключены обратные превращения дивинила в бутилены. 4А1 = — (К, +К, +К,)[А]+К,[В1+К4[С1, а 4В! (К1 + К5 + КЗ)[В1+ КЛА1+ Кб[С1. й~ [10) 4С1 — — (К4 + Кб+ К9)[С1+ Кз[А1+ К [В1. дг 401 — К7[А1+ Кз[В1+ К9[С1, ог где К; ([=1, 2, ..., 9) — константы скоростей псевдомономолекулярных реакций; [А1, [В1, [С1, [01 — концентрации реагентов в газовой фазе.

А, В, С, 0 — бутен-1; транс-бутен-2; цис-бутен-2, дивинил соответственно. При изучении кинетики реакций окислительного дегидрирования был установлен нулевой порядок по кислороду и первый порядок по бутиленам; специальными опытами было доказано отсутствие тормозящего влияния продуктов реакции. Кинетика реакций представлена системой дифференциальных уравнений: Таким образом, необходимо определить по экспериментальным данным константы скоростей девяти реакций и показать, что система (10) адекватно описывает кинетику реакций. Решение.

Обозначим начальные концентрации веществ А, В, С через к1, к2, лз соответственно. Рассматривая скорость накопления веществ как отклики у„, а начальные концентрации как независимые факторы, получим из (10) (11) у -~зо +~зз хз+~з~ х +~зз хз* а о„ Н '-О О ц=1, 2, 3,4. Переменные х„, х2, х, связаны с начальными концентрациями обычным линейным преобразованием о к — г. Х 3 3 3 Свободный член в уравнении (11) связан с линейными членами уо ~о ~о з с о «~ -н з~ .н 2 3 Кинетику дегидрирования бутенов исследовали при различных температурах. Для каждой температуры составлялась своя матрица планирования. Рассмотрим кинетический эксперимент при Т= 669 К.

Область исследования приведена в таблице. л — концентрация бутенов, мол. доли. Для каждой функции отклика (10) нужно найти только три коэффициента. Поэтому при планировании достаточно взять! 2 от ПФЭ 2'. Здесь, однако, использован ПФЭ типа 2 (табл.6), так как это позволя- ет точнее определить константы скоростей реакций и более подробно яуза.гп~й1.го/е-1 ~Ьгаг анализировать отклонения величин, предсказанных уравнениями от опытных данных. Для определения дисперсии воспроизводимости некоторые опыты дублировали. Таблица 6. Матрица планирования для исследования кинетики дегидрирования бутенов Скорость накопления веществ у (мол. доли/с) определили по форму- лам: [А!, — [А10 у~ 7 Р1, — [В!о У2 [01, — [Ф1о уд= Здесь индекс «О» относится к начальным концентрациям бутенов, а индекс т - к концентрациям, соответствующим времени т.

Таким обра- 25 мвимлтп1Ы.ги!е-ИЬгаг зом, начальные участки кинетических кривых были аппроксимированы прямой линией. При этом вычисляемые оценки скоростей накопления реагентов будут смещенными, причем смещения зависит от скорости реакции. Чем более активно в химическом отношении реагирующее вещество, тем большее смещение будут иметь соответствующие константы скорости.

В таблице приведены значения констант и их ошибок, полученные по данным планированного эксперимента и методом нелинейных оценок (МНО). Сравнение констант скоростей с их ошибками показывает, что ряд скоростей не выделяется на Фоне шума. Для уменьшения ошибок констант необходимо увеличить интервалы варьирования. Оценки полученных констант были уточнены методом нелинейных оценок (МНО). Согласно этому методу константы скоростей реакций должны быть подобраны таким образом, чтобы была минимальной сумма квадратов отклонений (8). Концентрации Е,, получены интегрированием системы (10) от 1=0 до 1=т при начальных условиях, приведенных выше. Суммирование поводилось по всем опытам, причем слагаемые входили с разными весами, так как было доказано, что ошибки воспроизводимости концентраций всех веществ однородны.

В качестве начального приближения были использованы константы, определенные по плану. Затем по критерию Фишера была проверена адекватность математической модели (10) эксперименту: 2 Ф (1С) 8.4 — 9 2Ь явим. пий1. гц~е-! ~ Ьгаг Сопоставление значения остаточной дисперсии с дисперсией воспроизводимости дает величину Р-отношения, близкую к единице, что указывает на адекватность системы (10). Литература 1. Семенов С.А.

Планирование эксперимента в химии и химической технологии. Учебно-методическое пособие. М,: ИПЦ МИТХТ, 2001 г.- 93 с. 2. Семенов С.А. Планирование эксперимента в химии и химической технологии. Часть 2. Учебно-методическое пособие. М.: ИПЦ МИТХТ, 2005 г.- 29 с. 3. Ахназарова С.Л., Кафаров В.В. Оптимизация эксперимента в химии и химической технологии. М.: Высш. школа, 1978. - 319 с. 4. Зайегйвайе Р.Е.