ИМИТАЦИОННАЯ МОДЕЛЬ ДЛЯ УПРАВЛЕНИЯ И ИЗУЧЕНИЯ ТЕМПЕРАТУРНОГО РЕЖИМА КОКСОВАНИЯ

Чистякова Т.Б., Бойкова О.Г., Бабина Е.В., Нестеров Д.М.

Санкт-Петербургский государственный технологический институт (Технический университет)

Источник:http://conf.bstu.ru/conf/docs/0030/0738.doc

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

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

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

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

Система 10 уравнений

где а - коэффициент температуропроводности; л - теплопроводность; k - коэффициент теплопередачи греющей стенки; t - время коксования; x - расстояние точки от оси камеры; L - ширина камеры; С_p - удельная теплоемкость при постоянном давлении; гамма_k - кажущийся удельный вес тела; гамма_kш - кажущийся удельный вес шихты; T_ш - температура в любой точке коксующейся шихты; F_ст - площадь греющей стенки; а_ог - коэффициент температуропроводности отопительного газа; л_ст - теплопроводность стенки; б_ст - толщина стенки; а_ш - коэффициент температуропроводности шихты; G_ог - расход отопительного газа; С_ог - теплоемкость отопительного газа; G_ш - расход шихты; Т_ст_ш - температура шихты у стенки; Т_оп_ш - температура шихты в осевой плоскости; Т_н_ш - начальная температура шихты; Т_н_ог - начальная температура отопительного газа; Т_к_ог - конечная температура отопительного газа; а_ст_ш - коэффициент температуропроводности шихты у стенки; С_ст_ш - теплоемкость шихты у стенки.

Решение (1) при условиях (2) - (6), позволяет определить время проведения процесса при заданных расходах шихты и отопительного газа и конструктивных характеристиках печи (ширина камеры, толщина стенки, материал, из которого изготовлена стенка) и температуру кокса при окончании процесса коксования, а также получить оценки тепловых свойств кокса по зависимостям (8)-(10).

Формализованное представление ММ может быть описано следующим набором:

Формализованное представление ММ

Вектор входных параметров:

Вектор входных параметров

где Х1 и Х2 – характеристики шихты и отопительного газа соответственно.

Характеристики шихты и отопительного газа

Вектор параметров и коэффициентов модели:

Вектор параметров и коэффициентов модели

где А1 и А2 – характеристики печи и шихты соответственно.

Вектор параметров и коэффициентов модели

Вектор выходных параметров:

Вектор выходных параметров

Вектором управляющего воздействия данной модели U_м является расход отопительного газа G_ог .

Вектор параметров и коэффициентов модели

Структура системы управления температурным режимом коксования с помощью разработанной имитационной модели представлена на рис.1.

Математическая модель (ММ) описывает поведение системы при различных условиях проведения эксперимента.

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

База знаний нештатных ситуаций (БЗ НС) содержит описания основных нештатных ситуаций, возможных при управлении объектом и рекомендации по их устранению.

Рис. 1 Структура системы управления, где ММ – математическая модель; БД ТП – база данных технологического процесса; БЗ НС – база знаний нештатных ситуаций; БЗ печей – база знаний печей; УСО – устройство связи с объектом; RNS – вектор рекомендаций по управлению в НС; Эмулятор НС – эмулятор нештатных ситуаций.

Рис. 1 Структура системы управления, где ММ – математическая модель; БД ТП – база данных технологического процесса; БЗ НС – база знаний нештатных ситуаций; БЗ печей – база знаний печей; УСО – устройство связи с объектом; RNS – вектор рекомендаций по управлению в НС; Эмулятор НС – эмулятор нештатных ситуаций.

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

Вектор рекомендаций по управлению в нештатных ситуациях (RNS) представляет собой рекомендации оператору коксовой батареи для действий при возникновении нештатных ситуаций. Формализованное представление нештатных ситуаций может быть представлено в виде вектора:

Вектор рекомендаций по управлению в нештатных ситуациях (RNS)

где Sit – описание НС, Prich{A1,A2,X1,X2} – причина ситуации, которая может быть задана изменением входных и настроечных переменных модели, Sov – совет оператору с рекомендациями о действиях в НС.

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

Нормирующий множитель для перехода к безразмерной температуре

где T0= T_н_ш – температуре загружаемой шихты; T_max – температура готового кокса. Tmin=T0. Tmax определяется из формулы:

Находим T_max

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

Далее вводится нормирующий множитель по X: X1=X/(L/2), где L – ширина камеры печи. На следующем шаге необходимо задать количество узлов сетки. Для этого задается время протекания процесса t в диапазоне от 13 до 27 часов, затем задается шаг по времени – и по х – h, что позволяет определить соответственно количество временных слоев – n и слоев шихты – k, для которых будет производиться расчет. После этого формируется система уравнений для n+1 временного слоя:

Система уравнений для n+1 временного слоя

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

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

Нештатные ситуации, связанные с температурным режимом коксования частично описаны в [3,4]. На расход отопительного газа в процессе коксования накладываются ограничения, связанные с особенностями и сроком эксплуатации теплотехнических агрегатов. Поэтому к разряду нештатных ситуаций в данной модели могут быть отнесены решения, оптимальные с точки зрения свойств кокса, но выходящие за рамки ограничений, накладываемых на G_ог . Кроме того, к нештатным ситуациям может быть отнесена ситуация, когда при заданных векторах Х,А недостижима возможность получения кокса, с заданными тепловыми свойствами Y0 .

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

Литература

1. Чистякова Т.Б. Интеллектуальные тренажерно-обучающие системы: Дисс. … док. техн. наук – СПб, 2000.

2. Чистякова Т.Б., Бойкова О.Г., Бабина Е.В. Моделирование тепловых свойств кокса // Математические методы в технике и технологиях – ММТТ-17: Сб. трудов Международ. науч. конф. – Кострома: Изд-во КГУ, 2004. - Т.1.

3. Бойкова О.Г. Гибридная экспертная система для управления процессами коксования: Дисс. … канд. техн. наук – СПб, 2000.

4. Блохина О.Ф. Оперативное управление коксохимической батареей. Дисс. ... канд. техн. наук – СПб, 2003.