Планирование эксперимента при оптимизации процессов химической техн
..pdfабсолютных значений коэффициентов В п и В22 |
измене |
|||
ние выходной переменной по осям |
Х\Х2 будет |
различ |
||
ным; |
в, г — один из коэффициентов |
близок или |
равен |
|
нулю; |
один из коэффициентов (В2 2 ) |
равен нулю, |
центр |
фигуры находится на бесконечности; поверхность от
клика представляет собой возрастающее |
возвышение — |
||||
« г р е б е н ь » — (г); |
возможен также |
случай, когда |
|||
В22 « 0 и центр фигуры находится |
в любой |
точке |
оси |
||
Х2) поверхность отклика представляет собой |
« с т а ц и о |
||||
н а р н о е в о з в ы ш е н и е » — (в). |
На |
практике |
слу |
||
чай (в) встречается редко. |
|
однако они |
по |
||
Указанные случаи |
идеализированы, |
могают исследователю ориентироваться в структуре по верхности отклика. Основные типы поверхностей откли ка для двух переменных приведены в табл. 2.9.
Т а б л и ц а 2.9. Типы поверхностей отклика канонического уравнения двух переменных
« |
S. |
|
Знаки коэффи |
|
|
|
н |
|
|
|
|||
5 |
° |
Коэффи |
циентов |
Геометрическая |
|
|
5 ° |
|
Тип кри |
|
|||
н |
X |
циенты |
|
вых в |
интерпретация |
Центр |
Q.& |
|
|||||
3 |
ш |
|
|
сечении |
|
|
О |
о |
|
В и |
в п |
|
|
X с |
|
|
|
|
|
1 В\\ = В22
2 |
Со |
II Со |
3Вц>В 22
4В\\>В22
5 В\\ = В22
6В\\ = В22
7В\ 1>7^22
8В22 = 0
9В22 = о
— |
— |
Окруж |
Круглая |
Максимум |
|
|
|
ность |
выпуклость |
Минимум |
|
+ |
+ |
То же |
Круглая |
||
|
|
Эллипс |
впадина |
Максимум |
|
|
|
Эллипсоид |
|||
|
|
|
ная |
выпук |
|
|
|
То же |
лость |
Минимум |
|
+ |
+ |
Эллипсоид |
|||
|
|
|
ная |
впадина |
Седловая |
+ |
— |
Гипер |
Симметрич |
||
— |
|
бола |
ное |
седло |
точка |
+ |
То же |
То |
же |
То же |
+» Вытянутое
|
|
|
седло |
|
Нет |
— |
0 |
Прямая |
Стационар |
||
|
|
Пара |
ный |
гребень |
На беско |
|
0 |
Возрастаю |
|||
|
|
бола |
щий |
гребень |
нечности |
Кратко изложим алгоритм канонического преобразо вания уравнения регрессии второго порядка в уравнение (2.103).
Начало координат переносят в новую точку фактор ного пространства 5. Для этого уравнение регрессии вто рого порядка дифференцируют по каждому фактору и приравнивают к нулю. Решая систему уравнений, нахо дят координаты нового центра, подставляют их в урав нение регрессии и получают значение выходной перемен ной Ys в центре факторного пространства 5. В уравнении регрессии исчезают члены первой степени и изменяется значение свободного члена
|
У— Ys = £ ibijxixj + £ lbiix* |
(2.104) |
|||
|
|
, 3= 1 |
1 = |
1 |
|
В новом |
центре |
поворачивают |
оси до |
совмещения |
|
с главными |
осями |
поверхности |
функции |
отклика. Эта |
операция осуществляется по правилам аналитической геометрии.
Например, |
для |
п = 2 характеристическое |
уравнение |
|
имеет вид |
|
|
|
|
№ = |
(Ьп — В) |
0,5ЪХ2 |
(2.105) |
|
|
= 0. |
|||
|
|
0,5^21 |
{р22 '— В) |
|
Решение (2.105) записывают уравнением |
|
|||
где |
|
В2 — а1В + а2 = 0 , |
(2.106) |
|
|
|
|
|
|
01 = |
(&п + |
Ьгг); |
° 2 = (^п^22 — 0,2561 2 )2. |
Два корня этого уравнения дают искомые значения ко эффициентов уравнения в канонической форме.
Для /1 = 3 характеристическая матрица (2.105) до полняется соответствующими строками и столбцами:
( Ь п - В ) |
0,56i2 |
0,56i3 |
|
НВ) = 0,562I |
|
(6 2 2 — В) |
0,5^23 = 0. (2.107) |
0,5631 |
0.5&32 |
(&зз — Я) |
|
Решение (2.107) записывают уравнением |
|||
В 3 — |
й \В 2 -}- (22В — а з = |
0 , |
|
где |
|
|
|
|
ai = |
2 ] Ьц\ |
|
|
|
i = 1 |
|
|
п |
|
|
02 = |
b i i b j j — 0,25 |
ьД; |
i, 3= 1 |
г<3 |
n |
n |
bij |
n |
Яз == JJ Ьц “|- 0,25 |
0,25 babiq. |
||
i= 1 |
\ 3= i |
i, 3= 1 |
|
|
|
|
гФз^Ч |
Правильность расчетов проверяют по формуле |
|||
г= 1 |
i—1 |
(2.108) |
|
|
|||
Пример. Произвести канонический анализ уравнения |
|||
регрессии |
|
|
|
^ = 85,14 + 3,43*1 — 1,32*2 + |
2,60*2— 1,19*2 + 3,00*!*2. |
Решение. Начало координат переносим в новую точ ку факторного пространства следующими действиями:
— ±- = 3,43 + 2-2,60*1 + 3,00*2 = 0;
дх\
Л
'ду
— = — 1,32 — 2-1,19*2 + 3,00*i = 0;
дх2
*is = —0,197; |
*2з = —0,802. |
Подставляем значения *is и x2s в уравнение регрес сии и получаем значение выходной переменной Ys=85,33 в центре S.
Для поворота осей в новом центре надо знать значе ния канонических коэффициентов В п и В22. Для вычис ления этих значений приравняем характеристический де терминант к нулю:
(2,60 — В) |
0,5-3,00 |
0,5-3,00 |
(— 1,19 — Б) |
или воспользуемся формулами (2.105)
щ = 2,60— 1,19 = — 1,41;
а2 = —2,60-1,19 — 0,25(3,00)2 = —5,34.
Получим уравнение
В2 — 1,41В — 5,34 = 0.
Корнями этого уравнения будут 5ц = 3,2, В22 = 1,71. Условие (2.108) выполняется:
2,60— 1,19 = 3,12— 1,71.
Уравнение регрессии в канонической форме примет вид
У — 85,33 = 3,12Х2— 1,71-Х2 .
1. 2
В соответствии с приведенной выше классификацией контурными ^кривыми в области оптимума являются ги перболы (.Бц < 0, В22 > 0).
Г л а вл 3
ОСОБЕННОСТИ ПРИМЕНЕНИЯ МЕТОДОВ ПЛАНИРОВАНИЯ ЭКСПЕРИМЕНТА ПРИ ИССЛЕДОВАНИИ СЛОЖНЫХ ОБЪЕКТОВ
ХИМИЧЕСКОЙ ТЕХНОЛОГИИ
Введение. В настоящее время очень редко приходит ся сталкиваться с необходимостью изучать обособлен ный технологический процесс. Чаще всего исследуемый объект представляет собой совокупность связанных.про цессов (говорят, что исследуемый объект обладает структурой). Очевидно, применение методов планирова ния эксперимента к таким объектам с одной стороны требует выработки новых рекомендаций и приемов, по зволяющих решить поставленные задачи моделирования и оптимизации, с другой стороны — для реализации экс перимента на сложном технологическом объекте (будем называть его системой), требуется своя система иссле дования, представляющая собой объединение технологи ческой установки, ЭВМ, исследователя и некоторых пра вил действий (или алгоритмов), обеспечивающих дости жение поставленных целей. На рис. 3.1 представлена схема такой системы исследования.
Глобальную задачу исследования как целенаправ ленного процесса можно сформулировать так: в соответ ствии с некоторым заданием, вызванным потребностью
Рис. 3.1. Схема системы исследования сложного объекта.
т
народного хозяйства в определенном продукте, необхо димо получить информацию о способе производства это го продукта (совокупности процессов), достаточную для проектирования химико-технологического комплекса за данной мощности. Решение глобальной задачи в значи тельной степени зависит от методов и средств, имеющих ся в распоряжении исследователя. Будем считать, что исследователь имеет возможность проводить экспери менты на технологических установках определенного масштаба, пользоваться методами теории эксперимента (а значит, и вычислительной машиной). Исследователь ограничен временными, физическими и материальными ресурсами.
Очевидно, решение глобальной задачи исследования с системных позиций — это решение сложной задачи оптимизации по многим критериям в условиях значи тельных ограничений и неопределенности. В настоящее время эта задача не формализована, тем не менее, она успешно решается исследователями с использованием проверенных практикой правил и принципов.
В этой главе рассмотрены особенности применения методов планирования эксперимента к сложным объек там химической технологии, а также правила, объеди ненные некоторой логической схемой достижения постав ленных целей.
3.1. Постановка задач исследования сложных объектов химической технологии
Введение. При исследовании сложных объектов хи мической технологии встречается ограниченное число типов соединения процессов в системы. В работе [38] рассмотрены пять типов: последовательное, последова тельно-обводное (с байпасом), параллельное, обратное (с рециклом) и перекрестное. Учитывая, что с точки зре ния использования методов теории эксперимента и опти мизации системы с байпасом и последовательным сое динением, а также системы с рециклом и перекрестные почти идентичны, ниже рассматриваются только три ти па соединения процессов — последовательное, парал лельное и с рециклом (рис. 3.2).
3.1.1. Исследование системы с последовательным соединением процессов. Исходные данные. Имеется си стема, процессы которой соединены последовательно (см. рис. 3.2, а) технологическими потоками, т. е. систе ма обладает структурой, которую можно выразить усло виями:
У\ = |
*4', |
|
У2= |
*?; |
|
в общем виде для i- и 6-го процессов |
|
|
yii = |
X i . |
(3.1) |
Необходимо найти математическую модель системы, используя методы планирования эксперимента, и решить
задачу оптимизации системы. |
Математическая |
модель |
|
ищется в виде регрессионного уравнения |
|
||
Л |
|
Ь0, Ь\,..., frjv+i), |
(3.2) |
yN = /(*1, _х2, .... * J V + I , |
|||
а задача оптимизации |
|
|
|
R = extr yNi |
(3.3) |
||
|
х |
|
|
где R — экстремальное значение функции цели. |
начиная |
||
План эксперимента. |
К каждому процессу, |
||
с первого, применяется |
соответствующий план |
экспери |
Рис. 3.2. Типы соединения процессов в системах исследования:
а — последовательное; б — параллельное; в — соединение с рециклом.
мента (первого или второго порядка) и определяются модели
уI |
= |
f\(xI, |
*2, |
-«з, б(|)); |
|
У2 = h(x4, |
|
|
3 |
||
х5, |
Хв, 6 (2>); |
||||
|
|
|
|
|
(3.4) |
y N |
--- |
{ N ( X N , |
JCJV + I . |
3 |
|
|
|
|
|
|
При этом найденное значение выходной переменной Ух (при постановке планируемого эксперимента на про цессе 1) является фактором х4 (при исследовании про цесса 2 (см. рис. 3.2, а)). Желательно при этом решить автономную задачу оптимизации для процесса /, т. е. найти
exlr у\ = Ru |
(3.5) |
Х|, х2, ж3
и в качестве нулевого уровня фактора х4 использовать
R\ = extr у\.
Решение задачи. Если не наметится противоречий в физико-химических связях уравнения
extr ук = Xi, |
(3.6) |
то, зафиксировав |
фактор X N = |
extr I J N - U м о ж н о решить |
|
глобальную задачу оптимизации (3.3) . Для |
этого доста |
||
точно найти |
|
|
|
extr ук = |
RN = R При xN = extr yN-i = |
const. |
|
Если физико-химические |
ограничения |
не позволят |
использовать условие (3.6), то для решения задачи (3.3) целесообразно использовать принцип динамического про
граммирования: |
задается |
ун и |
по |
модели |
у$ = |
= f(xN, xN+u |
определяются соответствующие XN и |
||||
xN+i. Ранее задается xN = |
ум-i и решается задача поис |
||||
ка соответствующих факторов для |
блока |
(N — 1). Рас |
|||
четная процедура поиска оптимальных факторов |
(навер |
||||
ное, обратная рассмотренной выше |
прямой процедуре) |
продолжается до тех пор, пока не будут найдены хи *2 , *3. Полученные результаты проверяются на эксперимен тальной установке и, если расхождения невелики, при нимается найденный технологический режим.
Замечание 1. Если в процессе оптимизации техноло гического режима на установке по каким-либо факторам будут достигнуты граничные условия, то они могут быть зафиксированы, и таким образом исключены из процеду ры планируемого эксперимента и поиска оптимальных значений выходных переменных.
Замечание 2. Если процессы в системе или сама си стема характеризуется больше чем одной переменной, то задача усложняется и переходит в область методов многокритериальной оптимизации.
3.1.2. Исследование системы с параллельным соеди нением процессов. Исходные данные. Исследуется систе ма, процессы в которой соединены параллельно общими технологическими потоками (см. рис. 3.2,6). Структура системы может характеризоваться условиями
*0 = |
2 *7; |
(3.7) |
'/о = |
2 > ь |
(з.8) |
|
г=1 |
|
Необходимо найти математическую модель системы, используя методы планирования эксперимента, и ре
шить задачу оптимизации системы. Математическая мо дель ищется в виде
= *(у. Ь«>), (3.9)
где i — 1,2, ..., Я — номер процесса.
Задачу оптимизации можно записать так:
R = extrу0 = |
y^Ri = |
extr V |
yi(x^\ |
j |
6 (0 ). |
(3.10) |
л: |
^ |
^ |
1 |
j |
|
Ялая эксперимента. Если выполняется условие (3.7) и возможно измерение, а также управление факторами *(|\ x(j[4 то задача получения математической
модели системы в виде (3.9) методами планирования эксперимента распадается на отдельные процедуры по числу блоков в исследуемой системе. Таким образом, математическая модель системы, состоящая из парал лельно соединенных процессов, представляет систему линейных и нелинейных регрессионных уравнений
У1 = Ш \ \ |
*(12>, |
.... |
*<!>, |
&(!>); |
|
|
У2 = |
/2 (х(2), |
*<*>, |
.... |
*<J>, .... 6(2)); |
(3.11) |
|
yN = |
fN(x(” )y *(*>, |
.... |
*<">, |
.... 6 W ). |
|
Решение задачи оптимизации. При выполнении усло вия (3.8) задача оптимизации системы также распада ется на автономные задачи поиска оптимума по каждо му процессу:
Ri = |
extr0 i(*(*>, |
*(i), .... |
6 <\)), |
(3.12) |
|
1 |
2 |
j |
|
а далее оптимум |
системы |
находится в |
соответствии |
с условием (3.10).
Замечание. Условие (3.8), как правило, выполняется для сложных объектов, у которых функция* цели у\ сов падает с выходной переменной, характеризующей техно логический режим (например, содержание некоторого продукта А в выходных потоках технологической систе мы). В тех случаях, когда у\ является показателем эф
фективности работы системы (надежность, устойчивость, рентабельность, стоимость и др.), условие (3.8) не вы полняется, и тогда используют взвешенные оценки р* для уг, получаемые экспертизами.
3.1.3. Исследование систем срециклом. Исходные данные. Значительно сложнее решаются задачи модели рования и оптимизации, когда в схеме имеется один или несколько рециклов. Структуру такойсистемы (см. рис. 3.2, в) можно выразить уравнениями:
*<»> = |
*о + *<*>; |
(3.13) |
у2 |
= х^К |
(3.14) |
|
2 |
|
Ставится задача: найти математическую модель системы методами теории эксперимента и определить оптимум у\\
R = |
extr у\ |
|
(3.15) |
|
при |
X |
|
|
|
|
|
|
|
|
У = fl(x0, * (1), |
*(1>, ЬО)); |
(3.16) |
||
|
2 |
3 |
j |
|
У2 = |
Х^, Ы2р . |
(3.17) |
||
Замечание 1. Возможен |
рецикл |
вида |
tj\ = ij2, т. е. |
выходной поток объекта /, характеризующийся перемен ной у 1, направляется на вход того же объекта.
План эксперимента. При наличии рецикла использо вание планов эксперимента для получения регрессион ных моделей невозможно вследствие коррелированное™ факторов между собой (х0 коррелировано с и др.).
Напомним, что одной из предпосылок регрессионного анализа является отсутствие корреляции между факто рами. При исследовании системы с рециклом возможны два случая.
1. Если существует возможность физически «разор вать» рецикл (между у2 и xf® поместить, например, ем
кость, накапливающую продукты выходного потока объ екта 2) у то моделирование сводится к соответствующей задаче последовательного соединения объектов в систе
му. Например для блока 1 (см. рис. |
3.2, в) |
факторами |
|
будут |
= *о + xf& при х№ = const, |
х^у |
а выход |
ной переменной у\.