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

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

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

Рис. 5.4. Блок-схема программы поиска оптимального режима методом симплексов.

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

Ионообмен как способ выделения веществ из раство­ ров, содержащихся в них, позволяет утилизировать цен­ ные продукты или очищать растворы от нежелательных примесей. Этот процесс осуществляется с помощью иони­ тов (катионитов или анионитов). Применительно к про­ цессу очистки анолита от ртути чаще всего в качестве ионитов (анионитов) используются сильноосновные ио-

шэобменные смолы потому, что в присутствии больших количеств хлорида натрия ртуть находится в виде хлоридных анионных комплексов.

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

Построение математической модели. Выбор выходно­ го параметра и факторов. Процесс ионообменного погло­ щения и концентрирования ртути из анолита на анионите АВ-17 определяется степенью очистки раствора от рту­ ти; эту количественную характеристику целесообразно принять в качестве выходной переменной у.

Кинетические исследования рассматриваемого про­ цесса показали, что на выходную переменную у сущест­ венно влияют следующие факторы: время контакта (Х{, ч), количество смолы, загружаемой в аппарат (Х2, г), расход анолита (Х3, л/ч).

Область изменения указанных факторов оценива­ лась, исходя из технологических условий проведения про­ цесса и исследований его кинетики. При этом приняты такие ограничения: 4 ,4 ^A j^ 8 ; 5 ^ X 2 ^ 15; 0 ,8 ^А з^2 .

Выбор центра плана и интервалов варьирования.

Данные предварительных опытов показали, что центр плана (нулевые уровни факторов) целесообразно помес­ тить в точку с координатами Х10 = 5 ч, Х2о=\1 г, Х30 — = 1,25 л/ч и принять следующие интервалы варьирова­ ния: A Ji=0,5 ч, ААГ2= 1,5 г, АХ3 = 0,25 л/ч. Эти условия, а также значения верхнего и нижнего факторов и коди­ рованные значения х\ сведены в табл. 5.6.

План эксперимента (матрица планирования), соот­ ветствующий ПФЭ типа 23, представлен в табл. 5.7, в ко­ торую также занесены экспериментальные значения вы­ ходного параметра в параллельных опытах ( т = 2): уиь

А

Уи2 и уи, построчные дисперсии s2u и значения уи, пред­

ставленные по уравнению регрессии со значимыми коэф­ фициентами.

Реализация матрицы планирования осуществлялась на опытной установке, на которой проводились предва­ рительные исследования; точность измерения уи и фак­ торов соответствовала предъявленным требованиям.

Т а б л и ц а 5.6.

Исходные данные

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

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

 

 

 

 

 

 

 

Факторы

 

Условия опыта

 

Обоппачс-

 

 

 

 

 

 

нис

*1.

Х2, г

 

Х3, л/ч

 

 

 

 

 

 

Нулевой уровень

 

 

 

 

5

11

 

 

1,25

Интервал

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

+1

АХ,-

 

0.5

1,5

 

0,25

Верхний

уровень

*,•„ =

X,-,.

 

5,5

12,5

 

1,5

Нижний

уровень

хщ =

—1

х,„

 

4,5

9,5

 

1,0

 

Т а б л и ц а 5.7. План и результаты эксперимента

 

 

 

План

 

Выходная переменная

 

Расчет

Номер

 

 

 

 

 

 

 

 

 

 

опыта

 

 

 

 

Уих

«и,

Уа

 

4

Уа

 

 

 

 

 

 

1

+ 1

+ i

+ i

+ i

93,18

92,62

92,9

'

0,16

80,94

2

+ 1

— 1

+ i

+ i

94,43

96,4

95,35

1,69

94,66

3

+ i

+ i

— 1

+ 1 .

80,99

77,61

79,3

 

5,72

80,36

4

+ i

— 1

— 1

+ i

84,20

82,8

83,5

 

1,00

84,10

5

+ 1

+ 1

+ 1

— 1

96,65

94,95

95,8

 

1,44

96,10

6

+ i

—1

+ i

— 1

99,07

97,93

98,5

 

0,64

99,84

7

+ i

+ i

— 1

— 1

89,25

92,75

86,0

 

21,16

85,60

8

+ i

— 1

— 1

— 1

92,91

90,09

91,5

 

4,00

89,30

Расчет коэффициентов уравнения регрессии и его статистический анализ проводились по алгоритму 1.5.1 и программе вычислений на ЦВМ «Мир-1». Результаты расчета приведены в табл. 5.8. В последних строчках этой таблицы записаны адекватные уравнения регрессии с переменными в кодированной форме и в натуральном масштабе.

Краткие выводы. Математическая модель — линей­ ное уравнение регрессии — пригодна для оценки влия­ ния каждого фактора на у (по величине и знаку коэффи­ циента), а также для определения оптимального режи­ ма процесса поглощения и концентрирования ртути из анолита на анионите АВ-17 (ионообменный способ).

5.1.5. Задача Б4. Извлечение ртути из отработанной серной кислоты. Информация о процессе. Согласно тех­ нологии производства хлора и едкого натра, получаемый

Т а б л и ц а 5.8. Результаты расчета коэффициентов и статистический анализ уравнения регрессии

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

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

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

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

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

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

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

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

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

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

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

10

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

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

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

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

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

Ь0

=

90,96;

Ь2 =

5,28;

Ь,2

=

0,57

* i,=

1,86;

63 =

2,59;

Ь13 =

0,19

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

&2з =

1,08

 

 

 

N

 

 

35,74

 

 

 

 

 

 

jH S* =

 

 

 

 

 

 

1

ю

05

 

 

 

 

 

 

 

Gр =

=

0,59

 

 

 

 

— —

 

 

 

 

 

 

35,74

 

 

 

 

 

 

 

GT=

0,68 (q =

5%,

fi = т — 1,

f2 =

N =

8)

12

 

 

 

 

 

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

 

 

 

 

 

su

однородны

 

So

=

-35,74 =

4,47

 

 

 

 

 

8

 

 

 

 

 

 

 

 

52 =

4,47

0,28

или

Sb =

0,53

---------=

bi

8,2

 

 

 

00

*

 

 

 

II

О

•о

II

0

 

 

и

5

 

 

 

 

 

л

 

 

 

 

 

К

=

3,5;

 

tbit -

1,08

 

 

tb2-

9,96

tbi3 =

0,36

 

 

tT=

 

 

К г

=

2-04

 

 

2;306

(q =

5%,

f =

N{m

-1) = 8)

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

кроме

tbi2 <

U,

h l3 <

< tr,

Кэ < К

612,

*13, h i

незна­

чимы

(из

уравнения регрессии ис­

ключаются)

 

 

 

 

 

2

6,48;

m — 2;

JV =

8;

/ ад =

5ад =

 

 

А

 

 

 

 

 

 

 

/V

вычисляется

по

урав­

= W— / = 4; у

нению

регрессии

со

значимыми ко­

эффициентами (/

=

4)

 

 

Fp

6,48

1,45

 

 

 

4,47

 

 

 

 

 

 

 

 

 

 

FT=

3,84 (q =

5%, f x = /ад =

 

= N — l =

45,

 

=

7V(m 1) = 8 )

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

уравнение

регрессии

адекватно процессу

 

 

 

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

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

п/п

15 Адекватное

уравнение

А

регрессии

(факторы в у = 90,36 — 1,86xi + 5,28х2 — 2,59хз

кодированной форме)

16Статистическая матема- А тическая модель (фак­ у = 53,19 — 3,72*! + 3,52Х2 — торы в натуральных — 10,36Хз единицах)

хлор-газ содержит влагу и ртуть в виде паров сулемы. Этот продукт перед поступлением потребителям подвер­ гается осушке орошением серной кислотой, которая раз­ бавляется и загрязняется ртутью (в кислоту переходит не только влага, но и сулема). Поэтому, отработанную разбавленную серную кислоту нужно обязательно очи­ щать от ртути, особенно в случае использования ее для получения удобрения (сульфата аммония). Многие опи­ санные в литературе способы очистки серной кислоты от ртути практически не применяются из-за существенных недостатков (низкая степень очистки — не более 90%, многостадийность, сложность аппаратурного оформле­ ния).

Экспериментальные исследования предложенного в 1973 г. сорбционного способа извлечения ртути из отра­ ботанной серной кислоты оксигидридом кремния пока­ зали, что разбавленная (15—20%-ная) серная кислота полностью очищается от ртути; с увеличением ее концен­ трации до 70% степень очистки уменьшается, а 70%-ная серная кислота практически не очищается от ртути этим способом. Из анализа результатов указанных исследова­ ний вытекает также, что на степень очистки разбавлен­ ной серной кислоты существенно влияет расход кислоты на входе в адсорбер (при изменении расхода 20%-ной кислоты от 5 до 12,5 мл/мин степень очистки составляет соответственно 98—92%) и количество адсорбента — оксигидрида кремния в аппарате (при увеличении со­ держания адсорбента от 3,5 до 13 г степень очистки воз­ растает от 88 до 99%). Изменение начальной концентра­ ции ртути от 2 до 16 мг/л приводит к незначительному колебанию степени очистки.

Построение математической модели. Выбор выходно­ го параметра и факторов. Анализ априорной информа­ ции и результаты кинетических исследований показыва­ ют, что в данном случае логично выбрать в качестве вы­ ходной переменной у — степень очистки серной кислоты от ртути, и факторами, существенно влияющими на про­ цесс очистки, расход отработанной кислоты в проточном адсорбере (Хи мл/мин), количество оксигидрида крем­ ния в аппарате (Х2, г) и концентрацию серной кислоты

№ , %).

Область изменения этих факторов можно определить на основе имеющейся информации о процессе; для реше­ ния рассматриваемой задачи выбраны следующие пре­ делы: 4 ^ X i^ l6 ; 14; 15^Х 3^ 4 0 .

Последующие этапы решения этой задачи выполня­ лись аналогично задачам Б1 и БЗ. Поэтому эти этапы

целесообразно

представить

соответствующими

табли­

цами.

 

 

 

 

 

 

 

 

Условия эксперимента приведены в табл. 5.9.

 

Т а б л и ц а 5.9.

Исходные

данные

для

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

 

 

 

 

 

 

 

Факторы

 

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

Обозна­

 

 

 

чение

 

л-2

*3

 

 

 

 

 

 

 

Нулевой уровень

 

 

*-о

10

8,25

25

Интервал варьирования

+1

АХ,

5

4,75

10

Верхний

уровень

хт =

Xt

15

13

35

Нижний

уровень

хш =

—1

x t.

5

3,5

15

Матрица планирования ПФЭ 23 приведена в табл. 5.10.

Результаты обработки эксперимента приведены в табл. 5.11.

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

Полученная статистическая математическая модель может быть использована как основа для решения зада-

 

Т а б л и ц а

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

 

 

 

 

План

 

Эксперимент

 

Расчет

Номер

 

 

 

 

 

 

 

 

 

опыта

 

Х

 

 

Уи,

Уиг

Уи

4

Уи

 

 

 

 

1

+ i

+ 1

+ i

+ i

77,18

78,92

78,05

1,514

77,883

2

+ i

— 1

+ i

+ i

82,23

82,71

82,47

0,115

83,013

3

+ i

+ i

—1

+ i

76,01

76,45

76,23

0,097

75,883

4

+ i

— 1

—1

+ i

80,94

81,14

80,94

0,08

80,963

5

+ i

+ i

+ i

— 1

92,33

92,77

92,55

0,097

93,073

6

+ i

- 1

+ i

— 1

98,88

99,32

99,10

0,097

98,203

7

+ i

+ i

— 1

— 1

90,69

91,27

90,98

0,168

91,023

8

+ i

— 1

- 1

— 1

95,54

96,10

95,82

0,157

96,153

чи оптимизации с целью определения параметров опти­ мального режима проведения процесса.

5. 1. 6. Задача Б5. Очистка сточных вод от ртути. Ин­ формация о процессе. В производстве хлора и едкого на­ тра ртутным электролизом весьма существенным кана­ лом потерь ртути являются потоки загрязненных рту­ тью сточных вод, которые образуются при промывке электролизеров и аппаратуры, при ремонтах и других технологических операциях. Сточные воды пополняются также за счет конденсата холодильников и стоков гидро­ затворов водорода и хлора. Спуск таких вод в канали­ зацию категорически запрещен.

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

С целью сбора информации о процессе проведены спе­ циальные кинетические исследования основные резуль­ таты которых заключаются в следующем:

а) высокая степень очистки (99,8%) достигнута при скорости подачи сточной воды в адсорбер 75 мл/мин, со­ держании ртути в ней на входе в аппарат 1 мг/л и коли-

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

15Адекватное уравнение регрессии (факторы в кодированной форме)

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

 

 

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

 

Ь0 =

87,018;

 

й2 =

1,025

 

й, =

—2,565;

 

й3 -

—7,595

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

 

 

N 2

 

 

 

 

 

 

2 ]

5U = 2,325

 

 

 

 

 

и—1

 

 

 

 

 

 

 

Gр

=

1,514

 

0,65

 

 

 

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

 

 

 

 

 

2.325

 

 

 

 

 

От =

0,68

(? =

5%,

/г =

N — 1 =

=

1,

h =

N = 8)

 

 

 

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

2

 

 

 

 

su— однородны

2

 

1

 

 

0,29

 

 

So =

----- 2,325 =

 

 

 

 

8

 

 

 

 

 

 

s2

=

0,29

 

0,018

или Sb = 0,135

---------=

bi

 

8-2

 

 

h 2 =

'

tb0 =

644,58;

 

7,59

tbi =

19,0;

 

 

h 3 =

59,22

tT=

2,306

(? =

5%,

f =

N(tn —

- 1 )

= 8 )

 

 

коэффициенты зна­

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

чимы

 

 

 

 

 

 

Saд == 0,836; tn — 2,

N -— 8;

е

 

W — / =

4;

А

вычисляется

/ад =

 

по уравнению регрессии со

значимы­

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

(/ =

4)

Fp =

0,836

 

2,88

 

 

 

— ------- =

 

 

 

 

 

0,29

 

 

 

 

 

FT=

3,84;

(q =

5%,

=

/ад =

= N — l =

4, f2 = W(m — 1) = 8)

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

уравнение

регрессии

адекватно процессу

 

 

 

Л

87,018 — 2,565*, + 1,025*2 —

у =

— 7,595*3

 

 

 

 

 

 

Л

188,216 — 0,513Xi +

0,216*2 —

у =

— 1,013*з

честве оксигидрида кремния 7 г; однако увеличение ис­ ходного содержания ртути до 5 мг/л привело к снижению степени очистки до 82% (при той же скорости и количе­ стве сорбента);

б) сточная вода очищается почти полностью при рас­ ходе воды 50 мл/мин, но при увеличении скорости пода­ чи степень очистки уменьшается и при 250 мл/мин со­ ставляет всего 80%;

в) степень очистки сточных вод растет от 82 до 98% при увеличении содержания оксигидрида кремния в ре­ акторе от 3 до 7 г; однако дальнейший рост количества сорбента не влияет на степень очистки.

Построение математической модели. Выбор исходно­ го параметра и факторов. Результаты предварительных исследований послужили основанием для выбора выход­ ным параметром у — степени очистки сточной воды от ртути, и факторами, существенно влияющими на процесс, исходную концентрацию ртути в сточной воде (Хь мг/л), ее расход (Х2, мл/ч) и количество адсорбента — оксиги­ дрида кремния (Яз, г).

Область изменения выбранных факторов также опре­ деляли исходя из результатов предварительных опытов:

0,5^ Z ,^ 5; ЗО^Хг^ЗОО; 2 < Х 3< 9 .

Т а б л и ц а

5.12. Исходные данные для

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

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

 

 

 

 

Факторы

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

Обозна­

 

 

 

чение

Xi

 

 

 

 

 

 

 

Нулевой уровень

* 0

2

150

5

Интервал

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

Xin

1,2

100

2

Верхний уровень

3,2

250

7

Нижний уровень

Xin

0,8

50

3

Последующие этапы построения математической мо­ дели, как и в задаче Б4, представлены в таблицах 5.12— 5.14.

Краткие выводы. Уравнение регрессии позволяет по значению и знаку коэффициентов регрессии (60> Ь\, &2> 6з) оценить интенсивность и направление влияния каж­ дой переменной на степень очистки сточных вод от рту­ ти оксигидридом кремния. Так, например, эффективно влияет на процесс очистки увеличение количества окси-

 

Т а б л и ц а

5.13.

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

 

Номер

 

 

План

 

Эксперимент

Расчет

*0

 

 

 

 

 

 

2

 

опыта

 

 

 

Уи2

Уи

 

 

 

 

 

 

su

Уи

1

+1

+ i

+ 1'

89.76

91,68

90,72

1,837

91.358

2

■f 1

- 1

+ 1

+ 1

95.95

97,31

96,63

0,925

95,104

3

+1

+ i

—1

+ 1

94,11

95,51

94,81

0,98

95,728

4

+ 1

- 1

—1

+1

99.98

99,02

99,5

0,461

99,474

5

+1

+ i

+1

- 1

82,94

81,56

82,25

0,952

81,763

6

+1

- 1

+ 1

—1

83,35'

84,91

84,13

1,220

85,509

7

+1

+ i

— 1

—1

86,43

87,97

87,2

1,186

86,133

8

+1

—1

—1

- 1

88,79

90,61

89,7

1,654

89,879

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

5.1.7. Задача Б6. Извлечение ртути из продукционно­ го едкого натра. Информация о процессе. В производ­ стве хлора и едкого натра ртутным электролизом унос ртути с потоком продукционного едкого натра неболь­ шой по сравнению с другими технологическими средами (составляет около 0,5% #в условном балансе ртути). Од­ нако очистка от ртути этого потока имеет значение, глав­ ным образом, с точки зрения обеспечения высокого каче­ ства едкого натра (особенно при использовании его в производстве искусственного волокна) и сохранения нор­ мального санитарного состояния производственных по­ мещений заводов-изготовителей и потребителей едкого натра. Указанные обстоятельства вынуждают приме­ нять наиболее эффективные способы извлечения ртути

из растворов едкого натра.

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