bbb

advertisement
РАЗРАБОТКА ПРИНЦИПОВ СОЗДАНИЯ СИСТЕМЫ УПРАВЛЕНИЯ
КОТЕЛЬНЫМ АГРЕГАТОМ ТГМЕ-464 В ПЕРЕХОДНЫХ РЕЖИМАХ С
ПОМОЩЬЮ ИСКУССТВЕННЫХ НЕЙРОННЫХ СЕТЕЙ
ЧЕРНОВ А.С., АБРУКОВ В.С., КОЩЕЕВ И.Г., ТРОЕШЕСТОВА Д.А.
ЧУВАШСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
ОГЛАВЛЕНИЕ
Введение …………………………………………………………………………. 2
1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПАРОГЕНЕРАТОРА ……………………2
2. АВТОМАТИЗИРОВАННАЯ СИСТЕМА УПРАВЛЕНИЯ
НА ОСНОВЕ ИНС …………………………………………………………….5
3. ЗАКЛЮЧЕНИЕ ………………………………………………………………..8
4. ЛИТЕРАТУРА …………………………………………………………………8
5 ОБОЗНАЧЕНИЯ ………………………………………………………………..8
6 ПРИЛОЖЕНИЕ
База данных для настройки нейроэмулятора ………………………….……10
ВВЕДЕНИЕ
Проблема создания автоматизированных систем управления в современной
теплоэнергетике является одной из наиболее актуальных. В частности, особое
внимание обращается на повышение уровня автоматизации, надежности,
безопасности и гибкости системы управления.
Данная работа посвящена разработке принципов создания системы управления
котельным агрегатом ТГМЕ-464 в переходных режимах (изменение нагрузки,
пуск и останов котельного агрегата и т.д.) на основе применения нейросетевых
технологий. Рассматривалась задача регулирования температуры пара в
пароперегревателе.
1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПАРОГЕНЕРАТОРА
Парогенерирующий канал в общем случае состоит из трех подсистем: оболочки и
двух сред, омывающих ее. В свою очередь каждая подсистема может быть
представлена в виде различных моделей, отличающихся набором допущений. При
аналитическом исследовании, прежде всего надо выбрать состав системы и
модели подсистем, входящих в нее. При этом каждый выбранный вариант дает
определенную модель парогенерирующего канала. Наличие трех подсистем и
большое количество частных моделей оболочки и потоков двух сред
обеспечивает большое разнообразие моделей парогенераторов. Все эти модели
можно разбить на три группы: распределенные, сосредоточенные и
комбинированные. Выбор модели определяется задачами исследования и
возможностями вычислительной техники.
Парогенератор обычно разбивается на участки, одни из которых рассматриваются
как системы с распределенными, а другие – с сосредоточенными параметрами.
Такое деление позволяет упростить решение. Потеря точности при этом может
быть небольшой, поскольку роль отдельных участков в формировании
динамических процессов неравноценна. Деление на участки позволяет также
использовать готовые решения. Расчетная, структурная схема математической
модели парогенератора зависит от конкретного вида технологической цепочки
элементов пароводяного и газового трактов и требований к степени полноты
контроля в динамике за состоянием параметров в различных точках тракта.
В общем случае, барабанный котельный агрегат можно представить в виде
последовательного соединения экономайзера, испарительной части, включающей
в себя барабан, опускные и подъемные трубы, и пароперегревателя.
Вычислительная модель каждого элемента аналогична схеме на рис.1, а
структурная схема объекта принимает вид рис.2.
Рис.1. Модель процессов в элементах котлоагрегата.
Рис.2. Структурная схема котлоагрегата.
Определим пароводяную смесь как двухфазный поток, каждая фаза которого
представляет собой вязкую, вообще говоря, сжимаемую жидкость. Основные
дифференциальные
уравнения,
полностью
определяющие
поведение
динамической системы в любой момент времени следующие [1].
Уравнение движения, удовлетворяющее закону Навье – Стокса, в произвольной
криволинейной системе координат
(1)
 u j  g ij u i  j u j   Fi +  j p ij .
Уравнение энергии


1
1
1
cv T  g ij u j  j =  Fi i ui  p ij  j ui +  i g ij  T  jT   qvi  .



(2)
Уравнение неразрывности
  g ij  j u j  0 ,
(3)
 p, v, T   0 .
(4)
где g ij - компоненты метрического тензора.
Уравнение состояния
Уравнения (1)-(4) вместе с краевыми условиями полностью описывают поведение
динамической системы в любой момент времени. Каждое из этих уравнений
может быть представлено как совокупность уравнений, описывающих движение
двухфазного потока с учетом турбулентного характера течения каждой фазы.
Рассмотрим замыкающее звено котельного агрегата – пароперегреватель,
выходные параметры которого являются выходными параметрами всего
теплового технологического цикла [2]. Здесь реализуется регулирование
температуры, а далее по паровому тракту – регулирование давления. Последнее
является существенным возмущающим воздействием и должно быть учтено в
математической модели, хотя расположение возмущения имеет место вне зоны
технологического цикла генерации и перегрева пара. Для дальнейшего изложения
определим основные абстракции, в рамках которых и будет создана
математическая модель пароперегревателя как динамического звена.
Конструктивно пароперегреватель представляет собой многорядный змеевик с
большим количеством изгибов. Создание математической модели полностью
отвечающей действительности, не представляется возможным, ввиду сложности
конструкции пароперегревателя, характера наружного обтекания поверхностей
нагрева продуктами сгорания (конвективный теплообмен) и наличия
радиационного теплообмена, характера движения нагреваемого теплоносителя
(турбулентность). Поэтому пароперегреватель будет представлен как
прямолинейный трубопровод, разделенный на три участка: входной – здесь поток
нагреваемого теплоносителя проходит тепловую и гидравлическую стабилизацию
(математически это выражается масштабными коэффициентами, учитывающими
тепловую и гидравлическую неравномерность на входе в реальный
пароперегреватель); управляющий – здесь стабилизированный поток
подвергается управляющему воздействию, в связи с этим параметры
теплоносителя изменяются согласно требуемому переходному режиму либо
регулируются до значений параметров предшествующего стационарного режима;
выходной – здесь поток нагреваемого теплоносителя вновь проходит тепловую и
гидравлическую стабилизацию.
Представим пароперегреватель как элемент с локально распределенными
параметрами: параметры распределены по времени – необходимое условие
моделирования динамического объекта, а также по дискретному пространству,
введение которого необходимо для получения из нелинейных дифференциальных
уравнений эволюционных уравнений динамического объекта.
2. АВТОМАТИЗИРОВАННАЯ СИСТЕМА УПРАВЛЕНИЯ НА ОСНОВЕ
ИНС
После перевода уравнений (1)-(3) в цилиндрическую систему координат,
дискретизации пространства и выделения в объеме пара некоторого конечного
элемента, была получена система конечно-разностных уравнений, например, для
уравнения массового расхода воды на инжекцию:
  v u  v
 M 
M    r ur   u
  uz
 v  r 



r
   ur

 uz


t
r  r

z
r 
 z  v 
 r

 M  M    r 1   ur  1   u
 1   uz


r


t
r 
r

z
  v u  v
 v  r 
  ,
 1    u r

 u z

r
r



z

 v 
(5)
где  - объемное паросодержание, т.е. доля поперечного сечения канала, занятая
в данный момент времени паром, u - скорость течения фазы, осредненная по всей
лощади, занимаемой ею в данный момент времени в данном поперечном сечении,
при этом одним штрихом обозначаются параметры воды, двумя – пара, конечноразностная форма имеет вид:
M i, j ,k ,t
 M 

t
ri
 i , j ,k 1, ui, j ,k 1,t
 ri 1 i 1, j ,k ,t ui1, j ,k ,t  i , j 1,k ,t ui, j 1,k ,t


r

i

r

z

 r
r 
1
  i , j ,k ,t ui, j ,k ,t  i 
 i  
  r   z 
vi1, j ,k ,t  vi, j ,k ,t ui , j 1,k ,t vi, j 1,k ,t  vi, j ,k ,t
vi, j ,k 1,t  vi, j ,k ,t

  i , j ,k ,t  ui1, j ,k ,t

 ui, j ,k 1,t
z
r
ri


M i, j ,k ,t 1  M i, j ,k ,t


t

 ri 


 vi, j ,k ,t 
1   i, j ,k 1,t ui, j ,k 1,t 
M i, j ,k ,t  ri 1 1   i 1, j ,k ,t ui1, j ,k ,t 1   i , j 1,k ,t ui, j 1,k ,t

 ri

ri 
r

z
 r
1
r 
 
 1   i , j ,k ,t ui , j ,k ,t  i 


r



z




v
v
u
v
 v

 1   i , j , k ,t  ui 1, j , k ,t i 1, j , k ,t i , j , k ,t  i , j 1, k ,t i , j 1, k ,t i , j , k ,t 
r
ri


 ui, j ,k 1,t
vi, j ,k 1,t  vi, j ,k ,t  ri 
 
 

   M    M  .

z
 vi , j ,k ,t 
(6)
Были выбраны управляемые и управляющие параметры [3]. Скорость изменения
температуры пара на выходе пароперегревателя была выбрана как управляемый
параметр, а массовый расход воды через инжектор был выбран в качестве
управляющего параметра. Аналитическая связь между ними была определено с
помощью математического аппарата производных Ли. При этом использовалось
модифицированное уравнение (6):
  4


 M
5 
 M 

 M   1  3  
 M   
 6  

t
 r
 t
 1   v 
  r 1   

 M  2    1  ,
(7)
 v
 *  T   T * , 1  L f *   0 * ,
       M        T     
L f   

  M    t
  T    t   





2  L f 1  11 ,
     

 
 t  
 

.
 t
(8)
Желаемое управление  может быть найдено из уравнения
 0.
(9)
 2 2  
2
Тогда
i
2

1   v   v   T  
 

M  5
1
*
 
   i ,

T

t
 

(10)
где  - некоторый оператор с действием, осуществляющим отображение вида
*i 
 i   n ,   0;1, n  1, 2, .

(11)
Уравнение (10) можно представить в конечных разностях следующим образом.
*1 
1    *3 M~  v k 1 2 c p,k 1
M 2  5 v k 1  vk  c p,k 1
Откуда
*2 
~
 A *2 M  c p,k 1 (T )
1
~
A M  cp, k 1 T 
так как при T   0
*1  1 .
~
Здесь принято M   M   138,8889 кг/с; M   3,1945 кг/с;
    0,5 ;
  1 с;
  0,95 ; 5  1 , A  2,2567  10 5 .
С помощью полученных аналитических связей была создана база данных для
обучения ИНС, часть которой показана ниже.
Работа эмулятора NNW предполагает
нормализацию
области
определения
*
функции 1 . Наиболее простым способом
нормализации
является
представление
входных данных в виде
*1  *1min
 *
.
1max  *1min
(12)
В этом случае можно реализовать значения
*1 min  *1  0 ,
  *1 max
но
при
распределение входных параметров будет
крайне неравномерным, что приведет к
ухудшению качества обучения. Поэтому в
подобных ситуациях, а также в случае,
когда значение входа лежит в диапазоне
0;  можно использовать нормализацию с
помощью функции вида

1
1  e u1
*
.
(13)
При моделировании системы управления
была
предопределена
следующая
конфигурация нейронной сети:
Tk1  ,
К/с
1    M  
M 2  5 
vk 1 2 cp,k 1

vk 1  vk cp ,k 1
*1 
0
0.100334672
0.202710035
0.309336249
0.422793218
0.546302489
0.684136808
0.84228838
1.029638557
1.138332713
1.260158218
1.398332713
1.557407725
1.74331531
1.964759657
2.234496949
2.572151622
3.009569674
3.602102448
1
0.732976
0.587377
0.490143
0.369166
0.296411
0.2294
0.187291
0.147342
0.131829
0.117151
0.103618
0.09186
0.08045
0.07002
0.059948
0.050464
0.041234
0.03207
в качестве входных полей использовались: Tk1  и (Tk1  )  1 ;
*
в качестве выходных полей использовались: *1 и  3 
1
:
~
A M  cp, k 1 T 
число входов: 2;
количество скрытых слоев: 2;
число нейронов в скрытых слоях: 6;
число выходов: 2;
для обучения использовалось: 80% записей;
скорость обучения: 0.1.
На момент 1400 эпох обучения было распознано 100% записей из базы данных.
ЗАКЛЮЧЕНИЕ
Была определена оптимальная архитектура ИНС, проведено обучение ИНС и
создана ИНС-модель системы управления. Был разработан прототип системы
автоматического управления температурой пара в пароперегревателе
котлоагрегата TGME-464. Он состоит из эмулятора ИНС, установленной на
персональном
компьютере,
промышленного
микроконтроллера
Philips
P89LPC935, кабельных линий и коммутационных устройств. Прототип может
быть использован как надстройка к основной системе управления. Полученный
прототип может работать, как в режиме уведомления (совета), так и в режиме
автоматического управления. Он может быть достаточно просто интегрирован в
существующие системы управления.
Особенности технологий ИНС управления [4]: свойство адаптируемости к новым
условиям и способность самообучению, обеспечивают простоту модернизации и
наращивания возможностей.
ЛИТЕРАТУРА
1. Седов Л.И. Механика сплошной среды. Т.1. М.: Наука, 1976. 536 с.
2. Серов Е.П., Корольков Б.П. Динамика парогенераторов. М.: Энергоиздат, 1981.
406 с.
3. Нейрокомпьютеры и их применение. Кн.8. Нейросетевые системы управления.
Учебное пособие для студентов ВУЗов /В.А. Терехов, Д.В. Ефимов, И.Ю.
Тюкин.; Под общ. ред. А.И. Галушкина.М.: ИПРЖР, 2002. 480 с.
4. Potential of artificial neural networks in power system operation. Damborg M.J., ElSharkawi M.A. and oth. IEEE Int. Symp. Circuits and Syst. La, May 1990. p. 29332937
ОБОЗНАЧЕНИЯ
T , T  температура воды и пара, К
M , M  массовый расход воды и пара,
кг/с
T , M  скорости изменения температуры
и расхода пара (обозначение
производной по времени), кг / с 2
c p (cv ) удельная изобарная (изохорная)
теплоемкость воды и пара,
кДж/(кг К)
v , v удельный объем воды и пара,
м 3 / кг
  ,  
  ,  
  ,  
плотность воды и пара, кг / м 3
теплопроводность воды и пара,
вт/(м∙К)
динамическая вязкость воды и
пара, Па∙с
u , u  скорость движения воды и пара,

 j , i
Ti,j ,k ,
1

i
i , i
Lf
м/с
объемное паросодержание
оператор ковариантного и
контравариантного
дифференцирования
значение температуры пара в
формализме конечных разностей,
порядок первых трех индексов
указывает на цилиндрические
координаты r , , z , индекс 
указывает на время
величина управления
величина возмущения
Макропеременные
постоянные или
слабоменяющиеся параметры
оператор Ли
Таблица 1. База данных для настройки нейроэмулятора.
Tk1  ,
(Tk1  )  1 ,
К/с
К/с

*
1
1   

~
u3*M 
M   
2
5
vk 1 2 cp, k 1

vk 1  vk cp , k 1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0.100334672
0.202710035
0.309336249
0.422793218
0.546302489
0.684136808
0.84228838
1.029638557
1.138332713
1.260158218
1.398332713
1.557407725
1.74331531
1.964759657
2.234496949
2.572151622
4.45522176
3.602102448
3.009569674
2.572151622
2.234496949
1.964759657
1.74331531
1.557407725
1.398332713
1.260158218
1.138332713
1.029638557
0.84228838
0.684136808
0.546302489
0.422793218
0.309336249
0.202710035
0.100334672
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
0.732976
0.587377
0.490143
0.369166
0.296411
0.2294
0.187291
0.147342
0.131829
0.117151
0.103618
0.09186
0.08045
0.07002
0.059948
0.050464

*3 
1
~
A M  c p, k 1 T 
0.119574963
0.119673637
0.119736514
0.11977697
0.119817453
0.119839956
0.119871474
0.119880482
0.119898502
0.119912021
0.119930051
0.119943577
0.119961616
0.119970638
0.119993198
0.120002225
0.120011253
0.120024797
0.120033828
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
3.009569674
3.602102448
4.45522176
0
0
0
0.041234
0.03207
0.022045
1
1
1
Download