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

книги из ГПНТБ / Автоматизация обогатительных фабрик

..pdf
Скачиваний:
50
Добавлен:
22.10.2023
Размер:
12.41 Mб
Скачать

 

 

 

до тех

пор,

пока значение другого

 

 

 

контролируемого параметра (напри­

 

 

 

мер, содержания

металла

в концен­

 

 

 

трате) не окажется меньше или боль­

 

 

 

ше допустимого.

 

 

 

 

 

 

После этого необходимо продол­

 

 

 

жить движение по градиенту содер­

 

 

 

жания до тех пор, пока

оно не вой­

 

 

 

дет в допустимую

область,

а затем;

 

 

 

снова двигаться по градиенту извле­

 

 

 

чения и т. д.

 

 

 

 

 

 

 

Такой метод движения обеспечи­

 

 

 

вает достижение

максимального из­

 

 

 

влечения при

заданном содержании

Содержание Рг05 S концентрате

металла

в концентрате.

Недостатка­

Рис. 117. Зависимость стоимости

ми его

являются необходимость по­

становки значительного числа экспе­

1 т фосфоритной

муки от содер­

жания

в ней Р20 5

риментов и относительная сложность

И с п о л ь з о в а н и е

движения.

 

 

 

к р и т е ­

э к о н о м и ч е с к о г о

 

р и я

п р и

и с с л е д о в а н и и

м е т о д о м

Б о к с а —

У и л с о н а.

Известно, что

использование в качестве оценки лю­

бого технологического параметра недостаточно для определения оптимальных условий процесса флотации. Необходимо прибегать к оценке процесса по нескольким показателям или отыскивать «компромиссные» условия. В то же время для быстрой и правильной оценки результатов флотации удобно пользоваться единым показателем, одним из объективных критериев оптимизации может быть принята экономическая эффективность [65].

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

= чрск,

где у — выход концентрата; |3 — содержание металла в концентрате; Ск — стоимость 1 т металла в концентрате.

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

Кэ = уСц.

Величина Ср упитывает стоимость 1 т готового концентрата в зависимости от процентного содержания в нем ценного материала. Ее определяют в соответствии с существующими тарифными расцен­ ками на выпускаемую продукцию.

180

Например, на комбинате «Фосфорит» стоимость (в руб.) 1 т фос­ форитной муки в зависимости от содержания (в %) в ней Р 20 5. изменяется следующим образом:

РргО„

19;

20

25

32

Ск

8,0;

9,7

12,5

16,0

Графически эта

зависимость

выражается

ломаной линией

(рис. 117). Для удобства расчетов ломаная кривая заменяется пря­ мой, проходящей через две крайние точки, соответствующие мини­ мальному и максимальному содержанию Р 20 5 в концентрате.

Зависимость Ср =

/ ф р г о „ )

найдем, используя уравнение прямой,,

проходящей через две точки. После подстановки координат X г = 19;.

= 8;

Х 2 — 32;

Y 2 = 16

и соответствующих преобразований:

получим

 

 

 

Ср = 0,62Х — 3,7, где Z = P ps0i.

Для определения оптимальных условий флотации фосфоритной руды был использован метод многофактор­ ного планирования (метод Бокса — Уилсона).

Для оценки оптимальности ис­ пользовали два критерия — извле­ чение Р 20 Бв концентрат и коэффи­ циент экономической эффективности. Расчет коэффициента экономической эффективности производили по фор­ муле К э = уСр = у (0,62 р — 3,7).

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

Врезультате проведенных иссле­ дований получены оптимальные ре­ жимы флотации: по I варианту — опыт№ 15 (табл. 17), по II варианту опыт № 13 (табл. 17).

ВI варианте при сравнительно высоком извлечении (е = 96,45%) получен концентрат с низким со­ держанием, равным 21,12%. При этом коэффициент экономической эф­ фективности составил 293,1.

Рис. 118. Диаграмма зависимости’ экономического показателя от из­ менения выхода концентрата, со­ держания и извлечения цинка:

1—4 — время измельчения

соответ­

ственно 11,4;

11,3;

11,2; 11,1

мин; 1

4 — добавка соды

соответственно 86,8;

79,6; 7,24;

65,2 г/т; 14 — добавка

талового мыла соответственно 766; 800; 838; 874 г/т; 1—4 — добавка солярово­ го масла соответственно 1600; 1700;

1800; 1900 г/т;

1—4 — добавка жидко­

го стекла соответственно

536;

542;

548; 554 г/т;

14 — время агитации

соответственно

5,4; 5,5; 5,7;

5,8

мин.

181

■182

Т а б л и ц а 1 7

*

 

 

 

Время

Расход

Расход

Расход

Расход

Время

 

 

 

 

 

 

 

измель­

талового

солярово­

жидкого

агитации

 

 

 

 

 

 

 

чения,

соды,

мыла,

го масла,

стекла,

с реаген­

 

 

 

 

Обозначение факторов

 

мин

г /т

г /т

г /т

г/т

тами, мин

е1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x t

X,

х

Xt

х .

X,

 

 

 

 

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

 

12

1300

550

1000

500

4,5

 

 

 

 

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

 

4

700

250

500

300

1,5

 

 

 

 

Рандомиза­

Порядковый

1

+

+

+

+

+

+

90.3

301.8

 

 

ция

номер

2

 

+

+

 

+

+

92.0

306.3

 

 

 

 

3

+

 

+

+

 

73.3

235.8

 

 

 

 

4

+

+

+

 

88,8

295.8

 

 

 

 

5

 

+

 

 

79.1

267.0

 

 

 

 

6

 

-1-

+

 

+

93.2

303,5

 

 

 

 

7

+

 

+

 

90,8

295.4

 

 

 

 

8

 

 

 

 

+

79.3

260.0

 

 

 

bt

 

-1,48

-2,48

2,80

4,93

0,25

2,25

 

 

 

 

 

bi * %i

 

5,92

1736

700

2465

75

3,38

85,85

2831,7

 

 

 

Л КВ

 

0,24

70

28

100

3

0,14

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11,70

1230

578

1100

503

4,64

Y

Р

е

Кэ

Н овый основной уровень

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

1

1

 

 

 

 

1

 

I вариант

Порядковый

9

11,52

1160

606

1200

506

4,78

 

мысленные

номер

10

11,28

1090

634

1300

509

4,92

опыты

 

И

11,04

1020

662

1400

512

5,06

_

_

 

 

12

10,80

950

690

1500

515

5,20

28,45

23,06

93,25

301,6

 

 

13

10,56

880

718

1600

518

•5,34

28,45

23,17

33,50

302,0

 

 

14

10,32

810

746

1700

521

5,48

29,85

22,55

95,9

306,9

 

 

15

10,08

740

774

1800

524

5,62

31,20

21,12

96,45

293,1

 

 

16

9,84

670

802

1900

527

5,76

29,35

22,94

94,70

308,7

-183

Ь[

 

 

-20,2

-82,1

113,5

159,5

17,0

77,0

 

 

 

 

-

 

 

 

 

 

 

 

 

 

 

 

 

ь, ■h

 

 

—80,8

-5 7 540

28 400

79 700

5100

115,5

 

 

 

 

 

 

 

 

 

 

 

 

 

У

Р

е

Кэ

^■кв

 

 

0,10

72

36

100

6,0

0,15

 

 

 

 

Новый основной уровень

11,9

1228

586

1100

506

4,65

 

 

 

 

II вариант Порядковый

9

11,8

1156

622

1200

512

4,80

 

 

 

мысленные

номер

10

11,7

1084

658

1300

518

4,95

опыты

 

И

11,6

. 1012

694

1400

524

5,10

_

_

 

 

12

11,5

940

730

1500

530

5,25

_

_

_

 

 

13

11,4

868

766

1600

536

5,40

28,0

25,50

95,2

340,0

 

 

14

11,3

796

802

1700

542

5,55

29,5

22,63

95,7

3050

 

 

15

11,2

724

838

1800

548

5,70

27,9

24,20

94,9

3150

 

 

16

11,1

652

874

1900

554

5,85

30,0

21,70

2920

 

 

 

По II варианту при более богатом

 

 

 

концентрате

(Р =

25,50%) извлечение

 

 

 

равно 95,2% и получен больший коэф­

 

 

 

фициент

экономической эффективности

 

 

 

К э =

340,0.

Оптимальные условия, по­

 

 

 

лученные по II варианту, являются бо­

 

 

 

лее приемлемыми

для условий данной

 

 

 

фабрики.

 

 

 

опы­

 

 

 

Сравнение лучших результатов

1 2

3

И

тов двух

вариантов позволяет сделать

П ар ам е т р ы

 

вывод о том, что поиск оптимальных

Рис. 119. Кривая значимости

условий с использованием коэффициен­

та экономической эффективности обес­

параметров процесса по

Сат-

печивает

такие

условия флотации и

терзвайту

 

 

предприятие

сможет

 

качество концентратов, при которых

получить

наибольший

экономический

эф­

фект.

118 показана

диаграмма

зависимости экономического

На рис.

показателя от изменения некоторых технологических параметров процесса. Диаграмма показывает, что экономическая эффективность является комплексным параметром, К э / (е, |5, Ср, а. . .), объек­ тивно учитывающим изменение качественных показателей флота­ ционного процесса. В определенных случаях она может быть исполь­ зована для оценки результатов процесса.

Метод стохастического равновесия. При исследовании сложных

многофакторных процессов

исследователю приходится иметь дело

с большим числом параметров,

хотя многие из них в дальнейшем

оказываются незначимыми.

В

то же время исследователь не мо­

жет произвести квалифицированного отбора наиболее существенных параметров, так как его выбор основан на интуиции.

В 1956 г. Саттерзвайтом для проведения отсеивающих экспери­ ментов предложен метод стохастического равновесия (случайного баланса) [92], в основе которого лежит предположение, что все па­ раметры исследуемого процесса по их значимости можно расположить в ряд (рис. 119).

Метод «случайного баланса» позволяет сделать выбор наиболее

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

Применение метода «случайного баланса» для выделения значи­ мых параметров рассмотрим на примере исследования процесса фло­ тации цинка.

 

Изучалось влияние на результаты процесса, т. е. на гщ, следу­

ющих параметров:

Х 3 — рас­

X

L— массы (веса) твердого; Х 2— расхода ксантата;

хода медного купороса; Х 4 — pH пульпы; Х 5 — степени

аэрации;

Х

в— времени флотации; X , — времени агитации пульпы с известью;

184

Т а б л и ц а I S

 

Параметр

 

X ,

Х 3

Х а

Лг опыта

 

 

Х а

X 5

X ,

X 7

X ,

X ,

Рандомизация

Блок 1

Блок 2

Блок 1

Блок 2

Блок 3

Блок 3

1

+

+

+

1

8

5

2

+

+

7

3

7

3

+

+

6

6

1

4

+

5

2

8

5

+

+

4

7

3

6

+

+

8

4

2

7

 

2

5

4

8

 

 

 

3

1

6

Х в — времени агитации пульпы с медным купоросом; Х д — расхода вспенивателя 1.

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

Для этого в нашем случае все параметры разбиты на три блока по три параметра в каждом. Для каждого блока составлена матрица планирования типа 23, и по таблице случайных чисел определена рандомизация опытов в каждом блоке (табл. 18).

Матрица планирования для отсеивающих экспериментов полу­ чена объединением из каждого блока строк с одинаковыми номерами рандомизации в одну строку (табл. 19).

 

 

 

 

 

 

 

 

Таблица 19

Лг опыта

X,

X.

Ха

Xt

Параметр

X,

X,

X,

X,

Ха

 

1

+

+

т1

_

_

_

+

__

+

2

+

+

+

3

+

+

+

+

+

4

+

+

+

5

+

+

+

+

+

6

+

+

+

+

7

+

+

+

_ 1

—.

+

+

8

 

+

 

+

+

-L.

 

 

+

 

 

i

 

 

1 Исследования проводились в лабораторных условиях на Мизурской обогати­ тельной фабрике.

185

В7=-95,2 + 91,75=-3,95

В8= -99,8 + 92,8=-2,0

В,= - 99,8 +91,9

= -Z,9

64 = /284= 3,22

Рис. 120. Диаграмма

рассеяния

результатов флотации для девяти

параметров

 

 

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

Диаграммы рассеяния представляют собой ряды распределения результатов опытов соответственно уровням варьирования исследу­ емых параметров.

Таблица 20

ез

6-

3

с

о

%

Количество твердого

Расход ксаитогената

Расход медного купороса

pH пульпы

Аэрация

Время флота­ ции

Время агита­ ции с СаО

Время агита­ ции о CuSOj

Расход вспсииватсля

X.

х.

X,

Xi

X,

X.

X,

■^8

X 9

е

е'

е"

е "' е " "

1

+

+

+

 

 

 

+

 

1

88,8

88,8

89,31

89,81

90,35

 

 

+

+

i

2

+

 

 

90,5

90,5

91,16

91,16

91,16

3

+

____

-1-

+

+

 

89,4

89,4

90,05

90,36

90,36

4

+

+

 

+

91,6

91,6

91,6

91,9

92,0

5

 

 

+

+

 

+

+

+■

96,4

89,95

91,17

91,07

92,11

6

+

____

+

+

+

 

_L

96,98

90,53

91,07

91,09

91,09

7

 

+

+

+

+

 

____

1

96,7

89,46

90,68

90,48

91,02

8

+

 

+

+

+

 

 

_i_

97,5

91,17

91,17

90,97

91,51

 

 

 

 

 

Д86

Для ряда обоих знаков каждого параметра строят медианы. Разность медиан параметра называют его вкладом. Выделение

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

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

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

Как видно на диаграмме, наибольшая разность медиан соответ­ ствует параметру Х 4 (pH пульпы). Для количественной оценки вклада, вносимого этим параметром, составляют таблицу с одним; входом:

 

-ЙГ4

96,40

88,8

96,98

89,4

96,70

91,6

97,50

90,5

Среднее + 96,8

Среднее— 90,35

Вклад Я4= 96,8 — 90,35 = 6,45, коэффициент

Процесс выделения существенных параметров с одновременным: получением приближенной статической модели процесса производят в несколько этапов.

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

Для исключения влияния выделенных параметров на остальные корректируют результаты эксперимента. Корректировку производят

следующим

образом: ко

всем

et-

прибавляют величину Ле; =

= — bj {tj +

1), где

tj — уровень,

на котором находится в данном

опыте выделенный

параметр

(tj

=

± 1); bj — коэффициент; j

номер параметра.

 

=

3,22,

следовательно, для всех опытов,

В нашем случае bj = Ь4

где параметр Х 4 находится на верхнем уровне (£4 = + 1 ), Ае4 = = —6,44, а в опытах, в которых Х 4 находится на нижнем уровне

(h = —1). Ае4 = 0.

После сложения е(. и Ае;- заполняют графу е' в табл. 20.

Для скорректированных результатов эксперимента снова строят диаграмму рассеяния (рис. 121), но для восьми параметров, так как Х 4 исключен.

187

£'

 

9

 

 

 

 

 

 

 

 

 

 

903

90’е

 

 

 

 

90,8•

90,5

.

89,9,| г -

M i

 

90,3

'.•W

.*

 

 

 

•89.7ШJ Г 89,8.

 

1 89,85 *

* 89,8

 

.89,8

 

• •

«

*.

 

. .

 

т

 

 

 

 

 

 

 

в

 

 

 

 

4-

 

 

 

 

- +

 

 

4-

 

 

 

-

-I-

-

+

- +

-

+

 

-

4*

 

Лi

яг

 

 

*s

хs

Ху

х3

Хд

Рис.

121. Диаграмма

рассеяния

для

скорректированных

результатов

флотации

 

 

 

 

 

 

 

 

 

 

 

 

Как видно из рис. 121, после исключения параметра Хл, вклады ■остальных параметров изменили свои значения. Наибольшими теперь являются вклады Х 3 н Х 8. Для их количественного определения построена таблица с двумя входами:

+

GO

- Д Г8

+ * з

 

 

- А 'з

89,95

Среднее

90,5

Среднее

89,46

89,7

89,4

89,9

88,75

Среднее

91,6

Среднее

90,53

89,64

91,17

91,38

I

.

89,7 + 89.64—89.9 — 91,38

=

—0,48

оз —------------------ -----------------

 

4

 

 

,

89,7 + 89,9 — 89.64-91.38

Л „„

Dg ----

------------------------------------- ^ --------------------------------------

U , O 0 .

Корректировку результатов выполняют аналогично корректи­ ровке, описанной выше, но теперь ко всем прибавляют Ае3 и

Ле8:

Де3 = - Ь3 (t3+ 1) = 0,48 ( ± 1 + 1),

Ае8 = — fc8 (f8 + 1) = 0,33 ( ± 1 + 1 ).

Таким образом, во всех опытах, в которых параметры Х 3 (расход медного купороса) и Х 3(время агитации пульпы с медным купоросом) находятся на верхних уровнях, к е(: добавляют Де3 = 0,96 и Де8 = 0,66.

Далее строят на базе вновь скорректированных результатов (графа е", табл. 20) диаграммы рассеяния для шести параметров (без Х 3 и Х 8) и т. д.

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

188

Значимость выделенных

факторов

проверяют

по

критерию

Стыодента. Фактор значим,

если bt

abt.

 

 

Значение оь известно из предыдущих исследований, проведенных

на той же руде, и равно 0,11.

t находим из приложения 1. Для числа

степеней свободы/ = N — 1

= 8 — 1 =

7 при р =

0,95,

t = 2,37.

Таким образом, при Ъ^ 0,11-2,37

0,26 фактор значим.

Если вы­

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

Последовательно, аналогично описанному выше способу, произ­ водят выделение всех существенных параметров.

В нашем примере методом «случайного баланса» из девяти иссле­ дуемых параметров только пять выделены в качестве значимых и включены в программу дальнейших исследований (Х 3, Х 4, Х 7, Ха, Х 0). Значения остальных параметров процесса застабилизированы на среднем уровне.

Результаты проведенного отсеивающего эксперимента предста­ влены в табл. 21.

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

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

К таким адаптационным методам оптимизации относится симп­ лексное планирование эксперимента.

Впервые этот метод был предложен Спиндлеем, Хецтом и Хнм-

суорзом в 1962 г. [93].

Исходная матрица планирования представляет собой правиль­

ный симплекс, т. в. совокупность п +

1 равноудаленных друг от

друга

точек

в

?г-мерном пространстве

(тг-количество параметров).

Например,

в одномерном

пространстве (п = 1)

симплексом

является отрезок прямой, в

двухмерном — равносторонний тре­

угольник и т.

д.

 

 

 

 

 

 

 

 

 

Таблица 21

Параметры в порядке убывания пх значимости

Обозначения

Ъ1

p H пул ьп ы

 

 

 

*4

3,22

Р а сх од

м едного

к уп ор оса

 

* 3

— 0,48

В рем я

агитации

пул ьп ы с медным к у п ор осом

■^8

— 0,33

В рем я

агитации пул ьп ы с известью

 

Xj

— 0,27

Р а сх од

вспеннвателя

 

* 9

— 0,24

189

x tZ X 2Z

Если

центр

симплекса

находится в на­

 

чале координат,

исходная

матрица имеет

 

следующий вид:

 

 

 

 

 

 

 

АУГ, .. •А п - А

 

 

 

RJU . ••A n-1A„

 

 

 

0 -

r 2 . ••Kn_JLn

 

 

 

0

 

0 . ■■Кп_гКп

 

 

 

0

 

0 . .. -Rn-JCn

 

 

 

0

 

0 .

• о- R n,

хйх2ч

 

 

 

 

 

 

 

Рис. 122. Иллюстрация

Rl

К

2(i +

l)

К‘

2 t ( i + l ) ’

зеркального отражения

 

i

= 1 ,

2, .

. п.

(6)

симплекса

 

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

После реализации исходной матрицы выявляется опыт, давший наихуднше результаты.

Для дальнейшего движения в область оптимума следующий опыт ставится в условиях, являющихся зеркальным отражением наихудшего опыта относительно противоположной грани симплекса.

Например (рис. 122), если из трех опытов в точках с координа­

тами Хг

X 2j худшей

оказалась точка Х 12Х 22, то ее зеркальным

отражением будет точка Х 14, Z 24-

Координаты этой точки опреде­

ляются

по формулам

 

 

 

 

Х 14 =

2Х1е- Х 12;

X 2i = 2X9c- X „

где Х и — координаты центра противоположной грани;

 

X 1с -

■У11 + -У1з

X,

х 21 + х 2

 

 

 

Х 12, Х 22 — координаты наихудшей точки (опыт № 2).

После реализации опыта в четвертой точке сопоставляются результаты 1, 3 и 4 опытов. Среди них снова определяется иаихудпгай, находится его зеркальное отражение и т. д.

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

симплекса Ур. Малое значение разности У — У 0 по сравнению

190

Соседние файлы в папке книги из ГПНТБ