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

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

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

п/п Наименование

1Коэффициенты регрес­ сии

2Построчные дисперсии

3Критерий Кохрена (расчетное значение)

4Критерий Кохрена (табличное значение)

5Условие однородности.

6Дисперсия воспроизво­ димости (ошибка опы­ та)

7Дисперсия коэффици­ ентов регрессии

8 Критерий Стьюдента

(расчетное значение)

 

 

Расчетные значения

 

 

 

Ьо =

 

90,618;

 

Ь2 =

—2,185

 

 

bi =

 

—1,873

 

Ь3 =

4,798;

—0,075

 

 

 

 

 

 

Ь12 =

 

 

 

 

 

 

bis =

—0,778

приведены в табл. 5.13

Ь23 =

0,445

 

 

 

 

*

 

2

 

 

 

 

 

 

 

 

 

= 9,215

 

 

 

 

 

 

и = I

1,837

 

 

 

 

 

 

Ор

 

 

 

 

 

 

 

 

9,215

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gr =

0,68

(q =

5%, fi =

т — 1

=

=

l,

f2 =

N =

8)

 

 

 

 

 

 

 

2

 

 

 

 

 

 

выполняется; su— однородны

 

 

2

 

1

 

 

1,152

 

 

 

 

So =

-----9,125 =

 

 

 

 

 

 

8

 

 

 

 

 

 

 

s2

=

1,152

 

 

Sb

= 0,268

---------0,072 или

\

 

8,2

 

 

 

l'

 

 

tb0

=

338;

tb2

=

8,15;

h l2 =

0,28

tbl =

6,99

h z =

17,9;

Ц 3 =

2,2

9 Критерий Стьюдента (табличное значение)

10 Условие значимости

11 Дисперсия адекват­

ности

12Критерий Фишера (расчетное значение)

13Критерий Фишера

(табличное значение) 14 Условие адекватности

<*„ = i.66

U = 2,306 (q = 5%, f — N(m —

-1) = 8)

выполняется,

кроме tbl2 <. t-t, h <

< tT и tb23 <

/T; b\2, bi3, b23 — незна­

чимы (из уравнения регрессии исклю-. чаются)

Яад =

3,445; m =

2,

iV =

8;

 

fan =

N — / =

4;

у

вычисляется

по

уравнению

регрессии

со

значимыми

коэффициентами

(/ =

4)

 

 

FP =

3,445

 

2,99

 

 

------ -------=

 

 

FT=

1,152

5%,

fi =

fa« =

 

3,84

(<7 =

 

= N —.l =

4,

f2

=

N(m — 1) =

8)

выполняется; уравнение регрессии адекватно процессу

Наименование

 

Расчетные значения

п/п

 

15

Адекватное

уравнение

А

9 0 ,6 1 8 —

1,873 x i — 2,1 8 5 х2 +

(/ =

 

регрессии

(факторы в

+ 4 ,7 9 8 х3

 

 

кодированной

форме)

 

 

 

16

Статистическая

мате­

у '=

85 ,0 2 4 —

1 ,5 6 1 *1 — 0 ,0 2 2 * 2 +

 

матическая

модель

+ 2 ,3 9 9 * з

 

(факторы в натураль­ ных единицах)

временно способствует протеканию цементационного процесса извлечения ионной ртути. Поэтому такой метод можно назвать экстракционно-цементационным спосо­

бом извлечения металлической и ионной форм ртути из растворй щелочи.

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

результаты этих исследований заключаются в следую­ щем.

яп ^ Повышение температуры раствора едкого натра от и до 120 С приводит к увеличению степени очистки. Ес-

ЛИ S r r b| ЧТ° ПРИ температуре 70°С вязкость экстрагента ЖСС существенно повышается и он может затверде­

вать (Гдл=68°С), то область изменения температуры сле­ дует ограничить интервалом 75—120°С.

2. Интенсивность перемешивания обеих жидких фаз исследуемой гетерогенной системы ЖСС — раствор ед­ кого натра мешалкой в значительной степени определяет кинетику процесса очистки. При вращении мешалки со скоростью около 400 об./мин исходная концентрация ртути в щелочи уменьшается незначительно, при

1800 об./мин — содержание ртути уменьшается более чем в 2,5 раза. Однако скорость вращения мешалки 2500 об./мин способствует ускоренному диспергированию сплава ЖСС в каустике и переходу его компонентов в раствор едкого натра (при скорости вращения мешалки более 3500 об./мин отмечено .превышение допустимых ГОСТом норм на содержание тяжелых металлов в гото­ вом продукте). Следовательно, число оборотов мешалки имеет ограничения, и изменение этого фактора можно определить интервалом 500—300 об./мин.

3.Количество экстрагента ЖСС существенно влияет на скорость и полноту извлечения ртути из раствора ка­ устика. Согласно проведенным исследованиям для лабо­ раторного реактора область изменения этого фактора следует ограничить интервалом 100—400 г.

4.Увеличение расхода щелочи интенсифицирует про­ цесс, но уменьшает степень очистки. Поэтому областью определения этого фактора будет расход 0,01—0,09 л/ч.

5.Опыты показали, что с ростом концентрации ртути

вэкстрагенте ЖСС существенно увеличивается время экстрагирования ее из раствора. При этом установлено, что эффективная очистка раствора едкого натра обеспе­ чивается лишь в том случае, если содержание ртути в ЖСС не превышает 1—2 мае. %; при большей концен­ трации ртути в ЖСС скорость и полнота извлечения ее из раствора заметно понижается. В лабораторном реак­ торе количество ртути в ЖСС практически не превышает 1 мае. %, поскольку объем потока мал и в сплав из ис­ ходного раствора каустика переходит небольшое коли­ чество ртути, которое не может составить заметной доли от общей массы экстрагента ЖСС. По этим соображе­ ниям в лабораторных условиях можно не учитывать вли­ яния рассматриваемого фактора. Однако в промышлен­ ных непрерывных проточных установках, когда перера­

батываются большие объемы щелочи, фактор накопления ртути нужно обязательно учитывать, а в технологичес­ кой схеме очистки каустика необходимо предусмотреть процесс регенерации экстрагента ЖСС.

6. Повышение концентрации ртути в исходном едком натре от 2 -10-4 до 12-10-4 мае. % мало сказывается на процессе очистки. Поскольку в продукционном едком на­ тре содержание ртути практически постоянно и ниже меньшего предела, то влияние на процесс этого фактора можно не учитывать.

Построение математической модели. Выбор выходно­ го параметра и факторов. В данном случае выходным параметром целесообразно выбрать степень очистки продукционного раствора едкого натра, а независимыми переменными — четыре фактора: температуру раствора едкого натра (Хи °С) интенсивность перемешивания (Х2, об./мин), количество экстрагента ЖСС (Хз, г), расход едкого натра (Х4, л/ч). При этом, ограничивая число фа­ кторов, преследовалась цель обойтись как можно мень­ шим количеством опытов в планируемом эксперименте (при /г= 4 N = 64, а при п '= 6 N '— 64).

Область изменения факторов имеет следующие огра­ ничения:

70 <

Х{ <

120;

100 <

Х3 < 300;

120 <

Х2 <

2600;

0,05 <

ХА< 0,15.

Последующие этапы построения СММ, аналогично предыдущим задачам, представлены в таблицах 5.15— 5.17.

Т а б л и ц а 5.15. Исходные

данные для

планирования

эксперимента

 

 

 

 

 

Факторы

 

Условия эксперимента

Обозначе­

 

 

 

 

 

ние

Xi

X2

*3

* 4

 

 

 

 

Основной уровень

 

Х 0

95

2000

173

0,09

Интервал

варьирования

 

AXi

15

400

25

0,01

Верхний

уровень

 

Xu

ПО

2400

200

0,10

Нижний

уровень

 

Xu

80

1600

150

0,08

Выполненные расчеты позволили записать искомую математическую модель процесса очистки от ртути про­ дукционного едкого натра в виде следующего полинома:

А

У = 6,29 • 10 -5 — 0,35 • 10 -5*! — 0,86 • 1 0 -5.V2 — 0,50 • 10”5х3 +

+ 1,05 • 10“5*4 — 0 , 1 0 - — 0,42 • 10-5лг2*40,14 • 10-5JC3A:4.

Краткие выводы. По полученному уравнению регрес­ сии можно рассчитать содержание ртути в каустике на выходе после очистки при любом сочетании независи­ мых переменных в изученном интервале их изменения и

 

Т а б л и ц а

5.16. План и результаты эксперимента

 

 

 

 

План

 

 

Эксперимент

 

 

Расчет

Номер

 

 

 

 

 

 

 

 

 

 

 

 

опыта

 

XХ

Д-2

 

г/Ц1-1П‘

Уи,-1°‘

уа- т

 

 

«и-10*

 

 

 

 

 

 

 

 

 

 

 

 

1

+ i

+ 1

+ 1

+ 1

+ 1

5,10

5,35

5,225

0,3125

4,9698

2

—1

+ 1

+ 1

+ 1

5,53

5,71

5,620

0,1620

5,8684

3

+ i

+ 1

—1

+ 1

+ 1

7,46

7,03

7,245

0,9245

7,5360

4

+ i

—1

—1

+ 1

4-1

8,60

8,83

8,720

0,2880

8,4346

5

+ i

+ 1

+ 1

—1

4-1

6,39

6,05

6,220

0,5780

6,2458

6

+ i

' —1

+ 1

—1

4-1

7,26

7,05

7,165

0,1805

7,1444

7

+ i

+ 1

—1

—1

4-1

8,73

9,02

8,875

0,4205

3,8120

8

+ i

—1

—1

— 1

4-1 10,05

9,26

9,685

3,1204

9,7106

9

+ i

+ 1

+ 1

+ 1

—1

4,31

3,80

4,055

1,3005

4,1798

10

+ i

—1

+ 1

+ 1

—1

4,56

4,48

4,520

0,0320

4,6748

11

+1

+ 1

—1

+ 1

—1

6,06

5,18

5,120

0,0720

5,0748

12

+ i

—1

—1

+ 1

—1

5,91

5,70

5,805

0,2205

5,5698

13

+ i

+ 1

+ 1

—1

—1

4,96

5,11

5,035

0,1125

4,9122

14

+ i

—1

+ 1

—1

— 1

5,25

5,61

5,565

0,0405

5,4072

15

+ i

+ 1

—1

—1

- 1

5,83

5,70

5,765

0,0845

5,3072

16

+ i

—1

—1

 

- 1

6,00

6,13

6,065

0,0845

6,3022

Т а б л и ц а

5.17.

Результаты

расчета

коэффициентов

и статисти­

 

 

ческий анализ уравнения регрессии

 

 

 

 

Наименование

 

 

 

Расчетные

значения

 

п/п

 

 

 

 

 

1

Коэффициенты регрес­

bQ= 6,29-10-5;

 

ЬХ2 =

0,06-10"5

 

сии

 

 

 

Ьх= —0,35-10-5;

Ь13 =

0,03 -Ю"5

 

 

 

 

 

Ь2 =

—0,86-10-5;

ЬХ4

=

—0,10* 10-5

 

 

 

 

 

Ь3 =

—0,50-10-5;

Ь23 =

—0,07- 10“ь

 

 

 

 

 

ЬА= 1,05-10-5;

 

Ь24 =

—0,42-10~5

2

 

 

 

 

 

 

 

 

&34 =

—0,14-10-5

Построчные дисперсии

приведены в табл. 5.16

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

2 ] Su =

7,93-10-11

 

 

 

3 * Критерий Кохрена (расчетное значение)

4Критерий Кохрена (табличное значение)

5Условие однородности

6Дисперсия воспроизво­ димости (ошибка опы­ та)

г*

_

3,12-10-11 — 0

 

Uп

--

-

 

U,OC7

 

1

II

7,93-10-11

т

II

нQ)

О 4*.

II СЛ о"

=

1 t f2 = N =

16)

II

 

 

 

выполняется;

2

 

 

su — однородны

 

2

 

1

 

 

 

s0

---------7,93- 10-п = 4,96-10-12

 

 

 

16-

 

 

 

п/п Наименование

7Дисперсия коэффици­ ентов регрессии

8 Критерий Стьюдента (расчетное значение)

9 Критерий Стьюдента (табличное значение)

10 Условие значимости

11Дисперсия адекватности

12Критерий Фишера (расчетное значение)

13Критерий Фишера (табличное значение)

14Условие адекватности

 

 

Расчетные значения

 

s2 =

4 ,96 - 1 0 - 12

= 0 ,1 5 5 - 1 0 - 12

1 6 -2

 

bi

 

 

 

 

 

..ли

^ = 0 . 3 9 - 1 0 - =

 

 

tbo =

159,8;

tt,t

= 26 ,66 ;

=

1,47

tb i =

8,85;

tbl2

=

1,44;

tb2A=

10,61

tb2 =

21 ,9 8 ;

h i3

=

0,74;

4 3i =

3 ,4 5

h 3 =

12,76;

tbu =

2 ,5 6

 

 

tT=

2,12

(<7

= 5 % ,

'

- 1 )

=

16)

кроме

U, tbi3 <

выполняется,

h i2 <

< и

и tb23 <

U\

bi2, b13, b23 незна­

чимы (из уравнения регрессии исклю­ чаются)

5ад = 3 ,6 7 - 1 0 - '2;

mд= 2, N = 16;

/ад == N / = 8;

у вычисляется

по уравнению регрессии со значимы­

ми коэффициентами

(/ =

8)

F p =

3 , 6 7 - 1 0 - 12

0,74

 

---------- — =

 

 

4 , 9 6 - 1 0 - 12

 

 

FT=

2 ,59

(g = 5%,

/, =

/а„ =

= N — 1 =

8 , fa = N(m 1) =

=1-6)

выполняется; уравнение регрессии адекватно процессу

получить количественную характеристику влияния раз­ личных факторов. Наличие эффектов взаимодействия между какими-либо параметрами процесса (например, Х { и Х2; Х2 и ХА\ Хг и ХА) свидетельствует об изменении интенсивности влияния одного из них в зависимости о^ значения другого.

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

Полученную адекватную математическую модель в виде уравнения регрессии можно использовать для реч

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

5.2. Планирование эксперимента при оптимизации трехстадийного процесса

извлечения скандия из многокомпонентного сырья

Введение. Данный раздел посвящен результатам ре­ шения задачи оптимизации многостадийного процесса с применением методов планирования эксперимента. Ста­ дии в процессе соединены последовательно. Поэтому по­ становка задач моделирования и оптимизации данного процесса соответствует постановке задачи в разделе 3.2. Применяемые в данной задаче методы отличаются от описанных в теоретической части пособия. Поэтому ни­ же даны пояснения по их использованию.

Описание процесса. Процесс извлечения скандия из бедного многокомпонентного сырья (содержание скан­ дия— 0,02% в расчете на оксид представляет собой по­ следовательное соединение трех стадий (рис. 5.5).

I стадия — фракционное осаждение гидроксидов, в результате которого отделяется основное количество же­ леза. Помимо скандиевого концентрата, получаются еще два товарных продукта — железо-титановый и маг­ ниевый концентраты. Независимые переменные (факто­ ры) : Х\ — значение pH; Х2 — продолжительность выдер­ живания осадка при данном pH, мин. Параметры опти­ мизации: ух— доля отделяемой примеси (железа), пере­ ходящей в концентрат, %; у2— доля скандия, переходя­ щего в концентрат, %; у3 — критерий экономической эф­ фективности стадии, характеризующий прибыль от обо­ гащения. Значение у2 определяется по формуле

Уз = mscCsc + 2 rrnCi — 2 L,

где т.— количество скандия; С, — стоимость целевого продукта (скандия) или товарных продуктов на стадии;

L— затраты-

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

Рис. 5.5. Схема технологического процесса извлечения скандия.

т. е. при получении нескольких товарных продуктов (по­ мимо скандиевого концентрата) их стоимостные функ­ ции складываются.

II стадия — щелочное отделение алюминия, в резуль­ тате которого выделяются скандиевый (целевой) и алю­ миниевый концентраты. Независимые переменные: X\

концентрация исходного раствора NaOH, %;-Х' — про­ должительность перемешивания, мин; Х'3— содержание

железа в исходном продукте, %. Параметрами оптими­ зации выбраны ранее принятые показатели: у'{— доля

отделяемой примеси (алюминия), переходящей в концен­ трат, %; у' — доля скандия, переходящего в концентрат,

%; У'3— критерий экономической эффективности стадии.

III стадия — селёктивная сорбция Иа фосфорсодер­ жащем сорбенте, в результате которой после дополни­

тельной очистки выделяется чистый оксид скандия — це­ левой продукт. Независимые переменные: Х'[ — концен­

трация исходного раствора H2S04, н; X'' — количество сорбента, мг; X" — содержание алюминия в исходном продукте, %. Параметры оптимизации: уг\ — содержа­ ние скандия в конечном продукте, %; у" — доля скандия,

переходящего в концентрат4, %.

Экономический критерий на последней стадии про­ порционален параметру оптимизации у", поэтому не рас­

сматривается.

Исследователи по результатам серии предваритель­ ных опытов убедились в том, что стадии процесса нахо­ дятся в условиях «почти стационарной области». Поэто­ му для каждой стадии процесса извлечения скандия ма­ тематическую модель целесообразно получать в виде ре­ грессионного уравнения второй степени по каждому уj.

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

Пояснение к планам Хартли. Планы Хартли являют­ ся центральными композиционными ортогональными планами второго порядка. Их можно выгодно использо­ вать вместо ЦКОП, когда априори известно, что часть коэффициентов bi и bij в модели отсутствует и имеются жесткие требования к числу опытов. Коэффициенты bi и bij рассчитываются так:

 

 

 

■Уа +У

(А)

 

 

 

~ 2^

 

X i u X j u y u

 

 

 

«=1_____________ — Уд + Уд

(Б)

 

2 п - р

 

 

 

 

 

 

 

S *2

х г

 

 

 

 

u=l iu

i u

 

 

где уа — значения выхода

в соответствующих звездных

точках;

2 п~р — количество

опытов

в центре композици­

онного

плана (обычно — это регулярная

реплика от

ПФЭ);

а — величина плеча звездных точек

(здесь а =

± 1).

Коэффициенты Ь0 и Ьц рассчитываются по формулам ЦКОП (1.162), (1.165).

 

 

Т а б л и ц а 5.18.

План и результаты эксперимента

 

 

 

Факторы

Расчетные столбцы фак­

 

Выходная

Наимено­

 

 

 

торов

 

переменная

вание

 

 

X а

XS

xxxt

 

 

 

 

 

 

**

У х

У я

У я

 

 

 

1

Я

Основной

3,50

45

 

 

 

 

 

 

уровень

 

 

 

 

 

 

Интервал

 

 

 

 

 

 

 

 

варьиро-

0,25

15

 

 

 

 

 

 

вания

 

 

 

 

 

 

 

 

 

План

 

 

 

 

 

 

Номер

X о

 

 

 

 

 

 

 

 

опыта

 

 

 

 

 

 

 

 

 

1

+ 1

0,588

0,809

0,3457

0,6544

0,4756

5,26

88,48

0,313

2

+ 1

0,951

-0,309

0,9044

0,0954

—0,2938

31,78

89,02

0,361

3

4-1

0

-1 .0 0

0

1

0

43.69

87,03

0,287

4

4-1

—0,951

-0,309

0,9044

0 0954

0.2938

61,08

96,25

0,293

5

4-1

-0,588

0,809

0,3457

0,6544

-0,4756

93,94

96,00

0,302

6

4-1

0

0

0

0

0

37,78

88,88

0,291

7

0

0

0

0

0

39,18

88,46

0,297

Пояснение к пентагональным планам. Пентагональные планы принадлежат к классу центральных компози­ ционных ротатабельных планов второго порядка, поэто­ му коэффициенты уравнения регрессии можно рассчи­ тать по алгоритму 1.6.2 (формула (1.103)). Дисперсия воспроизводимости рассчитывается по двум опытам в це­ нтре плана (один опыт входит в план, другой — ставится дополнительно).

5.2.1. Получение математической модели стадии фрак­ ционированного осаждения гидроксидов. Для получения математической модели I стадии процесса используется пентагональный план, обладающий, как уже отмечалось, одним преимуществом перед ЦКРП —: меньшим числом опытов. Модель ищем в виде регрессионного уравнения второго порядка.

Матрица планирования и результаты эксперимента представлены в табл. 5.18. Коэффициенты регрессии будем искать в матричной форме согласно выражению