УДК 621.321
Ресурсы стабилизации напряжения синхронного генератора АРВ сильного действия1
Д.В. Армеев, А.В. Чехонадских, Г.Б. Нестеренко ФГБОУВО «Новосибирский государственный технический университет» г. Новосибирск, Российская Федерация E-mail: alcheh@ngs.ru
Авторское резюме
Состояние вопроса: Распространение в энергетике новых типов генерирующих и распределительных систем, с одной стороны, и расширение технических возможностей регулирования с помощью устройств силовой электроники, с другой, стимулируют синтез алгоритмов управления, которые обеспечивали бы устойчивую работу электрогенератора в локальной сети с учетом ее особенностей и ограничений. Оптимальные настройки АРВ, реализующего ПИД или ПДД2-управление, остаются актуальным направлением теории регулирования. Наряду с этим синтез систем, оптимальных по степени устойчивости при числе настраиваемых параметров более трех, остается малоизученным в связи с крайней сложностью возникающих оптимизационных проблем. Материалы и методы: Нелинейная модель синхронного генератора с АРВ, реализующим ПИДД2-управление, построена в Matlab (Simulink) и линеаризована относительно действующих значений установившегося режима. Оптимальные по степени устойчивости расположения полюсов линейной САУ найдены алгебраическими средствами с помощью анализа критических корневых диаграмм, что позволяет определять оптимальные и субоптимальные расположения без ненадежных и обременительных численных процедур. Контроль переходных процессов и качественная оценка подавления типичного возмущения проводились в нелинейной модели. Результаты: Для четырехпараметрической ПИДД2-системы найдено экстремальное расположение полюсов с практически неограниченной степенью устойчивости, что достигается при чрезмерных значениях параметров регулирования. При фиксации интегрального либо пропорционального коэффициентов управления оптимизацией по трем параметрам получено удовлетворительное время подавления возмущения при отсутствии перерегулирования.
Выводы: настройка ПИДД2-управления позволяет достигать практически любой степени устойчивости линеаризованной модели генератора с АРВ за счет соответствующего увеличения коэффициентов регулирования. Целесообразным оказывается поиск оптимального регулятора с заданием приемлемого значения одного из коэффициентов. Это открывает перспективу стабилизации этим же методом двух- и трехмашинных систем АРВ сильного действия в локальных сетях электроснабжения.
Ключевые слова: синхронный генератор, стабилизация напряжения, установившийся режим, статическая устойчивость, полиномиальный синтез, ПИДД2-управление, оптимальное расположение полюсов, критические корневые диаграммы.
Resources of voltage stabilization of a synchronous generator by a multiparametric excitation controller
D.V. Armeyev, A.V. Chekhonadskikh, G.B. Nesterenko Novosibirsk State Technical University, Novosibirsk, Russian Federation E-mail: alcheh@ngs.ru
Abstract
Background: The wider use of new types of generation and distribution systems in the power industry, on the one hand, and the expansion of technical control options thanks to power electronics devices, on the other hand, stimulate the design of control algorithms which ensure stable operation of power generators in local networks with their specific conditions and limitations. The optimal settings of an automatic excitation controller, implementing PID or PDD2 control, remain the topical area of control theory. However, the problem of designing control systems of optimal stability with more than three control parameters is not studied well yet due to the extreme complexity of optimization problems.
Materials and methods: The nonlinear model of a synchronous generator with an automatic excitation controller, implementing PIDD2 control, was constructed in Matlab (Simulink) and linearized relative to the current values of the steady-state mode. The pole locations of the linear automatic control system ensuring the optimal stability were found algebraically by analyzing critical root diagrams, which allowed finding the optimal and suboptimal pole locations without
1
Работа выполнена при финансовой поддержке Министерства образования и науки РФ по государственному заданию № 2014/138, проект 1052.
unreliable and cumbersome numerical procedures. Transient control and qualitative estimation of typical perturbation suppression were carried out in the non-linear model.
Results: Extreme pole locations have been found for a four-parameter PIDD2 system with almost unlimited relative stability, which is achieved by setting excessive control parameter values. By fixing the proportional or the integral coefficient of optimization control by three control parameters the authors managed to achieve a satisfactory time of disturbance suppression without overshooting.
Conclusions: PIDD2-control setting enables reaching practically any stability degree of the linearized model of the generator with an automatic excitement controller through an appropriate control coefficient increase. It seems expedient to search for the optimal regulator by setting an acceptable value of one of the coefficients. This approach may be promising for stabilizing by the same method two- and three-generator systems with a multiparameter excitement controller in local power supply networks.
Key words: synchronous generator, voltage regulation, steady state mode, static stability, polynomial design, PIDD2 control, optimal pole location, critical root diagrams.
DOI: 10.17588/2072-2672.2017.1.024-032
Состояние вопроса. Одна из важнейших задач в производстве, передаче и распределении электроэнергии - поддержание заданного напряжения в узлах энергосистемы. Требования к уровню напряжения в сетях задаются, с одной стороны, условиями эксплуатации электрооборудования, а с другой - спецификой синхронной параллельной работы синхронных машин, в частности генераторов. Параллельная работа электроэнергетических систем, как правило, рассматривается в аспектах статической устойчивости нормального режима (способности системы возвращаться к исходному или близкому к нему режиму после малых возмущений) и динамической устойчивости (восстановление параметров режима после значительных возмущений, вызванных, в первую очередь, авариями, короткими замыканиями и т.п. [1]).
Каждый из современных силовых синхронных генераторов оснащается автоматическими регуляторами возбуждения (АРВ) сильного действия, действие которых необходимо и в нормальных, и в аварийных режимах. Их настройка находится в ведении специалистов ОАО «Системный оператор ЕЭС» и в сложившихся условиях не составляет трудной задачи2.
Однако общемировая тенденция в энергетике включает совершенствование электроэнергетических систем на основе инновационных технологий, среди которых «умные» сети типа SmartGrid или близкие к ним интеллектуальные электроэнергетические системы с активно-адаптивной сетью (ИЭС ААС), гибкие линии электропередачи (FACTS), широкомасштабные системы векторных измерений (WAMS) и близкие по смыслу системы мониторинга переходных режимов (СМПР), сети распределенной генерации, включающие нетрадиционные источники энергии, современные гибридные накопители энергии и др. [2, 3].
Такое обилие новых технологий и техники качественно меняет поведение энергосистем, что особенно проявляется в сетях типа SmartGrid
2 См.: Стандарт организации «Требования к системам возбуждения и автоматическим регуляторам возбуждения сильного действия синхронных генераторов». - М.: ОАО «СО ЕЭС», 2012. http://so-ups.ru/fileadmin/files/laws/standards/STO_ES_ AER_SA.pdf
с большим количеством и разнообразием источников электроэнергии. В этих условиях стандартные структуры регуляторов и способы их настройки оказываются недостаточными [2, 4, 5].
Мы рассматриваем задачу, связанную с регулированием напряжения в нормальном режиме и влиянием АРВ сильного действия на статическую устойчивость. Следует подчеркнуть, что традиционная для энергетики задача настройки АРВ с точки зрения теории регулирования осложняется тем, что АРВ воздействует на напряжение, а устойчивость энергосистемы контролируется, прежде всего, по углу электропередачи.
Основная функция АРВ в нормальном режиме - поддерживать напряжение на выводах генератора (или сборных шинах станции) с помощью регулирования напряжения возбуждения, в результате чего должны изменяться ток возбуждения, поле ротора машины, ЭДС и, в конечном итоге, регулируемое напряжение.
Входной величиной регулятора является отклонение напряжения от некоторого заданного значения (напряжения уставки иу). Критически важным для статической устойчивости оказывается эволюция другого параметра - взаимного угла 5 ротора генератора относительно вектора напряжения электрической системы ис, причем электромагнитный процесс в обмотке возбуждения описывается нелинейным дифференциальным уравнением второго порядка относительно угла 5.
В зависимости от системы возбуждения классифицируют АРВ двух типов. При электромашинном возбуждении используется регулятор пропорционального действия (ПД), что обычно обусловлено ограниченными ресурсами исполнительного устройства. При получении тока возбуждения с помощью силовой электроники применяется регулятор сильного действия (СД), в алгоритм которого включаются также производные отклонения напряжения и накопленное отклонение, т.е. интеграл, и др. [1].
Задача синтеза АРВ в последнем варианте сводится к выбору коэффициентов усиления при измеряемом отклонении напряжения и его производных в целях возможно быстрой стабилизации напряжения генератора иГ при сохранении устойчивости по углу 5 [1, 6].
Материалы и методы. Общепринятым подходом к нахождению коэффициентов регулятора системы автоматического управления (САУ) является линеаризация дифференциальных уравнений, описывающих переходные процессы, чтобы иметь возможность применения хорошо разработанной теории линейных систем и, в частности, перехода в частотную область. В этом случае основные критерии устойчивости связаны с расположением полюсов САУ в левой комплексной полуплоскости, что обеспечивает статическую устойчивость системы, а возможно больший сдвиг полюсов влево - лучшее время подавления возмущения. Для нахождения области в пространстве параметров, которая соответствует устойчивым значениям полюсов, используются различные приемы - от классического метода D-разбиения [7, 8] до современной теории робастных систем [7-10].
Однако метод D-разбиения эффективен для регуляторов с одним или двумя параметрами, а синтез трехпараметрических ПИД-регуляторов остается актуальным направлением ТАУ [10]. Настройка ПДД2-регуляторов, традиционно используемых в АРВ СД, осуществляется, главным образом, на основе методов с квадратичным функционалом качества, в ряде случаев приводящих итерационный процесс на границу устойчивости, либо эмпирически, что не гарантирует оптимальности. При синтезе электромеханических устройств, удовлетворяющих сложным критериям качества управления с применением наблюдателей, возникают разнообразные проблемы, новые подходы к решению которых предлагаются в [11]. Этим и рядом других факторов определяется практическая значимость создания новых подходов к настройке систем управления электрических машин и, в частности, АРВ.
Синтез оптимального четырехпараметри-ческого ПИДД2-регулятора осуществлен нами с помощью линеаризованной модели генератора в относительных единицах. Выбор действующих значений установившегося режима и переходных процессов замкнутой САУ проводился с помощью нелинейной модели, созданной стандартными средствами Matlab (Simulink). Для нахождения оптимальных по степени устойчивости замкнутой системы значений параметров используется метод оптимизации расположения полюсов САУ [12-14]. Это отличает метод, использованный нами от подхода [15], где оптимизация настроек ПДД2-регулятора возбуждения синхронного генератора решалась методом конечного градиента.
Уравнения генератора и его устойчивость в установившемся режиме. Основное уравнение, связывающее генератор и источник энергии мощностью РТ, выглядит следующим образом [1]:
При этом в самом генераторе активная и реактивная составляющие мощности связаны нелинейными уравнениями c напряжением возбуждения Eq и углом электропередачи 5:
Pr = E2qYi 1 sinа11 + EqUcY12 sin(5-a12);
Qr = E^Y-i cos a-, + EqUcYi2 cos(5-a12).
Связи напряжений (см. ниже) в обмотке возбуждения выражаются уравнениями:
(3)
dt2
-5 =
18000
Tj
(Pt - pr).
(1)
Eq - (xd - xd)UC Y12 cos(5 - a12)
^ 1 — (х^ — х^) Y i cos a^ i
d 1 ~HEq (Eqe - Eq)■
dt T d 0
Наконец, связь напряжения генератора с активной и реактивной составляющими мощности выражается соотношением
Ur = Ц—УиС2 + Qr xin )2 + (Pr xin )2. (4)
Значения параметров объекта в относительных единицах (о.е.) были приняты следующими (см. [16]):
• мощность первичного источника (например, турбины) РТ = 0,5 о.е.;
• постоянная механической инерции Tj = 1,703 c;
• собственная и взаимная проводимости ветвей генератора У11 = У12 = 0,4 о.е.;
• дополняющие фазные углы сопротивлений a11 = a12 = 0;
• напряжение системы UC = 1 о.е.;
• постоянная времени обмотки возбуждения при разомкнутой обмотке статора Td0 = 5 c;
• синхронное и переходное реактивное сопротивления по продольной оси xd = 2,3 о.е., x'd = 0,329 о.е. соответственно;
• внешнее сопротивление xin = 0,2 о.е.
Нелинейная модель такого объекта была
построена в пакете Matlab (Simulink). Рассчитанный при этих значениях установившийся режим разомкнутой системы (т.е. неуправляемого генератора) характеризуется следующими действующими значениями переменных:
• угол электропередачи 5 и 34,64°;
• активная мощность генератора РГ = 0,5 о.е.;
• реактивная мощность генератора Qr и 1,212 о.е.;
• напряжение генератора Ur и 1,246 о.е.;
• ЭДС холостого хода (напряжение возбуждения) и вынужденная ЭДС Eq = Eqe и 2,2 о.е.; переходная ЭДС Е q и 1,114 о.е.
Если одна из двух параллельных линий, связывающих генератор с энергосистемой, отключается, возникает типичный ступенчатый скачок напряжения ДЦ- и +0,04 о.е., который вводился в нелинейную модель в качестве возмущения. При подаче такого возмущения модель демонстрирует умеренную устойчивость по всем режимным переменным: время стабилизации
2
переходных процессов по напряжению Ur и углу 5 около 18 с; но она не обладает астатизмом, так как напряжение устанавливается на значении Ur » 1,357 о.е.
Стандартная процедура стабилизации ПД-регулятора (при K = K2 = 0) для нелинейной модели приводит к вполне удовлетворительной скорости стабилизации - время окончания переходных процессов чуть более 1 с, правда, при сохраняющейся статической ошибке: Ur и 1,253 о.е. (вместо Ur и 1,246 о.е.).
Ручная настройка ПИД-регулятора (в нашем случае при K2 = 0) может привести, например, к таким значениям параметров управления: K и 31,58; Кр и 26,94; K1 и -4,13; они обеспечивают астатизм системы и удовлетворительную скорость стабилизации (время завершения переходных процессов около 2,5 с); однако при этом возникает значительное перерегулирование (overshoot): возмущающий скачок AUr и +0,04 продолжается в начальной фазе процесса 0,0 - 0,2 c дальнейшим его нарастанием до AUr и +0,065.
Линеаризованная модель генератора и передаточная функция замкнутой системы. Уравнения генератора (1)-(4), линеаризованные относительно значений установившегося режима, в операторной форме принимают вид: s 2Д5 = -10569,6ДРг ; ДРГ = 0,2274ДЕ(? + 0,7240Д5; AQr = 1,4390ДЕд + 0,5252Д5;
(5)
(6)
(7)
ДЕ, = 4,7259ДЕ'q + 2,1179Д5 ; sAE'q= 0,2(AEqe - AEq); 1,246AUr = 0,2485AQ + 0,02 ДРГ. (8)
Исключая из (5)-(8) все переменные и выражая контролируемую переменную AUr через контролирующее воздействия AEqe, получаем передаточное соотношение генератора (plant):
м, NP' (s)AP
ДЦГ = 7Г^7 AEqe =
Dpi (s)
2(1,4288 s2 + 9559,2585)A E,
qe
10,4029s3 + 9,8326s2 + 1,32555s + 75243,968 ' Поскольку корни знаменателя попадают в правую полуплоскость, такой объект оказывается неустойчив к возмущениям в канале возбуждения; это не противоречит отмеченной выше устойчивости нерегулируемого генератора к возмущению на выводах. Техническое содержание и возможные причины подобной неустойчивости рассмотрены в [17, 18].
Для реализации АРВ сильного действия отрицательная обратная связь замыкается ПИДД2-регулятором (controller), все блоки которого (пропорциональный, интегральный и оба дифференциальных) хорошо отработаны на
практике . Поэтому его передаточная функция вместе с измерительным 1/^(5) и исполнительным ^2(5) устройствами (рис. 1), которые моделировались апериодическими звеньями с постоянными времени, принятыми равными 0,04 и 0,05 соответственно, принимает следующий вид:
AEqe = W2(s)
Ncon (s) Dcon (s)
W (s)AUr =
K, + Kps + Ks2 + K2s3 s(1 + 0,04s)(1 + 0,05s)
Д Ur.
Из Диофантова соотношения
f(s) = Dpi (s)Dcon(s) + Npi (s)Ncon (s) получается
приведенный характеристический многочлен замкнутой системы:
F (s) = s6 + (45,9452 + 137,3469K2 )s5 +
+(13248,64 + 137,347K )s4 +
+(5,8110 + 137,347Kp + 918903,56K2 )s3 + (9)
+(6696538,4 + 137,347K, + 918903,56^ )s2 + +(3616491,3 + 918903,56Kp )s + 918903,56K,.
Ресурсы стабилизации АРВ настройкой ПИДД2-управления. В самом сжатом виде идея отыскания параметров управления, обеспечивающих наибольшую степень устойчивости САУ с регулятором пониженного порядка, следующая. Рассмотрим наибольшую действительную часть полюсов системы (иными словами, правую границу расположения полюсов на комплексной плоскости) как целевую функцию, зависящую от параметров управления4. Минимизация этой функции, т.е. смещение влево правой границы, сначала оказывается сдвигом одного характеристического корня или комплексной пары; с какого-то момента на границе оказывается несколько действительных корней (кратный корень) или комплексных пар, и сдвиг должен осуществляться для всего ансамбля правых корней; пока таких корней немного, возможна дальнейшая минимизация [13]; когда же число связей между корнями, обусловленное их нахождением на правой границе, совпадает с числом параметров управления, дальнейшая минимизация может оказаться невозможной: в пространстве параметров возникает нульмерное критическое многообразие - одна или несколько точек.
Благодаря современным датчикам, быстродействие которых существенно превосходит характерное время развития переходных процессов, ДД2-регуляторы
убедительно реализуют теоретическую передаточную
функцию Ks + K2S2. Стоит отметить, что возможности современной силовой электроники значительно ослабляют требование правильности (properness) регулятора.
4 Понятие гурвицевой функции H(K) = max Re Zk как правой границы расположения полюсов в зависимости от
вектора K свободных параметров управления расширяет
общепринятое понятие степени устойчивости системы, допуская значения любого знака и тем самым включая в себя также «степень неустойчивости».
Рис. 1. Модель ПИДД2-регулятора в МАНАВ^тиНпк)
Вариантов подобных граничных расположений действительных корней и комплексных пар несколько, они зависят от числа параметров управления и называются критическими корневыми диаграммами [14]. В случае четырех свободных параметров возникает 13 диаграмм, которые и были нами изучены.
Часть из них не реализуется из-за недостаточной степени характеристического многочлена (9), другие не допускают устойчивых полюсов. Например, трехкратная комплексная пара полюсов г1-6 = х ± ¡у возникает только при х и 7,234; 30,924, т.е. при заведомо неустойчивой системе.
Единственное найденное устойчивое решение соответствует двукратному действительному корню х12 и -36,448 почти на одной вертикали с двойной комплексной парой, также почти кратной: г34 и -36,446 ± 66,786/ и
256 и-36,449 ± 66,787/ . Это решение соответствует значениям коэффициентов управления, два из которых неоправданно велики: К, и 48447,2 ; Кр и 3874,62 ; К, и 113,312 ;
К2 и 1,2577.
Этот глобальный минимум гурвицевой функции
Н(К,КР,К,К2) = тахРе^) и -36,45
соответствует чрезмерно высокой с практической точки зрения степени устойчивости ~ ехр(-36,45/); как замечено выше, он достигнут при не имеющих практического смысла значениях коэффициентов АРВ.
Из этого следует важный вывод: ПИДД2-регулятор позволяет достигать практически любой степени устойчивости замкнутой системы ценой увеличения значений коэффициентов регулятора, возможно, за рамками технической целесообразности.
Замечание. Моделирование переходных процессов в пакете МАНАВ^тиПпк) для экстремальных значений параметров К, ,Кр ,К1,К2
дает весьма высокую скорость подавления возмущения: начальный отрезок графика иГ©
практически вертикален, а завершение переходных процессов по напряжению происходит за 0,2 с; однако процессы по мощности и углу устанавливаются только за 1,2 с (рис. 2). Очевидно, что процессы в нелинейной модели не оправдывают степени устойчивости линейной модели по переменной иГ за счет сравнительно более медленного подавления возмущения переменных РГ и 5.
Тем самым вопрос об оптимальности ПИДД2-управления сводится к согласованию достаточной скорости стабилизации с приемлемыми значениями параметров регулятора К и Кр. Таким образом, остается установить, какие приемлемые значения коэффициентов регулятора обеспечивают удовлетворительное качество подавления возмущений. Это приводит к задачам поиска условного экстремума гурвицевой функции н(к/рд2) при фиксации одного из коэффициентов регулятора.
Ресурсы стабилизации АРВ настройкой ПДД2-управления. Первая возможность заключается в задании интегрального параметра: положим К, = 50.
Список критических корневых диаграмм в этом случае сокращается до восьми [14], из которых теоретически нереализуемыми оказываются две5.
Поиск реализующихся устойчивых критических расположений полюсов приводит трехкратному действительному корню на одной вертикали с простой комплексной парой, что выражается в следующих значениях полюсов: корни и-4,338,^ и-4,363±0,014/ представляют собой «почти трехкратный» действительный корень, который находится практически на той же вертикали, что и комплексная пара г45 и-4,354± 110,930/; шестой корень
г6 и -45,162 расположен значительно левее.
Это находящиеся на одной вертикали четыре простые комплексные пары либо три пары и действительный корень - в обоих случаях число корней превышает степень характеристического многочлена.
Рис. 2. Переходные процессы по мощности РГ (а), углу электропередачи 8 (б) и напряжению генератора иГ (в) в нелинейной модели генератора с ПИДД2-регулятором при значениях коэффициентов, соответствующих максимальной степени устойчивости линеаризованной модели
Это расположение полюсов возникает при параметрах управления Кр«31,657,
К «1,409 , К - 0,153 (как условлено выше, К= 50).
Время стабилизации при этом около 2,5 с; однако при моделировании процессов с такими параметрами управления в нелинейной модели возникает заметное перерегулирование: если начальное возмущение ДЦ- = 0,040 о.е., то в переходном процессе через 0,2 с возникает всплеск диГ(0,2) - 0,052 о.е.
Замечание. Критическое корневое расположение другого типа: двукратный действительный корень 2 «-0,05 на одной вертикали с
простыми парами
' 3,4
¡-0,05 ± 81,71/
5,6
-0,05 ±1647,,5 / находится практически на
границе устойчивости и за счет большой колебательности второй из них теряет смысл; слишком велики и значения двух коэффициентов: Кр «1981,6 и К « 19705,7 (К2 «-0,33).
Ресурсы стабилизации АРВ настройкой ИДД2-управления. Фиксируем коэффициент пропорционального усиления Кр = 50. Список возможных критических диаграмм тот же, и наибольшая устойчивость достигается при аналогичном расположении корней: «почти трехкратный действительный» корень ^ «-5,379, г23 «-5,396±0,010/ находится приблизительно на той же вертикали, что и про-
стая пара г45 «-5,391 ± 110,831/; шестой корень г6 « -44,306 находится существенно левее.
Такое расположение полюсов достигается при параметрах управления К, «93,0 ,
К « 3,522 , К2 « 0,184. Переходные процессы при таких значениях параметров в нелинейной системе также дают несколько более убедительный результат (рис. 3). Время подавления возмущения чуть больше 2 с, перерегулирования по напряжению не возникает.
Нелинейная модель с таким регулятором демонстрирует также удовлетворительную ро-бастность, сохраняя форму и численные значения переходных процессов при изменении параметров управления на ±5 %.
Замечание. Еще одно решение с меньшей степенью устойчивости реализуется при
параметрах
управления
К. «15,62,
К1 « 54,80, К2 « 23,21. Степень устойчивости Н(К/12) «-0,869 достигается при аналогичном
расположении полюсов: «почти трехкратный действительный» корень
^ «-0,8695, г23 «-0,869 ± 0,005/ находится на
одной вертикали с комплексной парой г45 «-0,8690±82,300/, а шестой корень
г6 « -3229,4 - далеко слева.
а)
б)
в)
Рис. 3. Переходные процессы по мощности РГ (а), углу электропередачи 5 (б) и напряжению генератора иГ (в) в нелинейной модели генератора с ПИДД2-регулятором при значениях коэффициентов регулятора К и 93,0, К и 93,0, Кр = 50, К и 3,522, К2 и 0,184
Однако при моделировании управления с такими параметрами регулятора нелинейная система оказывается неустойчивой. По-видимому, это объясняется сравнительно большими значениями дифференциальных коэффициентов К1,2 и резким всплеском управляющего воздействия, возникающим при численной аппроксимации неограниченных первой и второй производных ступенчатого возмущения. Технически это означало бы недопустимо высокую чувствительность к шумам.
Результаты. Моделирование синхронного генератора с ПИДД2-регулятором возбуждения осуществлялось двумя способами: поиск оптимального регулятора осуществлялся с помощью модели, линеаризованной для действующих значений установившегося режима, а подбор самого установившегося режима и контроль переходных процессов для оптимальных регуляторов выполнялись в нелинейной модели, простроенной стандартными средствами пакета МаАаЬ.
Хотя переходные процессы в нерегулируемом генераторе оказались устойчивы к скачкообразному возмущению на шинах, принятые параметры установившегося режима в линеаризованной модели привели к неустойчивым полюсам передаточной функции объекта (от вынужденной ЭДС Бце в качестве контролирующего воздействия к напряжению на зажимах иГ в качестве контролируемой величины) и тем самым к содержательной задаче его стабилизации.
Оптимальные регуляторы для линеаризованной модели находились алгебраическими
средствами с помощью критических корневых диаграмм и корневых многочленов по той же схеме, как это было сделано в [12-14].
Четырехпараметрический ПИДД2-регуля-тор обеспечивает практически не ограниченную степень устойчивости линеаризованной системы за счет чрезмерно больших значений коэффициентов усиления; однако при моделировании такого регулятора в нелинейной системе подавление возмущения происходит быстро, но немгновенно (на этом сказывается более медленное протекание процессов стабилизации угла электропередачи и мощности). Это приводит к задаче поиска оптимального регулятора в рамках принятых ограничений на параметры управления.
Среди рассмотренных возможностей наилучший результат был получен при фиксированном пропорциональном коэффициенте Кр = 50; удовлетворительное время подавления возмущения и отсутствие перерегулирования достигнуты при умеренных значениях и невысокой чувствительности к изменениям других коэффициентов усиления.
Выводы. Полученные результаты позволяют рассчитывать на успешное применение этого метода (линеаризации нелинейной модели относительно значений установившегося режима, поиска оптимальных настроек ПИДД2-регулятора возбуждения линейной модели методом критических корневых диаграмм и контроля переходных процессов в нелинейной модели при найденных параметрах управления) для эффективной стабилизации двух- и трехмашинных систем генерации электроэнер-
гии в локальных сетях. Отдельного изучения заслуживает вопрос о связи степени устойчивости линеаризованной модели со временем подавления возмущения в исходной нелинейной модели.
Список литературы
1. Веников В.А. Переходные электромеханические процессы в электрических системах. - М.: Высш. шк., 1985. - 536 с.
2. Momoh J. SMART GRID Fundamentals of Design and Analysis // IEEE PRESS. - Hoboken, New Jersey: John Wiley & Sons, 2012.
3. Kassakian J., Schmalensee R. The future of the electric grid: An interdisciplinary MIT study. Technical report / Massachusetts Institute of Technology, 2011.
4. Demeo A., Peterson M.L. Community Smart Grid Utilizing Dynamic Demand Response and Tidal Power for Grid Stabilization // Smart Grid and Renewable Energy. - 2013. -No 4. - P. 465-472.
5. Sharaf A.M., Gandoman F.H. FACTS Based Stabilization for Smart Grid Applications // International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering. - 2014. - Vol. 8, No 11. - P. 1807-1811.
6. Lakshma Naik V. A Comparison Scheme of Dynamic Voltage for Smart Electric Grid Stabilization and Efficient Utilization using FACTS // International Journal of Latest Trends in Engineering and Technology. - 2015. - Vol. 6, Issue 2. -P. 353-364.
7. Жданов П.С. Вопросы устойчивости электрических систем. - М.: Энергия, 1979.
8. Поляк Б.Т., Щербаков П.С. Трудные задачи линейной теории управления. Некоторые подходы к их решению // Автоматика и телемеханика. - 2005. - № 5. -С. 7-46.
9. Polyak B.T., Gryazina E.N. Stability domain in the parameter space: D-decomposition revisited // Automatica. -2006. - Issue 42. - P. 13-26.
10. Datta A., Ho M.-T., Bhattacharyya S.R. Structure and Synthesis of PID Controllers. - N.Y.: Springer, 2000.
11. Анисимов А.А., Тарарыкин С.В., Аполон-ский В.В. Параметрическая оптимизация электромеханических систем с регуляторами и наблюдателями состояния // Вестник ИГЭУ. - 2016. - Вып. 2. - С. 21-26.
12. Чехонадских А.В. О ступенчато-дифференциальной оптимизации корней характеристического многочлена САУ // Науч. вестник НГТУ. - 2008. - № 4(33). -С. 205-208.
13. Корюкин А.Н., Чехонадских А.В. Предел устойчивости трехмассовой системы с регулятором 3-го порядка. Ч.1 // Сб. науч. тр. НГТУ. - 2011. - № 4(66). - С. 3-22.
14. Чехонадских А.В. Экстремальные расположения полюсов систем автоматического управления с регулятором пониженного порядка // Автоматика и телемеханика. -2014. - № 10. - С. 6-24.
15. Armeev D.V., Chekhonadskikh A.V., Voevoda A.A. Modal optimization of AVR for synchronous generator using the finite gradient // International Siberian conference on control and communications (SIBC0N-2015): proc., Omsk, 21-23 May, 2015. - Omsk: IEEE, 2015.
16. Шойко В.П. Автоматическое регулирование в электрических системах. - Новосибирск: НГТУ, 2012.
17. Хрущев Ю.В., Заподовников К.И., Юшков А.Ю. Электромеханические переходные процессы в электроэнергетических системах. - Томск: Изд-во ТПУ, 2012.
18. Москвин И.А. Колебательная статическая устойчивость электроэнергетической системы // Вестник ИГЭУ. -2013. - Вып. 5.- С. 46-50.
References
1. Venikov, V.A. Perekhodnye elektromekhanicheskie protsessy v elektricheskikh sistemakh [Electromechanical
transients in electric systems]. Moscow, Vysshaya shkola, 1985. 536 p.
2. Momoh, J. SMART GRID Fundamentals of Design and Analysis. IEEE PRESS. Hoboken, New Jersey, John Wiley & Sons, 2012.
3. Kassakian, J., Schmalensee, R. The Future of the Electric Grid: An interdisciplinary MIT Study. Technical Report. Massachusetts Institute of Technology, 2011.
4. Demeo, A., Peterson, M.L. Community Smart Grid Utilizing Dynamic Demand Response and Tidal Power for Grid Stabilization. Smart Grid and Renewable Energy, 2013, no. 4, pp. 465-472.
5. Sharaf, A.M., Gandoman, F.H. FACTS Based Stabilization for Smart Grid Applications. International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering, 2014, vol. 8, no. 11, pp. 1807-11.
6. Lakshma Naik V. A Comparison Scheme of Dynamic Voltage for Smart Electric Grid Stabilization and Efficient Utilization using FACTS. International Journal of Latest Trends in Engineering and Technology, 2015, vol. 6, issue 2, pp. 353-364.
7. Zhdanov, P.S. Voprosy ustoychivosti elektricheskikh sistem [Stability problems of electric systems]. Moscow, Ener-giya, 1979.
8. Polyak, B.T., Shcherbakov, P.S. Trudnye zadachi li-neynoy teorii upravleniya. Nekotorye podkhody k ikh resheniyu [Difficult Problems of Linear Control Theory: Possible Approaches to their Solution]. Avtomatika i telemekhanika, 2005, vol. 66, no. 5, pp. 681-718.
9. Polyak, B.T., Gryazina, E.N. Stability domain in the parameter space: D-decomposition revisited. Automatica, 2006, issue 42, pp.13-26.
10. Datta, A., Ho, M.-T., Bhattacharyya, S.R. Structure and Synthesis of PID Controllers. N.Y., Springer, 2000.
11. Anisimov, A.A., Tararykin, S.V., Apolonsky, V.V. Pa-rametricheskaya optimizatsiya elektromekhanicheskikh sistem s regulyatorami i nablyudatelyami sostoyaniya [Parametrical optimization of regulators and state observers in electromechanical systems]. Vestnik IGEU, 2016, issue 2, pp. 21-26.
12. Chekhonadskikh, A.V. O stupenchato-differentsial'noy optimizatsii korney kharakteristicheskogo mnogochlena SAU [On stage-differential optimization of control system characteristic roots]. Nauchnyy vestnik NGTU, 2008, issue 4(33), pp. 205-208.
13. Koryukin, A.N., Chekhonadskikh, A.V. Predel ustoychivosti trekhmassovoy sistemy s regulyatorom 3-go poryadka, ch. 1 [Stability limit of a 3-mass system with a third-order controller, p. 1]. Sbornik nauchnykh trudovNGTU, 2011, no. 4(66), pp. 3-22.
14. Chekhonadskikh, A.V. Ekstremal'nye raspolozheniya polyusov sistem avtomaticheskogo upravleniya s regulyatorom ponizhennogo poryadka [Extreme pole locations in automatic control systems with a reduced-order controller]. Avtomatika i telemekhanika, 2014, no. 10, pp. 6-24.
15. Armeev, D.V., Chekhonadskikh, A.V., Voevoda, A.A. Modal optimization of AVR for synchronous generator using the finite gradient. Proc. International Siberian conference on control and communications (SIBC0N-2015), Omsk, 21-23 May, 2015. Omsk, 2015.
16. Shoiko, V.P. Avtomaticheskoe regulirovanie v elek-tricheskikh sistemakh [Automatic control in electric systems]. Novosibirsk, NGTU, 2012, 195 p.
17. Khrushchev, Yu.V., Zapodovnikov, K.I., Yushkov, A.Yu. Elektromekhanicheskie perekhodnye protsessy v elektroener-geticheskikh sistemakh [Electromechanical transients in electrical power systems]. Tomsk, Izdatel'stvo TPU, 2012.
18. Moskvin, I.A. Kolebatel'naya staticheskaya ustoychi-vost' elektroenergeticheskoy sistemy [Oscillatory steady state stability of an electrical power system with an interconnection containing controlled series capacitors]. Vestnik IGEU, 2013, issue 5, pp. 46-50.
Aрмеев Денис Владимирович,
ФГБОУВО «Новосибирский государственный технический университет»,
кандидат технических наук, доцент кафедры автоматизированных электроэнергетических систем, e-mail: denarpost@yandex.ru
Armeyev Denis Vladimirovich, Novosibirsk State Technical University,
Associate Professor of the Department of Automated Electric Power Systems, e-mail: denarpost@yandex.ru
Чехонадских Александр Васильевич,
ФГБОУВО «Новосибирский государственный технический университет», доктор технических наук, профессор кафедры алгебры и математической логики, e-mail: alcheh@ngs.ru
Chekhonadskikh Aleksandr Vasilyevich, Novosibirsk State Technical University,
Professor of the Department of Algebra and Mathematical Logic, e-mail: alcheh@ngs.ru
Нестеренко Глеб Борисович,
ФГБОУВО «Новосибирский государственный технический университет», магистрант кафедры автоматизированных электроэнергетических систем, e-mail: nesterenkogbn@mail.ru.
Nesterenko Gleb Borisovich, Novosibirsk State Technical University,
Master Course Student of the Department of Automated Electric Power Systems, e-mail: nesterenkogbn@mail.ru