Адаптивная система управления содержанием оксида кремния в шлаках при переработке медно-никелевых руд

УДК 004.8+669+519.7 Дата подачи статьи: 27.05.2014

АДАПТИВНАЯ СИСТЕМА УПРАВЛЕНИЯ СОДЕРЖАНИЕМ ОКСИДА КРЕМНИЯ В ШЛАКАХ ПРИ ПЕРЕРАБОТКЕ МЕДНО-НИКЕЛЕВЫХ РУД

С.Г. Фомичева, к.т.н., профессор, зав. кафедрой; А.В. Конев, аспирант

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

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

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

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

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

концентрате, SiO2 в концентрате, а также корректно учитывать при прогнозировании ряд ретроспективных параметров: в штейне (истор.), Fe в штейне (истор.), SiO2 в штейне (истор.), в шлаке (истор.), Fe в шлаке (ис-тор.).

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

Received 27.05.2014

ADAPTIVE CONTROL SYSTEM FOR SILICON OXIDE CONCENTRATION IN SLAGS AT PROCESSING COPPER-NICKEL ORES

Fomicheva S. G., Ph.D. (Engineering), Professor, Head of Chair; Konev A. V., Postgraduate Student (NorilskIndustrial Institute, 50 let Oktyabrya St. 7, Norilsk, 663310, Russian Federation, ist@norvuz.ru) Аbstract. The article presents the results of designing self-adaptive regulators based on the fuzzy-neural networks. Such networks are used to balance the content of silicon oxide in slags when processing copper-nickel ores to copper and nickel cathodes. The paper shows that these metallurgical processes are highly inertial and demand two-level regulation system: rough and exact adjustment. Thus, an exact adjustment can adequately be carried out in completely automatic mode reducing operational loading for actions of a technologist-dispatcher. Operating influences of an independent adaptive regulator provide higher balancing accuracy comparing to manual management. There is a structure of a hybrid fuzzy-neural network which is a base for an adaptive regulator development. The paper shows revealed fuzzy tendencies for dependences of the silicon oxide content in slags from operated parameters of high orders of nonlinearity. They are experimentally proved.

The configuration of an author's hybrid network structure includes five interlocked fuzzy-neural subnets with feedback coupling and delay elements. Delay elements simulate the time lag displaying processes iteration depth. The lag in turn also is a fuzzy variable. The given configuration allowed realizing adequately soft operating influences on such parameters as the general sand content; a dust flow; a blasting enrichment; Ni+Cu in a concentrate; S in a concentrate; Fe in a concentrate; SiO2 in a concentrate. Also it is correct to consider a number of retrospective parameters at forecasting: Ni+Cu in matte (history); Fe in matte (history); SiO2 in matte (history); Ni+Cu in slag (history); Fe in slag (history).

Keywords: metallurgical processes automation, self-adaptive regulators, fuzzy-neural networks, the intelligent analysis of multivariate data, copper-nickel ores.

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

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

13

— FeS +--Fe3O4 ^

10 10 34

^ FeO + ilo SO2 - AH298 50,58 кДж/(гхмоль),

FeO+^SiO2 ^ 22

^ ^Fe2SiO4 +AH°29S5,43 кДж/(гхмоль),

ZnO + ^SiO2 2 2

^ ^Zn2SiO4 +ДЯ2°9814,63 кДж/(гхмоль),

2

1

PbO + - SiO2 22

^-РЪ+ДЯ2°9814,63 кДж/(гхмоль),

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

Оксид железа, образующийся при окислении концентрата, флюсуется с использованием SiO2 для получения жидкого шлака. Количество SiO2 выбирается таким образом, чтобы обеспечить низкую растворимость № и ^ в шлаке и текучесть шлака, достаточную для его простого выпуска и четкого разделения штейн/шлак [2]. Увеличение содержания оксида кремния в шлаках снижает их плотность, следовательно, увеличивается разность в плотностях штейна и шлака. Такое же влияние

на плотность шлака оказывает и повышение температуры. В среднем с повышением температуры шлака на 100 °С плотность его прямолинейно уменьшается на 0,2-0,3 г/см3.

Следует отметить, что содержание SiO2 влияет и на вязкость шлака. При содержании оксида кремния выше 36 % начинает снижаться растворимость цветных металлов в шлаке. Оптимальным считается содержание оксида кремния на уровне 32-35 %.

Задача поддержания оксида кремния в обозначенных границах лежит на сменном технологе и в большей степени зависит от его квалификации. Проведенный авторами анализ современного состояния вопроса построения математических моделей с целью обеспечения управления, в том числе автоматизированного, сложными многомерными металлургическими процессами, и ПВП в частности, показывает, что проводившиеся до последнего времени работы в основном сводились к разработке термодинамических, гидродинамических и балансовых прогнозных моделей [3-5]. Данные модели служили, как правило, для определения границ безаварийной работы технологических агрегатов, их конструктивных особенностей, параметров основных технологических режимов, то есть носили статический характер. Соответственно, из-за отсутствия в данных моделях динамических составляющих создание на их базе управляющих моделей невозможно в принципе. При построении АСУ для ПВП ограничиваются, как правило, созданием контуров стабилизации входных материальных потоков, а всю ответственность за безаварийность оперативного управления и качество выходного продукта возлагают на оператора-технолога.

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

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

В данной статье основное внимание уделено структуре экспертного модуля, на базе которого

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

Постановка задачи. Формализация входных параметров модели

Классический подход к построению систем управления объектами предлагает макроописание системы «Внешняя среда - Объект управления -Система управления», представленное на рисунке 1, где объект управления (ОУ) рассматривается отдельно от системы управления (СУ), блока датчиков, системы распознавания (СР) и исполнительных устройств (ИУ).

Возмущающее Среда ¡воздействие

ОУ Блок — СР Целевая

датчиков V* функция

ИУ СУ

Рис. 1. Классическое представление системы «Внешняя среда - Объект управления - Система управления»

Fig. 1. Classical representation of "Environment -Controlled object - Controlling system " system

В действительности ОУ (в данном случае печь взвешенной плавки), как правило, является частью среды Ж, в которой среда ( обусловлена человеческим фактором (рис. 2). Следовательно, с точки зрения управляющей системы среда Ж вовлечена в цикл управляемого воздействия.

Тем самым в определенном смысле для СУ стирается грань между собственно ОУ и средой, поскольку СУ управляет системой «Внешняя среда - ОУ - СУ» [6]. Это требует предусмотреть и определить не один, а несколько контуров управления как маршруты в системе, по которым могут проходить воздействия и информация. В частности, минуя среду ( (сменного технолога), предусмотреть прямые связи между блоком датчиков и исполнителем.

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

Среда U

Среда W

Среда Q

ОУ

Блок датчиков

ИО

СУ

Формирование и распознавание

образов

Аппарат оценки целевой функции

Формирование базы знаний

Принятие решений

Определение времени принятия решений

Рис. 2. Структура функциональных подсистем СУ Fig. 2. A structure of functional subsystems CS

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

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

6

р = 1 --

-Е (R - R )2

п(п -1 )(п + 1) ,=1 где Я* - ранг наблюдения х, в ряду х; К/ - ранг наблюдения у, в ряду у.

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

Таблица 1

Значения коэффициента взаимной корреляции

Table 1

The values of a cross-correlation coefficient

Анализируемый фактор Коэффициент взаимной корреляции с SiO2 в шлаке

Общий расход песка 0,211771

Расход шихты 0,049199

Обогащение дутья -0,243223

Скорость КВС 0,144039

Температура в отстойнике 0,144141

№+Си в концентрате 0,063082

S в концентрате 0,322642

Fe в концентрате 0,291398

SiO2 в концентрате -0,150673

№+Си в штейне (истор.) -0,167156

Fe в штейне (истор.) 0,207017

SiO2 в штейне (истор.) 0,048660

Ni+Cu в шлаке (истор.) -0,120593

Fe в шлаке (истор.) -0,018739

Для дальнейшего анализа и последующего обучения нейронной сети факторы были нормализованы путем наложения на шкалу [0, 1] по выра-

жению x = —;-;——

i max mm

a — a

где x, - i-й фактор из мно-

жества нормализованных входных и выходных параметров а,. Для проверки взаимного влияния данных факторов (параметров) на итоговое значение У (нормализованное содержание SiO2 в шлаке) в пакете программ Statistica'2011 были проведены

многофакторный и многомерный анализы и построены соответствующие уравнения регрессии. Наименьшие абсолютные и среднеквадратические ошибки на обучающей выборке продемонстрировало уравнение регрессии, соответствующее выражению (6). Объем исходной выборки для каждого xi составлял 8 707 значений (ежечасные показания значений наблюдения в течение года), уровень значимости коэффициентов уравнения 0,005. Незначимые взаимодействия были элиминированы по /-критерию Стьюдента:

7=0,857839+0,200159*Х1+0,03488*Х3--0,205173*Х4-0,093226*Х6+0,62798*Х7--0,321175*Х8-0,225513*Х9-0,07388*Х13--0,358075*Х14-0,572854*Х15--0Д4791*Х2*Х12-0,033176*Х5*Х11+ +0,123209*Х2*Х5*Х11,

где Х1 - общий расход песка; Х2 - расход шихты; Х3 - расход пыли; Х4 - обогащение дутья; Х5 -температура в отстойнике; Х6 - №+Си в концентрате; Х7 - S в концентрате; Х8 - Fe в концентрате; Х9 - SiO2 в концентрате; Х10 - А1203 в концентрате; Х11 - №+Си в штейне (истор.); Х12 - Fe в штейне (истор.); Х13 - SiO2 в штейне (истор.); Х14 - №+Си в шлаке (истор.); Х15 - Fe в шлаке (истор.).

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

В частности, коэффициент взаимной корреляции Спирмена (5) между действительным и прогнозным значением 7 (содержание оксида кремния), полученным данным уравнением регрессии, составил 0,523778, хотя список входных парамет-

ров экспертами, имеющими многолетнюю практику управления данным процессом, сомнению не подвергался и дополнительными факторами расширен не был. Следовательно, значения коэффициентов в полученных уравнениях регрессии динамически меняются, к тому же высококвалифицированные сменные технологи опытным путем научились их контролировать. Поскольку при аппроксимации нелинейных динамических систем нейро-нечеткие регуляторы хорошо зарекомендовали себя во многих промышленных средах [5-7], было принято решение использовать аппарат гибридных сетей для модуля экспертной системы, который прогнозирует содержание SiO2 в шлаке на несколько часов вперед, чтобы была реализована функция «профессионального совета» для не очень опытных сменных технологов в условиях чрезвычайно высокой ротации кадров в районах Крайнего Севера - это один из критичных факторов. А также разработать автономный адаптивный регулятор тонкой настройки, который оптимизирует управляющие воздействия в квазистационарных периодах процесса и тем самым позволяет частично снять проблему человеческого фактора в управлении металлургическим процессом. Задача управления при этом сводится к задаче нечеткой оптимизации, где целевой функцией Е является минимизация среднеквадратического отклонения регулируемого выхода 7 ^Ю2 в шлаке) от нечеткого эталона 7э «около 34 %» с нечеткими границами «не более ±2 %» на выборке объемом /:

1 t /- Л2

E=7 5 - 7 ' )

^ min.

Отметим, что среднеквадратическая невязка 7 от нечеткого нормализованного эталона (около 0,5) вследствие регулирующих воздействий высо-

0 001 0 001 0 000 0 000 0 000

mmiDiDr^r^cocot^t^oo

mmiDiDr^r^cocot^t^oO'

IT IT ID ID

Рис. 3. Прогноз содержания оксида кремния в шлаке ПВП с использованием регрессионной модели на час вперед

Fig. 3. The results offorecasting a content of silicon oxide in slag flash smelting furnace using regression model

on one hour ahead

0 001

0 001

0 001

0 001

0 000

0 000

коквалифицированных сменных технологов на наблюдаемой (годовой) выборке составила 0,005, а прогноз на основе (6) привел к значению сред-неквадратической невязки, равной 0,013.

Архитектура предлагаемой нейро-нечеткой сети

Для выявления существенных и несущественных входов нейро-нечеткой сети был использован метод средних нечетких кривых, предложенный Лином и Каннингэмом [8], преимущество которого в исследовании влияния каждого входа х, на выход У по отдельности, что дает возможность визуализировать результаты с помощью двумерных кривых, а также позволяет эксперту выполнять сравнительную оценку значимости входов и оценивать тенденцию влияния заданного х, на выход У. Помимо этого, с помощью метода средних нечетких кривых можно определять структуру сети и выбирать начальные значения ее параметров, обеспечивающих ее быструю настройку.

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

x ) = exp

f » л

Л

b

ширина окрестности

точки х,* определяется на основании параметра Ь,. Согласно Лину и Каннингэму, рекомендуемое значение Ь, составляет 20 % от ширины Ах, диапазона изменения соответствующей входной величины х,. Среднюю нечеткую кривую получаем, вычислив утест (х*) для произвольного входного х,* в соответствии с выражением

S xtk ) • Уt

y ( x. ) = -

s mean^ I ' m

S Xt)

где I - номер входа; k - номер измерения.

Средние нечеткие кривые для рассматриваемой системы сведены в таблице 2. Чем больше ве-

личина

^Уmean () |Уmean ()max Уmean ()п

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

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

В таблице 2 жирным шрифтом выделены наиболее значимые с точки зрения влияния на У входные параметры, при этом значимость истори-

ческих параметров подтверждает сильную инерционность рассматриваемого процесса.

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

Временной лаг также является нечеткой переменной, диапазон изменения которой лежит в границах от 2 до 8 часов для всех входных параметров, кроме Х5 - температуры в отстойнике ПВП. Значение временного лага для Х5 изменяется от 10 до 15 минут.

Описываемая нейро-нечеткая система, состав параметров которой (табл. 1) был подвергнут экспертной оценке и ранжирован, представлена в виде каскада адаптивных систем нечеткого вывода АОТБ (рис. 4).

Рис. 4. Иерархическая структура разработанной адаптивной нейро-нечеткой сети

Fig. 4. An hierarchical structure of the developed adaptive fUzzy-neural network

На основе экспертных оценок средних нечетких кривых были сформированы термы нечетких тенденций и построены правила, которые составляют ядро системы нечеткого вывода типа Сугено (рис. 4). A-Net - система нечеткого вывода для агрегирования данных параметров X1, X3, X4 и X6, содержит 81 правило; B-Net - по X7, X8 и X9, содержит 27 правил; C-Net агрегирует данные более

t=i

и

Входные параметры адаптивной нейро-нечеткой модели Input parameters of an adaptive fuzzy-neural model

Таблица 2 Table 2

№

Группа

Переменная

Aymean

Средняя нечеткая кривая

Технологические

Общий расход песка

0,041103

Расход шихты

0,037705

Расход пыли

0,109245

Обогащение дутья

0,091745

Температура в отстойнике ПВП

0,090034

0.2 0.4 0.6 0JB

Концентрат

Ni+Cu в концентрате

0,122571

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

0,128472

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

0,093292

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

0,132276

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

0,063013

[< 02 04 Об

02 0.4 0.6

Q-Э 04 0.6

Q-Э 0 А 0-6

Исторические данные по штейну

Ni+Cu в штейне

0,095809

Fe в штейне

0,066399

SiO2 в штейне

0,047669

-3 04  Е

[< 02 04 Об

Исторические данные по шлаку

Ni+Cu в шлаке

0,372507

Fe в шлаке

0,522363

SiO2 в шлаке

0,551556

0 ОД 0,4 0,6 0,8

1

2

3

4

высокого порядка нелинейности X2, X5, X11, X12 и содержит 243 правила; D-Net - система нечеткого вывода для агрегирования данных параметров X13, X14 и X15, содержит 27 правил.

F-Net - система нечеткого вывода для агрегирования векторов предыдущих сетей (нечетких

тенденций) и переменной Х16, содержит 5 входов и 108 правил. Элементы задержки Ь отмечают подаваемые на вход ретроспективные данные.

В иерархической структуре каждая нейро-нечеткая сеть типа АОТК имеет пять слоев: термы входных переменных, антецеденты нечетких

правил, нормализация степеней выполнения правил, заключения правил, агрегирование результата. В процессе обучения производилась настройка wij - весов абсцисс носителей Ь„ Уу - весов абсцисс высот х,* термов входных лингвистических переменных (/' = 1,16, ] = 1, Nterms), сит+1)к - коэффициентов линейных заключений Сугено (где т -количество И-связок в антецедентах), а также zik -

весов правил нечеткого вывода (/ = 1,16,

к = 1, Ыги^) для каждой сети методом обратного распространения ошибки в соответствии со следующими выражениями:

Wj (t +1) = Wj (t) -n

dE

dWj

dE

^

Ci(m+1)k (t + 1) = Ci(m+1)* (t) -П

dE

i(m+1)k v

dCi

(i +1) = Zk « -n

dE dzik

где п - скорость обучения сети, ^у(-)е[0, 1],

1 ' / Л2

уу(-)е[0, 1], г,к(-)е[0, 1], Е = - £ (У - У') .

I '=1

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

IX ( +1) > (1 -в),

Wj (t +1) =

Wj (t + 1), Wj (t + 1) <8,

11, V, (t + 1) > (1 -s),

v .. (t +1) = \ jV j ( ) |Vj (t + 1), V,.,. (t + 1) <s,

C,(m+1)k (t + 1) =

f0, cKm+i)k (t +1) <e,

I c

i(m+1)k v

Zik (t +1) =

0, Z,k (t + 1) <8,

z k (t + 1), z k (t + 1) >8,

где е - априорная малая пороговая константа.

В таблице 3 представлены функции принадлежности термов нечетких тенденций, являющихся входами системы нечеткого вывода F-Net после обучения алгоритмом обратного распространения

Функции принадлежности термов нечетких тенденций Membership functions of fuzzy trends terms

Таблица 3

Table 3

Таблица 4

Нормализованные значения содержания оксида кремния в шлаке на сутки

Table 4

The normalized values of forecasting a content of silicon oxide in slag per day

Час SiO2 в шлаке Прогноз Абсолютная ошибка

SiO2, нейроннаясеть SiO2, регрессия нейронной сети регрессии

1 0,534 0,520069927 0,607590648 0,014 -0,073

2 0,568 0,544541968 0,596700079 0,024 -0,029

3 0,557 0,576808962 0,594877295 -0,020 -0,038

4 0,534 0,565227411 0,601351714 -0,031 -0,067

5 0,523 0,541833426 0,609263965 -0,019 -0,087

6 0,500 0,530684357 0,583347127 -0,031 -0,083

7 0,500 0,507250527 0,573836196 -0,007 -0,074

8 0,500 0,507280949 0,581416924 -0,007 -0,081

9 0,489 0,507584816 0,590933649 -0,019 -0,102

10 0,489 0,496862291 0,592433513 -0,008 -0,104

11 0,489 0,497778245 0,600927239 -0,009 -0,112

12 0,489 0,498942001 0,609953821 -0,010 -0,121

13 0,443 0,499806522 0,618931646 -0,057 -0,176

14 0,500 0,454961897 0,636495236 0,045 -0,136

15 0,511 0,513263847 0,673975467 -0,002 -0,163

16 0,534 0,524270693 0,665766827 0,010 -0,132

17 0,545 0,547468183 0,657618939 -0,002 -0,112

18 0,545 0,557096557 0,649706965 -0,012 -0,104

19 0,534 0,556472453 0,641281792 -0,022 -0,107

20 0,534 0,544258208 0,626940119 -0,010 -0,093

21 0,523 0,543461532 0,6142034 -0,021 -0,091

22 0,523 0,533506036 0,663458543 -0,011 -0,141

23 0,511 0,532862532 0,658574478 -0,021 -0,147

24 0,511 0,520268262 0,655440789 -0,009 -0,144

ошибки. Пики термов «норма» функций принадлежности выходов сетей 71 и 72 показывают, что зоны без регулировок параметров Х1, Х3, Х4, Х6, Х7, Х8 и Х9 очень малы, управляющие воздействия на них при значениях «ниже нормы» более интенсивные, нежели при значениях «выше нормы».

Пики термов «выше среднего» для выхода 73 и «ниже нормы» для 74 указывают зоны тонкой регулировки параметров Х2 и Х5, поскольку на па-

раметры Х11 и Х12 управляющими воздействиями повлиять уже невозможно. Термы Х16 подтверждают практически линейную зависимость между текущим и предыдущим значениями содержания SiO2 в шлаке.

Результат обучения сети представлен на рисунке 5. (Для сравнения: на рисунке 3 приведен прогноз, полученный до проектирования нейро-нечеткой сети с помощью уравнения (6).) В таблице 4 показаны нормализованные значения про-

0 001 0 001 0 001 0 001 0 000 0 000 0 000 0 000

^ai^ai^wmwmwrsir^rsir^rsiiD^iD^iDOLnoLno^ai^ai^wmwmwrsir^rsir^rsiiD^HiD^HiDOLnoLno^rai^rai^rw ^^rsirsimm^^LnLniDiDr^r^wwaiaioo^^rsirsirsimm^^LnLniDiDr^r^wwaiaioo^^rsirsimm^r^r^rLnLniDiD

Рис. 5. Результаты прогнозирования содержания оксида кремния в шлаке авторской нейро-нечеткой сетью Fig. 5. The results offorecasting a content of silicon oxide in slag by the authors 'fuzzy-neural network

0 000

гноза оксида кремния в шлаке. Обучение сетей A-Net, B-Net, C-Net, D-Net составило 2 000 эпох. Обучение сети F-Net составило 5 000 эпох. Ошибка обучения всей гибридной сети - 0,026246, ошибка тестирования - 0,027246, ошибка на контрольной выборке составила 0,033756. Время обучения сети не превышает 4 минут.

Коэффициент взаимной корреляции между прогнозным (на один час вперед) и реальным значениями составил 0,953887263, среднеквадратиче-ское отклонение равно 0,041850444 (для уравнения регрессии (6) среднеквадратическое отклонение равно 0,081505196), среднеквадратическая невязка от нечеткого эталона равна 0,004435377 (что соответствует уровню высокопрофессиональных технологов), на два часа вперед - коэффициент взаимной корреляции равен 0,880619, среднеквадратическая невязка от нечеткого эталона равна 0,005782496.

С целью проверки обобщающих способностей разработанной адаптивной сети авторы имитировали отключение ряда датчиков на короткие (от 2 до 4 часов) интервалы, а также локально случайно искажали параметры Х3, Х7, X12. Сеть не выходит из устойчивого состояния, что указывает на хорошую степень обобщения разработанной сети. При этом коэффициент взаимной корреляции прогнозного и экспериментального содержания оксида кремния не опускался ниже значения 0,607939.

В рабочем режиме разработанная адаптивная нейро-нечеткая система прогнозирует содержание оксида кремния на 1, 2, 3 и 4 часа вперед, параллельно ее копия работает в режиме обучения. При снижении коэффициента взаимной корреляции между прогнозным и фактическими значениями ниже 0,85 управление передается переобученной копии, в свою очередь в оригинальную систему перегружаются адаптированные параметры настройки сети, а сама система переключается в режим фонового обучения.

В заключение отметим, что разработанная авторами нейро-нечеткая система способна обеспечить качество балансировки содержания оксида кремния в шлаках при переработке медно-никелевых руд в медные и никелевые катоды не ниже, чем это обеспечивается высокопрофессиональными сменными технологами. На программную реализацию данной системы получено свидетельство о государственной регистрации программ для ЭВМ № 2014611976 (автор Конев А.В.). Реализация SCADA-модуля, интегрированного с АСУТП верхнего уровня, на основе описанной модели позволит с точностью ~ 88 % прогнозировать процент содержания оксида кремния в шлаке на будущий период (не менее 4 часов), что позволит сократить процент потерь цветных металлов со шлаком. На Надеждинском металлургическом заводе комбината «Норильский никель» прово-

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

Литература

1. Производство металлов за Полярным кругом: технологическое пособие; [под ред. Н.Г. Кайтмазова]. Норильск: Антей лимитед, 2007. 296 с.

2. Дэвенпорт У.Г., Джоунс Д.М., Кинг М.Дж., Партелпо-ег Е.Г. Взвешенная плавка: контроль, анализ и оптимизация; [пер. с англ. под ред. Р.В. Старых]. М.: Изд-во МИСИС, 2006. 400 с.

3. Конев А.В., Фомичева С.Г. Разработка экспертной системы для управления взвешенной плавкой никелевого концентрата // Перспективы развития информационных технологий. Новосибирск: Изд-во ЦРНС, 2012. № 9. С. 85-88.

4. Фомичева С.Г., Филенко М.В., Братухин Р.И. Интеграция распределенных информационных систем горнодобывающей корпорации // Организатор производства. 2008. № 1. С. 76-83.

5. Спесивцев А.В. Металлургический процесс как объект изучения: новые концепции, системность, практика. СПб: Изд-во политех. ун-та, 2004. 306 с.

6. Жданов А.А. Автономный искусственный интеллект. М.: Бином. Лаборатория знаний, 2012. 359 с.

7. Юнусов Т.Р., Ярушкина Н.Г., Афанасьева Т.В. Моделирование трафика терминал-сервера на основе анализа нечетких тенденций временных рядов // Программные продукты и системы. 2007. № 4. С. 15-19.

8. Lin Y., Cunningham G.A. A new approach to fuzzy-neural system modeling: IEEE Transactions on Fuzzy Systems, 1995, vol. 3, no. 2, pp. 190-198.

9. Пегат А. Нечеткое моделирование и управление. М.: Бином. Лаборатория знаний, 2011. 798 с.

References

1. Kaytmazov N.G. (Ed.) Proizvodstvo metallov za Polyarnym krugom: tekhnologicheskoe posobie [Metal production behind the Polar circle: technological guide]. Norilsk, Antey Ltd. Publ., 2007, 296 p.

2. Davenport W.G., Jones D.M., King M.J., Partelpoeg E.H. Flash smelting: analysis, control and optimization. Warrendale, Pennsylvania, The Minerals, Metal, Materials Society, 2001, 321 p. (Russ. ed.: Vzveshennaya plavka: kontrol, analiz i optimizatsiya. Moscow, 2006, 400 p.).

3. Konev A.V., Fomicheva S.G. Developing expert system to control flash smelting of nickel concentration. Perspektivy razvitiya informatsionnykh tekhnologiy [IT prospects]. Novosibirsk, 2012, no. 9, pp. 85-88 (in Russ.).

4. Fomicheva S.G., Filenko M.V., Bratukhin R.I. Integration of distributed information systems in mining corporation. Organizator proizvodstva [Proguction organizer]. 2008, no. 1, pp. 76-83 (in Russ.).

5. Spesivtsev A.V. Metallurgicheskiy protsess kak obyekt izucheniya: novye kontseptsii, sistemnost, praktika [Metallurgical process as a study object: new concepts, systemacity, practice]. St. Petersburg, 2004, 306 p.

6. Zhdanov A.A. Avtonomny iskusstvenny intellekt [Autonomous artificial intelligence]. Moscow, 2012, 359 p.

7. Yunusov T.R., Yarushkina N.G., Afanasyeva T.V. Modeling terminal-server traffic based on timing series fuzzy trends analysis. Programmnye produkty i sistemy [Software & Systems]. 2007, no. 4. pp. 15-19.

8. Lin Y., Cunningham G.A. A new approach to fuzzy-neural system modeling: IEEE Transactions on Fuzzy Systems. 1995, vol. 3, no. 2, pp. 190-198.

9. Pegat A. Nechetkoe modelirovanie i upravlenie [Fuzzy modeling and control]. Moscow, Binom Publ., 2011, 798 p.