УДК 621.18
ПРИМЕНЕНИЕ ОБОБЩЕННОГО ТЕРМОДИНАМИЧЕСКОГО АНАЛИЗА В ЗАДАЧЕ ОПРЕДЕЛЕНИЯ КООРДИНАТ ТЕХНОЛОГИЧЕСКИХ ОБЪЕКТОВ УПРАВЛЕНИЯ
Д.Ю. ТВЕРСКОЙ, канд. техн. наук
Рассмотрена задача определения управляемых координат пылесистем прямого вдувания котлов тепловых электростанций как сложных технологических объектов управления. Дано обоснование обобщенных термодинамических координат пылесистем, применение которых позволяет оценивать неконтролируемые технологические параметры и решать широкий класс задач управления и диагностирования.
Ключевые слова: обобщенный термодинамический анализ, управляемые координаты, пылесистемы прямого вдувания.
USING GENERIC THERMODYNAMIC ANALYSIS IN DETERMINATION TASK OF CONTROL TECHNOLOGICAL OBJECT COORDINATES
D.Yu. TVERSKOY, Candidate of Engineering
The author considers the determination task of controllable coordinates for coal-pulverization systems with direct injection of thermal power plant boilers as complex technological control units. The author describes the generic thermodynamic coordinates of coal-pulverization systems. The usage of coordinates allows to estimate uncontrolled technological parameters and to solve a wide range of control and diagnostics tasks.
Keywords: generic thermodynamic analysis, controllable coordinates, coal-pulverization systems with direct injection.
мах пылеприготовления, сложностью протекания мельничных процессов размола, сушки, перемещения пылепотоков и невозможностью непосредственного контроля основных параметров. При этом на физические переменные ТОУ (технологические параметры: расход су-шильно-вентилирующего агента, загрузка мельницы топливом, температура аэросмеси, влажность пыли и др.) накладываются расчетно-конструкторские ограничения, при выполнении которых имеет место устойчивый массопо-ток топливовоздушной смеси и управляемость изделия в заданном диапазоне нагрузок. Основные допущения для нелинейной динамической модели рассматриваемого класса объектов сведены в таблицу [4, 6, 7, 8, 9] (табл. 1).
Ниже задача определения координат сложного ТОУ рассмотрена на основе развития обобщенного термодинамического анализа [10-12].
Обобщенные термодинамические потенциалы и координаты систем пылеприготовления по схеме прямого вдувания. Для того чтобы провести обобщенный термодинамический анализ эффективности той или иной сложной системы, необходимо одной группе переменных ТОУ приписать содержательный смысл обобщенного потенциала (обобщенных сил) X, а другой группе переменных - смысл обобщенных термодинамических координат xj. При этом Xjdxj должно иметь смысл элементарной работы. Обобщенная физическая величина X является количественной мерой интенсивности процесса обмена. Обобщенная физическая величина Xj рассматривается как объект переноса, иными словами, является координатой состояния системы. Кроме того, одной из переменных необходимо приписать смысл абсолютной температуры, или энтропии.
Введение. Задача получения необходимой информации о координатах технологических объектах управления (ТОУ) решается, как правило, из опыта или на основе экспериментальных исследований [1-3].
Идентификация тОу связана с определенными рисками при проведении экспериментальных исследований и с ограниченными возможностями непосредственного контроля технологических параметров [4, 5]. Например, для пылеугольных котлов с пылесистемами прямого вдувания задача идентификации осложнена спецификой технологических участков, в которых формируются потоки топливовоздушных смесей, направляемые в топку котла (см. рисунок).
мельницам
Принципиальная схема пылесистем прямого вдувания котлов ТЭС: 1 - бункер сырого угля; 2 - питатель сырого угля с электроприводом; 3 - электропривод мельницы; 4 - мельница; 5 - сепаратор; 6 - канал подачи готового топлива в горелки; 7 - горелка; 8 - топка котла; 9, 10 - регулирующие органы подачи первичного и вторичного воздуха
Эти особенности связаны с физикой формирования топливовоздушных потоков в систе-
Таблица 1. Основные ограничения (допущения)
Наименование фактора Исходные посылки
Расчетно-конструкторские ограничения на значения технологических параметров заданы типоразмером мельницы - топливовоздушный поток устойчивый, режим расчетный (номинальный), установившийся
Распределение пыли внутри объема мельницы и в потоке мельничного продукта неравномерное (неизвестное)
Состояние потоков топливовоздушной смеси однозначно определяются характеристиками пылепотоков (концентрация, скорость движения, температура)
Химические реакции между компонентами пылевоздушной смеси отсутствуют
Сушка топлива заканчивается в объеме мельницы
Расход сушильно-вентилирующего агента по пространственным координатам не изменяется
«Элементарный объем» размеры не определены
Скорость центра масс «элементарного объема» для всех компонентов многофракционного потока топливовоздушной смеси V- = м/с с Ш
Оценка плотности р3 «элементарного объема» потока мельничного продукта (пылевоздушной смеси в объеме установки) - м 3 р3 = —, кг/м И3 V
Оценка плотности р2 «элементарного объема» потока готовой топливовоздушной смеси - мп 3 р2 = ——, кг/м3 V
Масса готовой пыли определяется через ее концентрацию р*2 в потоке готовой первичной топливовоздушной смеси мп = Р2 • V, кг
Концентрация р*2 пыли в потоке готовой первичной топливовоздушной смеси —* В2 <3 Р2 = -=£-, кг/м вс
Концентрация р*3 мельничного продукта —* В3 3 Р3 = , кг/м3 вс
Для устойчивого массопотока (расчетный режим, см. п.1) р3 > р* , кг/м3
Полагая, что рз = р3 ММ с открытым ротором (встроенный гравитационный сепаратор) м ^ В3 = — мвс вс V 3 V с
Полагая, что В3 = кхВ^ , получаем исходную посылку (допущение) математической модели П. Профоса кх - экспериментальный коэффициент кх В2 =м , кх = м?с. хвс V х VB2
Полагая, что Р3 > р3 ММ с закрытым ротором (вынесенный сепаратор инерционный, центробежный) В3 = а — Мвс 3 V с
Рассмотрим переменные системы пыле-приготовления прямого вдувания как ТОУ в указанном содержательном смысле, выполняя требование сохранения размерности, Дж: бА, = Х,бх,, (1)
где X - обобщенный потенциал; бх, - изменение обобщенной термодинамической координаты состояния системы.
Надо отметить, что в выборе координат существует известная свобода, а следовательно, существует известная свобода в выборе и сопряженных с ними потенциалов. Правильность выбора потенциалов по имеющимся зарядам проверяют по формуле
X, =
дА_
дх і
(2)
где А = ^А/ - сумма элементарных работ.
Рассмотрим следующие виды работ, совершаемых в системе пылеприготовления: работу перемещения потока мельничного продукта в пространстве мельничной системы; работу перемещения потока готовой пыли; работу по преодолению сопротивления установки; работу
по преодолению гравитационных сил; технические работы, определяемые работой ротора; тепловую работу в системе.
Гидродинамические работы. Работа перемещения топливовоздушного потока в объеме системы мельница-сепаратор-пылепровод равна изменению кинетической энергии системы:
бА = б
= Му6у + — бМ, 2
(3)
где V = чс - скорость, м/с; М - масса перемещаемого топлива в потоке, кг.
Учитывая правила определенной свободы выбора зарядов, нетрудно видеть, что количественная мера интенсивности процесса переноса (обобщенный потенциал) есть удельный объем системы (объем системы, отнесенный к сечению ротора), который определяется строго конструкцией установки, а термодинамическими координатами, характеризующими состояние системы, являются сложные комплексиро-ванные параметры. При этом координата х-| может быть определена по косвенным показа-
телям и использована при оценке состояния процесса.
Аэродинамическая работа, совершаемая потоком топливовоздушной смеси, направлена на преодоление сопротивления установки:
СА3 = цСт, (4)
где ц - гидродинамический потенциал, Дж/кг; т - масса перемещаемого продукта, кг, бт = С(рV) = Vdр, где V- объем системы, м3; р -3 Р
плотность потока, кг/м ; ц =—, Р - давление в
р
системе, Па.
В результате аэродинамический (гидродинамический) потенциал (обобщенная сила системы) определяется давлением в системе, Па:
*3 = Р,
а обобщенная координата (аэродинамический заряд) - объемом системы, м3:
х3 = V.
Очевидно, что увеличение объема системы (изменение типоразмера мельницы) ведет к увеличению аэродинамической работы, совершаемой потоком. Рассматривая процесс в действующей системе приготовления топливовоздушной смеси, в которой объем V величина фиксированная, а работа совершается при соответствующей потере напора, и учитывая известную свободу в выборе X, и х, , можно записать
*3 = V;
х3 = Р.
Работа преодоления гравитационных сил определяется перемещением потока мельничного продукта (в том числе, готовой пыли) в поле сил тяжести (по отношению к мельничному продукту гравитационные силы выполняют своего рода роль ступени сепарации):
СА4 = С(вЛ), (5)
где в = Мд - вес тела, Н; Л - элементарное перемещение по вертикали, которое определяется приведенной высотой установки, м; М - масса перемещаемого продукта, кг.
В результате гравитационный потенциал системы (обобщенная гравитационная сила) определяется приведенной высотой Л (для молотковых мельниц - отношение объема к сечению ротора), а координатой состояния системы является вес переносимой субстанции, кг-м/с2 = Н: х4 = Мд . (6)
Полученный результат согласуется с известным расчетом шахтных (гравитационных) сепараторов молотковых мельниц [6]. Технические работы Работа вращения ротора. Ротор мельницы представляет собой вращающуюся систему размольной камеры, которая обладает кинетической энергией вращения. Эта энергия определяется по формуле
(7)
где J - момент инерции системы, Нм с2; ю - угловая скорость вращения системы, 1/с.
Определим составляющие (7). Известно, что кинетическая энергия вращения находится пу-
2
т.1/.
тем суммирования кинетических энергий
2
всех частиц системы, обладающих массой т, и скоростью V,. Скорость каждой частицы находится по формуле
= ю • 0,
где г, - расстояние от рассматриваемой частицы до оси вращения, м.
Суммирование дает
,2-2
Ет/ Х""' т/О) г со V"' 2 «Ю
-2~ = ^^Т~=Т ^ тг =-^.
Таким образом, момент инерции системы
равен
и=Х т г2,
/'
угловая скорость, р/с,
= 2пп
Ю= "б^
где п - число оборотов ротора мельницы, об/мин.
Окружная скорость определяется скоростью бил:
жйп
V, = иб =---,
1 б 60
а ^т I - массой Мр вращающихся частей ро/'
тора и массой М топлива, циркулирующего в системе:
^ т I = тр + М = М,
В результате получаем
( _____2
dAs = d
z-
= d I u6Mx
(8)
v i
Потенциал ротора (обобщенная сила) мельницы определяется квадратом окружной скорости бил (м2/с2), а координатой состояния системы является масса Mx системы. Полагая, что mp = const (режимные периоды времени), то обобщенная координата будет определяться только массой топлива.
Работа механического (дисперсионного) превращения в камере размола. Известно, что при размоле топлива изменяется масса m, отдельных составляющих массовой загрузки. Например, масса исходного продукта, поступающего в мельницу, с плотностью р1 превращается в мельничный продукт с плотностью р3, а затем в сепараторе - в готовую топливовоздушную смесь с плотностью р2. Чтобы учесть эту особенность процесса дисперсионного превращения, формулу для работы представим в обобщенном виде, Дж:
САб = ЦбСтб, (9)
где ц6 - механический потенциал системы, Дж/кг; т6 - «дефект» массы, кг.
Очевидно, ц6 равен обобщенному потен-
22
циалу ротора мельницы, м /с , или, приводя к единой размерности, получаем, м,
Xб = M-б =
где ускорение свободного падения д можно интерпретировать как нормирующий множитель.
«Дефект» массы т6 определяется разностью количества М топлива, циркулирующего в системе, и массы Мв возврата мельничного продукта, кг: т6 = М - Мв,
т.е. т6 характеризуется массой Мп готового продукта, покидающего установку. Тогда
dAб = 2 u2d(M
(10)
где координата состояния, характеризующая механическое превращение, с учетом нормирующего множителя приобретает вид
х6 = вп = в2 Vg.
вс
Поскольку работа, затрачиваемая на измельчение, есть доля полной работы вращения ротора, то полезный механический потенциал камеры размола не больше, очевидно, обобщенного потенциала ротора мельницы. Тогда
СА6 = ±и2 -пСМ, (11)
где п - доля обобщенного потенциала, затрачиваемая на размол и доведение топлива до заданного дисперсионного состояния (другими словами, обобщенный КПД пылесистемы). Исходя из равенства работ А6 = А*, получаем, что отношение полезно затраченной работы к полной равно
п =
M - Mg =1 ^ .
M
M
B
Учитывая, что M - Mg = — V , находим
п = -
Go
M
p2.
p3
(12)
Полученный результат (14-15 %) хорошо согласуется с экспериментальными оценками [7, 8], а также отражает физику реальных процессов. Например, при износе бил загрузка М растет при одновременном снижении производительности В2 мельничной системы, т.е. п уменьшается.
Термическая работа в системе определяется количеством переданной теплоты, Дж:
бА8 = бО = ГбБ, (13)
где Т - абсолютная температура, К; 5 - энтропия, Дж/К;
dQ = d (cmMT + cM mMT),
где Ст, См - теплоемкость соответствующего топлива и металла, Дж/кг-К; тм - масса металла, аккумулирующая тепло. Будем полагать, что теплоемкость и масса металла не изменяются (Ст, См = const). Тогда при T = const для конкретной марки топлива находим dQ = d (cmMT),
и термическая работа определяется как
dA8 = Td(cmM), (14)
где термический потенциал, К,
*8 = Т,
а термической координатой является параметр, имеющий размерность энтропии S, Дж/К:
x8 = cmM.
Таким образом, термический потенциал системы определяется расчетной температурой сушильного агента, которая ограничена условиями нормальной (взрывобезопасной) эксплуатации пылесистемы.
Так как теплоемкость твердого топлива величина достаточно устойчивая к изменениям (Ст = const), то термическое состояние системы при заданном потенциале полностью определяется ее массовой загрузкой. Иными словами, обеспечив с помощью регулятора подачи топлива массовую загрузку мельницы постоянной, мы стабилизируем энтропию системы.
Обсуждение результата. Анализ координат пылесистемы как ТОУ позволяет выделить четыре группы параметров (табл. 2). Первая группа параметров определяется массовой загрузкой М мельницы топливом и характеризует общее состояние установки. Вторая группа определяется перепадом давлений на установке Р, работающей под наддувом, и характеризует нарушение устойчивости массопотока. Третья группа параметров определяется вектором координаты переносимой субстанции (мельничный продукт, готовое топливо и возврат мельничного продукта на домол). Эта группа параметров может служить теоретической основой для формирования комплексиро-ванных сигналов для решения информационно-управляющих задач. Четвертая группа параметров определяется дефектом массы металла ротора Мх, что связано с износом бил, и может служить для решения информационнодиагностических задач. Обобщенные координаты характеризуют режим работы установки и представляют собой показатели, которые в условиях АСУТП на базе ПТК могут быть сформированы в виде сигналов-параметров и в виде сложных комплексированных сигналов. Эти показатели могут служить теоретическим обоснованием параметров переопределенного объекта рассматриваемого класса и их комплексирования как способа определения неконтролируемых технологических параметров и совершенствования АСУТП электростанций [13].
o
Таблица 2. Сводная таблица обобщенных термодинамических потенциалов и координат системы пылеприготовления
Вид работы А Потенциал Координата Примечание
X разм. Хі разм.
Работа перемещения потока мельничного продукта Аі 1 V 2 їй м § ^ Н в Мвс Б3вс В3 = с х = 3 с 3 V 1 /.й
Работа перемещения потока готового топлива А2 1 V 2 їй м ^0 їй Н
Работа аэродинамических сил Аз а) V V м Р Па а,б - варианты представления
б) їй м Рій Н
Работа гравитационных сил А4 V їй м Мд Н
Работа вращения ротора А5 1 2 а) ^ иб 1 иб б ) 2 д в) й2к2пр д м2 С2 м м Мх Мхд Ох кг Н Н а,б,в - варианты представления
Работа дисперсионного превращения в камере размола А6 12 а) 2 иб м2 с2 м В2 V Оо кг Н а,б - варианты представления
1 иб б) —-) 2 д о2 vg Оо
А'в 1 2 а) 2 иб 1 и- б ) 6 ) 2 д м2 с2 м М о кг Н а,б - варианты представления
Работа дисперсионного превращения в сепараторе А7 1 V а ) 2 /.й 1 V б ) д 2 їй* м м2 с2 о 0возв Н Н Овозв~Мвозв' д
Термическая работа Ав т К СтМ Дж/К
В частности, для рассматриваемого класса ТОУ обобщенными термодинамическими координатами, характеризующими состояние пылесистемы как термодинамической системы, служат режимные параметры пылесистемы: количество топлива, циркулирующего в системе мельница-сепаратор; гидродинамический перепад давлений на установке и некая новая обобщенная координата переносимой субстанции, характеризующая кинетическую энергию потока мельничного продукта.
Обобщенными потенциалами являются расчетно-конструктивные показатели, которые характеризуют типоразмер мельницы и могут служить основой оценки эффективности топливоприготовительного оборудования и его конструктивного совершенствования.
Таким образом, применение обобщенного термодинамического анализа позволяет дать строгое решение задачи определения управляемых координат сложного объекта, в том
числе оценивать неконтролируемые технологические параметры
Список литературы
1. Автоматизация настройки систем управления / В.Я. Ротач, В.Ф. Кузищин, А.С. Клюев и др.; под ред. В.Я. Ро-тача. - М.: Энергоатомиздат, 1984.
2. Разработка и внедрение систем регулирования основных параметров котла в составе АСУТП энергоблока 500 МВт Рефтинской ГРЭС / В.А. Биленко, Н.Н. Деркач, Э.Э. Микушевич, Д.Ю. Никольский // Теплоэнергетика. -1999. - № 10. - С. 2-9.
3. АСУТП энергоблока ПГУ-450 Калининградской ТЭЦ-2 / А.А. Умрихин, А.С. Ладохин, Ю.С. Тверской и др. // Технология АСУТП электростанций / под ред. Ю.С. Тверского. -Иваново, 2005. - С. 57-60.
4. Тверской Ю.С. Автоматизация котлов с пылесисте-мами прямого вдувания. - М.: Энергоатомиздат, 1996.
5. Агафонова Н.А., Таламанов С.А., Твеpской Ю.С. Анализ промышленных методик идентификации на основе критерия минимума дисперсии частотных характеристик // Автоматика и телемеханика. - 1998. - № 6. - С.117-129.
6. Расчет и проектирование пылеприготовительных установок котельных агрегатов (нормативные материалы). - Л.: ЦКТИ, 1971. - 309 с.
7. Осокин В.П. Молотковые мельницы. - М.: Энергия, 1980.
8. Левит Г.Т. Пылеприготовление на тепловых электростанциях. - М.: Энергоатомиздат, 1991.
9. Тверской Д.Ю. Особенности построения нелинейной динамической модели пылесистем прямого вдувания котлов // Вестник ИГЭУ. - 2008. - Вып. 3. - С. 83-89.
10. Вейник А.И. Термодинамика. - Минск: Высш. шк., 1965.
11. Вейник А.И. Новая система термодинамики обратимых и необратимых процессов. - Минск: Высш. шк., 1966.
12. Бурдаков В.П., Данилов Ю.И. Физические проблемы космической тяговой энергетики. - М.: Атомиздат, 1969.
13. Тверской Д.Ю., Тверской Ю.С. Задачи и проблемы совершенствования АСУТП энергоблоков в направлении их интеллектуализации / Управление и информационные технологии: сб. докл. IV Всерос. науч. конф., 10-12 октября 2006 г. - СПб., 2006. - С. 230-236.
Тверской Дмитрий Юрьевич,
ООО «Научно-производственное предприятие ИНОТЭКС», кандидат технических наук, исполнительный директор, телефон (493б) бб-97-57, e-mail: kafsu@su.ispu.ru