CC BY
16
Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. Т. 61, № 5 (2018), с. 451-462 Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. V. 61, No 5 (2018), pp. 451-462 451
https://doi.org/10.21122/1029-7448-2018-61-5-451-462 УДК 62.5
Параметрическая оптимизация системы автоматического управления мощностью энергоблоков 300 МВт в режиме постоянного давления пара перед турбиной
Г. Т. Кулаков1*, А. Т. Кулаков1*, К. И. Артёменко1*
'^Белорусский национальный технический университет (г. Минск, Республика Беларусь)
© Белорусский национальный технический университет, 2018 Belarasian National Technical University, 2018
Реферат. Рассмотрен вопрос параметрической оптимизации систем автоматического управления мощностью энергоблоков (САУМБ) 300 МВт Лукомльской ГРЭС в режиме постоянного давления перегретого пара перед турбиной. В течение 1974-1979 гг. на восьми энергоблоках Лукомльской ГРЭС были внедрены САУМБ с ведущим котельным регулятором мощности, которые будут вынуждены работать в широком диапазоне изменения нагрузок в связи с планируемым вводом двух энергоблоков Белорусской АЭС: первого - в 2019 г., второго - в 2020 г. Суммарная мощность Белорусской АЭС составит 2400 МВт. Кроме того, ужесточились современные требования к регулированию частоты и перетоков активной электрической мощности в энергосистеме: время достижения половинного значения необходимого изменения мощности должно составлять 10 с в пределах нормального и аварийного резервов; время достижения полного значения необходимого изменения мощности должно составлять 30 с в пределах нормального резерва и 2 мин - в пределах аварийного. В связи с этим актуальной становится задача повышения эффективности работы энергоблоков, работающих в переменной части графика электрических нагрузок, за счет использования современных методов расчета параметров динамической настройки регулирующих устройств САУМБ. Приведена методика параметрической оптимизации типовой САУМБ, позволяющая улучшить качество регулирования мощности и давления перегретого пара перед турбиной. Описанная методика иллюстрируется результатами компьютерного моделирования переходных процессов в системе при отработке задания и внутренних возмущений, которые подтверждают правильность предлагаемой методики по сравнению с известными методами оптимизации типовых САУМБ.
Ключевые слова: автоматическое управление, электрическая мощность, режим постоянного давления пара, параметрическая оптимизация системы
Для цитирования: Кулаков, Г. Т. Параметрическая оптимизация системы автоматического управления мощностью энергоблоков 300 МВт в режиме постоянного давления пара перед турбиной / Г. Т. Кулаков, А. Т. Кулаков, К. И. Артёменко // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2018. Т. 61, № 5. С. 451-462. https://doi.org/10.21122/ 1029-7448-2018-61-5-451-462
Адрес для переписки
Кулаков Геннадий Тихонович Белорусский национальный технический просп. Независимости, 65/2, 220013, г. Минск, Республика Беларусь Тел.: +375 17 293-91-45 tes_bntu@tut.by
Address for correspondence
Kulakov Gennady T.
Belarusian National Technical University
65/2 Nezavisimosty Ave.,
220013, Minsk, Republic of Belarus
Tel.: +375 17 293-91-45
tes_bntu@tut.by
Parametric Optimization for Automatic Control System of 300 MW Power Units at a Constant Steam Pressure Upstream of the Turbine
G. T. Kulakov1*, A. T. Kulakov1*, K. I. Artsiomenka1*
1)Belаrusian National Technical University (Minsk, Republic of Belarus)
Abstract. The problem of parametric optimization of automatic control systems of power units (ACSPU) of 300 MW at the Lukoml'skaya GRES (Lukoml Condensing Power Plant) at a constant superheated steam pressure upstream of the turbine is under consideration. During 1974-1979, eight units of the Lukoml'skaya GRES implemented ACSPU with a leading boiler power control, which will be forced to work in a wide range of loads due to the planned commissioning of two power units of the Belarusian NPP: the first one - in 2019, and the second one - in 2020. The total capacity of the Belarusian Nuclear Power Plant will amount to 2400 MW. In addition, modern requirements for the regulation of the frequency and flows of active electric power in the power system have been tightened: the time to reach the half value of the required power alteration should be 10 seconds within the normal and emergency reserves; the time to reach the full value of the required power alteration should be 30 seconds within the range of the normal reserve and 2 minutes - within the range of the emergency reserve. With this regard, increasing the efficiency of the units operating in the variable part of the electric load schеdule by the use of modern methods of calculation of parameters of the dynamic adjustment of the regulating devices of ACSPU becomes an urgent problem. The methodology of parametric optimization of the typical ACSPU that make it possible to improve the quality of power regulation and the pressure of the superheated steam upstream of the turbine is presented. The described technique is illustrated by the results of computer simulation of transient processes in the system when the task and internal disturbances are worked out, which confirm the correctness of the proposed technique in comparison with the known methods of optimization of the typical ACSPU.
Keywords: automatic control, electric power, constant steam pressure regime, system's parametric optimization
For citation: Kulakov G. T., Kulakov A. T., Artsiomenka K. I. (2018) Parametric Optimization for Automatic Control System of 300 MW Power Units at a Constant Steam Pressure Upstream of the Turbine. Епе^еИка. Proс. GS Higher Educ. Inst. аnd Power Eng. Assoc. 61 (5) 451-462. https://doi.org/10.21122/1029-7448-2018-61-5-451-462 (in Russian)
Введение
Вопросам автоматического регулирования перетоков активной мощности и частоты энергообъединения посвящено большое количество публикаций [1-14], исследующих системы автоматического управления мощностью энергоблоков [1, 4, 7-11], системы общестанционной части автоматического управления мощностью станции [5], способы управления энергоблоком в аварийных режимах энергосистемы [6], алгоритмы работы систем автоматического регулирования перетоков мощности энергосистемы [2, 3], анализ динамических характеристик энергоблоков в широком диапазоне изменения нагрузок энергоблоков [9, 11, 12].
В настоящее время существенно ужесточились требования к качеству поддержания мощности и частоты энергообъединения: время достижения половинного значения необходимого изменения мощности должно составлять 10 с в пределах нормального и аварийного резервов; время вхождения в зону ±1 % номинальной мощности энергоблока Ыном требуемого значения нагрузки: 30 с - в пределах нормального; 2 мин - в пределах аварийного резерва [12].
Вместе с тем повышение качества поддержания требуемой мощности энергоблоков обычно достигается увеличением степени форсировки котла в 1,5 раза, для дубль-блока мощностью 300 МВт с котлом ТГМП-114 и для моноблока с котлом ТГМП-314 она возрастает до 1,6 раза [12]. Решение этого противоречия приводит к существенному усложнению схем управления и методов их оптимизации в широком диапазоне изменения нагрузок энергоблока. Особенно актуальна эта проблема для Белорусской энергосистемы в связи с планируемым вводом двух энергоблоков Белорусской АЭС суммарной мощностью 2400 МВт, которые будут работать в базовой части графика электрических нагрузок энергосистемы.
Особенности структурно-параметрической оптимизации предлагаемых систем управления мощностью энергоблоков
Структурная схема моделирования переходных процессов типовой системы автоматического управления мощностью энергоблоков с форсированием сигнала задания по мощности и ведущим котельным регулятором при постоянном давлении пара перед турбиной приведена на рис. 1.
Рис. 1. Структурная схема моделирования переходных процессов типовой системы автоматического управления мощностью энергоблоков с форсированием сигнала задания по мощности
Fig. 1. The block-diagram of a typical ACSPU transient simulation with power signal forcing
На рис. 1 использованы следующие обозначения: Ызд - заданное значение мощности энергоблока; Щ - фактическая электрическая мощность энергоблока; - внутреннее возмущение; хзд1 - задание регулятору топлива; а - доля ошибки регулирования по мощности с выхода измерительного блока котельного регулятора мощности (КРМ) на вход турбинного регулятора мощности (ТРМ); Вт - изменение расхода топлива; р0д, р0 - соответственно заданное и текущее значения давления перегретого пара перед турбиной; крк - перемещение регулирующих клапанов турбины.
Исходные данные для расчета параметров динамической настройки типовой САУМБ при моделировании переходных процессов энергоблока 300 МВт Лукомльской ГРЭС приведены в [14].
Передаточная функция опережающего участка по расходу топлива (газа)
К 1,2
Ж (р) = Жв , (р) =-25-=-, (1)
оп Вт ,Лро (Топ Р +1)(° оп Р +1) (7р +1)(0, 7р +1)
где Коп - коэффициент передачи; Нро - перемещение регулирующего органа подачи топлива; Топ, соп - большая и меньшая постоянные времени передаточной функции опережающего участка по расходу топлива (газа), с.
Передаточная функция по электрической мощности при скачкообразном изменении задания котельному регулятору нагрузки
КеТ1Р 1,6е-11,2 р
К в (Р) =-1-=-, (2)
Мф, Вт (Т1 р +1)^ р +1) (101р +1)(19 р +1)
где К1 - коэффициент передачи; Т1, ст1 - большая и меньшая постоянные времени передаточной функции по мощности при изменении расхода топлива, с; т1 - время запаздывания по каналу регулирующего воздействия, с.
Передаточная функция по давлению перегретого пара перед турбиной при скачкообразном изменении задания котельному регулятору нагрузки
К2в~11Р 1,6е"19р
ЖР в (Р) =-2-=-, (3)
Ро' Вт (Т2 р +1)(с 2 р +1) (56р +1)(14 р +1)
где К2 - коэффициент передачи; Т2, с2 - большая и меньшая постоянные времени передаточной функции по давлению пара перед турбиной при изменении расхода топлива, с; т2 - время запаздывания по каналу регулирующего воздействия, с.
Передаточная функция по электрической мощности при скачкообразном перемещении регулирующих клапанов турбины
То р (Т3 р +1) 169,8 р (2,3 р + 1)
Жмфн (р)=-=-, (4)
ф, рк (Т4р +1)(Т5 р +1) (168р + 1)(48р + 1)
где Т0, Т3, Т4, Т5 - постоянные времени передаточной функции по мощности при возмущении регулирующими клапанами турбины, с.
Передаточная функция по давлению перегретого пара перед турбиной при скачкообразном перемещении регулирующих клапанов турбины
К3(Г6 p +1) -0,25(4,6 р +1)
К л (р) =-=-' (5)
> Т7 р +1 90 р +1 (5)
где К3 - коэффициент передачи; Т6, Т7 - постоянные времени передаточной функции по давлению перегретого пара перед турбиной при возмущении регулирующими клапанами турбины, с.
Котельный регулятор мощности и регулятор топлива (РТ) представляют собой каскадную систему автоматического регулирования, в которой стабилизирующий ПИ-регулятор топлива настраивают на оптимальную отработку внутреннего возмущения, а корректирующий котельный ПИ-регулятор - на оптимальную отработку задающего сигнала.
Регулятор топлива (рис. 1, передаточная функция Кр1) настраивают по передаточной функции опережающего участка (1) по методу частичной компенсации [13]. Вначале рассчитывают относительное значение коэффициента передачи стабилизирующего регулятора К
К = КрКоп = 0,7395Т^1 + 1] -1, (6)
где Т - относительная постоянная времени передаточной функции опережающего участка (1), равная
Т
Т = -<*■; (7)
аоп
Кр1 - коэффициент передачи стабилизирующего регулятора. Затем с учетом (6) и (7) находят абсолютное значение Кр1
К
Кр1 = — • (8)
оп
Относительное значение времени интегрирования стабилизирующего регулятора I рассчитывают с учетом (6) и (7) по формуле
1 = ^ 6,3К . (9)
а°п т[1+^
Затем находят абсолютное значение времени интегрирования Ти1
Ти1 = 1 ° оп. (10)
Котельный ПИ-регулятор мощности с передаточной функцией Жр2 настраиваем по передаточной функции объекта (2) при возмущении расходом топлива, используя метод полной компенсации в частном виде [13].
Время интегрирования КРМ
Ти2 = Т1 + СТЬ (11)
а коэффициент передачи
Т +С1
Кр2 = '2 1 , (12)
р2 4Е,2 К1Т1
где Е - коэффициент демпфирования, равный единице, позволяющий устранить перерегулирование при отработке скачка задания.
Структуру ТРМ формируем на основе передаточной функции оптимального регулятора с использованием передаточной функции (5) по давлению пара перед турбиной при возмущении перемещением регулирующих клапанов турбины [14]
жР3 ( Р) = Жтрм (р) = [Жро, ИрК (р) ]-1 г/;1("=1) (р), (13)
где Жзрдс(п=1) (р) - передаточная функция разомкнутой системы в виде идеального интегрирующего звена,
Жзр;1("=1)(р) = (14)
Тзд1 р
Тзд1 - единственный расчетный параметр динамической настройки ТРМ.
С учетом передаточных функций (5), (13) и (14) передаточная функция ТРМ примет вид реального ПИ-регулятора
Т р +1
ЖТРМ (р) = (Т ,1УГ-■ (15)
К 3(Т6 р +1)Тзд1 Р
Численное значение Тзд1 рассчитывают по формуле
Тзд1 = 0,618Т7. (16)
Расчет параметров динамической настройки дифференциатора инвариантности осуществляют с учетом передаточных функций (4) и (5) из условия
Жинв ( р) = КдиТдир = Ж4(р) = Т0 р(Т3 р + 1)(Т7 р + 1)
д Тдир +1 Т3 (р) (Т4р + 1)(Т5р +1)К3 (Т6р +1),
откуда численное значение времени дифференцирования
ти = Т4 + Т5 + Тб - Т7 - Т3, (18)
а коэффициент передачи Ки принимают равным единице.
Передаточную функцию формирователя сигнала задания (ФСЗ) представляют в виде звена быстрого реагирования
ЖФСЗ (p) = 1 + Гд = 1 +
p (i+кдд )т; p+1 ^ p+i
тзд 11 Тд p +1
тзд I 1 Тд p +1
тзд I 1
T p +1
(19)
Постоянную времени передаточной функции (19) с учетом (4) рассчитывают по формуле
тдд = 3,09 (t4 + т5).
(20)
Численное значение постоянной времени числителя передаточной функции ФСЗ (19) определяют также с учетом ряда чисел правила «золотого сечения», приняв за целое три
Тзд = 3,382 (Т4 + Т5). (21)
Анализ результатов переходных процессов предлагаемых систем
Графики переходных процессов САУМБ с использованием пакета 81-mulink MatLab при отработке скачка задания по мощности для трех вариантов оптимизации параметров динамической настройки систем приведены на рис. 2-5.
1,0
у
I g, 0,8
л о
а й
S 6- ПЛ
™ И 0,6
К Р
s
U -8 I
0,4
0.2
0
— 1 -2 -3
50
150
200
100 Время, с
Рис. 2. Фактическая мощность при отработке скачка задания Fig. 2. Available power in the process of response to input step
_ 2,0
ffi № 3 I i 1
1,5
и H
I * 1.0
Щ о
0,5
h I О I
— 1 — 2 — 3
0
500
1500
2000
1000 Время, с
Рис. 3. Расход топлива при отработке скачка задания Fig. 3. Flow of fuel in the process of response to input step
— 1 — 2 — 3 -
500
1000 Вре:,м. с
1500
:oos
Рис. 4. Перемещение регулирующих клапанов турбины при отработке скачка задания Fig. 4. Turbine regulation valve motion in the process of response to input step
0,20 0,15 0,10 0,05 0
-0.05 -0,10 -0.15
-i — 1 — 2 -3
0
200
400 600 Время, с
800
1000
Рис. 5. Давление перегретого пара перед турбиной при отработке скачка задания Fig. 5. Over-heated steam pressure upstream of the turbine in the process of response to input step
Графики переходных процессов при внутреннем возмущении изображены на рис. 6-9.
15
х10-
и я
а я
а Р,
м о
(U fS
s э Э &
и
и
Я
10
и
И
3
х t О
-5
.................. .................. — i -2 -3
\
(c\
200
400 600 Время, с
800
1000
Рис. 6. Фактическая мощность при внутреннем возмущении Fig. 6. Available power in the process of response to internal disturbance
Рис. 7. Расход топлива при внутреннем возмущении Fig. 7. Flow of fuel in the process of response to internal disturbance
Время, с
Рис. 8. Перемещение регулирующих клапанов турбины при внутреннем возмущении Fig. 8. Turbine regulation valve motion in the process of response to internal disturbance
Рис. 9. Давление перегретого пара перед турбиной при внутреннем возмущении
Fig. 9. Over-heated steam pressure upstream of the turbine in the process of response to internal disturbance
На рис. 2-9 кривые 1 соответствуют схеме САУМБ (рис. 1, вариант I), оптимизированной по описанной выше методике, кривые 2 - варианту II и отличаются тем, что регулятор топлива настроен с учетом передаточной функции (1) по методу симметричного оптимума [15] с уменьшенным в два раза коэффициентом передачи:
Tri = 4ооп; (22)
■ (23)
1 оп
Кривые 3 соответствуют варианту III и отличаются от варианта I тем, что, во-первых, дифференциатор инвариантности превращен в пропорциональное звено с коэффициентом передачи, равным единице, во-вторых, на входе КРМ, оптимизированного по методу полной компенсации в частном виде, установлен ограничитель (±0,618), который уменьшает корректирующий сигнал на входе в РТ, что приводит к уменьшению максимального изменения расхода топлива при отработке задания.
Прямые показатели качества при отработке задания Ызд и внутреннего возмущения f1 сведены в табл. 1.
Таблица 1
Прямые показатели качества переходных процессов сравниваемых вариантов систем автоматического управления мощностью энергоблоков
Direct indicators of the quality of transients of the compared variants of systems of automatic control of power units
Вариант Вид возмущения tp, с AßT1 AhPK Aho
I ^зд 25 3,500 1,000 +0,1630 -
fl 900 0,125 0,007 +0,0065 0,00150
II ^зд 25 2,130 0,950 +0,1750 -
fl 900 0,280 0,050 0,0510 0,01250
III ^зд 25 1,100 1,000 0,1200 -
fl 250 0,122 0,009 0,0065 0,00135
Обозначения в табл. 1: tp - полное время регулирования; AßT1 - относительное максимальное изменение расхода топлива; AhpK - относительное максимальное изменение положения регулирующих клапанов турбины; Ah0 - максимальное отклонение давления перегретого пара перед турбиной; AN^ - максимальное отклонение фактической мощности при внутреннем возмущении.
Из анализа прямых показателей качества сравнительных вариантов САУМБ следует, что лучшие характеристики при отработке скачка задания и внутреннего возмущения обеспечивает вариант III системы с ограничением корректирующего сигнала на выходе из КРМ и подающий на вход последнего сигнал по давлению перегретого пара перед турбиной в диапазоне от 100 до 70 % номинальной мощности энергоблока.
ВЫВОДЫ
1. Предложена система автоматического управления мощностью энергоблоков (вариант III), отличающаяся от типовой наличием ограничите -ля корректирующего сигнала на выходе котельного регулятора мощности и заменой дифференциатора инвариантности на сигнал по давлению пере-
гретого пара перед турбиной, подаваемого на вход котельного регулятора мощности в режиме постоянного давления пара перед турбиной, при повышенном качестве управления с минимальными отклонениями расхода топлива и перемещения регулирующих клапанов турбины.
2. Результаты моделирования переходных процессов предлагаемой системы автоматического управления мощностью энергоблоков по сравнению с такими же типовыми системами, настроенными по современным методам параметрической оптимизации, показали существенное улучшение качества управления: при отработке скачка задания Л^зд полное время регулирования 1р во всех вариантах составляет 25 с. Однако при этом максимальное отклонение давления перегретого пара перед турбиной в третьем варианте уменьшилось на 35,8 % по сравнению с первым вариантом и на 45,8 % - по сравнению со вторым. Кроме того, в предлагаемом варианте системы автоматического управления мощностью энергоблоков максимальное изменение расхода топлива сократилось по сравнению со вторым в 1,94 раза, а по сравнению с первым - в 3,18 раза. Причем как максимальное изменение расхода топлива, так и перемещение регулирующих клапанов турбины оказались меньше лучших зарубежных аналогов [11].
3. При отработке внутреннего возмущения /1 время регулирования сократилось в 3,6 раза по сравнению с двумя первыми вариантами при примерно одинаковых других показателях качества.
ЛИТЕРАТУРА
1. Система регулирования энергоблока: а. с. 657179 СССР: МПК Е 01 К 13/02 / Г. Т. Кулаков, М. Н. Терешко, В. И. Литвинец [и др.]; дата публ. 15.04.1979.
2. Способ автоматического регулирования перетоков активной мощности и частоты энергообъединения: а. с. 864427 СССР: МПК Н 02 I 3/06 / Ю. М. Терезов, Г. Т. Кулаков, В. В. Гуренко; дата публ. 15.09.1981.
3. Устройство для автоматического регулирования перетока мощности между двумя частями энергосистемы: а. с. 864426 СССР: МПК Н 02 I 3/06 / Ю. М. Терезов, Г. Т. Кулаков, Л. М. Плахов, В. И. Литвинец; дата публ. 15.09.1981.
4. Система автоматического регулирования мощностью энергоблока: а. с. 877091 СССР: МПК Е 01 К 13/02 / Г. Т. Кулаков, А. А. Москаленко, И. Ю. Костив, Н. В. Молотков; дата публ. 30.10.1981.
5. Система группового автоматического управления мощностью энергоблоков тепловой электростанции: а. с. 1053221 СССР: МПК Н 02 I 3/46 / Г. Т. Кулаков, А. А. Москаленко,
A. Д. Качан; дата публ. 07.11.1983.
6. Способ управления энергоблоком в аварийных режимах энергосистемы: а. с. 881356 СССР: МПК Е 01 К 13/02 / Г. Т. Кулаков, И. Ю. Костив; дата публ. 15.11.1981.
7. Система автоматического регулирования частоты и мощности блочных ТЭС с прямоточными котлами / Н. И. Давыдов [и др.] // Теплоэнергетика. 1979. № 8. С. 2-6.
8. Система автоматического управления энергоблоком сверхкритических параметров /
B. И. Литвинец [и др.] // Энергетика. Изв. высш. учеб. заведений. 1980. № 8. С. 40-45.
9. Экспериментальное определение и анализ динамических свойств энергоблоков мощностью 300 МВт / Г. Т. Кулаков [и др.] // Энергомашиностроение. 1984. № 1. С. 9-11.
10. Система автоматического регулирования мощностью энергоблока паровой котел - турбина: пат. 2315871 Рос. Федерации: МПК Е 01 К 13/02 / Н. И. Давыдов; дата публ. 27.01.2008.
10. Система автоматического регулирования мощностью энергоблока паровой котел - турбина: пат. 2315871 Рос. Федерации: МПК Е 01 К 13/02 / Н. И. Давыдов; заявитель и патентообладатель Москва, ВТИ. № 2006121529/06; заявл. 20.06.2006; опубл. 27.01.2008, Бюл. № 3. 7 с.
11. Зорченко, Н. В. Исследование влияния формирующих сигналов в системе автоматического управления мощностью энергоблока на его приемистость / Н. В. Зорченко, Н. И. Давыдов, А. А. Григоренко // Теплоэнергетика. 2006. Т. 53, № 10. С. 810-818.
12. Разработка и внедрение САРЧМ крупных энергоблоков / В. А. Биленко [и др.] // Теплоэнергетика. 2008. Т. 55, № 10. С. 824-837.
13. Кулаков, Г. Т. Инженерные экспресс-методы расчета промышленных систем / Г. Т. Кулаков. Минск: Вышэйш. шк., 1984. 192 с.
14. Теория автоматического управления теплоэнергетическими процессами / Г. Т. Кулаков [и др.]; под общ. ред. Г. Т. Кулакова. Минск: Вышэйш. шк., 2017. 238 с.
15. Фрер, Ф. Введение в электронную технику регулирования / Ф. Фрер, Ф. Орттенбургер. М.: Энергия, 1973. 192 с.
Поступила 12.04.2018 Подписана в печать 27.06.2018 Опубликована онлайн 28.09.2018
REFERENCES
1. Kulakov G. T., Tereshko M. N., Litvinets V. I., Volnyanko M. G., Pirogov V. G., Kurilin L. A., Manuilov V. K. (1979) Power Unit Regulating System. USSR Author's Certificate No 657179 (in Russian).
2. Terezov Yu. M., Kulakov G. T., Gurenko V. V. (1981) Method for Automatic Control of Active Power Transfers and Power Pool Frequency. USSR Author's Certificate No 864427 (in Russian).
3. Terezov Yu. M., Kulakov G. T., Plakhov L. M., Litvinets V. I. (1981) Device for Automatic Control of Power Transfers between Two Parts of Power System. USSR Author's Certificate No 864426 (in Russian).
4. Kulakov G. T., Moskalenko A. A., Kostiv I. Yu., Molotkov N. V. (1981) The Automatic Control System for a Power Unit Capacity Regulation. USSR Author's Certificate No 877091 (in Russian).
5. Kulakov G. T., Moskalenko A. A., Kachan A. D. (1983) System of Automatic Regulation of Thermal Power Units' Group Capacity. USSR Author's Certificate No 1053221 (in Russian).
6. Kulakov G. T., Kostiv I. Yu. (1981) Method of Power Unit Control in Emergency Modes of Power System. USSR Author's Certificate No 881356 (in Russian).
7. Davydov N. I., Malamed A. D., Trakhtenberg M. D., Fotin L. P. (1979) The System of Automatic Regulation of Frequency and Power of Block-Type Thermal Power Plants with Straight-Through Boilers. Teploenergetika = Thermal Engineering, (8), 2-6 (in Russian).
8. Litvinets V. I., Kulakov G. T., Molev V. V., Red'ko I. I., Polishchuk S. I. (1980) Automatic Control System of Power Unit with Supercritical Parameters. Izvestiya Vysshikh Uchebnykh Zavedenii. Energetika [Proceedings of CIS Higher Education Institutions. Energetika], (8), 40-45 (in Russian).
9. Kulakov G. T. (1984) Experimental Determination and Analysis of Dynamic Properties for 300 MW Power Units. Energomashinostroenie [Power Plant Engineering], (1), 9-11 (in Russian).
10. Davydov N. I. (2006) System of Automatic Control of Power of Steam Boiler-Turbine Power Unit. Russian Federation Patent No 2315871 (in Russian).
11. Zorchenko N. V., Davydov N. I., Grigorenko A. A. (2006). Investigating the Effect of Forcing Signals in a System for Automatically Controlling the Output of a Power Unit on its Pickup Capability. Thermal Engineering, 53 (10), 810-818. https://doi.org/10.1134/s0040601506100107.
12. Bilenko V. A., Melamed A. D., Mikushevich E. E., Nikol'skii D. Y., Rogachev R. L., Romanov N. A. (2008) Development and Application of Automatic Frequency and Power Control Systems for Large Power Units. Thermal Engineering, 55 (10), 824-837. https://doi.org/10.1134/s00 40601508100030.
13. Kulakov G. T. (1984) Engineering Proximate Methods of Design Calculation for Industrial Systems. Minsk, Vysheishaya Shkola Publ. 192 (in Russian).
14. Kulakov G. T., Kulakov A. T., Kravchenko V. V., Kuchorenko A. N., Artsiomenka K. I., Kov-rigo Yu. M., Golinko I. M., Bagan T.G., Bunke A. S. (2017) Automatic Control Theory for Thermal Power Processes. Minsk, Vysheishaya Shkola Publ. 238 (in Russian).
15. Frer F., Orttenburger F. (1973) Introduction to Regulating Electronic Technology. Moscow: Energiya Publ. 192 (in Russian).
Recеived: 12 April 2018
Accepted: 27 June 2018
Published online: 28 September 2018
