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

На правах рукописи

Кравцова Ольга Александровна

Алгоритмизация систем диагностики состояния элементов конструкций металлургических агрегатов

Специальность 05.13.06 - Автоматизация и управление технологическими процессами и производствами (промышленность)

Автореферат диссертации на соискание ученой степени кандидата технических наук

□03472947

Томск - 2009

003472947

Работа выполнена на кафедре автоматизации производственных процессов Кузбасской государственной педагогической академии (г. Новокузнецк)

Научный руководитель

Официальные оппоненты:

доктор технических наук профессор Кошелев Александр Евдокимович

доктор технических наук профессор Дмитриев Вячеслав Михайлович (Томский государственный университет систем управления и радиоэлектроники)

кандидат технических наук доцент Буинцев Владимир Николаевич (Сибирский государственный индустриальный университет, г. Новокузнецк)

Ведущая организация - Кемеровский государственный университет

Защита состоится « 18 » июня 2009 года в 15 час. 15 мин. на заседании диссертационного совета Д 212.268.02 при Томском государственном университете систем управления и радиоэлектроники (г. Томск, пр.Ленина, 40).

С диссертацией можно ознакомиться в библиотеке Томского государственного университета систем управления и радиоэлектроники (г. Томск, ул. Вершинина, 74).

Автореферат разослан « 18 » мая 2009 года.

Ученый секретарь диссертационного совета

Цц/

Мещеряков Р. В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ ВВЕДЕНИЕ

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

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

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

Объектом исследования. Высокотемпературные металлургические агрегаты и элементы их конструкции.

Предметом исследования являются системы диагностики технического состояния элементов конструкций металлургических агрегатов.

Цели и задачи диссертации. Основной целью работы является разработка методики и алгоритмов диагностирования технического состояния элементов конструкций металлургических агрегатов с применением методов обнаружения локальных особенностей контролируемых сигналов. В соответствии с основной целью решаются следующие задачи: ■ обобщение и дополнение теоретических и практических разработок по технической диагностике элементов конструкций металлургических агрегатов;

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

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

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

■ конкретизация алгоритмов оценивания технического состояния элементов конструкций металлургических агрегатов.

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

Научная новизна

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

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

■ Алгоритм рекуррентного определения вейвлет-свертки, уменьшающий вычислительные затраты на 19-23%, по сравнению с классическим алгоритмом расчета вейвлет-свертки.

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

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

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

Практическая значимость

■ Предложенные методика и алгоритмы позволяют:

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

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

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

- повысить информативность краевых зон вейвлет-спектрограмм и обеспечить их построение в режиме реального времени;

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

■ Разработанные методика и алгоритмы применимы к любым диагностируемым технологическим объектам.

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

Предмет защиты и личный вклад автора. На защиту выносятся:

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

■ алгоритмы повышения информативности краевых зон вейвлет-спектрограмм и рекуррентного расчета вейвлет-свертки;

■ структура распараллеленного алгоритма вейвлет-анализа анализируемых сигналов с использованием рекуррентного определения вейвлет-свертки;

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

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

Реализация результатов работы. Разработанные алгоритмы обнаружения дефектов технического состояния элементов конструкций металлургических агрегатов непрерывного и циклического действия, включающие алгоритм определения возможного будущего технического состояния элементов конструкций объектов управления, приняты к использованию на ОАО "Новокузнецкий металлургический комбинат" (ОАО НКМК). Разработанные методы и алгоритмы используются в учебном процессе Сибирского государственного индустриального университета (г. Новокузнецк) и Кузбасской государственной педагогической академии (г. Новокузнецк).

Апробация работы. Основные положения диссертации освещались на 8 конференциях: 5-я Всероссийская научно-практическая конференция "Системы автоматизации в образовании, науке и производстве" (Новокузнецк, 2005); 2-я

Всероссийская научно-практическая конференция "Моделирование, программное обеспечение и наукоемкие технологии в металлургии" (Новокузнецк, 2006); 3-я Международная научно-практическая конференция "Технолого-экономическое образование в XXI веке" (Новокузнецк, 2004, 2005, 2006); Всероссийская научно-практическая конференция "Металлургия: новые технологии, управление, инновации и качество" (Новокузнецк, 2006, 2008); 4-я Международная научно-практическая конференция "Электронные средства и системы управления: опыт инновационного развития " (Томск, 2007).

Публикации. По теме диссертации опубликовано 16 работ, из них 2 статьи в изданиях, рекомендованных ВАК, 2 патента РФ на изобретение и 14 работ - материалы научно-технических и научно-практических конференций.

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения, списка литературы и содержит 179 страниц основного текста, включая 59 рисунков и 9 таблиц. Библиографический указатель содержит 134 наименования.

СОДЕРЖАНИЕ РАБОТЫ

Глава 1. Основы построения систем технической диагностики элементов конструкций металлургических агрегатов

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

В настоящее время развитие как отечественных, так и зарубежных систем диагностики элементов конструкций металлургических агрегатов осуществляется по нескольким направлениям: 1) использование математических моделей нестационарного теплообмена (ОАО "Черметавтоматика", ОАО "Новолипецкий металлургический комбинат", ОАО "Магнитогорский металлургический комбинат", компания "Даниэли Корус", завод фирмы Stahlwerke Bremen GmbH); 2) ультразвуковой контроль диагностируемых элементов конструкции (Днепропетровский металлургический институт).

Большой вклад в совершенствование, например, систем диагностики состояния футеровки горна доменной печи внесли авторы: Гаврилов Е. Е., Спирин Н. А., Серов Ю. В., Кудинов Г. А, М. Шульте, Васильев П. Г, и Голосков Б. В. Суть большинства их разработок заключается в оперативном определении остаточной толщины футеровки горна доменной печи с опорой на решения прямой и обратной задач нестационарного теплообмена. Отличительной особенностью некоторых систем диагностики состояния футеровки горна доменной печи является наличие характерных вариантов ее разгара, с одним из которых отождествляется текущее состояние огнеупорной кладки горна.

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

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

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

Рис. 1. Обобщенная структура системы технической диагностики элементов конструкций металлургических агрегатов

Впервые задача обнаружения изменения свойств временных последовательностей поставлена Е. Пейджом. Колмогоровым А. Н. и Ширяевым А. Н. предложен метод определения моментов изменения свойств в рассматриваемых временных рядах данных при заданном уровне ложных тревог, применение которого на практике требует решения сложных нелинейных дифференциальных уравнений. Телькснис Л. А. предложил метод обнаружения наиболее вероятных изменений свойств коррелированных случайных процессов, описываемых стохастическими дифференциальными

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

Глава 2. Совершенствование алгоритмов систем диагностики элементов технологических объектов

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

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

Проведенный сравнительный анализ эффективности рассматриваемых методов показал следующее:

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

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

- вейвлет-анализ характеризуется стабильностью выявления моментов изменения структурных свойств рассматриваемых временных последовательностей даже при малых их изменениях;

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

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

- вейвлет-анализ выявляет без ошибки по времени все моменты изменений локальных особенностей помехоискаженных временных

последовательностей данных.

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

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

О 0.063 0,25 1.53 Дисперсия

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

Повышение информативности краевых зон вейвлет-спектрограмм. Решение данной задачи направлено на устранение ограниченных возможностей применения математического аппарата вейвлет-анализа в системах оперативной диагностики.

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

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

■ выбор типовой функции в зависимости от специфики объекта управления и контролируемой величины;

■ формирование объема экстраполируемых данных: х"(1 + \),х"(1 + 2),...,х"(1 + ма), где г"(/)-данные модельной экстраполируемой последовательности; Лм- количество весовых коэффициентов в вейвлет-базисе на максимальном (А-ом) уровне детализации сигнала (в случае неверно сформированного объема данных модельного сигнала имеющиеся

изменения, находящиеся в краевых зонах, остаются недоступными для анализа);

■ корректировка начального значения модельной реализации, которая предполагает анализ последнего фактического значения исследуемой реализации и изменение начальной координаты модельного сигнала относительно последнего исходного данного: Х"(1 + 1) = Д-Х"(/), где Xм(1 +1)- первое данное модельной последовательности; Х"(1) - последнее данное исходной временной последовательности; Д-корректирующий коэффициент.

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

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

■ экстраполяционный подход: например, при использовании аппроксимирующего полинома второй степени: х""'"(/) = а,1 + агР, число математических операций, затраченных на экстраполяционное продолжение сигнала для проведения вейвлет-анализа, равно х = 4ЛГА ■ При применении

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

• предлагаемый подход: число математических операций, затраченных на генерацию экстраполируемой последовательности для дальнейшего проведения вейвлет-анализа, равно х = ^А +1, что значительно меньше, чем

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

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

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

Алгоритм рекуррентного определения вейвлет-свертки. Традиционно для вычисления значений вейвлет-свертки используется следующее математическое выражение:

(I -Ь\

Ь)=\а\~у2 -где - значения вейвлет-свертки

(коэффициенты вейвлет-преобразования); где / = 1,2..Х - исследуемая реализация; па = 1 ...А'а - число весовых коэффициентов вей влет-базиса на выбранном уровне детализации сигнала. В работе предлагается в системах диагностики элементов конструкций объектов управления использовать для расчета значений вейвлет-свертки разработанную в ходе выполнения исследований рекуррентную форму: 1/2

\У(1,а,Ь) = Ш(1-\,а,Ь) + \а~' где п = \...Ыа - количество весовых

коэффициентов вейвлет-базиса на выбранном уровне детализации сигнала; а - параметр масштабирования; Ъ - параметр смещения по оси времени.

Блок-схема предлагаемого варианта рекуррентного расчета значений вейвлет-свертки представлена на рис. 3. Оценка работоспособности предлагаемого алгоритма расчета значений вейвлет-свертки проводилась на модельных и натурных данных.

с:

Выбор начальных условий

Расчет на а-ом уровне детализации сигнала начального значения

вейвлет-свертки: (1,а, Ь) =\ а\']12

Расчет приращения на текущий момент времени:

Дх,

Расчет значений вейвлет-свертки в последующие моменты времени: ¡У(1,а,Ь) = №(1-\,а,Ь) + У~ ^6прДг(р),гдеп = 1...№-число

р=1-\а*\

весовых коэффициентов вейвлет-базиса на выбранном уровне детализации сигнала; а — параметр масштабирования; Ь — параметр смещения по оси времени.

Сохранение значения вейвлет - свертки IV^ (/, ¿7, Ь) и текущих

приращений АХ/.

Проверка условия: 1>Ь

Смена уровня детализации сигнала: Д/=а,+1, где а, - текущий уровень детализации сигнала;

а¡- новый уровень детализации сигнала. Определение коэффициентов вейвлет-базиса на а-ом уровне детализации сигнала.

Проверка условия: а,>А, где А - максимально заданный уровень детализации сигнала.

Представление результатов вейвлет - анализа:

1) в виде вейвлет - спектрограммы;

2) в виде временных последовательностей для каждого 6/ уровня.

Рис. 3. Блок-схема алгоритма вейвлет-анализа с применением рекуррентной формы расчета вейвлет-свертки

Диаграмма, характеризующая временные затраты на реализацию разработанного алгоритма расчета значений вейвлет-свертки и алгоритма, базирующегося на традиционном варианте, представлена на рис. 4. Анализ временных затрат проводился на модельных последовательностях данных в 50, 100, 200, 500 отсчетов для 100 уровней детализации сигнала, использовалась машина РепйипЯУ 2.4ГНг/768МЬ-0011.

Рис. 4. Временные затраты на определение значений вейвлет-сверток при применении разных методов их расчета

Применение разработанного алгоритма рекуррентного оценивания вейвлет-свертки позволяет снизить временные затраты на получение значений вейвлет-свертки на 19-23%.

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

В процессе функционирования предлагаемой структуры на каждом уровне детализации (в каждом канале) выполняются следующие функции:

• формирование начальных условий;

• расчет весовых коэффициентов вейвлет-базиса;

• определение вейвлет-свертки на всех заданных уровнях детализации сигнала;

• выявление структурных локальных изменений измерительных сигналов.

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

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

х(1) - фактические значения анализируемого сигнала в / - текущий момент времени; у/ - тип вейвлет-базиса; 1У(1,а,Ь) - значения вейвлет-свертки; апа2,...,ам — уровни детализации сигнала; М - заданное число уровней детализации сигнала; Ь - параметр смещения по оси времени.

Рис. 5. Многоканальная структура распараллеленного алгоритма вейвлет-анализа измерительных сигналов

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

Последовательная форма построения вейвлет-спектрограммы

|ООООЕН> РОП-ООО^ [СНООООО'

Параллельная форма построения

вейвлет-спектрограммы <1

-ООО-ООО- , Ч>СНОООО I -СНОП-ООО | ООО-ООО I

□ Последовательная форма построения вейвлет-спектрограммы аа построения вейвлет-спектрограммы

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

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

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

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

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

Для проверки работоспособности предлагаемого алгоритма получения возможных оценок будущего состояния диагностируемого объекта были использованы натурные данные о температуре футеровки на горизонтах № 1.2, №3 и №5 горна доменной печи №5 ОАО НКМК за сентябрь 2004 года. Некоторые из реализаций натурных данных представлены на рис. 7.

В состав рассматриваемого алгоритма входят следующие функциональные блоки:

■ накопления начального объема анализируемых данных: х0, ____х„ (п -

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

■ формирования эффективного объема выборки (детально рассмотрено в задаче "Повышение информативности краевых зон вейвлет-спектрограмм");

а - горизонт №3 сечение №2, б - горизонт №3 сечение №4; 1-8 - номера структурно однородных участков анализируемых реализаций. Рис. 7. Реализации фактической температуры (Т°С) футеровки горна доменной печи от различных термодатчиков

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

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

с 1-1/"*

5 =-, где п - число данных, оставшихся в усеченном ряде;

1-1/л

п - число данных в исходной реализации.

■ полиномиальной аппроксимации и экстраполяции последовательности, отображающей динамику показателей сложности (рис. 8);

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

I %85 _

5 *

Р §0.75 -я и

" 0.7 -

0 2 4 6 8

Рис. 8. Фактические - модельные — 8М(/) и экстраполированное значения— 8(/+1) показателей сложности для анализируемого ряда данных (Т°С) (горизонт №3, сечение №4)

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

|(/+|)»0.775

5„(/) = 0.0049/ - 0.0549/ + 0.9204 _¿.

Номер участка

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

а - участок №7 (горизонт №3 сечение №4); б - участок №8 (горизонт №3 сечение №4); Т'(1), т"(1) - варианты прогнозируемых рядов данных с одинаковыми показателями сложности.

Рис. 9. Фактические и вариантные экстраполированные траектории динамики температуры футеровки горна доменной печи

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

Табл. 1. Значения среднемодульного отклонения (СМО).

Горизонт, сечение № участка СМО при экстраполировании с использованием

полиномиальной экстраполяции по одному смежному структурно-однородному участку полиномиальной экстраполяции по всей реализации показателя сложности "S"

Гориз. 3 сеч. 2 7 1.21 10.25 0.30483

8 1.42 21.2 1.63618

Гориз. 3 сеч. 4 7 2.28 9.2 0.38744

8 5.46 28.8 0.50817

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

Глава 3. Совершенствование систем диагностики элементов конструкций металлургических агрегатов

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

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

Модернизация системы диагностики состояния футеровки горна доменной печи

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

В качестве прототипа выбрана система контроля разгара горна доменной печи №5 "Футеровка" ОАО "Новокузнецкий металлургический комбинат" (НКМК). В данной системе измерение температур в футеровке производят при помощи датчиков температуры, месторасположение которых определяется номерами поясов (и), радиусов (г), на которые условно разбита футеровка металлоприемника печи и номером (т) установки термодатчика по направлению г-го радиуса. На основании показаний группы датчиков решается обратная задача теплопроводности, однако, точное ее решение для условий горна доменной печи невозможно из-за неопределенности и существенных упрощений относительно ряда технологических характеристик. Далее на экран монитора выводится топология распределения тепловых полей по всему радиусу горна на текущий момент времени. Однако для правильного оценивания развивающейся ситуации технологическому персоналу необходимо наблюдать топологию разгара горна как в предшествующие, так и в будущие моменты времени.

Отличительной особенностью предлагаемого алгоритма является то, что информация с каждого функционирующего термодатчика подвергается проверке на устойчивость считывания и на достоверность, из достоверных данных формируются эффективные временные последовательности, которые затем подвергают структурному анализу с целью выявления локальных изменений с последующим запоминанием их координат. Далее проводится анализ приращений температуры за текущие сутки по /и, датчику Дс (/) = *(/)-*(/-23), величина которого не должна превышать порог в 2°С, и за

7

текущую неделю по т,• датчику: ДЛ(/) = ^АС(/), максимальное изменение

которого не должно превышать 15°С (пороговые величины в 2°С и 15°С приняты из расчета максимально допустимого изменения температуры футеровки горна доменной печи на 50°С за месяц), по результатам которого осуществляется представление информации технологическому персоналу о

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

Рис. 10. Отображение изменения температурного поля футеровки горна

доменной печи

Предлагаемый подход позволяет обнаружить локальные изменения рабочего пространства (разгар футеровки, "зарастание" горна) металлоприемника доменной печи на ранней стадии и предоставлять информацию об изменениях в температурном режиме горна в динамике. Например, по данным о температуре футеровки горна доменной печи №5 ОАО НКМК за август 2006 года, начинающийся интенсивный разгар огнеупорной кладки с использованием предлагаемого подхода обнаружен на 710 дней раньше, чем фактически зафиксированный разгар.

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

Структура модернизированной системы технической диагностики состояния продувочных устройств на АКОС представлена на рис. 11.

Рабочая станция: Celeron CPU 3.2 GH. 2 Гб ОЗУ

Рис. 11. Структура предлагаемой системы технической диагностики состояния продувочного устройства на агрегате комплексной обработки стали

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

Для апробации разработанного алгоритма использовались натурные данные ОАО "Новокузнецкий металлургический комбинат". Результаты диагностирования продувочного устройства представлены на рис. 12. Так, например, в ситуации, когда наблюдается рост значений вейвлет-свертки анализируемых сигналов о расходе газа и положении фурмы, а значения вейвлет-свертки сигнала о давлении газа изменяются незначительно и не превышают установленный диапазон, можно говорить об образовании в продувочной фурме свища.

Рис. 12. Реализации фактических значений контролируемых величин процесса выплавки стали с возможным дефектом фурмы.

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

Заключения и выводы

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

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

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

■ Выявлены предпосылки, достоинства и недостатки классического математического аппарата вейвлет-анализа.

■ Предложен метод коррекции процедуры вейвлет-анализа устраняющий неопределенность концевых зон вейвлет-спектрограмм с использованием типовых функций, сокращающий объем вычислений при построении вейвлет - спектрограмм в 2-4 раза, а так же предложена рекуррентная форма вычисления вейвлет-свертки, уменьшающая временные затраты на 19-23% по сравнению с классической формой расчета вейвлет-свертки.

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

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

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

Труды по теме диссертации

Научные журналы и издания, рекомендованные ВАК.

1. Кошелев, А.Е. Подсистема диагностики состояния элементов конструкций металлургических агрегатов в составе АСУ ТП / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова // Доклады Томского государственного университета систем управления и радиоэлектроники. - 2008. - 2(18), ч. 2 -С. 110-116.

2. Кошелев, А.Е. Структурный анализ измерительных сигналов в доменном и сталеплавильном производствах / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова // Известия вузов. Черная металлургия. — 2006. - №12. - С. 43-47.

Материалы международных и всероссийских научно-технических и научно-

практических конференций.

3. Криволапова, Л.И. О применении вейвлет - анализа в автоматизированной системе контроля состояния футеровки горна доменной печи / Л.И. Криволапова, O.A. Кравцова // Системы автоматизации в образовании, науке и производстве: Материалы 4-ой Всероссийской научно-практической конференции. - Новокузнецк: СибГИУ, 2003. - С 377-379.

4. Кошелев, А.Е. Анализ остаточной толщины футеровки горна доменной печи / А.Е. Кошелев, O.A. Кравцова // Технолого-экономическое образование в XXI веке: Материалы Всероссийской конференции -

Новокузнецк: КузГПА, 2004. - С 115-119.

5. Кошелев, А.Е. Применение вейвлет-анализа для диагностики состояния продувочной фурмы при внепечной обработке стали / А.Е. Кошелев, O.A. Кравцова // Технолого-экономическое образование в XXI веке: Материалы Всероссийской конференции - Новокузнецк: КузГПА, 2004. - С 90-95.

6. Кошелев, А.Е. Прогнозирование нестационарных рядов данных с применением показателя их сложности / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова // Технолого-экономическое образование в XXI веке: Материалы 2-й Международной научно-практической конференции в 4-х томах. Том 4. - Новокузнецк: Изд-во КузГПА, 2005. - С. 77-83.

7. Кошелев, А.Е. Применение вейвлет-анализа в задаче диагностики металлургических агрегатов / А.Е. Кошелев, O.A. Кравцова // Системы автоматизации в образовании, науке и производстве: Материалы 5-й Всероссийской научно-практической конференции - Новокузнецк: СибГИУ, 2005.-С 212-213.

8. Криволапова, Л.И. Применение показателя сложности измерительных сигналов в системах управления технологическими процессами / Л.И. Криволапова, O.A. Кравцова // Моделирование, программное обеспечение и наукоемкие технологии в металлургии: Материалы 2-й Всероссийской научно-практической конференции - Новокузнецк: СибГИУ, 2006. - С 370-374.

9. Кошелев, А. Е. Построение вейвлет - спектрограммы с использованием параллельных вычислений / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова // Технолого-экономическое образование в XXI веке: Материалы 3-й Международной научно-практической конференции в 4-х томах. Том 2. - Новокузнецк: Изд-во КузГПА, 2006. - С. 74-78.

10. Криволапова, Л.И. Обучающая система по управлению тепловым состоянием доменной печи / Л.И. Криволапова, М.В. Ляховец, O.A. Кравцова, A.B. Овечкин // Технолого-экономическое образование в XXI веке: Материалы 3-й Международной научно-практической конференции в 4-х томах. Том 2. - Новокузнецк: Изд-во КузГПА, 2006. -С. 91-95.

11. Криволапова, Л.И. Модернизация автоматизированной системы диагностики состояния футеровки горна доменной печи / Л.И. Криволапова, O.A. Кравцова // Металлургия: новые технологии, управление инновации и качество: Материалы Всероссийской научно-практической конференции. - Новокузнецк: Изд-во СибГИУ, 2006. - С. 214-218.

12. Кошелев, А. Е. Обобщенное представление технической диагностики состояния элементов конструкций металлургических агрегатов / А. Е. Кошелев, Л. И. Криволапова, О. А. Кравцова // Электронные средства и системы управления: опыт инновационного развития. Материалы 4-й международной научно-практической конференции в 2-х частях. Часть 2. -Томск: Изд-во В-Спектр, 2007. - С. 50-53.

13. Кошелев, А.Е. Использование оценки сложности измерительных сигналов в системах диагностики состояния объектов управления / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова // Электронные средства и системы управления: опыт инновационного развития. Материалы 4-й международной научно-практической конференции в 2-х частях. Часть 2. -Томск: Изд-во В-Спектр, 2007. - С. 53-56.

14. Кошелев, А.Е. Модернизация АСУ внепечной обработки стали / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова // Металлургия: новые технологии, управление инновации и качество. Материалы Всероссийской научно-практической конференции. - Новокузнецк: Изд-во СибГИУ, 2008. -С. 284-289.

Патенты на изобретение Российской Федерации.

15. Патент РФ №2299910 С2. МПК С 21 В 7/24. Способ контроля состояния футеровки горна доменной печи / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова. - Опубл. 27.05.2007, Бюл. изобр. №15. - С. 554-555.

16. Патент РФ №2299914 С2. МПК С 21 С 7/072, С 21 С 5/30 Способ контроля состояния фурмы при продувке расплава в ковше / А.Е. Кошелев, Л.И. Криволапова, O.A. Кравцова. - Опубл. 27.05.2007, Бюл. изобр. №15. -С. 557.

Кравцова Ольга Александровна

Алгоритмизация систем диагностики состояния элементов конструкций металлургических агрегатов

Автореферат диссертации на соискание ученой степени кандидата технических наук

Подписано в печать 15.05.2009 г. Формат бумаги 60*84 1/16. Бумага писчая. Печать офетная. Усл. печ. Л 1,0. Уч.-изд.л. 1,1. Тираж 100 экз. Заказ № 346.

ГОУ ВПО "Сибирский государственный индустриальный университет" 654007, г. Новокузнецк, ул. Кирова, 42 Издательский центр ГОУ ВПО "СибГИУ"

Введение.

Глава 1. Основы построения систем технической диагностики элементов конструкций металлургических агрегатов.

1.1 Обобщенное представление технической диагностики.

1.2 Аналитический обзор действующих систем диагностики элементов конструкций технологических объектов.

1.2.1 Системы контроля состояния футеровки горна доменной печи.

1.2.2 Система контроля состояния процесса продувки при внепечной обработке стали.

1.2.3 Методы обнаружения изменений свойств анализируемых процессов.

1.3 Основы вейвлет-преобразования.

Выводы по первой главе.

Глава 2. Совершенствование алгоритмов систем диагностики элементов технологических объектов.

2.1 Сопоставление методов структурного анализа временных рядов данных

2.2 -Задача повышения эффективности вейвлет-преобразования.

2.3 Применение показателей "сложностного" подхода в задаче оценивания возможного будущего состояния диагностируемого объекта.

2.3.1 Понятие сложности в системах управления.

2.3.2 Сравнительный анализ различных показателей сложности.

2.3.3 Задача прогназирования будущего состояния диагностируемого объекта.

Выводы по второй главе.

Глава 3. Совершенствование систем диагностики элементов технологических объектов.

3.1 Модернизация системы диагностики состояния футеровки горна доменной печи.

3.1.1 Описание действующей системы контроля разгара горна ОАО НКМК.

3.1.2 Модернизированная система диагностики состояния футеровки горна доменной печи.

3.2 Модернизация системы диагностики состояния продувочной фурмы на установке продувки стали газом.

3.2.1 Описание действующей системы диагностики состояния продувочной фурмы АКОС.

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

Выводы по третьей главе.

Заключения и выводы.

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

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

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

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

Предметом исследования являются системы диагностики технического состояния элементов конструкций металлургических агрегатов.

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

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

Научная новизна

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

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

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

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

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

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

Практическая значимость

Предложенные методика и алгоритмы позволяют:

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

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

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

- повысить информативность краевых зон вейвлет-спектрограмм и обеспечить их построение в режиме реального времени;

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

Разработанные методика и алгоритмы применимы к любым диагностируемым технологическим объектам.

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

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

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

Реализация результатов работы. Разработанные алгоритмы обнаружения дефектов технического состояния элементов конструкций металлургических агрегатов непрерывного и циклического действия, включающие алгоритм определения возможного будущего технического состояния элементов конструкций объектов управления, приняты к использованию на ОАО "Новокузнецкий металлургический комбинат" (ОАО НКМК). Разработанные методы и алгоритмы используются в учебном процессе Сибирского государственного индустриального университета (г. Новокузнецк) и Кузбасской государственной педагогической академии (г. Новокузнецк).

Апробация работы. Основные положения диссертации освещались на 8 конференциях: 5-я Всероссийская научно-практическая конференция "Системы автоматизации в образовании, науке и производстве" (Новокузнецк, 2005); 2-я Всероссийская научно-практическая конференция "Моделирование, программное обеспечение и наукоемкие технологии в металлургии" (Новокузнецк, 2006); 3-я Международная научно-практическая конференция "Технолого-экономическое образование в XXI веке" (Новокузнецк, 2004, 2005, 2006); Всероссийская научно-практическая конференция "Металлургия: новые технологии, управление, инновации и качество" (Новокузнецк, 2006, 2008); 4-я Международная научно-практическая конференция "Электронные средства и системы управления: опыт инновационного развития " (Томск, 2007).

Публикации. По теме диссертации опубликовано 16 работ, из них 2 статьи в изданиях, рекомендованных ВАК, 2 патента РФ на изобретение и 14 работ - материалы научно-технических и научно-практических конференций.

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения, списка литературы и содержит 179 страниц основного текста, включая 59 рисунков и 9 таблиц. Библиографический указатель содержит 134 наименования.

ЗАКЛЮЧЕНИЯ И ВЫВОДЫ

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

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

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

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

Предложен метод коррекции процедуры вейвлет-анализа, устраняющий неопределенность концевых зон вейвлет-спектрограмм с использованием типовых функций, сокращающий объем вычислений при построении вейвлет - спектрограмм в 2-4 раза, а так же предложена рекуррентная форма вычисления вейвлет-свертки, уменьшающая временные затраты на 19-23% по сравнению с классической формой расчета вейвлет-свертки.

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

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

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

1. Авдеев В.П. Динамические показатели распределенности и нестабильности процессов / В.П. Авдеев, JI.A. Сульман, Б.И. Ашпин и др. // Известия ВУЗов. Черная металлургия—1980. — №8. С. 96-100.

2. Авторское свидетельство СССР 1305179 А1. МКИ С 21 С 7/064. Способ внепечной обработки стали / В.М. Тарасов, C.B. Лепорский, И.М. Рожков и др. Опубл. 23.08.87, Бюл. изобр. №15.

3. Авторское свидетельство СССР 1397487 А1. МКИ С 21 В 7/24, G 01 F 23/28. Способ измерения толщины футеровки доменной печи / В.Е. Ленченко, П.Г. Васильев, Н.Е. Алпаев и др. Опубл. 23.05.88, Бюл. изобр. №19.

4. Алексеев Д.И. К вопросу закупоривания продувочных фурм для внепечной обработки / Д.И. Алексеев, C.B. Казаков, А.Г. Свяжин // Известия ВУЗов. Черная металлургия.- 1993. — №3. — С. 31—35.

5. Аснис А. Е. Повышение прочности сварных конструкций / А.Е. Аснис, Г.А. Иващенко. Киев: Наукова думка, 1978. - 193 с.

6. Астафьева Н.М. Вейвлет-анализ: основы теории и примеры применения //Успехи физический наук.-1996.- Т.166. №11. С. 1145-1170.

7. Байшев Ю.П. Доменные печи и воздухонагреватели (конструкции, эксплуатуционные воздействия, свойства материалов, расчеты) / Ю.П. Байшев. Екатеринбург, 1996. - 987 с.

8. Баххофен Х.Ю. Условия течения и износ в горне доменной печи / Х.Ю. Баххофен, Х.П. Рютер, М. Петере и др. // Черные металлы.-2000. -№ 2.-С. 24-31.

9. Бестужев Н.И. Внепечная обработка сплавов / Н.И. Бестужев // Литейное производство. 1999. № 3. - С. 17 - 18.

10. Бородкин Л.И. Алгоритм обнаружения моментов изменения параметров уравнения случайного процесса / Л.И. Бородкин, В.В. Моттль // Автоматика и телемеханика.-1976. № 6. - С. 23—32.

11. Бурков В.Н. Локальный структурный анализ- реализаций данных / В.Н. Бурков, Т.В. Киселева // Автоматика и телемеханика.-2004. №12. -С. 153-158. ^

12. Васильев П.Г. Ультразвуковой контроль износа футеровки шахтыдоменной печи / П.Г. Васильев, В.Е. Левченко, Н.Е. Алапаев,

13. A.B. Шульга // Металлургическая и горная промышленность.-1992. -№3.-С. 3-5. I

14. Вихлевщук В.А. Ковшовая доводка стали / В.А. Вихлевщук,

15. B.C. Харахулах, С.С.Бродский. — Днепропетровск: Системные ^ технологии, 2000. 190 с.

16. Вихлевщук В.В. Установки доводки стали в ковшах большой емкости / В.В. Вихлевщук, A.C. Стороженко, C.B. Лепорский, В.А. Синельников, В.А. Сахно // Сталь.-1989. №11. - С. 34-37.

17. Вихлевщук В.А. Создание современных процессов внепечной обработки стали для ресурсосбережения и освоения новых видов продукции / В.А. Вихлевщук // Металлургическая и горнорудная промышленность.-1999.-№ 5. С. 33 - 35.

18. Вишняков А.Н. Адаптивный метод обнаружения нарушений закономерностей по наблюдаемым данным / А.Н. Вишняков, Я.З. Цыпкин // Автоматика и телемеханика.-1991. — №6. — С. 127—135.

19. Вишняков А.Н. Адаптивный метод обнаружения нарушений закономерностей по наблюдаемым данным при наличии помех / А.Н. Вишняков, Я.З. Цыпкин // Автоматика и телемеханика.-1991. — №12. -С. 128-137.

20. Володин А.Ф. Продувка стали азотом в ковше / А.Ф. Володин, Н.М. Блащук, В.Н. Мачикин и др. // Черная металлургия: Бюл. ин-та "Черметинформация".-1985. № 23. С. 41.

21. Воробьев В.И. Теория и практика вейвлет-преобразования / В.И. Воробьев, В.Г. Грибунин СПб: Военный университет связи.-1999. -204с.: ил.

22. Гаврилов Е.Е. Диагностирование состояния футеровки металлоприемника доменной печи по тепловым нагрузкам / Е.Е. Гаврилов, А.Г. Ульянов, В.В. Канаев и др. // Сталь.-1987. №1. - С 13-15.

23. Галемин И.М. Углеродистая футеровка доменных печей / И.М. Галемин,

24. A.B. Горох -М: Металлургия.-1964. 100 е.: ил.

25. Гер Т. Наблюдение за агрегатами в черной металлургии / Т. Гер, Г. Люльф // Черные металлы.-2001. № 6. — С. 36-40.

26. Гладышев В.В. Количественные критерии влияния цинка на стойкость футеровки горна и работу доменной печи / В.В. Гладышев,

27. B.В. Филипов, B.C. Рудин // Сталь. 2001. - №1. - С. 6-10.

28. Голосков Б.В. Новые методы исследования доменного процесса: Труды научно-технического общества черная металлургия / Б.В. Голосков. -М.: Металлургиздат, 1958. Т. 17. - С. 566-578.

29. ГОСТ 20911-89 "Техническая диагностика. Термины и определения". -М.: Изд-во стандартов, 1990.- 13 с.

30. ГОСТ Р 22.0.05-94 "Безопасность в чрезвычайных ситуациях. Технологические чрезвычайные ситуации. Термины и определения". -М.: Изд-во стандартов, 1995.

31. Гребенюк Е.А. Обнаружение изменений свойств нестационарных случайных процессов / Е.А. Гребенюк // Автоматика и телемеханика.-2003.-№12.-С. 44-59.

32. Дарховский Б.С. Апостериорное обнаружение момента "разладки" случайной последовательности / Б.С. Дарховский, Б.Е. Бродский // Теория вероятностей и ее применения.-1980 — Т. 25. Вып. 3. — с. 476 -489.

33. Дарховский Б.С. Непараметрический метод для апостериорного обнаружения момента "разладки" последовательности независимых случайных величия / Б.С. Дарховский // Теория вероятностей и ее применения.-1976. Т. 21. - Вып. 1.-е. 180-134.

34. Дарховский Б.С. О подходе к оценке сложности кривых / Б.С. Дарховский, А.Я. Каплан, СЛ. Шишкин // Автоматика и телемеханика-2002. №3. - с. 134-140.

35. Даценко Н.М. Вейвлетный анализ временных рядов и динамика атмосферы / Н.М. Даценко, Д.М. Сонечкин // Известия ВУЗов. Прикладная нелинейная динамика-1993. Т.1, - № 1,2.-е 9-14.

36. Демидов Е.Е. Нелинейный корреляционный анализ / Е.Е. Демидов, Ю.В. Даревская, O.A. Моренков. и др. // Обозрение прикладной и промышленной математики. — М.: ТВП, 1999. — Т.6. — Вып. 1.-е. 4-57.

37. Дробышевский О.В. Выбор оптимальной конфигурации сопла фурмы для продувки металла в ковше / О.В. Дробышевский, A.B. Явойский, В.Г. Востриков, А.Н. Резяпкин // Известия ВУЗов. Черная металлургия.-1994. №1. - С. 23-25.

38. Дубовик В.Я. Об автоматическом контроле состояния горна доменной печи / В.Я. Дубовик, A.A. Плышевский, Н.И. Валов // Известия ВУЗов. Черная металлургия.-1982. №4. — С. 4-6.

39. Дьяконов В.П. Вейвлеты. От теории к практике / В.П. Дьяконов. М.: СОЛОН-Р, 2002. - 448 е.: ил.

40. Жуков М.В. Выбор методов и средств диагностирования и восстановления с учетом экономической целесообразности /М.В. Жуков, В.В. Карибский // Приборы и системы. Управление, контроль, диагностики.-2001. №3. - С. 58-62.

41. Ивахненко А.Г. Предсказание случайных процессов / Г.А. Ивахненко, В.Г. Лапа. Киев: Наукова думка, 1971. - 416с.: ил.

42. Инструкция по футеровке и сушке доменных печей и их вспомогательных устройств. Министерство черной металлургии СССР.-М., 1985.- 127 с.

43. Иодковский С.А. Состояние и перспективы развития внепечной обработки стали / С.А. Иодковский // Труды четвёртого конгресса сталеплавильщиков. М., 1997. - С. 237 - 243.

44. Иодковский С.А. Состояние и перспективы развития внепечной обработки стали / С.А. Иодковский // Черная металлургия.-1999. № 9 -10.-С. 21 -25.

45. Камерон J1.B. Разработка новой системы перемешивания металла аргоном / J1.B. Камерон, С.Б. Макин, М.А. Роберте // Труды четвёртого конгресса сталеплавильщиков. М., 1997. - С. 247-251.

46. Каминская В.А. Обнаружение изменений параметров процесса авторегрессии / В.А. Каминская, Д.А. Шипените // Труды АН ЛитССР.-1975. сер Б, т. 4(89). - С. 143-147.

47. Касти Д. Большие системы: Связность, сложность и катастрофы / Д. Касти; пер. с англ.; под ред. Ю.П. Гупало, A.A. Пионтковского. М.: Мир, 1982.-216 е.: ил.

48. Клигене Н.И. Оценка момента изменения параметров распределения случайных последовательностей / Н.И. Клигене // Теория вероятностей и ее применения.-1973. — Т. 18 — Вып. 3 С. 677-678.

49. Клигене Н.И. Точное распределение оценки максимального правдоподобия момента изменения параметров авторегрессии / Н.И. Клигене // Статистические проблемы управления. Вильнюс: Изд-во Ин-т математики и кибернетики АН ЛитССР.-1978. - №. 31. - С. 929.

50. Клигене Н.И. Исследование точности оценки максимального правдоподобия момента изменения параметров уравнения авторегрессии/ Н.И. Клигене // Статистические проблемы управления.

51. Вильнюс: Изд-во Ин—т физики и математики АН ЛитССР, 1975. — №. 12.- С. 42-70.

52. Клигене Н. Методы обнаружения моментов изменения свойств случайных процессов / Н. Клигене, Л. Телькснис // Автоматика и телемеханика.-1983. №10. - С. 5-56.

53. Ковальски В. Стойкость доменных печей: современный уровень, развития и мероприятий по продлению компании / В. Ковальски, Х.Б. Люнген, К.П. Штреккер // Черные металлы.-1999. №9. - С. 26-35.

54. Коган А.Е. Внепечные и ковшевые процессы: учебное пособие / А.Е. Коган.-Новокузнецк: Изд-во Сибирского металлургического института, 1990. С. 99.

55. Коган А.Е. Внепечные и ковшевые процессы: Обработка стали при пониженных давлениях: учебное пособие / А.Е. Коган. — Новокузнецк: Изд-во Сибирского металлургического института, 1992. С. 85.

56. Коган А.Е. Особенности гидродинамической обстановки у погружаемой фурмы при продувке стали в ковше нейтральным газом / А.Е. Коган // Изв.ВУЗов. Черная металлургия.-1998. №4. - С. 10 - 12.

57. Короновский A.A. Введение в непрерывный вейвлет-анализ для специалистов в области нелинейной динамики / A.A. Короновский, А.Е. Храмов // Известия ВУЗов. Прикладная нелинейная динамика.-2001.- Т.9, № 4, 5 с. 7-49.

58. Кошелев А.Е. Применение вейвлет-анализа для диагностики состояния продувочной фурмы при внепечной обработке стили / А.Е. Кошелев, O.A. Кравцова // Технолого-экономическое образование в XXI веке.

59. Материалы Всероссийской конференции: В 2 томах. Том 2. — Новокузнецк: Изд-во КузГПА, 2004. С. 90-95.

60. Кривандина В.А. Теплотехника металлургического производства. Теоретические основы: учебное пособие для вузов: В 2Т. Т1 / В.А. Кривандина, В.А. Арутюнов, В.В. Белоусов и др. М.:МИССИС, 2002. - 608 с.

61. Кривенко О.В. Типовые методы оценки контролируемых показателей при решении задач управления / О.В. Кривенко // Приборы и системы. Управление, контроль, диагностики.-2003. №11. - С. 50-54.

62. Кудинов Г.А. Компьютерная диагностика разгара огнеупорной кладки горна и лещади доменных печей / Г.А. Кудинов, В.А. Кришталь, Е.Е. Лысенко // Сталь.-1997. №10. - С. 10-12.

63. Кудинов Г.А. Охлаждение современных доменных печей / Г.А. Кудинов М.: Металлургия, 1988. - 256 с.

64. Куземко Р.Д. Анализ режимов работы для интенсивной продувки стали в 160-т ковше / Р.Д. Куземко, В.А. Наумов, М.А. Поживанов и др. // Сталь.-1997 №2. - С 15-19.

65. Лаар Р. Техническое обслуживание горна доменной печи / Р. Лаар Э. Стейн Калленфельс, М. Геердес // Сталь.-2005. №10. - С. 18-23.

66. Лавров В.В. Разработка и испытание систем диагностики огнеупорной кладки и контроля теплового режима фурменной зоны доменной печи. Автореферат дисс. канд. техн. наук. Сп. 05.16.02 Металлургия черных металлов / В.В. Лавров - Екатеринбург, 1996. - 23 с.

67. Ладонкий А.Н. Диагностика аварийных ситуаций доменной печи с использованием экспертной системы / А.Н. Ладонкий, Б.Н. Жеребин, В.П. Чистов, Пареньков и др. // Металлург. 1997. - №12. — С.46.

68. Левкович-Маслюк Л. Дайджест вейвлет-анализа / Л. Левкович—Маслюк //Компьютера.-1998-№8 (236)-С. 31-37.

69. Левченко В.Е. Оценка влияния некоторых факторов на тепловые потери с водой, охлаждающей шахту доменной печи / В.Е. Левченко, П.Г. Васильев, Н.Е. Алпаев и др. // Известия ВУЗов. Черная металлургия.-1985. №9 - С. 39^3.

70. Липейка А.К. Об определении моментов изменения свойств авторегрессионной последовательности / А.К. Липейка // Статистические проблемы управления — Вильнюс: Изд-во Ин—т математики и кибернетики АН ЛитССР, 1979. № 39. - С. 9 - 23.

71. Липейка А.К. Определение изменений свойств последовательностей авторегрессии / А.К. Липейка // Статистические проблемы управления -Вильнюс: Ин т физики и математики АН ЛитССР, 1975.-№ 12,-с. 27-41.

72. Липейка А.К. Определение моментов изменения свойств авторегрессионной последовательности / А.К. Липейка // Статистические проблемы управления. Вильнюс: Изд-во Ин-т математики и кибернетики АН ЛитССР, 1977. - № 24. - С. 27 — 71.

73. Липейка А.К. Оценка точности момента изменения свойств случайной последовательности / А.К. Липейка // Статистические проблемы управления. Вильнюс: Изд-во Ин - т физики и математики АН ЛитССР, 1973.-№7.-с. 76-87.•'.'•■ 172

74. Лукашин Ю.П. Адаптивные методы краткосрочного прогнозирования временных рядов / Ю.П. Лукашин. — М-: Финансы- и статистика; 2003. — 415 е.: ил.

75. Мизин Д.А. Оценка энтропии динамической системы; / Д.А. Мизин // Автоматика и телемеханика —2002. — №11. — С. 183—189!,

76. Монтвилас А.М. Определение изменения свойств случайных сигналов при неизвестных параметрах этих сигналов / А.М. Монтвилас // Статистические проблемы управления. Вильнюс: Изд-во Ин—т физики и математики АН ЛитССР; 1973. - № 7. -С.8-20.

77. Мясникова Н.В. Аппроксимативный способ вейвлет-анализа /. Н.В. Мясникова, М.П. Берестень // Датчики и системы.-2003. №1. - С. 17-20.

78. Никифоров* И;В. Модификация и исследование процедуры кумулятивных сумм /. И.В. Никифоров // Автоматика и телемеханика.-1980.-№9.-С. 74-80.

79. Никифоров И.В. Применение кумулятивных сумм для обнаружения и изменения характеристик случайного процесса / И.В. Никифоров // Автоматика и телемеханика-1979.;-№. 2. С. 48-58.

80. Никифоров И.В. Применение последовательного анализа к процессам авторегрессии / И.В. Никифоров // Автоматика и телемеханика.-1975. -№8.-С. 174-177.

81. Охотский В.Б. Массообмен газов при продувке в ковше / В.Б. Охотский // Изв. ВУЗов. Черная металлургия.-1996. № 12. - С.8-12.

82. Патент РФ 1838743" A3. МКИ F 27 D 19/00, 21/04. Способ контроля работы печи / Н. А. Спирин, В. С.Новиков, Ю. В. Федулов и др. Опубл. 30.08.93, Бюл. изобр. №32.

83. Патент РФ №2101366 С1. ^ЩИ С21 С7/072. Способ оперативной оценки состояния фурмы при продувке расплава в ковше / В.М. Свекров, В.И. Веревкин, В.К. Буторин и др.-Опубл. 10.01.98, Бюл. изобр. №1.

84. Патент РФ №2113507 СГ МКИ С 21 С 7/072. Способ определения наличия свищей при продувке расплава газом в ковше / В.И. Веревкин, В. К. Буторин, Е.Е. Кошелев и др. Опубл. 20.06.98, Бюл. изобр. №17.

85. Пилюшенко B.JI. Способы вдувания газа при внепечной обработке стали / B.JI. Пилюшенко, А.Н. Смирнов, С.П. Еронько // Черная металлургия: Бюл. института "Черметинформация".-1988. №20.-С 2—13.

86. Пирожников В.Е. Автоматизация сталеплавильного производства / В.Е. Пирожников. М.: Металлургия, 1985. - 184 с.

87. Погоржельский В.И. Контролируемая прокатка / В.И. Погоржельский, Д.А. Литвиненко, Ю.Д. Матросов и др. — М.: Металлургия, 1979.-184 с.

88. Проталинский О.М. Система диагностики предаварийных ситуаций / О.М. Проталинский // Приборы и системы. Управление, контроль, диагностики.-2003. №12. - С. 40^13.

89. Ракитин В.И. Практическое руководство по методам вычислений с приложением программ для персональных компьютеров / В.И. Ракитин, В.Е. Первушин. — М.: Высшая школа, 1998. — 384 с.

90. Редль 3. Моделирование процессов ковшовой металлургии / 3. Редль // Черные металлы.-1995. № 4-5. - С. 31 - 34.

91. Рингель Д. Новая футеровка доменной печи II на заводе фирмы STAHLWERKE BREMEN. / Д. Рингель, Й. Янц, Г. Треккер // Черные металлы.-2000. № 11. - С. 18-22.

92. Сенкус А. Об оценке момента изменения параметров авторегрессионной последовательности / А. Сенкус // Статистические проблемы управления. Вильнюс: Изд-во Ин-т физики и математики АН ЛитССР, 1973.-№7.-С. 19-29.

93. Сергии М.Ю. Структурная оптимизация процесса наблюдения в системе управления / М.Ю. Сергин // Приборы и.системы. Управление, контроль, диагностики.-2001. №3.- С. 1-5.

94. Серов Ю.В. Новые информационные технологии контроля работы горна доменных печей / Ю.В. Серов, В.Г. Макиенко, В.Н. Бражко и др. // Сталь.-1997. № 10. - С. 4-9.

95. Сивирюк В.П. Магнитный контроль напряженного состояния кожуха доменной печи / В.П. Сивирюк, В.П. Горбачев, C.B. Серебряков и др. // Сталь.-1999. №5. - С. 24-25.

96. Смирнов H.A. Внепечная обработка чугуна и стали / H.A. Смирнов // Сталь-1995. N 10. - С. 34 - 36.

97. Соколовский П.И. Малоуглеродистые и низколегированные стали / П.И. Соколовский. М.: Металлургия, 1996. - 216 с.

98. Спирин H.A. Диагностика состояния футеровки доменных печей по температурному полю кладки / H.A. Спирин, Ю.В. Федулов, B.C. Новиков и др. // Сталь—1997. №10. - С. 13-16.

99. Стариков B.C. Огнеупоры и футеровки в ковшевой металлургии: учебное пособие для вузов / B.C. Стариков, М.В. Темлянцев, В.В. Стариков.-М.: МИССИС, 2003.-238 с.

100. Старков Е.Ф. Вейвлет-анализ электрофизических сигналов / Е.Ф. Старков // Датчики и системы—2005. №2. - С. 10-13.

101. Телькснис JI.A. Изменение параметров распределений случайных процессов методом Бейеса / JT.A. Телькснис // Теория и практика измерений статистических (вероятностных) характеристик. Материалы 1-й Всесоюзной конференции. Ленинград, 1973. - С. 87-96.

102. Телькснис Л.А. О применении оптимального бейесова алгоритма обучения при определении моментов времени изменения свойств случайных сигналов / Л.А. Телькснис // Автоматика и телемеханика.— 1969.-№6.-С. 52-58.

103. Телькснис Л.А. Определение изменений свойств случайных последовательностей при неполных априорных данных / Л.А. Телькснис // Статистические проблемы управления. Вильнюс: Изд-во Ин-т физики и математики АНЛитССР, 1975. - № 12 - С. 12 - 26.

104. Телькснис Л.А. Определение наиболее вероятных моментов времени изменения свойств случайных сигналов / Л.А. Телькснис // Автоматика и вычислительная техника.-1970. № 1, - С.24—27.

105. Телькснис Л.А. Определение наиболее вероятного момента времени изменения характера случайного процесса / Л.А. Телькснис,

106. B.Ю. Черняускас // Нелинейные и оптимальные системы. Материалы 1- • го Всесоюзного симпозиума по статистическим проблемам в технической кибернетике. М.: Наука, 1971. - С. 223-228.N

107. Удовенко В.Т. Исследование и разработка технологии продувки конвертерной стали азотом в ковше / В.Т. Удовенко, А.Д. Шевченко,

108. C. А. Донской и др. // Сталь—1979. № 1. - С. 26 - 27.

109. Фильтрация и стохастическое управление в динамических системах / под ред. К.Т. Леондеса. М.: Мир, 1980 - 408 с.

110. ПО.Хаванова H.A. Методы анализа временных рядов: учеб. пособие. / H.A. Хаванова, И.А. Хаванов. Саратов: Изд-во ГосУНЦ "Колледж", 2001.-120с.: ил.

111. Ш.Хайнке Р. Автоматическая система измерения и ремонта огнеупорной футеровки в электропечах фирмы BUDERUS / Р. Хайнке, К. Шеферс, К. Васмут, Б. и др. // Черные металлы.-2003. № 2. - С. 32—35.

112. Шакиров И.В. Методы и алгоритмо-программное обеспечение вейвлет-анализа статических изображений: Дис. . канд. техн. наук. Сп. 05.13.18 Математическое моделирование, численные методы и комплексы программ / И.В. Шакиров. — Томск, 2003. - С 161.

113. Шалтяните В. Вычисление момента изменения авторегрессионной последовательности измеряемой с помехами / В. Шалтяните // Статистические проблемы управления. Вильнюс: Изд-во Ин-т математики и кибернетики АН ЛитССР, 1978. - № 31. - С. 31-42.

114. Шамовский В.Э. Структурный анализ сигналов при алгоритмизации доменного производства: Дисс. . канд. техн. наук. Сп. 05.13.14

115. Автоматическое управление и регулирование в металлургической промышлешюсти / В.Э. Шамовский.-Новокузнецк, 1973. С 222.

116. Швагер Д. Технический анализ. Полный курс / Д. Швагер. М.: Альпина, 2001.-805с.

117. Швадкий B.C. Диагностики футеровки доменной печи с использованием математической модели ее нестационарного теплообмена / B.C. Швадкий, Ю.В. Федулов, H.A. Спирин, В.В. Лавров // Известия ВУЗов. Черная металлургия.-1997. №8. - С. 40-43.

118. Шеффер Х.-Ю. Состояние диагностики машин в черной металлургии / Х.-Ю. Шеффер // Черные металлы.-2001. № 6. - С. 28-32.

119. Ширяев А.Н. Задача скорейшего обнаружения нарушения стационарного режима / А.Н. Ширяев.-Докл. АН СССР. 1961. Т.138. №5. - С. 10391042.

120. Ширяев А.Н. Обнаружение спонтанно возникающих эффектов / А.Н. Ширяев.-Докл. АН СССР, 1961. Т. 138. - №4. - С. 799-801.

121. Шнееров Я.А. Полустойкая сталь / Я.А. Шнееров, В.А. Вихлевщук М.: Металлургия, 1973. - 368 с.

122. Штайгер А.Ф. Разработка и применение моделей и алгоритмов для автоматизированной доводки плавки в агрегатах внепечной обработки стали: Дисс. .канд. техн. наук/ А.Ф. Штайгер-Новокузнецк, 1998—155 с.

123. Шульте М. Контроль износа горна доменной печи с помощью датчиков теплового потока / М. Шульте, Р. Клима, Д. Рингель, М. Фосс // Черные металлы—1998. №4. - С. 17-22.

124. Эппле У. Автоматическое генерирование диагностической информации/ У. Эппле // Черные металлы-2004. № 4. - С. 57-62.

125. Эрлих А. Технический анализ товарных и финансовых рынков. Прикладное пособие / А. Эрлих. -М.: ИНФРА-М, 1996. 176с.ч ,

126. Akaike H. A new look at the statistical model identification. / H. Akaike // IEEE Trans. Automat. Control, 1974. V. 19, N 6, pp. 716-722.

127. Akaike H. procedure for the modeling of non stationary time series. / H. Akaike, G. A Ketagava // Ann. Inst. Statist, 1978. V 30, N 3, pp. 351-363.

128. Johnson R. A. The effect of serial correlation on the performance of CUSUM tests I. / R. A. Johnson, M. L. Bagshaw // Technometrics, 1974. V. 16, №1, pp. 73-80.

129. Johnson R. A. The effect of serial correlation on the performance of CUSUM tests II. / R. A. Johnson, M. L. Bagshaw // Technometrics, 1975. V. 17, №1, pp. 103-112.

130. Jones R. H. A method for detecting change in a time series applied to newborn EEG. / R. H. Jones, D. H. Crowell, L. E. Kapuniai // Electroencephalography and Clinical Neurophysiology, 1969. V. 27, pp. 436440.

131. Jones R. H. An adaptive method for testing for change in digitized cardiotachometer data. / R. H. Jones, D. H. Crowell, L. E. Kapuniai // IEEE, Trans, on Bio-Medical Engineering, 1971, V 5, pp. 360-365.

132. Jones R. H; Crowell D. H., Kapuniai L. E. Change detection model for serialy correlated multivariate data. / R. H. Jones, D. H. Crowell, L. E. Kapuniai // Biometrics, 1970, V. 26, N 2, pp. 269-280.

133. Newbold P. M. Detection of changes in the characteristics of a Gauss. / P. M. Newbold, Yu Chi Ho. // IEEE Trans. Aerospace and Electronic Systems, 1968, V.AES-№5, pp. 707 718.

134. Page E. Continuous inspection schemes. / E. Page — Biometrika, 1954, V. 41, N2, pp. 100-114.

135. Segen J. Detecting change in time series. / J. Segen, A. Sanderson // IEEE Trans. Inform. Theory, 1980, V. IT-26, N2, pp. 250-355.