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

Изобретение относится к способу автоматического управления процессом плавки медно-никелевого сульфидного сырья в печи Ванюкова при переработке сульфидных шихт на штейн. Способ включает постоянный контроль параметров процесса, корректировку управляющих параметров для стабилизации содержания меди, причем в качестве основного параметра Q_i выбирают удельный расход кислорода на тонну только металлосодержащих в м³/т, разбивку общего диапазона Q_i на ряд областей, соответствующих: нормальному течению процесса Q₀ в интервале 150-250 м³/т, склонности к холодному течению процесса Q_-1 в интервале 150 м³/т и ниже, склонности к горячему течению процесса Q₊₁ в интервале 250 м³/т и выше, дополнительно учитывают степень согласованности процесса, выход значений по основному или какому-либо из параметров процесса за границы допустимых режимов области согласованности интерпретируют как конфликт K_i и по корректировочным моделям ведут поиск вариантов выхода процесса в область согласованности, а в качестве конфликтов выбраны: K₁ - переокисление металлосодержащих (МС), K₂ - недоокисление МС, K₃ - избыток флюсующих, K₄ - недостаток флюсующих, K₅ - горячий ход печи, K₆ - холодный ход печи, K₇ - невозможность прямого определения температуры расплава в реакционной зоне печи, K₈ - невозможность прямого определения физического объема поступающей шихты, при этом проверяют уровень степени согласованности величин и качества загрузки, дутьевых режимов с заданным содержанием меди в штейне по полиномиальным моделям. Обеспечиваются оперативное управление процессом плавки, визуализация процесса, стабилизация качества продуктов плавки, получение неизмеряемых или плохо измеряемых параметров технологического процесса и оценивание состояния агрегата косвенными методами, снижение энергоемкости процесса переработки шихты,- стабилизация температурного режима процесса при сохранении плановых заданий и целей и создание программного продукта. 3 ил., 2 пр.

 

Заявляемое техническое решение относится к области металлургии, а именно к процессу плавки медно-никелевого сульфидного сырья в печи Ванюкова.

Известно техническое решение по патенту №2048534, опубл. 31.03.1992, RU. Способ автоматического управления конвертерной плавкой включает контроль состояния шлака во время продувки по косвенному параметру S(t), сравнение параметра S(t) с заданными по технологии пороговыми значениями и изменение положения фурмы Н, расхода I кислорода через фурму, подачу сыпучих в конвертер в зависимости от изменения параметра S(t). В способе автоматического управления конвертерной плавкой дополнительно определяют диапазон изменения параметра S(t) во время продувки, разбивают его на ряд областей S_i, соответствующих нормальному S₀ состоянию шлака, S₁ склонности шлака к сворачиванию, S₂ - предполагаемому сворачиванию шлака, S₃ сворачиванию шлака, S_-1 склонности шлака к выбросам, S_-2 предполагаемым выбросам шлака, S_-3 выбросам шлака, во время продувки определяют область нахождения параметра S(t) и продолжительность t пребывания его в этой области и управление ведут в зависимости от того, в какой области находится основной параметр.

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

Известно техническое решение по патенту №2204616, опубл. 10.11.2000, RU. Способ автоматического управления процессом обжига никелевого концентрата в печи кипящего слоя включает разбивку диапазона основного параметра на ряд областей, установление принадлежности основного параметра к одной из областей и в зависимости от того, в какой области находится основной параметр, изменение загрузки сыпучих. В качестве основного параметра берут температуру в реакционной зоне печи, а разбивку ее диапазона производят на области: T_-1 (800-900)°C, T₀ (850-1150)°C, T₊₁ (1130-1200)°C. Дополнительно устанавливают градиент температуры и его направление и в зависимости от сочетаний к диапазону температуры и градиента температуры выбирают один из трех режимов управления печью кипящего слоя: при T (800-850)°C и gradT любой, T (850-900)°C и gradT равен 0, процесс относят к области T_-1 и процесс ведут с зажженными горелками до перехода его в область Т₀; при Т (850-900)°C и gradT>0, T (900-1130)°C и gradT любой T, (1130-1150)°C и gradT 0 процесс относят к области Т₀ и ведут в автоматическом режиме с расчетом управляющего воздействия величины загрузки сыпучих по аналитическим полиномам; при Т (1130-1150)°C и gradT>0, Т (1150-1200)°C и gradT любой процесс относят к области T₊₁ и при времени нахождения в такой ситуации 6-8 мин прекращают автоматическое управление загрузкой сыпучих, информацию о текущих параметрах обжига отображают на мониторе оператора, а процесс ведут с повышением скорости подачи пыли до перехода его в область T₀.

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

Известно техническое решение по авт.св. №1625015, опубл. 10.06.1999, С22В 5/02, SU, Способ управления автогенной плавкой медного сульфидного сырья, включающий подачу реагентов в плавильный агрегат с расплавлением шихты и образованием штейно-шлаковой эмульсии с последующим расслаиванием ее на шлак и штейн, непрерывный раздельный выпуск продуктов из печи, отбор проб исходных материалов и продуктов, их анализ, постоянный контроль параметров процесса по показаниям приборов и по данным анализа, внесение корректировок в управляющие параметры, при этом проводят отбор не менее двух проб штейно-шлаковой эмульсии, по результатам их анализа строят прямую зависимости содержания меди в эмульсии от содержания серы в эмульсии, описываемую уравнением вида: Cu_э=a+b·S_э, и определяют содержание меди в донном штейне по уравнению Cu_шт=a+b·S_шт, где Cu_э, S_э - содержание меди и серы в эмульсии; S_шт - содержание серы в штейне, определяемое по эмпирическим данным; а, b - коэффициенты уравнения прямой, и вносят изменения в управляющие параметры процесса с учетом полученного значения Cu_шт.

Однако способ не обеспечивает оперативное управление процессом плавки в зависимости от состава подаваемой шихты и достижения требуемого качества готовых продуктов. При этом для получения значимых коэффициентов регрессии a и b по законам статистики требуется значительное количество точек, что связано с существенным интервалом времени течения процесса, а также неясен смысл коэффициентов, поскольку уравнений два и для обоих одинаковые значения а, b.

Известно техническое решение по патенту №2456353 от 09.08.2010, МПК С22В 1/10, RU. Способ автоматического управления содержанием меди в штейне, включающий постоянный контроль параметров процесса, корректировку управляющих параметров для стабилизации содержания меди в штейне, в качестве основного параметра выбирают общий расход шихтовых материалов, производят разбивку его диапазона на области: G_-1 (менее 60 т/ч), G₀ (от 60 до 180 т/ч), G₊₁ (более 180 т/ч), и в зависимости от того, в какой области находится основной параметр, осуществляют изменение общего расхода шихтовых материалов и технического кислорода до достижения области G₀, дополнительно определяют соотношение расхода технического кислорода дутья на тонну шихтовых материалов и при градиенте его изменения свыше 10% от регламентируемого корректируют общий расход шихтовых материалов и технического кислорода в зависимости от того, в какой области находится основной параметр, до достижения области G₀.

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

Недостатком способа является то, что достижение указанной в патенте области G₀ (от 60 до 180 т/ч) по сути не является гарантом управления содержанием меди в штейне, поскольку без количественных определений в режиме реального времени соотношения расхода технического кислорода дутья на тонну шихтовых как показателя физико-химического состояния процесса в целом утверждения носят декларативный характер. Кроме того, практика ведения процесса Ванюкова показывает необходимость и возможность автоматического управления в области G_-1 (менее 60 т/ч) как основного режимного параметра в силу вынужденных обстоятельств, например отсутствия достаточного количества шихты.

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

Технический результат достигается тем, что способ автоматического управления качеством конечных продуктов процесса Ванюкова при переработке сульфидных шихт на штейн включает получение, анализ и обработку данных АСУ ПВ, контроль состояния процесса по косвенному параметру Q_i, разбивку диапазона основного параметра на ряд областей и установление принадлежности основного параметра к одной из областей, сравнение с заданными по технологии пороговыми значениями расхода кислорода, изменения состава и количества шихтовых материалов в зависимости от изменения параметра Q_i, причем в качестве основного параметра Q_i выбирают удельный расход кислорода на тонну металлосодержащих в м³/т, разбивку общего диапазона Q_i на ряд областей, соответствующих: нормальному течению процесса Q₀ в интервале 150-250 м³/т, склонности к холодному течению процесса Q_-1 в интервале 150 м³/т и ниже, склонности к горячему течению процесса Q₊₁ в интервале 250 м³/т и выше, дополнительно учитывают степень согласованности процесса, выход значений по основному или какому-либо из параметров процесса за границы допустимых режимов области согласованности интерпретируют как конфликт K₁ и по корректировочным моделям ведут поиск вариантов выхода процесса в область согласованности,

а в качестве конфликтов выбраны:

К₁ - переокисление металлосодержащих (МС);

К₂ - недоокисление МС;

К₃ - избыток флюсующих;

К₄ - недостаток флюсующих;

К₅ - горячий ход печи;

К₆ - холодный ход печи;

К₇ - невозможность прямого определения температуры расплава в реакционной зоне печи;

К₈ - невозможность прямого определения физического объема поступающей шихты,

при этом проверяют уровень степени согласованности величин и качества загрузки, дутьевых режимов с заданным содержанием меди в штейне по полиномиальной модели (1):

где в кодированном виде представлены:

x₁=(X₁-5)/5, где X₁ - отклонение загрузки флюсов от заданной нормы, т/ч;

x₂=(X₂-6)/6, где X₂ - нормированное отклонение удельного расхода кислорода на тонну МС, б/р;

x₃=(X₃-5)/5, где X₃ - отклонение расчетного содержания меди в штейне от нормируемого, %;

x₄=(X₁-0,6)/0,2, где X₄ - расчетное содержание меди в шлаке, %;

x₅=(X₅-1340)/40, где X₅ - температура расплава, °C;

Y₁ - степени согласованности, б/р,

корректировку параметров процесса с визуализацией на АРМ оператора показаний основных параметров в виде графиков и/или на индикаторах по моделям (2) - (5):

- загрузка металлосодержащих (мс)

где в кодированном виде представлены:

x₆=(X₆-15000)/10000, где X₆ - расход кислорода, м³/ч;

x₇=(X₇-2400)/1400, где X₇ - расход природного газа, м³/ч;

x₈=(X₈-0,05)/0,05, где X₈ - доля лежалых в МС, б/р;

x₉=(X₉-17)/8, где X₉ - перепад температуры воды на сливе 1-2 ряда кессонов, °C;

Y₂ - уставка скорости загрузки металлосодержащих, т/ч;

- загрузка флюсов

где в кодированном виде представлены:

x₁₀=(X₁₀-60)/40, где X₁₀ - загрузка МС, т/ч;

x₁₁=(X₁₁-0,125)/0,125, где X₁₁ - доля техногенных от МС, б/р;

x₁₂=(X₁₂-1)/1, где X₁₂ - отношение обороты/лежалые;

x₁₃=(X₁₃-200)/50, где X₁₃ - удельный расход кислорода, м³/т;

Y₃ - уставка скорости загрузки флюсов, т/ч,

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

- содержание меди в штейне

где в кодированном виде представлены:

x₁₃=(X₁₃-200)/50, где X₁₃ - удельный расход кислорода, м³/т;

x₁₄=(X₁₄-0,125)/0,125, где X₁₄ - доля богатых оборотов от МС, б/р;

x₁₅=(X₁₅-1340)/40, где X₁₅ - температура расплава, °С;

x₁₆=(X₁₆-60)/5, где X₁₆ - содержание кислорода в КВС, %;

x₇=(X₇-2400)/1400, где X₇ - расход природного газа, м³/ч;

x₁₁=(X₁₁-0,125)/0,125, где X₁₁ - доля техногенных от МС, б/р;

Y₄ - содержание меди в штейне, %;

- содержание меди в шлаке

где в кодированном виде представлены:

x₁₀=(X₁₀-60)/40, где X₁₀ - загрузка МС, т/ч;

x₁₇=(X₁₇-28)/4, где X₁₇ - содержание диоксида кремния в шлаке, %;

x₁₈=(X₁₈-0,6)/0,3, где X₁₈ - качество загрузки, б/р;

x₁₃=(X₁₃-200)/50, где X₁₃ - удельный расход кислорода, м³/т;

x₁₄=(X₁₄-0,125)/0,125, где X₁₄ - доля богатых оборотов от МС, б/р,

Y₅ - содержание меди в шлаке, %,

- содержание диоксида кремния в шлаке

x₂₁=(X₂₁-9)/3, где X₂₁ - загрузка флюсов измеренная или рассчитанная по (3), т/ч;

x₂₂=(X₂₂-5)/5, где X₂₂ - загрузка руды, т/ч;

x₂₃=(X₂₃-3)/3, где X₂₃ - загрузка лежалых измеренная или рассчитанная как 0,05*Y₂, т/ч;

x₂₄=(р-1)/1, где р - количество залитых ковшей конвертерного шлака;

x₂₅=(X₂₅-60)15, где X₂₅ - содержание кислорода в КВС, %;

Y₆ - содержание диоксида кремния в шлаке, %,

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

где в кодированном виде представлены:

x₁₁=(X₁₁-0,125)/0,125, где X₁₁ - доля техногенных от МС, б/р;

x₁₂=(X₁₂-1)/1, где X₁₂ - отношение обороты/лежалые, б/р;

x₁=(X₁-5)/5, где X₁ - отклонение загрузки флюсов от заданной нормы, т/ч;

x₂₁=(X₂₁-15)/5, где X₂₁ - влажность шихты, %;

x₂₂ - содержание меди в оборотах, качественная переменная (низкое «-1», высокое «+1»);

Y₇ - качество загрузки, б/р, в интервале от 0,3 до 0,9 при средней 0,6,

и при сочетаниях Q_i∈Q₀ и уровне согласованности выше или равному регламентируемому 0,45 процесс продолжают в основном режиме, а при ниже регламентируемого - идентифицируют конфликт и корректируют изменением составляющих загружаемой шихты путем изменения входных переменных до достижения уровня согласованности по (1) не ниже 0,45 при требуемом уровне содержания меди в штейне, и при расхождении расчетного содержания меди в штейне от заданного на величину 1-1,5% пересчитывают величины управляющих воздействий по (2) и (3) до достижения заданного диапазона по меди в штейне по модели (4) с отслеживанием расчетных значений меди в шлаке по модели (5) ниже критического 0,8;

при Q_i∈Q₊₁ идентифицируют конфликт K₁, и/или K₄, и/или K₅,

при ∈ идентифицируют конфликт K₂, и/или K₃, и/или K₆,

анализируют разрешения конфликтов по (1)-(7) и по результатам корректируют процесс изменением составляющих загружаемой шихты путем изменения значений входных переменных до достижения области Q₀;

при возникновении конфликта K₇ значение температуры расплава в реакционной зоне печи аппроксимируют полиномиальной моделью (8)

где в кодированном виде представлены:

x₁₀=(X₁₀-60)/40, где X₁₀ - загрузка МС, т/ч;

x₁₃=(X₁₃-200)/50, где X₁₃ - удельный расход кислорода, м³/т;

x₁₉=(X₁₉-0,15)/0,15, где X₁₉ - доля инертных, б/р;

x₇=(X₇-2400)/1400, где X₇ - расход природного газа, м³/ч;

x₂₀=(X₂₀-10)/10, где X₂₀ - влажность шихты, %;

Y₈ - температура расплава в реакционной зоне, °C,

при возникновении конфликта K₈

физический объем поступающей шихты на плавку рассчитывается с учетом его влажности для каждого питателя отдельно по формуле (8) с последующим суммированием по всем питателям

где γ объемный вес продукта, т/м³;

V=36V_max, м/ч;

V_max=ω*r - максимальная скорость движения ленты питателя, устанавливаемая экспериментально или по паспорту устройства, м/с;

ω=πn/30 - угловая скорость вращения барабана, с^-1;

r - радиус барабана, мм.

n - число оборотов барабана;

S - площадь сечения призмы шихтового материала на ленте питателя, м²;

36 - коэффициент перевода размерности скорости из «м/с» в «м/ч».

Изобретение иллюстрируется чертежами Фиг.1, Фиг.2 и Фиг.3.

На представленных чертежах:

Фиг.1 - блок-схема автоматического управления, построенная на базе интеллектуальной системы автоматического управления процессом Ванюкова (ИАСУ ПВ), где

1 - печь Ванюкова, в которой реализуется процесс Ванюкова;

2 - течки для подачи шихты в реакционную зону печи;

3 - сборочный конвейер шихтовых материалов из бункеров с информационными каналами

3.1 - насыпного веса в т/час,

3.2 - скорости движения ленты;

4 - сборочный конвейер шихтовых материалов из бункеров с информационными каналами

4.1 - насыпного веса в т/час,

4.2 - скорости движения ленты в м/с;

5 - бункеры шихтовых материалов:

5.1-5.4 - металлосодержащие (МС),

5.5-5.6 - флюсы,

5.7 - обороты,

5.8 - уголь;

с информационными каналами

5.1.1-5.8.1 - насыпного веса в т/час,

5.1.2-5.8.2 - скорости движения ленты в м/с;

6 - информационный канал подачи кислорода в смеситель (не показан) в м³/час;

7 - информационный канал подачи воздуха в смеситель в м³/час;

8 - подача кислородно-воздушной смеси (КВС) из смесителя с информационными каналами

8.1 - подачи расхода КВС в реакционную зону печи в м³/час,

8.2 - содержания кислорода в КВС в %;

9 - подачи природного газа в реакционную зону печи в м³/час;

10 - информационный канал заливки конвертерного шлака в реакционную зону печи в количестве ковшей;

11 - содержание меди в штейне с информационными каналами

11.1 - расчетное по модели в оперативном режиме в %,

11.2 - по данным химических анализов периодического опробования штейна в %;

12 - содержание меди в шлаке с информационными каналами

12.1 - расчетное по модели в оперативном режиме в %,

12.2 - по данным химических анализов периодического опробования шлака в %;

13 - перепад температур воды на выходе и входе кессонов 1-2 рядов в °C;

14 - температура расплава в °C;

15 - блок распознавания и предварительной обработки входного потока информации о состоянии процесса;

16 - блок переключения системы на ручное управление;

17 - блок расчетов по предлагаемому алгоритму при работе в автоматическом режиме;

18 - АРМ оператора;

19 - блок управления подачей природного газа;

20 - блок управления подачей МС;

21 - блок управления подачей флюсов;

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

Фиг.2 - схема к обоснованию возникновения конфликтов, где

зона 1 - конфликты K₃ или K₄ из-за избытка или недостатка по отношению к металлосодержащим флюсующих, обусловливающих низкое качество шихты;

зона 2 - конфликты K₁ или K₂ при неблагоприятных сочетаниях качества загрузки с подачей дутья и/или избытком подачи природного газа;

зона 3 - конфликты K₁, K₂, K₅, K₆ при недостатке или избытке дутья;

зона 4 - конфликты K₁-K₆ при неблагоприятных сочетаниях загрузки и дутья;

зона 5 - конфликты K₁-K₄ при недостатке или избытке загружаемых шихтовых материалов;

зона 6 - конфликты K₃-K₆ при неблагоприятных сочетаниях загрузки и ее качества, которые могут возникать при необходимости переработки нестандартных материалов.

Фиг.3 - интерфейс ИАСУ ПВ на АРМ оператора, где в реальном масштабе времени отражена основная информация о течении процесса:

поле 1 - визуализирует значения расчетов по основной переменной «Удельный расход кислорода на тонну металлосодержащих» в м³/т;

поле 2 - отражает прогностическую способность ИАСУ ПВ-3 по содержанию меди в штейне, рассчитываемой по модели (4) в %;

поле 3 - расчетное содержание диоксида кремния в шлаке по модели (6) в %;

индикатор 4 - содержание меди в шлаке, рассчитываемое по модели (5) в %;

индикатор 5 - процент флюсов от металлосодержащих в загрузке;

индикатор 6 - расчетное значение качества загрузки по модели (7), б/р;

индикатор 7 - температура расплава, измеренная или расчетная по (8) в °C;

поле 8 - почасовые значения расхода шихтовых материалов по бункерам, рассчитанные по (9) с учетом их влажности в тоннах;

поле 9 - отражает названия конфликтов, имеющих место в текущий момент времени;

поле 10 - радиокнопки выбора режимов управления по переменной «Удельный расход кислорода на тонну металлосодержащих» в м³/т;

поле 11 - радиокнопка режима заливки конвертерного шлака.

Сущность заявляемого технического решения заключается в следующем.

Схема для осуществления способа автоматического управления качеством конечных продуктов процесса плавки в печи Ванюкова при переработке сульфидной шихты, представленная на Фиг.1, включает печь Ванюкова (ПВ) 1, в своде которой размещены течки 2 для подачи шихты в реакционную зону. Загрузочный тракт загрузки в печь 1 состоит из сборочных конвейеров 3 и 4, бункеров шихтовых материалов 5.1-5.8 с питателями. Измерительные каналы 3.1 и 4.1 связаны с конвейерными весами (на схеме не показаны), отображающими текущие веса шихтовых материалов на конвейерах 3 и 4, связаны с тиристорными преобразователями (не показан), питающими электродвигатели 3.2 и 4.2 для определения текущих скоростей сборочных конвейеров.

Измерительные каналы 5.1.1-5.8.1 связаны с конвейерными весами (на схеме не показаны) питателей бункеров 5.1-5.8 для отображения текущих весов шихтовых материалов, связаны с тиристорными преобразователями (не показаны), питающими электродвигатели 5.1.2-5.8.2, для определения текущих скоростей питателей.

Измерительные каналы 6 и 7 связаны с расходомерами (на схеме не показаны) подачи в расплав кислорода и воздуха в виде кислородно-воздушной смеси (КВС) 8 и природного газа 9 соответственно. При этом измерительные каналы 8.1 и 8.2 отображают общий расход КВС в м³/час и содержание кислорода в КВС в %. Измерительный канал 9, связанный с расходомером (на схеме не показан), отображает общий расход природного газа в м³/час в плавильную зону.

Измерительный канал 10 связан с заливкой конвертерного шлака в ванну расплава. Факт заливки фиксируется вручную. Измерительные каналы 11.1 и 12.1 предназначены для определения расчетных содержаний меди в штейне и шлаке в % по моделям (4) и (5) соответственно, 11.2 и 12.2 - путем производства химического анализа периодически отбираемых продуктов для корректировки процесса. Измерительные каналы 13 и 14 предназначены для определения теплового состояния печи путем определения перепада температур в °C на выходе и входе 1-го и 2-го ряда кессонов, а 14 - для определения температуры расплава в реакционной зоне печи в °C.

Все измерительные каналы 3.1, 3.2, 4.1, 4.2, 5.1.1-5.8.1, 5.1.2-5.8.2, 6, 7, 81, 8.2, 9, 11.1, 12.1, 13, 14 предназначены для получения информации о мгновенных значениях соответствующих параметров, имеют прямой выход на модуль 15 и вместе с периодически определяемыми 11.2, 12.2 служат для распознавания входного потока после предварительной обработки измеренных основных параметров, определения принадлежности к одной из областей по косвенному параметру Q_i и/или установления наличия конфликта. Модуль 15 связан с переключающим модулем 16, реализующим возможность включения либо режима автоматического управления с помощью модуля 17 расчетов по предлагаемому способу, либо его выключение и перевод всей информации на монитор 18 АРМ оператора для ручного управления. Модуль расчетов 17 связан с устройствами 19, 20 и 21 выработки управляющего воздействия дутья и загрузки шихтовых материалов соответственно. В свою очередь, они связаны с электродвигателями 5.1.2-5.8.2 питателей шихтовых материалов для изменения скоростей питателей и задвижек (не показаны) подачи кислорода 6, воздуха 7 к смесителю (не показан), общего расхода КВС 8.1 и содержания кислорода в нем 8.2. При этом оператор при переходе на ручное управление по переключающему модулю 16 также может воздействовать по каналам связи 22 на скорости загрузки шихтовых материалов 5.1.2-5.8.2 и параметров КВС 6,7, 8.1, 8.2 и природного газа 9 непосредственно, минуя блоки управления 19, 20 и 21.

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

В печи 1, Фиг.1, с помощью технологических операций, включающих заливку расплава выше уровня фурм, подачи дутья 6, 7, 8.1, 8.2 и 9 через фурмы, загрузку шихты из бункеров 5.1-5.8 через течки 2 конвейерами 3 и 4, устанавливают режим, при котором начинается интенсивное плавление сульфидной шихты с выделением значительного тепла за счет протекания экзотермических реакций. Посредством измерительных каналов 3.1, 3.2, 4.1, 4.2, 5.1.1-5.8.1, 5.1.2-5.8.2, 6, 7, 8.1, 8.2, 9, 11.1, 12.1, 13, 14 определяют скорости загрузки шихты, количество и процент обогащения кислородом КВС, количество потребного природного газа и доводят процесс плавления до выпуска через сифоны готовой продукции в виде штейна 11 и шлака 12 с последующим отбором проб для определения химического состава. Данные о состоянии процесса плавления поступают на вход модуля предварительной обработки входной информации 15 и определения текущей ситуации. Если ситуация не требует переключения в ручной режим работы печи, переключающий модуль 16 переводит работу на вход модуля 17, в котором производятся вычисления по прогностическим полиномам согласно предлагаемому способу.

Пример идентификации конфликтов показан на фиг.2.

Множество фактически встречающихся режимов при эксплуатации процесса в промышленных условиях гораздо шире «области согласованности» и полностью покрывает ее, как схематически показано на Фиг.2. При этом, находясь в различных зонах допускаемых фактических режимов, но вне «области согласованности», процесс «испытывает» постоянные «конфликты» при стремлении достичь этой зоны. Так, при идентификации состояния процесса в зоне 1 понятие «конфликт» может выражать низкое качество шихты по металлосодержащим из-за избытка - конфликт K₃ или недостатка - K₄ флюсующих, в зоне 2 - конфликты K₁ или K₂ при неблагоприятных сочетаниях качества загрузки с подачей дутья и/или избытком подачи природного газа в зону плавления, в зоне 3 - конфликты K₁, K₂, K₅, K₆ при недостатке или избытке дутья, в зоне 4 - конфликты K₁-K₆ при неблагоприятных сочетаниях загрузки и дутья, в зоне 5 - конфликты K₁-K₄ при недостатке или избытке загружаемых шихтовых материалов, в зоне 6 - конфликты K₃-K₆ при неблагоприятных сочетаниях загрузки и ее качества, которые могут возникать при необходимости переработки нестандартных материалов.

На Фиг.3 приведен пример, отражающий в реальном масштабе времени основную информацию о течении процесса на интерфейсе И АСУ ПВ. Графики окон 1-3 представляют изменение переменных на двухчасовом интервале, который приурочен к дискретности отбора проб и проведения химанализа продуктов плавки.

В поле 1 визуализированы графики изменения переменной «Удельный расход кислорода на тонну металлосодержащих»:

А - при ведении процесса рядовым оператором;

Б - значения расчетов управления по предлагаемому способу.

В поле указаны также нижний 150 и верхний 250 пределы удельного расхода кислорода в м³/т, выход за которые приводит к конфликтам K₂, K₆ и K₁, K₅ соответственно. Так, до 17:15 значения переменной были существенно ниже 180, доходя кратковременно даже до 150 м³/т. Визуализация управления позволила оператору принять управляющие решения для поддержания более высоких значений удельного расхода, не выходящих за ограничительные пределы. Такое ведение процесса вызывает ответную реакцию по содержанию меди в штейне, рассчитываемой по модели (4), как показано на графике А в окне 2.

При управлении по предлагаемому способу процесс не испытывает таких скачков ни по удельному расходу (поле 1 график Б), ни по содержанию меди в штейне (поле 2 график Б).

Наблюдаемые пики в интервале времени от 17:51 до 18:03 по обеим переменным обусловлены технологической операцией фурмования - прочистки устья фурм при их «затягивании» затвердевшим расплавом.

К 18:40 значения графиков А и Б в поле 1 сблизились, что повлекло за собой сближение графиков А и Б в поле 2 по содержанию меди в штейне. Фактическое значение меди в штейне 56%, график B в поле 2, является отражением информации по каналу 11.2 при получении информации о результатах химанализа отобранной пробы штейна. Пониженное содержание меди в штейне против прогнозируемого 60% является результатом ведения процесса по удельному расходу (окно 1, график А) в заниженном режиме по сравнению с установленным 200 м³/т, как показано в поле 10. Химанализ отражает в данном случае осредненную за два часа величину содержания меди в штейне.

Поле 3 визуализирует прогнозные значения содержания диоксида кремния в шлаке по модели (6) в %. По данной информации оператор отслеживает качество шлака в реальном режиме времени.

В качестве индикаторов на экраны выведены:

поле 4 - содержание меди в шлаке, рассчитываемое по модели (5) в %;

поле 5 - процент флюсов от металлосодержащих в загрузке;

поле 6 - расчетное значение качества загрузки по модели (7), б/р;

поле 7 - температура расплава, измеренная или расчетная по (8) в °C;

поле 9 - отражает названия конфликтов, так в текущий момент времени идентифицирован конфликт K₃ - избыток флюсующих;

поле 10 - радиокнопки режимов управления по удельному расходу кислорода на тонну металлосодержащих в м³/т позволяют менять тактику расчетов по моделям в зависимости от производственного задания, так выставленный режим 200 м³/т соответствует содержанию меди в штейне 55-60%;

поле 11 - радиокнопка режима заливки конвертерного шлака, содержание диоксида кремния в котором существенно ниже, реализует действие информационного канала 10, влияющего на достоверность определения содержания диоксида кремния в отвальном шлаке ПВ.

Особенность поля 8 состоит в определении почасовых значений расхода шихтовых материалов по бункерам или по (9), если весовое хозяйство некорректно отражает загрузку печи. Это предоставляет возможность повышения точности и достоверности балансов материальных потоков шихтовых материалов посменно.

Заявляемое техническое решение «Способ автоматического управления качеством конечных продуктов процесса плавки в печи Ванюкова при переработке сульфидной шихты на штейн» реализовано в виде программного обеспечения Интеллектуальной АСУ для действующей печи Ванюкова (ИАСУ ПВ) на Медном заводе Заполярного филиала ОАО «ГМК «Норильский никель».

Выводы

1. Ввод в промышленную эксплуатацию на Медном заводе ЗФ ОАО «ГМК «Норильский никель» интеллектуальной автоматизированной системы для управления процессом плавки в печи Ванюкова ИАСУ ПВ по переработке сульфидного медно-никелевого сырья, использующей заявляемый способ, подтвердил, что в условиях нечеткости исходной информации оправданным является подход создания экспертной системы ИАСУ ПВ на основе синтеза двух уникальных методологий: построение сложных многофакторных моделей в виде семантической сети на специфичном ограниченном онтологическом словаре и мультимодельной системы аналитических выражений, полученных с использованием неявных экспертных знаний о процессе.

2. Заявляемый способ и построенная на его основе ИАСУ ПВ обладают рядом специфических преимуществ, главными из которых являются следующие:

- возможность оценивания состояния процесса плавки, реализуемого в агрегате ПВ, по многим критериям с идентификацией конфликтов при их возникновении и разрешения каждого конфликта индивидуально в режиме реального времени;

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

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

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

Способ автоматического управления процессом плавки медно-никелевого сульфидного сырья в печи Ванюкова при переработке сульфидной шихты на штейн, включающий постоянный контроль параметров процесса, корректировку управляющих параметров для стабилизации содержания меди, причем в качестве параметра выбирают расход шихты, производят разбивку его диапазона на области и в зависимости от того, в какой области находится параметр, осуществляют изменение расхода шихты и технического кислорода до достижения соответствующей области, дополнительно определяют соотношение расхода технического кислорода дутья на тонну шихты, отличающийся тем, что в качестве основного параметра Q_i выбирают удельный расход кислорода на тонну только металлосодержащих в шихте в м³/т, производят разбивку общего диапазона Q_i на ряд областей, соответствующих нормальному течению процесса Q₀ в интервале 150-250 м³/т, склонности к холодному течению процесса Q_-1 в интервале 150 м³/т и ниже и склонности к горячему течению процесса Q₊₁ в интервале 250 м³/т и выше, дополнительно учитывают степень согласованности процесса, причем выход значений параметров процесса за границы допустимых режимов области согласованности интерпретируют как конфликт K_i и по корректировочным полиномиальным моделям ведут поиск вариантов выхода процесса в область согласованности, а в качестве конфликтов выбраны:
K1 - переокисление металлосодержащих;
K2 - недоокисление металлосодержащих;
K3 - избыток флюсов;
K4 - недостаток флюсов;
K5 - горячий ход печи;
K6 - холодный ход печи;
K7 - невозможность прямого определения температуры расплава в реакционной зоне печи;
K8 - невозможность прямого определения физического объема поступающей шихты, при этом проверяют степень согласованности величины загрузки и дутьевых режимов с заданным содержанием меди в штейне по полиномиальной модели (1):

Y1=0,4383-0,0352*x1-0,1008*x2-0,0914*x3-0,0352*x4-0,0164*x5
+0,0164*x1*x2+0,0258*x1*x3+0,0539*x2*x3+0,0164*x2*x4
+0,0164*x2*x5+0,0258*x3*x4+0,0164*x3*x5-0,0164*x2*x3*x5,	(1)

где в кодированном виде представлены:
x₁=(X₁-5)/5, где X₁ - отклонение загрузки флюсов от заданной нормы, т/ч;
х₂=(Х₂-6)/6, где Х₂ - нормированное отклонение удельного расхода кислорода на тонну МС, б/р;
x₃=(Х₃-5)/5, где Х₃ - отклонение расчетного содержания меди в штейне от нормируемого, %;
x₄=(X₁-0,6)/0,2, где Х₄ - расчетное содержание меди в шлаке, %;
x₅=(Х₅-1340)/40, где Х₅ - температура расплава, °С;
Y₁ - степени согласованности, б/р,
корректировку параметров процесса с визуализацией на АРМ оператора показаний основных параметров в виде графиков и/или на индикаторах по полиномиальным моделям (2)-(5):
- загрузка металлосодержащих (МС)

Y2=49,875+26,375*x6-4,875*х7+3,875*х8+2,625*х9-1,375*x6*х7
+2,375*x6*x8+3,625*x6*x9-1,375*х7*x8+l,125*x8*x9
-1,125*x6*x7*x9-1,125*x7*x8*x9,	(2)

где в кодированном виде представлены:
x₆=(Х₆-15000)/10000, где Х₆ - расход кислорода, м³/ч;
x₇=(Х₇-2400)/1400, где Х₇ - расход природного газа, м³/ч;
x₈=(Х₈-0,05)/0,05, где X₈ - доля лежалых материалов в МС, б/р;
x₉=(Х₉-17)/8, где Х₉ - перепад температуры воды на сливе 1-2 ряда кессонов, °C;
Y₂ - уставка скорости загрузки металлосодержащих, т/ч;
- загрузка флюсов

Y3=9,25+5,25*х10-0,75*х12+1,25*x13-0,75*x10*x12+l,25*x10*x13
-0,5*x11*x13-0,25*x12*x13-0,5*x10*x11*x13-0,25*x10*x13*x14,	(3)

где в кодированном виде представлены:
x₁₀=(Х₁₀-60)/40, где Х₁₀ - загрузка МС, т/ч;
x₁₁=(Х₁₁-0,125)/0,125, где Х₁₁ - доля техногенных материалов от МС, б/р;
x₁₂=(X₁₂-1)/1, где Х₁₂ - отношение оборотные/лежалые материалы;
х₁₃=(Х₁₃-200)/50, где Х₁₃ - удельный расход кислорода, м³/т;
Y₃ - уставка скорости загрузки флюсов, т/ч,
текущие расчетные значения содержаний меди в конечных продуктах плавки определяют
по полиномиальным моделям:
- содержание меди в штейне

Y4=57,820+5,391*x13+0,875*x14+0,820*x15-0,273*x16-0,547*х7
+0,547*х11-0,383*х13*x7-0,164*x13*x11-0,164*x14*x7-0,164*x14*x11
-0,164*x15*x16-0,219*x15*x11-0,766*x16*x7+0,273*x13*x14*x16
+0,711*x13*x16*x7+0,164*x13*x16*х11-0,164*x14*x15*x11
-0,164*x14*x16*x11-0,219*x15*x16*x11,	(4)

где в кодированном виде представлены:
х₁₃=(Х₁₃-200)/50, где Х₁₃ - удельный расход кислорода, м³/т;
х₁₄=(Х₁₄-0,125)/0,125, где Х₁₄ - доля богатых оборотных материалов от МС, б/р;
x₁₅=(Х₁₅-1340)/40, где Х₁₅ - температура расплава, °C;
x₁₆=(X₁₆-60)/5, где X₁₆ - содержание кислорода в КВС, %;
x₇=(Х₇-2400)/1400, где Х₇ - расход природного газа, м³/ч;
х₁₁=(Х₁₁-0,125)/0,125, где Х₁₁ - доля техногенных материалов от МС, б/р;
Y₄ - содержание меди в штейне, %;
- содержание меди в шлаке

Y5=0,7088+0,0056*х10-0,0619*х17-0,0225*x18-0,0169*х13+0,0225*х14
+0,0056*x10*x18-0,0338*x10*x13+0,0056*x17*x18+0,0056*x18*x13
-0,0113*x10*x17*x18-0,0169*x10*x17*x13+0,01*x17*x18*x13,	(5)

где в кодированном виде представлены:
х₁₀=(Х₁₀-60)/40, где Х₁₀ - загрузка МС, т/ч;
х₁₇=(X₁₇-28)/4, где Х₁₇ - содержание диоксида кремния в шлаке, %;
x₁₈=(X₁₈-0,6)/0,3, где X₁₈ - загрузка, б/р;
х₁₃=(Х₁₃-200)/50, где Х₁₃ - удельный расход кислорода, м³/т;
x₁₄=(X₁₄-0,125)/0,125, где Х₁₄ - доля богатых оборотных материалов от МС, б/р,
Y₅ - содержание меди в шлаке, %,
- содержание диоксида кремния в шлаке

Y6=27,44+1,31*x21+0,63*х22-0,19*х23-1,0*х24-1,06*х25
+0,19*x21*х23-0,25*x21*х24+0,13*х22*х23+0,13*х22*Х5
+0,25*х23*х24+0,19*х23*х25+0,38*x24*x25-0,19*x21*x23*x25	(6)

x₂₁=(X₂₁-9)/3, где X₂₁ - загрузка флюсов, измеренная или рассчитанная по (3), т/ч;
х₂₂=(Х₂₂-5)/5, где Х₂₂ - загрузка руды, т/ч;
х₂₃=(Х₂₃-3)/3, где Х₂₃ - загрузка лежалых материалов, измеренная или рассчитанная как 0,05*Y₂, т/ч;
х₂₄=(р-1)/1, где р - количество залитых ковшей конвертерного шлака;
х₂₅=(Х₂₅-60)/5, где Х₂₅ - содержание кислорода в КВС, %;
Y₆ - содержание диоксида кремния в шлаке, %,
и при сочетаниях Q_i∈Q₀ и уровне согласованности выше или равному регламентируемому 0,45 процесс продолжают в основном режиме, а при ниже регламентируемого - идентифицируют конфликт и корректируют изменением составляющих загружаемой шихты путем изменения входных переменных до достижения уровня согласованности по (1) не ниже 0,45 при требуемом уровне содержания меди в штейне, и при расхождении расчетного содержания меди в штейне от заданного на величину 1-1,5% пересчитывают величины управляющих воздействий по (2) и (3) до достижения заданного диапазона по меди в штейне по полиномиальной модели (4) с отслеживанием расчетных значений меди в шлаке по полиномиальной модели (5) ниже критического 0,8;
при Q_i∈Q₊₁ идентифицируют конфликт K₁, и/или K₄, и/или K₅,
при Q_i∈Q_-1 идентифицируют конфликт K₂, и/или K₃, и/или K₆,
анализируют разрешения конфликтов по (1)-(6) и по результатам корректируют процесс изменением составляющих загружаемой шихты путем изменения значений входных переменных до достижения области Q₀;
при возникновении конфликта K₇ значение температуры расплава в реакционной зоне печи аппроксимируют полиномиальной моделью (7)

Y7=1322,47+2,84*Х1+13,59*Х2-5,91*Х3+7,34*Х4
+5,09*Х1*Х2+5,09*Х1*Х4-1,91*Х3*Х4-1,97*Х3*Х5
+1,97*Х1*Х2*Х3-1,91*Х1*Х2*Х4,	(7)

где в кодированном виде представлены:
х₁₀=(Х₁₀-60)/40, где Х₁₀ - загрузка МС, т/ч;
x₁₃=(X₁₃-200)/50, где Х₁₃ - удельный расход кислорода, м³/т;
x₁₉=(X₁₉-0,15)/0,15, где Х₁₉ - доля инертных, б/р;
Х₇=(Х₇-2400)/1400, где Х₇ - расход природного газа, м³/ч;
х₂₀=(Х₂₀-10)/10, где Х₂₀ - влажность шихты, %;
Y₇ - температура расплава в реакционной зоне, °С,
при возникновении конфликта K₈
физический объем поступающей шихты на плавку рассчитывают с учетом его влажности для каждого питателя отдельно по формуле (8) с последующим суммированием по всем питателям

F=γVS,

(8)

где γ - объемный вес продукта, т/м³;
V=36 V_max, м/ч;
V_max=ω*r - максимальная скорость движения ленты питателя, м/с;
ω=πn/30 - угловая скорость вращения барабана, с^-1;
r - радиус барабана, мм;
n - число оборотов барабана;
S - площадь сечения призмы шихтового материала на ленте питателя, м²;
36 - коэффициент перевода размерности скорости из м/с в м/ч.