CC BY
744
СИСТЕМА АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ ТЕМПЕРАТУРЫ ПЕРЕГРЕТОГО ПАРА БАРАБАННОГО КОТЛА.
Журавлев А.А., Шит М.Л., Попонова О.Б., Шит Б.М., Зубатый А.Л.(Институт
энергетики АНМ)
AUTOMATIC CONTROL SYSTEM OF THE DRUM BOILER SUPERHEATED
STEAM TEMPERATURE.
The control system of the temperature of the superheated steam of the drum boiler is examined. Main features of the system are the PI-controller in the external control loop and introduction of the functional component of the error signal of the external control loop with the negative feedback of the error signal between the prescribed value of steam flowrate and the signal of the steam flowrate in the exit of the boiler in the internal control loop.
Key words: steam boilers, automatic control of superheaters.
SISTEM DE DIRIJARE A TEMPERATURII ABURULUI SUPRAANCALZIT
CAZANULUI CU TAMBUR.
Este examinat sistem de tip cascad cu PI-regulator contur extern cu introducerea componentei functionale a semnalului de eroare intre contur extern §i bucla negativa a semnalului de eroare intre sarcina de debit a aburului §i semnalului de debit de abur la ie§ire din cazan in contur intern.
Cuvinte cheie: cazane de abur cu tambur, dirijarea supraancalzitoarelor.
СИСТЕМА АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ ТЕМПЕРАТУРЫ ПЕРЕГРЕТОГО ПАРА БАРАБАННОГО КОТЛА
Рассмотрена каскадная система с ПИ-регулятором во внешнем контуре, с введением функциональной составляющей сигнала рассогласования внешнего контура и отрицательной обратной связью по сигналу рассогласования между заданием по расходу пара и сигналом расхода пара на выходе котла во внутреннем контуре.
Ключевые слова: паровые котлы, автоматическое управление пароперегревателями.
1. Объект управления. Рассматривается пароперегреватель котла в виде двух теплообменников типа «газ-пар» конвективного типа (рис.1). В необогреваемую рассечку подведен трубопровод для впрыска собственного конденсата котла, представляющий собой пароохладитель. В пароохладителе происходит смешивание пара и собственного конденсата и охлаждение пара. Для контроля температуры пара установлено две термопары: после пароохладителя и после пароперегревателя. Схематично, пароперегреватель выполнен следующим образом. В верхнюю секцию пароперегревателя поступает пар из котла. После этой секции идет необогреваемый участок, где происходит смешивание собственного конденсата котла и пара и нагрев пара газами, которые уже нагрели вторую по ходу пара секцию. Далее пар поступает во вторую секцию, связанную с общим паропроводом, расположенную ниже первой секции, где происходит начальный нагрев пара. Эта секция обогревается дымовыми газами в первую очередь. За этой секцией установлена термопара, измеряющая температуру перегретого пара (см. рис.2).
Поэтому, когда изменяется температура, расход дымовых газов, расход пара через пароперегреватель, происходит изменение как температуры в рассечке, так и температуры на выходе ПП. Причем, коэффициенты усиления по каналам возмущения температуры в рассечке и температуры после 1111 различны, что объясняется тем, что
на температуру в рассечке действует поток пара, охлажденный в первой по ходу дымовых газов секции ПП.
Секция 2 пароперегревателя является регулируемой. Измеряют температуры пара на входе и выходе из секции, (иногда, расход пара), и по этим параметрам ведут процесс.
Из этого рисунка видно, что при изменении расхода пара, температуры и расхода дымовых газов изменяются синхронно как температура пара в месте впрыска конденсата, так и температура пара после пароперегревателя (ПП).
Рис. 1. Схема пароперегревателя.
2. Математическая модель пароперегревателя как объекта управления температурой пара.
Объектом регулирования является участок котла между точкой ввода регулирующего воздействия и точкой контроля регулируемой температуры, включающей в себя радиационные, полурадиационные и конвективные поверхности нагрева. На температуру пара оказывают влияние большое количество факторов: из них наиболее существенными оказываются: нагрузка котла, загрязнение поверхностей нагрева, ошлаковывание топки, температура питательной воды, избытки воздуха.
При снижении нагрузки удельное тепловосприятие радиационных поверхностей нагрева повышается, а конвективных - уменьшается. Кроме того, на температуру перегрева оказывают большое влияние такие возмущения, как случайные возмущения подачи топлива при заданной нагрузке, изменение тягодутьевого режима и т.п.
В последнее время САР ТПП проектируются с вводом дополнительного воздействия по расходу пара. Встречаются САР с различного рода связями типа «люфт» между отдельными схемами регулирования расхода пара на участках парового тракта котла, ПИ и ПИД - регуляторы, САР ТПП с корректирующим и стабилизирующим регуляторами, которые, в некоторых случаях, могут иметь значительные преимущества перед типовыми с исчезающим импульсом из промежуточной точки. Для регулирования впрыска широко применяются шиберные клапаны.
Для исследуемой системы характерно большое транспортное запаздывание в объекте управления, нестационарность и априорная неопределенность математической модели ПП, что приводит, в конечном счете, к ручной перенастройке коэффициентов регулятора и необходимости проектировать систему из условий робастности.
Передаточная функция объекта регулирования включает в себя:
WTH - передаточную функцию участка трубопровода от места впрыска до места установки термопары (передаточная функция опережающего участка);
WT — передаточную функцию 1111 между местами установки термопар 1 и 2, ее называют передаточной функцией инерционного участка;
WmWT - передаточную функцию главного участка;
WD - передаточную функцию 1111 при возмущении расходом пара;
Wp — передаточную функцию 1111 при топочных возмущениях по теплоте дымовых газов.
Wm(p) = К
vpr
ТзР + 1. Т4р + Ґ
(1)
WT (р) = кте
к Р+
Ткр
Тр+1
(2)
0,2Тр +1
WD(p) = -kL
1 - ехр-
т0р +
Ткр Тр + 1
Тр + 1
т0р +
Ткр Тр + 1
1-е
t0p+
Ткр
Тр+1
(3)
Передаточная функция Жр — по теплоте газов, уходящих из топки, имеет вид
Wp=-
1 -ехр-
т0р +
Ткр Тр + 1
Тр + 1
т0р +
Ткр Тр + 1
1-е
t0p +
Ткр
Тр+1
(4)
Характерные параметры участка пароперегревателя, состоящего из одного пакета, определяют по формулам
ос F
£ _ ^IC1 INC .
DCST ’
Г = 3600 GmCm ; сс F
^1C1 INC
w
(5)
(6) (7)
где а1С — условный коэффициент теплоотдачи от стенки змеевика к пару на внутренней поверхности змеевиков, кДж 1(м2 • ч • °С);
С5Т - средняя изобарная теплоемкость пара, проходящего через пакет,
кДж/(кг ■ °С);
1'].., 1’оит - соответственно внутренняя и наружная поверхности змеевиков пакета,
омываемые, соответственно, паром и дымовыми газами, м2;
См - масса металла змеевиков пакета, кг;
/ - длина пути пара от входа до выхода из обогреваемой поверхности нагрева, м; м> — средняя скорость пара, м/с;
- - ^ • (8)
а2 - коэффициент теплоотдачи от стенки к пару на внутренней поверхности змеевиков, кДж/(м2 - ч- °С);
Лм - коэффициент теплопроводности металла змеевиков при его средней температуре, кДж/(м ■ ч ■ °С);
^ш^оит - внутренний и внешний диаметры змеевика, м.
Коэффициент усиления по впрыску, °С /т-ч1;
и -/
к..=
vpr
1 vpr
D С ■
D0CSTIN ;
Сшы - средняя изобарная теплоемкость пара на входе в пакет, кДж/(кг ■ С) ;
/, - энтальпия пара перед пакетом за местом впрыска, кДж / кг ;
1ург - энтальпия впрыскиваемой воды (впрыскиваемого конденсата), кДж / кг .
Коэффициент усиления по расходу пара, °С / т • ч-1 для участков пароперегревателя с конвективным теплообменом рассчитывается по формуле:
, 1 + 0,2 т
кп =
кт
( кт \
\-е 1+т
V У
D С
О STOUT
где,
С’stout ~ средняя изобарная теплоемкость пара на выходе из пакета, кДж/(кг ■ °С) ;
surf ~ приращение энтальпии в пароперегревателе при стационарном состоянии до
возмущения, кДж / кг .
Коэффициент усиления по температуре пара, 0С /°С, для участков с конвективным теплообменом
s-i кт
кт=±шш-е~^п. ^О)
С STOUT
ОС F
m =-+J2L- (11)
oc2F1N
ах — наружный коэффициент теплоотдачи от дымовых газов к стенке змеевика,
кДж/(м2 - ч- °С).
Коэффициент усиления при топочных возмущениях для пароперегревателей с конвективным теплообменом при возмущении температурой газов:
s-f | ( кт \
_ ^STIN STOUT Л с1 ---
2 С
stout V. У
1-е
1+т
(12)
Выражения (1) - (12) заимствованы из [13].
Упростим выражение: А = ехр-
А Ткр л Тр + 1
для чего воспользуемся методикой [15].
-Ткр Тр + \
= -к
1-
V
Ткр
Тр + \
-к+-
= -к + -
Тр + 1
W(p) = е Tp+l = е Tp+l = ekeTp+l =
k\ekem
-рч
(13)
(Тр + \)(Т2р + \) (Тр + \)(Т2р + \У к\ = ек ',к2 = к\-е~к =1;
Т2 = Ту/2к-1, тг = т\ к-1 -72Та
Запишем упрощенное выражение для передаточной функции Жв (Р), предварительно введем следующее обозначение:
т — т 4-f
WD{p) = Где,
t Т
Т — 0
кт - Т22- 0,2Тр +1 ТХМТ22-Р1 +Т2МТ22-р + \
(TdP +1)
{Тр + \){Т2р + \) 1 + 0,5рг0
+ Тк
-к,
2 (tJ + ТкУ
Т\\ = Т-Т2-т0;
Т2\ = т0 ■ Т + Т2 +2-Г2-Г; Т22 = т0+2-(Т+ Т2) + к2т0 =2(Т+ Т2 + т0)
3. Синтез системы управления пароперегревателем.
На рис.2. приведена функциональная схема управления пароперегревателем. На этой схеме приняты следующие условные обозначения: ТЕ—термопары (установлены после пароохладителя и после пароперегревателя); FE- расходомеры (измеряют расход пара и собственного конденсата на впрыск); SPH1, SPH2 — первая и вторая секции пароперегревателя по ходу пара; DSH - пароохладитель; COL - коллектор.
к
k
а
Рис.2. Функциональная схема пароперегревателя как объекта управления.
На рис.3 приведена структурная схема пароперегревателя как объекта управления в системе регулирования температуры.
На этом рисунке приняты следующие обозначения:
Жс (р) - передаточная функция основного регулятора САР, р — оператор дифференцирования,
Ж2-передаточная функция компенсирующего звена расхода пара,
ЖЗ - передаточная функция компенсирующего звена температуры дымовых газов, ЖГ1,ЖГ2 ~ передаточные функции термопар,
^а'2 — передаточные функции звеньев, компенсирующих инерционность термопар, 1¥8_ты - передаточная функция, связывающая расход топлива и температуру пара после №>:_Тс - передаточная функция, связывающая расход топлива и температуру пара после пароперегревателя,
Ж,,Л,2 - передаточная функция, связывающая расход газа и расход пара в канале возмущения на температуры пара после пароохладителя,
ЖР81 - передаточная функция, связывающая расход газа и расход пара в канале возмущения на температуры пара после пароперегревателя,
Ж/_7(н - передаточная функция, связывающая температуру дымовых газов и температуру пара после пароохладителя,
И/т_Гс - передаточная функция, связывающая температуру дымовых газов и температуру пара после пароперегревателя,
Т5Т-Р¥ - задание температуры пара после пароперегревателя, пароохладителя,
Рис.3. Структурная схема системы управления температурой перегретого пара.
Ор - приращение расхода топлива, С5 - приращение расхода пара, Т8 - температура дымовых газов.
Как видно из рассмотрения структурной схемы системы управления пароперегревателем, инерционные звенья, включенные в цепи сигналов главной обратной связи и сигнала опережения по производной сигнала из промежуточной точки, существенно замедляют реакцию системы. Поэтому включим последовательно с этими звеньями, звенья, компенсирующие инерционность термопар.
Передаточная функция компенсирующего звена имеет вид:
Р^К2(КЗ) ~ кт (7г2(Т3)Р ■*” 1) (14)
В этом уравнении постоянные времени зависят от ряда факторов, определяемых как конструкцией узлов измерения температур, так и скоростью и температурой и давлением перегретого пара. Определение этих постоянных времени представляет собой решение классической задачи параметрической идентификации и не входит в задачу данного исследования.
4. Моделирование системы управления пароперегревателем при детерминированных воздействиях.
На рис.4. приведен результат моделирования системы при задании по температуре пара относительно базового значения -300С и возмущении по расходу топлива равному 10% от базового значения и возмущении коэффициентом избытка воздуха.
I FUEL FLOW RATE
g_____________I___________I____________i___________I___________I____________I___________I____________
0 1000 2000 3000 4000 5000 6000 7000 8000
0 1000 2000 3000 4000 5000 6000 7000 8000
50----------------1--------------1-------------1--------------1 , ■ ■ ■ ------,
| STEAM TEMPERATURE ERROR |
I I I i
0 1000 2000 3000 4000 5000 6000 7000 8000
I i i
0 1000 2000 3000 4000 5000 6000 7000 8000
STEAM VALVE STROKE
0 1000 2000 3000 4000 5000 6000 7000 8000
Рис.4. Графики работы модернизированной системы при положительном знаке изменения расхода пара и температуры дымовых газов.
На верхнем графике приведена зависимость расхода топлива от времени (так называемый, единичный скачок). На втором сверху графике приведено изменение температуры дымовых газов при увеличении коэффициента избытка воздуха по
отношению к базовому значению. На третьем сверху графике приведена зависимость ошибки регулирования от времени. На четвертом сверху графике зависимость
выходной температуры от времени. На нижнем графике приведена зависимость управления от времени. Время переходного процесса составляет для этого варианта 445,2 с., а значение интегрального квадратичного критерия 1.1874e+004.
На рис. 5. приведены графики работы для системы без коррекции при тех же начальных условиях и возмущениях, что в первом примере.
Рис. 5. Графики работы типовой системы при ненулевых начальных условиях и постоянно действующих возмущениях.
Для этого варианта расчета время переходного процесса составляет 475,35 с., что больше времени переходного процесса для скорректированной системы на 7%, а значение критерия качества составляет 1.2802e+004, что больше чем у модернизированной системы на 7,8%. Рассмотрим вариант расчета, когда начальные условия отрицательны.
Рис.6. Графики работы модернизированной системы при положительном знаке изменения расхода пара и температуры дымовых газов и отрицательном задании по
температуре пара.
Для этого варианта расчета время переходного процесса составляет 639,5 с., значение квадратичного критерия качества составляет 1.7075е+004.
Рис.7. Графики работы типовой системы при отрицательном знаке изменения расхода пара и температуры дымовых газов.
Для этого варианта время переходного процесса составляет 1187,1 с, а значение квадратичного критерии я качества - 4.1230е+004, что соответственно на 85,6% и 241,5%: больше, чем у модернизированной системы. Рассмотрим работу
модернизированной системы управления при тех же возмущающих воздействиях по расходу пара и температуре дымовых газах, что и в предыдущих примерах, но при нулевых начальных условиях. На рис. 8 приведены графики работы для модернизированной системы, а на рис.9 графики работы системы для типовой системы.
FUEL FLOW RATE
FUEL GAS TEMPERATURE
• | STEAM TEMPERATURE ERROR |
i i i i i
STE :AM OUTPUT TEMPERATURE |
STEAM VALVE STROKE |
1000 2000 3000 4000 5000 6000 7000
Рис.8. Графики работы модернизированной системы при положительном знаке изменения расхода пара и температуры дымовых газов и нулевых начальных условиях
по температуре пара.
Время переходного процесса (время достижения отклонением пара от заданного (базового) значения температуры равного 20С) составляет для этого случая 568 с., а значение квадратичного критерия качества (интеграла квадрата рассогласования) равно 1.5158е+004. Рассмотрим тот же вариант расчета, но для типовой системы (рис. 9).
Время переходного процесса (время достижения отклонением пара от заданного (базового) значения температуры равного 20С) составляет для этого случая 983,4 с., а значение квадратичного критерия качества (интеграла квадрата рассогласования) 2.9140е+004.
Как следует из изложенного, время переходного процесса для модернизированной системы сократилось на 73%, а значение квадратичного критерия качества сократилось на 92%. Как известно [20], снижение температуры пара на 100С эквивалентно перерасходу топлива на 0,2%. Так как происходит уменьшение величины критерия качества в модернизированной системе, по сравнению с типовой, следовательно, при сбросах нагрузки происходит экономия тепловой энергии (за счет снижения отрезка времени, в течение которого снижается температура пара). При набросах нагрузки за счет снижения времени превышения температурой пара заданного значения повышается долговечность пароперегревателя.
FUEL FLOW RATE
1000
2000
3000
4000
5000
6000
7000
8000
f
FUEL GAS TEMPERATURE -
8000
8000
! | STEAM OUTPUT TEMPERATURE |
i 1 i i i
1000
2000
3000
4000
5000
6000
7000
8000
I STEAM VALVE STROKE
-
1000
2000
3000
4000
5000
6000
7000
8000
Рис.9. Графики работы типовой системы при положительном знаке изменения расхода пара и температуры дымовых газов и нулевых начальных условиях по
температуре пара.
Рассмотрим влияние на температуру пара изменения температуры дымовых газов, возникающей при работе подсистемы регулирования «топливо- воздух».
Обозначим:
(I- расход природного газа на котел;
Ол - расход воздуха на котел;
ке - отношение массового расхода воздуха к расходу газа, необходимое для обеспечения стехиометрического режима сгорания газа;
Та,ТА — температуры газа и воздуха, поступающие на горение, соответственно;
са,сА - изобарные теплоемкости газа и воздуха, соответственно;
Ва - удельная теплота сгорания природного газа.
Запишем уравнение материального баланса между теплом, уходящим с дымовыми газами и теплом, приходящим с воздухом и газом, поступающим на горение.
&оТаса + ^АТАсА + г/0оВо = (Оа + Ол )Тарсар (15)
Пренебрежем изменением КПД котла при изменении его нагрузки, а также изменением температур воздуха и газа, поступающими на горение, присосами воздуха по тракту котла, потерями тепла с химическим и механическим недожогом, и определим зависимость температуры газов в зависимости от изменения соотношения
воздух - газ.
Обозначим коэффициент избытка воздуха - СС .
Тогда
^о^аса с’паК1'лс. , + пОаВа (Сг0 + акеСа)Тарсар
(16)
Сократим выражение (16) на . В результате получим:
Таса + акеТАсА + гіВа - (1 + аке)Торсор или Таса + акеТАсА + г]Ва _
т = ° °
GF
(17)
(1 + а*е)се
Определим зависимость между приращением температуры дымовых газов и коэффициентом избытка воздуха, для чего продифференцируем второе выражение в (17) по ОС.
В результате получим:
АТОР =
ґ5Та^ V да ;0
_ Таса (1 + а?оке)сар сар Таса + а0кеТАсА + г/Ва _ £Дег
(1 + «оке)сор 2 ’ (18)
Индекс «О» при означает, что мы рассматриваем динамику процессов в относительной близости от некоторого установившегося режима работы котла, принятого за базовый.
Так как процессы в следящих системах «топливо воздух», в основном, имеют линейный характер, то оператор, связывающий изменение рассогласования коэффициента избытка воздуха и температуру дымовых газов Тра, является линейным, т.е.,
Щти а=^^ = —-—; (19)
ТРО~а Аа(р) Тоар +1 К '
Величина Т определяется инерционностью топки.
Так как канал, по которому движутся дымовые газы, имеет сложную форму, присутствуют пульсации давления в газовоздушном тракте, то измерение температуры газов в одной точке, является непредставительным, а измерение температуры в нескольких точках с последующим осреднением может потребовать большого количества точек измерения и дорого. В связи с этим, вычисление коэффициента избытка воздуха и последующее вычисление отклонения температуры газов от базового значения и введение соответствующего корректирующего воздействия на впрыск является рациональным способом компенсации влияния температуры дымовых газов на температуру пара.
Составим условие полной инвариантности выходного сигнала относительно температуры дымовых газов:
^ТР0-а^Т-Те = (20)
Из этого выражения можно вычислить Ж3 .
W2 - WTFG-aWT-Te _ kK^ T4P + l к2(Тр + 1)(Т2р + 1)(2 + рТ0)
WTHWSH kvpr TDGp + \ TPG_Tep +1 T3p +1 (2-рт0У (21)
Выражение (21) является сложно реализуемым, поэтому ограничимся следующим представлением W3.
W3 « £3, (22)
где ^ 3 является функцией расхода пара и конструктивных параметров пароперегревателя и должно определяться в результате идентификационных экспериментов.
Корректирующие связи по расходу пара и температуре дымовых газов обеспечивают снижение времени переходного процесса регулирования температуры пара и повышение качества переходного процесса, величина интегрального квадратического критерия снижается от 7 до 37%).
ЗАКЛЮЧЕНИЕ
Разработан новый закон регулирования температуры перегретого пара на базе которого синтезирована система управления, позволяющая повысить точность поддержания температуры пара на выходе пароперегревателя до 2 гр.С (при норме до 4,5 гр.С), сократить до 2 раз длительность переходных процессов, повысить качество
управления котлом, что позволит снизить расход топлива при работе котла за счет
уменьшения периодов времени при которых происходит снижение температуры пара ниже допустимого значения.
Разработан новый метод компенсации инерционности сигналов термопар и инерционности канала влияния расхода пара на температуру пара после пароперегревателя (функциональной зависимости температуры пара на выходе пароперегревателя от расхода пара). Использование методов компенсации:
инерционности сигналов термопар в пароперегревателе; инерционности сигнала расхода пара, позволяет сократить до 2 раз длительность переходных процессов.
ЛИТЕРАТУРА
1. Сидельковский Л.Н., Юренев В.Н. Котельные установки промышленных предприятий: Учебник для вузов. - 3-е изд. - М.: Энергоатомиздат, 1988.-528с.
2. Субботина Н.П. Водный режим и химический контроль на ТЭС. М.:
Энергоатомиздат, 1985. - 312 с.
3. Energy Efficiency Handbook /Alliance to Save Energy, Council of Industrial Boiler Operators, U.S. DOE Office of Industrial Technologies. - 1998. - 64 p.
4. Пухов Г.Е., Жук К.Д. Синтез многосвязных систем управления по методу обратных операторов. «Наукова Думка», Киев, 1966.
5. Бессекерский В.А. Цифровые автоматические системы. Наука, М., 1976, с.576.
6. Давыдов Н.И., Тюпина Т.Г. Исследование системы регулирования температуры пара с двумя опережающими скоростными сигналами. Теплоэнергетика, №10, 2002, с. 17-21.
7. Холщев В.В. Регулирование температуры перегретого пара на барабанном котле. Теплоэнергетика, №8, 2003, с.39-43.
8. Пикина Г.А., Чикунова А.М. Математические модели противоточного конвективного пароперегревателя котла в системе регулирования температуры. Теплоэнергетика, №8, 2002, с.25-33.
9. Пикина Г.А. Сравнительный анализ линейных моделей противоточного конвективного пароперегревателя котла в системе регулирования температуры. Теплоэнергетика, №10, 2002, с.22-25.
10. Ефимов Н.Н. Основные закономерности изменения количества аккумулированного тепла в пароперегревателях. Теплоэнергетика, №1, 1999, с.36-39.
11. Sang Hyuk Lee, Jaesop Kong, Jin H.Seo Observers for bilinear systems with unknown inputs and application to superheater temperature control. Control Engineering Practice.1997, Vol.5, N4, pp.493-506.
12. Sima W. Numerical modeling of steam superheaters. Energy 26(2001), pp. 1175-1184.
13. Клюев А.С. Наладка систем автоматического регулирования барабанных котлов. - М.: Энергоатомиздат, 1985.-280с.
14. Клюев А.С., Товарнов А.Т. Наладка систем автоматического регулирования котлоагрегатов. М; «Энергия», 1970, 280 с.
15. Шевяков А.А., Яковлева Р.В. Инженерные методы расчета динамики теплообменных аппаратов. Машиностроение, М., 1968., 320 с.
16. Давыдов Н.И. Исследование системы регулирования температуры пара с двумя опережающими скоростными сигналами. Теплоэнергетика, №10, 2002, с. 17-21.
