Оптимизация эксперимента в химии и химической технологии
..pdfТабличное значение критерия Стьюдента /0,05 (12) =2,18
_ 2 I9i8 6 _ |
2- 8,86 |
|
t\ — |
— 3,о > 2,18, |
= 3,1 > 2,18 . |
|
5,6 |
5,6 |
Таким образом, оба коэффициента значимы.
Далее снимают эффекты значимых факторов для выявления более слабых эффектов. Для этого вычитают ^значения 2bj из всех значений у, для которых фактор х, находился на уровне +1. Так, при снятии эффекта фактора вычли значение 2-9,86 из у ъ У а , Уб , у 8, t/ю, у \ 2, г/н. У \б , а при снятии эффекта фактора Х в вычли значе ние 2-8,86 из уа , У 5, Уб, y i , У % у ю , У 1з, У 16- Полученные результаты значений у1 приведены в табл. 57. Затем всю обработку опытных данных повторяют, начиная с построения диаграммы рассеяниядля откорректированных значений у. Критерием для окончания отсева служит дисперсионное отношение
где б-2Ср— дисперсия результатов относительно среднего арифмети ческого этих результатов.
В рассматриваемом примере пришлось повторять весь цикл семь раз, выделенные эффекты и их оценки приведены в табл. 60.
Т а б л и ц а 60 Сводная таблица результатов отсеивающего эксперимента
Номер |
Выделен |
Оценка |
Дисперсия |
Номер |
|
Выделен |
Оценка |
Дисперсия |
|
коэффици |
|
коэффици |
|||||||
этапа |
ные |
ентов рег |
ср |
этапа |
|
ные |
ентов ре-г- |
s2 |
|
|
эффекты |
рессии |
|
|
эффекты |
ресии |
ср |
||
I |
*1 |
9.86 |
214,6 |
V |
|
|
—0,48 |
|
|
|
8.86 |
*13 |
2,84 |
||||||
|
X6 |
|
|
|
+0,55 |
||||
И |
х\ Хз |
3,14 |
|
|
|
*2 |
Х5 |
|
|
27,35 |
VI |
|
|
|
0,55 |
|
|||
|
х 6 х7 |
2,53 |
|
Хз |
|
2,08 |
|||
|
|
|
|
||||||
III |
XQ Х 5 |
—2,31 |
10,45 |
VII |
|
X 7 |
—0,25 |
|
|
IV |
X6 Х ю |
1,0 |
3,89 |
|
1,73 |
||||
|
|
хд |
|
- 0 ,4 |
|||||
|
|
|
|
|
|
|
|
||
Дисперсия воспроизводимости 5 2воспр—0,48, /в о сп р —-5. Вычислен |
|||||||||
ное значение /•’-критерия |
/=• = 1,73/0,48 = 3,62 |
|
|
|
|
||||
|
|
|
|
|
|
|
|||
меньше табличного /'0,95 (15,5) =4,6. |
|
дисперсиями <S2Cp |
|
||||||
Следовательно, различие между |
И S2uocnp |
||||||||
незначимо. |
|
|
|
|
|
|
|
|
|
Ручная обработка результатов в методе случайного баланса |
|||||||||
чрезвычайно трудоемка. Предложен |
алгоритм |
обработки |
резуль |
татов случайного баланса на ЦВМ, так называемая «ветвящаяся стратегия» [29]. Разработан алгоритм для выделения наибольших эффектов по диаграммам расс.еяния. Этот этап не вносит ничего но вого по сравнению с ручной обработкой. Для количественной оцен
Если зависимость константы скорости от температуры подчиня ется закону Аррениуса, то имеем
К = К0ехр( — - J r ) , |
(V. 162) |
где Ко — предэкспоненциальный множитель; |
Е — энергия |
актива |
|||
ции; R — газовая постоянная; Т — абсолютная температура. |
|
||||
Логарифмирование уравнения |
(V. 161) |
с учетом (V.162) дает |
|||
|
|
TTL |
|
|
|
I n l T = i n K o - - J r |
+ |
^ ^ l n |
[ A |
/]. |
(V.163) |
Обозначив переменные следующим образом: |
|
||||
In W = y , |
|
In K Q |
|
|
|
In 1/7'=Z 1, |
|
— E j R = a b |
|
||
In [A ]= Z 2) |
|
Pl = a2 > |
|
||
In[An] —Ал+Ь |
|
Pm |
|
|
|
получим линейное уравнение |
|
|
|
|
|
/72 +1 |
|
|
|
||
У = ао+ |
2 |
aigi- |
|
|
(V.165) |
/ - 1
Для определения коэффициентов уравнения (V.165) методом планирования экспериментов можно использовать линейные орто гональные планы с числом опытов N ^ m + 2 .
Пример 12 [33]. Исследовалась кинетика процесса сополимеризации а, (3, p-трифторстирола с метакриловой кислотой при небольших степенях превращения исходных мономеров в полимер. Уравнение формальной кинетики процесса рас сматривалось в следующем виде:
|
|
W=K[Ci + C2]ni [И]п\ |
|
|
( V . 166) |
||
где [Ci] |
и [С2] — концентрации исходных мономеров, |
моль/л; |
[И] — концентра |
||||
ции инициатора (перекиси бензоила), моль/л; щ |
и п 2 — порядок реакции по сум |
||||||
марной |
концентрации |
мономеров |
и инициатору; |
Е — энергия |
активации, |
||
кДж/моль; К — суммарная константа скорости реакции. |
|
|
относительно |
||||
Уравнение (V. 166) |
позволяет |
определить |
порядок реакции |
||||
суммарной концентрации мономеров |
и инициатора, энергию |
активации и кон |
станту скорости реакции. Константа скорости реакции К зависит от температуры по уравнению (V.162). Б эксперименте суммарная концентрация исходных моно меров изменялась в диапазоне 8,08+-11,18 моль/л, концентрация инициатора — в диапазоне 0,0504-4-0,1512 моль/л, температура реакции 60+-800 С.
Р еш е н и е . Логарифмируя уравнения (V.166) с учетом (V.162), получим
Ig W = !g /Со + щ lg [Cl + С2] + Л2 Ig [И] — j \ g e j - |
(V. 167) |
|
|
(* ? )Ч |
t { z l f s l |
0\2_2 |
|
|
|
f (*°) S*s |
|||
|
|
sl + |
|
Д^2 |
A*S |
|
|
A«f |
|||
|
|
(0 ,977)2 |
( — 1,0591)2 |
(0 ,00292)2 |
|
|
«1,35 j / " 1 + (0 ,0703)2 + |
(0 ,238)2 |
(0 ,00008)2 = 0 ,53. |
||
Из |
(V. 171) следует, |
что ошибки |
констант зависят от интервалов варьиро |
||
вания |
факторов. За счет |
расширения интервалов варьирования можно умень |
|||
шить ошибки определения констант. |
|
|
|
Рассмотренный метод линеаризации кинетических уравнений приводит к нарушению условий максимума правдоподобия и полу чению смещенных оценок для констант, так как константы опреде ляются из условия минимума квадратичной формы:
®i = S {У1 - h i ) 2 = 2 (ig w t - Ig^ )2 ,
i-l /=1
а не из
®2 = | ] /=1
Кроме того, зависимость (V.161) достаточно формальна. Опи раясь на стадийный механизм реакций, был© предложено опреде лять константы скоростей каждой стадии из условия минимума квадратов отклонений [34]:
тп
Фз = 2 |
2 " Т " icjl ~ ®л)2- |
(v -172> |
)Z\ |
Sc) |
|
где п — число опытов; т — число компонентов; sC]— ошибка воспро
изводимости определения концентрации /-го компонента. Установлено, что при определении концентраций веществ без
систематической ошибки оценки констант, минимизирующие квад ратичную форму Фз, будут несмещенными. Вычисление концентра
ций Сц производится или на основе интегральной формы кинетиче ского уравнения, или численным интегрированием системы кине тических уравнений.
Минимум квадратичной формы Фз определяют- каким-либо ме тодом нелинейного программирования. Сходимость решения зави сит от начального приближения.
Для реакций, кинетика которых описывается системой обыкно венных дифференциальных уравнений с правыми частями в виде полиномов по концентрациям реагирующих веществ, предлагается для получения оценок констант пользоваться дифференциальными уравнениями, а не их интегральными формами [35]. Например, для последовательно-параллельной реакции вида
С
(V. 173)
Е
система кинетических уравнений имеет вид:
d [В]
W2 = - ^ L = k l [ A ] 2 - k 3 [A] [В],
(V. 174)
^ | = £ ^ - = Лз[А] [в1—
Из (V.174) следует, что уравнения кинетики для рассматривае мого класса реакций в общем виде соответствуют квадратичному полиному, обычно применяющемуся для описания почти стацио нарной области:
У.= h + 2 ъ)х ) + 2 buJx ux ) + |
2 ь11х)> |
(V.175) |
|
} |
и<) |
) |
|
в котором независимыми |
факторами х |
являются |
концентрации |
реагирующих веществ, а откликом у — скорость образования того или иного компонента. Константы скоростей реакций первого по рядка интерпретируются при таком подходе как линейные эффекты bj, константы скоростей смешанного второго порядка — как эффек ты взаимодействия bUj и константы скоростей реакции второго по
рядка— как квадратичные эффекты — |
Ьц. Если |
схема реакций |
заранее неизвестна, но известны все те |
вещества, |
которые могут |
участвовать в реакции, то план эксперимента составляется таким образом, чтобы по результатам опытов можно было найти незави симые оценки всех коэффициентов полиномов второго порядка. Далее проводится проверка значимости коэффициентов и опреде ляется схема реакции.
При протекании реакции концентрации действующих веществ непрерывно меняются, поэтому при планировании эксперимента в качестве факторов берут начальные концентрации реагентов, а в качестве отклика — начальные скорости реакции. Для определения последних обычно проводят графическое или численное дифферец-
цирование начальных участков кинетических кривых или обраба тывают данные, полученные на проточно-циркуляционных уста новках.
Пример 13 [35]. Изучалась кинетика реакций каталитического окислительногодегидрирования бутенов в дивинил. На основании предварительных исследований был принят механизм реакций, в котором исключены обратные превращения дивинила в бутилены,
А, В, С, D — бутен-1; транс-бутен-2; цис-бутен-2, дивинил соответственно.
При изучении кинетики реакций окислительного дегидрирования был уста новлен нулевой порядок по кислороду и первый порядок по бутиленам; специ альными опытами было доказано отсутствие тормозящего влияния продуктовреакций. Кинетика реакций представлена системой дифференциальных уравнений:
|
|
at |
= |
— (fix + |
+ |
£7) [А] + |
^2 [В] + |
£4 [С], |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
d [В] |
= ' — (fii + |
£5 + |
^в) И |
+ |
^1 [А] + |
^6 [С], |
|
||||||
|
|
“ |
|
|
|||||||||||
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
(V. 176) |
|
|
|
d ГС] |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
= |
— (^4 + |
+ |
£9) [С] + |
k3 [А] + |
ks [В], |
|
|||||||
|
|
~ |
|
|
|||||||||||
|
|
a t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- ^ |
' = А |
7 [А] + |
А8 [В ]+ А 9 [С]> |
|
|
|
|
|||||||
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
где ki (t= l, 2, |
. . . , |
9) — константы скоростей |
псевдомономолекулярных |
реакций; |
|||||||||||
[А], [В], [С], |
[D] — концентрации реагентов в газовой фазе. |
|
|
||||||||||||
Таким образом, необходимо определить по экспериментальным данным кон |
|||||||||||||||
станты скоростей |
девяти |
реакций |
и показать, |
что система |
(V.176) адекватна |
||||||||||
описывает кинетику реакций. |
начальные |
концентрации веществ А, В, С через |
|||||||||||||
Р е ш е н и е . |
Обозначим |
||||||||||||||
?i, z2, z3 соответственно. Рассматривая |
скорости накопления веществ как откли |
||||||||||||||
ки уUt а начальные концентрации как независимые факторы, получим из |
(V.176) |
||||||||||||||
|
Ум = |
^ои “Ь ^1их \ "Ь ^2их 2 + |
и-*з» |
|
ц = 1, |
2, |
3, 4. |
(V. 177) |
|||||||
Переменные хи х2, Хз связаны с начальными концентрациями обычным ли |
|||||||||||||||
нейным преобразованием |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
*) = (*]—tylbzj. |
|
|
|
|
|||||
Свободный |
член в уравнении |
(V.177) |
связан |
с линейными членами |
|
||||||||||
|
|
|
ьои —^1и |
Д*1 + ^0 |
4 |
+ |
зи Д*з |
|
(V .178) |
||||||
|
|
|
Дг2 |
|
|||||||||||
|
|
|
|
|
Кинетику дегидрирования бутенов исследовали при различных температу рах. Для каждой температуры составлялась своя матрица планирования. Рас смотрим кинетический эксперимент при Г=669 К. Область исследования приве дена в таблице.
т
*5 |
0,01500 |
0,01158 |
0,00993 |
|
0,00500 |
0,00386 |
0,00331 |
||
Azj |
■z— концентрация бутенов, мол. доли.
Для каждой функции отклика (V.176) нужно найти только три коэффици ента. Поэтому при планировании достаточно взять V2 от ПФЭ 23. Здесь, однако,
яспользован ПФЭ типа 23 (табл. 62), так как это позволяет точнее определить константы скоростей реакций и более подробно анализировать отклонения вели чин, предсказанных уравнениями от опытных данных.
Для определения дисперсии воспроизводимости некоторые опыты дубли ровали.
|
|
|
|
|
|
|
Т а б л и ц а 62 |
|
|
Матрица планирования для исследования кинетики |
|
||||
|
|
|
дегидрирования бутенов |
|
|
||
Номер |
X t |
*2 |
Хг |
Ух |
У* |
У8 |
У*. |
опыта |
|||||||
1 |
—1 |
—1 |
—1 |
—0,288 |
0,046 |
0,003 |
0,272 |
2’ |
—1 |
—1 |
+ i |
—0,263 |
0,077 |
—0,102 |
0,255 |
3. |
—1 |
+1 |
—1 |
—0,286 |
—0,124 |
0,041 |
0,257 |
4 |
+ i |
—1 |
- л |
—0,737 |
0,068 |
0,033 |
0,370 |
5 |
+ i |
+ 1 |
+1 |
—0,678 |
—0,063 |
—0,060 |
0,426 |
6 |
+ i |
+ 1 |
—1 |
—0,725 |
- 0,043 |
0,058 |
0,342 |
7 |
+ i |
—1 |
+ i |
—0,705 |
0,124 |
—0,130 |
0,356 |
8 |
—1 |
+ 1 |
+1 |
—0,276 |
—0,139 |
—0,180 |
0,262 |
Скорость |
накопления |
веществ |
у (мол. |
доли/с) определяли по формулам: |
У1 = |
[А]Т- [ А ] 0 |
[B ],-[B J 0 |
(V. 179) |
------- 1-------- . |
У2 = ----------------- |
||
|
тз |
тз |
|
Уз = |
[C]t - [С]0 |
[D]x— [D]0 |
|
|
^ 4 = ------ :----- * |
|
Здесь индекс «О» относится к начальным концентрациям бутенов, а индекс т — к концентрациям, соответствующим времени т. Таким образом, начальные участки кинетических кривых были аппроксимированы прямой линией. При этом вычис ляемые оценки скоростей накопления реагентов будут смещенными, причем сме щения зависят от скорости реакции. Чем более активно в химическом отношении реагирующее вещество, тем большее смещение будут иметь соответствующие константы скорости. В таблице приведены значения констант и их ошибок, полу ченные по данным планированного эксперимента и методом нелинейных оценок
(МНО).
|
|
|
|
К о н с т а н т ы с к о р о с т е й , с -1 |
|
|
|
||
|
кг |
л2 |
кг |
к* |
кь |
кг |
*7 |
кг |
к9 |
По плану |
7,128 |
0,74 |
6,62 |
5,60 |
3,47 |
2,57 |
10,80 |
1,08 |
2,90 |
По МНО |
9,07 |
0,75 |
9,38 |
6,51 |
3,64 |
2,58 |
29,59 |
. 1,11 |
3,10 |
s(kj |
1,8 |
ч 2,6 |
1,8 |
3,5 |
2,6 |
3,5 |
1,8 |
2,6 |
3,5 |
Сравнение констант скоростей с их ошибками показывает, что ряд констант не выделяется на фоне шума. Для уменьшения ошибок констант необходимо увеличить интервалы варьирования. Оценки полученных констант были уточнены методом нелинейных оценок (МНО). Согласно этому методу константы скоро стей реакций должны быть подобраны таким образом, чтобы была минимальной
сумма квадратов отклонений (V.172). Концентрации Cji получены интегрирова нием системы (V.176) от t= 0 до t.=т при начальных условиях (см. таблицы на с. 248). Суммирование проводилось по всем опытам, причем слагаемые входили с равными весами, так как было доказано, что ошибки воспроизводимости концент раций всех веществ однородны. В качестве начального приближения были исполь зованы константы, определенные по плану. Затем по критерию Фишера была проведена адекватность математической модели (V.176) эксперименту:
2 ___ Ф 3 ( ^ ) m ln
S°CT- 8-4 -9 ■
Сопоставление значения остаточной дисперсии с дисперсией воспроизводи мости дает величину ^-отношения, близкую к единице, что указывает на адек ватность системы (V.176).
ГЛАВА VI
ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТА ПРИ ИЗУЧЕНИИ ДИАГРАММ СОСТАВ — СВОЙСТВО
1. Метод симплексных решеток. При изучении свойств смеси, зависящих только от соотношений компонентов, факторное прост ранство представляет собой правильный (q—1)-мерный симплекс. Для систем выполняется соотношение
|
я |
* |
|
2 * / = !, |
(VI. 1) |
|
/-1 |
|
где Хг^О — концентрация |
компонента; |
q — количество компо |
нентов. |
|
|
Для двухкомпонентных систем симплекс — прямая: содержание |
||
компонентов определяется |
соотношением |
отрезков. При q= 3 пра |
вильный симплекс — равносторонний треугольник. Каждая точка треугольника отвечает одному определенному составу тройной системы и, наоборот, каждый состав представляется одной опреде ленной точкой. Состав может быть выражен в мольных, весовых или объемных долях или процентах. Вершины треугольника соот ветствуют чистым веществам, стороны — двойным системам. Опус-