Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Планирование эксперимента при оптимизации процессов химической техн

..pdf
Скачиваний:
15
Добавлен:
15.11.2022
Размер:
9.83 Mб
Скачать

абсолютных значений коэффициентов В п и В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,

х^у

а выход­

ной переменной у\.