СПОСОБ УПРАВЛЕНИЯ АГЛОМЕРАЦИОННЫМ ПРОЦЕССОМ
CONTROL METHOD OF AGGLOMERATION PROCESS

Изобретение относится к черной металлургии и может быть использовано при подготовке сырья к доменной плавке, в частности для управления агломерационным процессом.
Известен способ спекания агломерата, при котором последовательно подают на движущиеся спекательные тележки два слоя агломерационной шихты различной толщины и с различным содержанием в каждом из них твердого топлива и влаги, после чего зажигают шихту и создают под спекательными тележками разрежение - см. RU 2164252, С22В 1/16, 20.03.2001 /1/. В процессе спекания скорость движения спекательных тележек устанавливают по зависимости:
жc
вc
нс
нс
нс
вc
вc
V=K*Q*T*W*H*P/T* W*H*B,
жc
вс
нс
вс
нс
вс
нс
-7
где V - скорость движения спекательных тележек, м/мин; Q - расход железорудной составляющей шихты, кг/т агломерационной шихты; Т, Т - расход твердого топлива соответственно в верхнем и нижнем слоях, %; W, W - содержание влаги соответственно в верхнем и нижнем слоях, %; Н, Н - толщина слоев шихты, соответственно верхнего и нижнего, мм; Р - разрежение под спекательными тележками, мм вод. ст.; В - необходимая основность готового агломерата, равная 0,5-0,7, безразмерная; К - эмпирический коэффициент, характеризующий скорость спекания шихты в обоих слоях, равный (8,5-80)*10, м·т/мин·кг (мм. вод. ст.), при этом толщину нижнего слоя шихты устанавливают в пределах 0,2-0,8 от суммарной толщины слоя шихты. В качестве железорудной составляющей шихты используют шлам, смесь железорудных концентратов, отсев агломерата, окалину. Расход каждого компонента устанавливают соответственно в пределах (0,05-0,13), (0,7-0,8), (0,11-0,12) и (0,03-0,065) от расхода железорудной составляющей шихты.
Известен способ управления агломерационным процессом, раскрытый в RU 2114190, С22В 1/20, 27.06.1998 /2/, в котором регулируют скорость перемещения спекательных тележек, определяемую по следующей зависимости:
V=K*T*m*H*g/t*M*h*Q,
где V - скорость движения спекательных тележек, м/мин; Т - содержание железорудного концентрата в шихте, %; t - содержание твердого топлива в шихте, %; М, m - содержание твердого топлива в шихте соответственно в нижнем и верхнем слоях, %; Н, h - толщина нижнего и верхнего слоев шихты, мм; Q, g - содержание влаги в шихте соответственно в нижнем и верхнем слоях, %; К - эмпирический коэффициент, характеризующий скорость спекания шихты в обоих слоях, равный 0,1-0,2 м/мин. При этом отношение массы железорудного концентрата к массе твердого топлива в шихте поддерживают в пределах 12-20.
В качестве наиболее близкого аналога (прототипа) принят способ управления процессом спекания шихты на агломашине, согласно которому осуществляют регулирование скорости движения спекательных тележек с учетом технологических параметров процесса - см. RU 2037540 С1, 19.06.1995 /3/. В данном способе контролируют массу производимого концентрата на конец текущей смены, сравнивают его с заданным и устанавливают толщину слоя шихты на агломашине в зависимости от прогнозируемого отклонения от заданного на смену значения. Способ направлен на повышение прочности агломерата при выполнении планового задания по объему его получения.
Приведенные способы регулирования агломерационного процесса путем расчета скорости движения спекательных тележек учитывают определенный набор технологических параметров, влияющих на процесс агломерации. Однако ни в одном из приведенных способов не учитываются такие важные параметры, как температура и давление газов в коллекторе, температура в последних вакуум-камерах, температура шихты, температура зажигания (по зонам горна) и др. Поэтому для расчета скорости движения спекательных тележек используется эмпирический коэффициент, имеющий широкий диапазон изменений (как, например, в способе /1/), заданы значительно ограниченные диапазоны расходов каждого компонента аглошихты. Кроме того, в приведенных способах не учитывается изменение во времени состояния оборудования, характеристики сырья и т.п.
Важной характеристикой для управления агломерационным процессом является относительный выход агломерата, который определяется отношением массы спека после грохочения (выделения годного агломерата) к массе сырой шихты, поступающей на агломашину.
Задачей предлагаемого изобретения является создание способа управления процессом, учитывающего максимальное количество параметров, влияющих на массу готовой продукции, с учетом адаптации во времени, характеристик сырья и оборудования.
Технический результат, полученный при решении указанной задачи, заключается в обеспечении получения заданного выхода агломерата, при его высоком качестве, за счет регулирования скорости движения спекательных тележек.
Приведенный технический результат достигается за счет того, что в способе управления агломерационным процессом, при котором осуществляют регулирование скорости движения спекательных тележек с учетом технологических параметров процесса, скорость движения спекательных тележек регулируют для обеспечения заданной величины относительного выхода агломерата, прогнозное значение которого определяют методом регрессионного анализа по уравнению связи, включающему в качестве случайных величин выбранные из массива исходных технологических параметров методом корреляционного анализа значимые технологические параметры процесса за последний, предварительно заданный период работы агломашины, величину которого определяют из условия получения нормального закона распределения упомянутых случайных величин, при этом скорость движения спекательных тележек изменяют в случае установления прогнозного значения относительного выхода агломерата менее заданного.
Для формирования массива исходных технологических параметров используют следующие обозначения:
- относительный выход агломерата,
А
Ш
1
2
3
4
5
к
6
к
7
8
L-1
9
L
10
11
1
12
2
13
N
где М - масса спека агломерата после грохочения, т/ч; М - масса сырой шихты, поступающей на агломерацию, т/ч; x(В) - расход собственного возврата, относительные единицы; х(О) - расход доменного возврата, относительные единицы; х(К) - расход коксовой мелочи, относительные единицы; x (V) - скорость движения спекательных тележек, м/мин;х (Р) - давление газа в коллекторе, мм вод. ст.; х (Т) - температура газа в коллекторе, °С; х (W) - влажность шихты, загружаемой на агломашину, %; х (T) - температура в предпоследней вакуум-камере, °С; х (T) - температура в последней вакуум-камере, °С; х (H) - скорость вращения окомкователя, об/мин; х (ЖР) - расход концентрата 1, относительные единицы; x (ЖР) - расход концентрата 2, относительные единицы; х (ЖР)- расход концентрата N, относительные единицы.
Величину относительного выхода готового агломерата определяют с учетом приведенных выше обозначений по уравнению связи:
0
1
1
2
2
3
3
13
13
0
1
2
3
13
у=A+Ах+Ах+Ах+…+Aх, где А - свободный член; A, A, А,…, А - коэффициенты уравнения связи.
Технологический параметр процесса относят к значимому, если коэффициент парной корреляции между ним и относительным выходом агломерата составляет не менее 0,2.
Предлагаемый способ управления агломерационным процессом основан на изменении скорости движения спекательных тележек, определенной в соответствии с расчетом и с учетом всех параметров, которые прямо или косвенно могут влиять на относительный выход агломерата. Кроме того, введена адаптация параметров во времени, поскольку учитывается состояние основного технологического оборудования, время транспортного запаздывания, изменения характеристик сырья.
На фиг.1 приведена блок-схема, показывающая принцип работы системы управления.
1
2
3
11
12
13
ш
4
5
6
7
8
9
10
а
1
2
3
11
12
13
В блок 1 ввода исходных данных вводят информацию о характеристиках аглошихты и технологических параметрах процесса: расход собственного возврата, т/ч (x); доменного возврата, т/ч (х); коксовой мелочи, т/ч (х); железорудного концентрата 1, 2, N, т/ч (х, x, х); массе шихты, т/ч (М); скорости движения спекательных тележек, м/мин (х); давлении газа в коллекторе, мм вод. ст. (х); температуре газа в коллекторе, °С (х); влажности шихты, загружаемой на агломашину, %(х); температуре в предпоследней вакуум-камере, °С (х); температуре в последней вакуум-камере, °С (х); скорости вращения окомкователя, об/мин (х); массе годного агломерата, т/ч (М). Далее в блоке 2 рассчитывается в относительных единицах расходы собственного возврата (x), доменного возврата (x), коксовой мелочи (х), железорудного концентрата 1 (х), 2 (x), N (х). Затем в блоке 3 рассчитывается выход агломерата, в относительных единицах (у), путем деления массы годного агломерата на массу шихты, поступающей на агломашину. Расчет выхода агломерата производится с учетом транспортного запаздывания между временем поступающей на агломашину шихты и временем получения готового агломерата. В блоке 4 формируется матрица технологических параметров. В блоке 5 выполняется расчет и ранжирование коэффициентов линейной парной корреляции по убыванию. Составляется следующая матрица - блок 6, в которую включены значимые параметры. На основании данных матрицы значимых параметров в блоке 7 выполняется расчет уравнения связи. В блоке 8 рассчитывается прогнозное значение выхода агломерата. В блоке 9 происходит сравнение прогнозного и заданного значения выхода агломерата и, если прогнозное значение выхода агломерата больше или равно заданному, то продолжается вся предыдущая обработка информации. В случае если прогнозное значение меньше заданного, то изменяют скорость движения спекательных тележек таким образом, чтобы относительный выход агломерата был не меньше заданного. Новое значение скорости движения спекательных тележек рассчитывается по уравнению связи, в которое подставляется заданное значение относительного выхода агломерата. Управление величиной скорости спекательных тележек выполняется в блоке 12.
Технический эффект при реализации способа заключается в обеспечении такой скорости движения спекательных тележек, при которой относительный выход агломерата будет не менее заданного. Для чего фиксируют все измеряемые параметры аглошихты и технологические режимы работы агломашины (часовые значения), за последние часы работы, с учетом времени транспортного запаздывания.
Значение периода работы, в течение которого фиксируют измеряемые параметры и технологические режимы, принимают таким, чтобы закон распределения случайной величины (технологический параметр) был нормальным. Как показала практика таких расчетов, обычно при шестидесяти значениях каждого параметра закон распределения его величины становится нормальным. Причем информация о работе агломашины ежечасно обновляется, что обеспечивает адаптацию полученных результатов анализа во времени.
Методом корреляционного анализа выбирают параметры и характеристики технологического процесса, влияющие на величину относительного выхода агломерата (значимые параметры). Коэффициент линейной парной корреляции определяется по формуле [1]:
где x - текущее значение аргумента;  - среднее значение аргумента; у - текущее значение функции;  - среднее значение функции. Коэффициенты корреляции рассчитывают между относительным выходом агломерата (у) и каждым параметром
1
2
3
13
(x, х, х,…, х). Параметры (столбцы матрицы), имеющие коэффициенты корреляции менее 0,2, исключают из первичной матрицы, как малозначимые. Установленная величина коэффициента корреляции обусловлена тем, чтобы доверительный уровень вероятности величины относительного выхода был 95%.
Полученная вторичная матрица значимых технологических параметров используется для расчета уравнения связи, имеющего вид:
0
1
1
2
2
3
3
13
13
у=А+А·x+А·х+А·х+…+А·х,
0
где А - свободный член;
1
2
n
A, A,…, А - коэффициенты уравнения связи;
1
2,
n
x, x…, х - аргументы уравнения связи.
Значения коэффициентов уравнения связи находят [2] из решения системы нормальных уравнений по способу наименьших квадратов. Прежде всего, все переменные зависимости выражаем в стандартизированном масштабе с помощью формул перехода
где у - значение функции;
у
σ - среднее квадратическое отклонение функции,
i
где x - значение аргумента (параметра) в натуральном масштабе;
- среднее значение аргумента (параметра);
xI
t - значение параметра в стандартизированном масштабе.
Уравнение связи в стандартизированном масштабе будет иметь вид:
1
x1
2
x2
13
13
t=β·t+β·t+…+β·t,
x1
x2
x13
1
2
13
где t, t,…, t - стандартизированные значения переменных х, х,…,x;
1
2
13
β, β,…, β - стандартизированные коэффициенты множественной регрессии.
0
0
A·b=0.
Используя метод наименьших квадратов, исследуем на минимум функцию:
1
2
13
Запишем систему уравнений для определения параметров зависимости у от х, х,…,х
yx1
1
x1x2
2
x1x3
3
x1x13
13
r=b+rb+rb+…+rb,
yx2
x2x1
1
2
x2x3
3
x2x13
13
r=rb+b+rb+…+rb,
yx3
x3x1
1
x3x2
2
3
x3x13
13
r=rb+rb+b+…+rb,
• • • • • • • • • • • •
yx13
x13x1
1
x13x2
2
x13x3
3
13
r=rb+rb+rb+…+b,
yx1
yx2
x1x2
x12x13
где r, r, r, …, r - соответствующие коэффициенты линейной парной корреляции.
1
2
3
13
Из приведенной системы уравнений находим b, b, b,…, b.
где Δ - основной определитель;
i
Δ - дополнение определителя по i.
Рассчитанный относительный выход агломерата используют в качестве прогнозного значения. В случае если прогнозное значение выхода агломерата меньше заданного, регулируют скорость спекательных тележек.
Анализ работы агломашины показывает, что с течением времени коэффициенты параметров, влияющих на величину относительного выхода агломерата, изменяются, становится другим вид зависимости, кроме того, параметры, влияющие на относительный выход агломерата, могут исчезать из уравнения связи и вместо них могут появляться другие параметры. Таким образом, уравнение связи может видоизменяться.
Расчеты зависимостей относительного выхода агломерата, выполненные для агломашины АКМ-312 ОАО «Северсталь» в различные периоды ее работы, подтвердили это. Результаты приведены в таблицах.
Так, например, в первый период (табл.1) уравнение связи имеет вид
8
9
6
5
4
13
у=-0,203+0,00295х+0,00565х-0,004х-0,00054х-0,0385х-2,67х.
Во второй период (табл.2) уравнение связи имеет вид
8
9
4
6
5
у=0,57+0,015х+0,0014х-0,072х-0,0005х-0,00038х.
В третий период (табл.3)
1
2
3
4
5
y=0,17-0,0036x-0,0034х+0,415х-0,083х-0,0001х.
В приведенных примерах меняется как вид зависимости, так и значимые параметры уравнения связи.