Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,087

SYNTHESIS OF CONTROL SYSTEMS PRIME MOVERSSYNCHRONOUS GENERATORS

Melnikov D.V. 1 Kyaw Thu Aung 1 Oukar Min 1 Min Chzho Tu 1
1 Bauman Moscow State Technical University (Kaluga Branch)
В статье на базе современных проекционно-матричных методов представлен алгоритм параметрического синтеза системы регулирования первичного двигателя электрических синхронных генераторов, которым служит паровая турбина. Алгоритм основан на представлении процессов в форме разложения по базисной системе функций, использует аппарат структурных преобразований и методы математического программирования. Отличительной особенностью алгоритма является его ориентация на класс нелинейных систем, причем количество нелинейных элементов не имеет принципиальных ограничений и существенно не отражается на быстродействии вычислительных алгоритмов. Это стало возможно благодаря разработке «быстрых» алгоритмов замены нелинейных элементов соответствующими эквивалентными матричными операторами. В случае детерминированной постановки нелинейный элемент заменяется эквивалентным линейным элементом с переменным коэффициентом передачи, что эквивалентно матричному оператору умножения. Критерий эквивалентности – равенство выходных сигналов нелинейного и эквивалентного ему линейного элемента, заданного матричным оператором, при отработке конкретного входного сигнала. Алгоритм приводится на примере конкретной энергетической турбины, математическая модель которой описывается структурной схемой.
The article on the basis of the modern projection-matrix methods the algorithm parametrize-tion of synthesis of system of regulation of the primary motor electric synchronous generators, which is a steam turbine. The algorithm is based on processes in the form of decomposition on basic system functions, uses the machine structural reforms and methods of mathematical-sky programming. A distinctive feature of the algorithm is its focus on a class of nonlinear systems, and the number of nonlinear elements is not fundamental limitations, and not significantly effect on performance of the computational algorithms. It became possible thanks to the development of the «fast» algorithms replacing the nonlinear elements corresponding equivalent matrix operators. In the case of a deterministic formulation of the nonlinear element is replaced by the equivalent linear element with a variable transfer coefficient, which is equivalent to a matrix multiplication operator. Equivalence criterion is the equality of output signals of nonlinear and equivalent linear element defined matrix operator, when developing a specific input signal. The algorithm provides for example, a particular energy include turbine mathematical model which describes the structural scheme.
power turbine
generator
matrix operator
orthonormal basis
the regulatory system
the synthesis controller
algorithm
a mathematical model
1. Aung Ch.S., Makarenkov A.M., Mo P.S. Identifikacija sluchajnyh parametrov matematiche-skoj modeli jelektrogidravlicheskogo sledjashhego privoda // Fundamentalnye issledova-nija. 2016. no. 2–2. рр. 231–235.
2. Gusev A.S., Svechkarev S.V., Plodistyj I.L. Matematicheskaja model pervichnyh dvigatelej sinhronnyh generatorov // Izvestija Tomskogo politehnicheskogo universiteta. 2005. T. 308. no. 7. рр. 216–221.
3. Kornjushin Ju.P, Melnikov D.V., Egupov N.D., Kornjushin P.Ju. Issledovanie i raschet pa-rametrov jelementov sistemy regulirovanija chastoty vrashhenija rotora turbiny s uchetom pa-rametricheskoj neopredelennosti matematicheskoj modeli // Vestnik MGTU im. N.Je. Bauma-na. Ser. Estestvennye nauki. no. 1.рр. 78–93.
4. Melnikov D.V., Egupov N.D. Sintez sistem regulirovanija jenergeticheskih turbin v uslovijah parametricheskoj neopredelennosti // Izvestija Tulskogo gosudarstvennogo universiteta. Tehnicheskie nauki. 2011. no. 5–1. рр. 108–113.
5. Melnikov D.V., Kornjushin Ju.P., Min Chzho Tu, Chzho Tu Aung, Okar Min. Algoritm sinteza sistemy regulirovanija chastoty vrashhenija rotora jenergeticheskoj turbiny // Nauchnoe obo-zrenie. 2015. no. 20. рр. 138–143.
6. Melnikov D.V., Kornjushin P.Ju., Min Ch.T., Chzho T.A., Okar M. Proekcionno-matrichnyj podhod k analizu i sintezu sistem upravlenija jelektrojenergeticheskih sistem // Nauchnoe obozrenie. 2015. no. 2. рр. 88–97.
7. Melnikov D.V. Metod avtomatizirovannogo issledovanija sistem regulirovanija jenergeti-cheskih turbin pri sluchajnyh vozmushhenijah: dis. ... kand. tehn. nauk. Kaluga, 2002.
8. Melnikov D.V., Min Chzho Tu. Matematicheskaja model sistem regulirovanija jenergeti-cheskih turbin s reguliruemymi otborami para // Inzhenernyj zhurnal: nauka i innovacii. 2015. no. 2 (38). рр. 13. http://engjournal.ru/catalog/pmce/tctp/ 1366.html.
9. Melnikov D.V., Okar Min, Chzho Tu Aung. Matematicheskoe modelirovanie sistem reguliro-vanija jenergeticheskih turbin // Nauchnoe obozrenie. 2015. no. 20. рр. 144–149.
10. Melnikov D.V. Proekcionno-matrichnyj metod sinteza kontura regulirovanija chastoty vrashhenija rotora parovoj turbiny // Vestnik Moskovskogo gosudarstvennogo tehnicheskogo universiteta im. N.Je. Baumana. Serija: Mashinostroenie. 2013. no. 4 (93). рр. 43–53.
11. Melnikov D.V., Fisher M.R. Matematicheskaja model kontura regulirovanija chastoty vra-shhenija rotora parovoj turbiny K-800-130 // Vestnik Moskovskogo gosudarstvennogo tehni-cheskogo universiteta im. N.Je. Baumana. Serija: Mashinostroenie. 2011. no. 7. рр. 197–215.
12. Normy uchastija jenergoblokov teplovyh jelektrostancij v normirovannom pervichnom regu-lirovanii chastoty i avtomaticheskom vtorichnom regulirovanii chastoty i peretokov aktivnoj moshhnosti // Standart organizacii. Otkrytoe akcionernoe obshhestvo «Sistemnyj operator edinoj jenergeticheskoj sistemy». Izdanie oficialnoe. M., 2013. http://so-ups.ru/fileadmin/files/laws/standards/sto_002-2013_freq_regulation.pdf.
13. Pupkov K.A., Egupov N.D., Kolesnikov L.V., Lukashenko Ju.L., Melnikov D.V., V.M. Ry-bin, Trofimov A.I. Matrichnye metody rascheta i proektirovanija slozhnyh sistem avtomaticheskogo upravlenija dlja inzhenerov / pod red. K.A. Pupkova i N.D. Egupova. M.: Izdatelstvo MGTU im. N.Je. Baumana, 2007. 664 р.
14. Okar Min, Melnikov D.V. Algoritm rascheta nelinejnyh sistem upravlenija proek-cionno-matrichnym metodom // Inzhenernyj zhurnal: nauka i innovacii. 2014. no. 12 (36). рр. 17. http://engjournal.ru/catalog/it/asu/1268.html.
15. Chzho Tu Aung, Melnikov D.V. Algoritm issledovanija nelinejnyh sistem avtomaticheskogo upravlenija v stohasticheskih rezhimah // Inzhenernyj zhurnal: nauka i innovacii. 2014. no. 4 (28). рр. 3. http://engjournal.ru/catalog/it/asu/1270.html.

Одним из наиболее сложных элементов современных энергетических систем является первичный двигатель электрических синхронных генераторов, которым обычно служит паровая турбина со своими системами регулирования. Основным способом получения информации, необходимой для проектирования, эксплуатации и развития электроэнергетических систем, является математическое моделирование [2]. Динамика первичных двигателей существенно влияет на нормальные и аварийные процессы производства, распределения и потребления электроэнергии. Поэтому требование к уровню адекватности модели первичного двигателя постоянно возрастает. А это означает, что на сегодняшний день описывать модели систем регулирования приходится в классе нелинейных систем.

Наиболее наглядной и исходной формой описания сложной динамической модели энергетической системы является структурная (операторная) схема, которая не только иллюстрирует содержание модели, но и показывает структуру моделируемых объектов, что позволяет легко ориентироваться в ней и адаптировать ее для конкретных целей. Поэтому актуальной является задача разработки алгоритмов синтеза систем регулирования, заданных структурной схемой.

Решение поставленной задачи осложнено тем, что в настоящее время в электроэнергетических системах находится в эксплуатации множество различных типов турбин и котлоагрегатов с разнообразными системами и законами регулирования. Хотя передаточные функции отдельных элементов и звеньев известны, разработка полной математической модели для конкретного типа оборудования и компоновки первичного двигателя, а также расчета системы регулирования является уникальной [2, 8, 9, 11]. Поэтому возникает задача разработки универсального алгоритма синтеза, который может использоваться для стационарных и нестационарных систем, как для линейных, так и нелинейных систем управления. Такие алгоритмы могут быть построены на использовании метода матричных операторов расчета и проектирования сложных систем [6, 13]. Для иллюстрации таких алгоритмов в работе приводится синтез системы регулирования турбины производства Калужского турбинного завода. Хотя модель не является громоздкой, но на ее примере довольно подробно иллюстрируются все особенности алгоритма, что не снижает его общности по отношению к системам высокого порядка [10].

Общий алгоритм синтеза системы регулирования

Рассмотрим систему регулирования турбины ПТ-25/30-90/10М производства Калужского турбинного завода в классе нелинейных систем (рис. 1) [7]. Эта система предназначена для поддержания частоты вращения ротора турбины на требуемом уровне.

На рис. 1 обозначено: φ – относительное отклонение угловой частоты вращения ротора от номинального значения ωНОМ, Melnikov01.wmf; λ – относительное отклонение электрической нагрузки генератора; ξ – относительное изменение расхода пара в турбине; μ – относительное отклонение поршня сервомотора; z – относительное отклонение поршня золотника; δ – степень неравномерности системы регулирования; Tz – постоянная времени золотника; Tμ – постоянная времени сервомотора; Tξ – постоянная времени камеры эквивалентного объема (время полного заполнения камеры при определенном расходе пара); Tφ – постоянная времени ротора; θ – эквивалентный коэффициент самовыравнивания.

pic_22.wmf

Рис. 1. Структурная схема нелинейной системы регулирования энергетической турбины ПТ-25/30-90/10М (ПТ-12/15-35/10М)

Нелинейность F(z), определяющая форму кромок отсечного золотника, описывается следующей зависимостью:

Melnikov02.wmf

где Δz0 – относительные высоты треугольного профиля кромок поршня отсечного золотника; ε0 – величина, характеризующая степень нечувствительности сервомотора.

Требуется определить параметры системы регулирования δ, Tz, Tμ таким образом, чтобы φЭ(t) и φ(t) (реальный переходной процесс) были близки, например, в среднеквадратическом смысле:

Melnikov03.wmf

Представим модель системы регулирования турбины с нелинейностью в операторной форме (в матричном виде) (рис. 2).

На схеме обозначено:

Melnikov04.wmf

Melnikov05.wmf Melnikov06.wmf

Melnikov07.wmf Melnikov08.wmf

где AИ – матричный оператор интегрирования. Из рис. 2 имеем следующие векторно-матричные соотношения:

Melnikov09.wmf Melnikov10.wmf

Melnikov11.wmf Melnikov12.wmf Melnikov13.wmf

Melnikov14.wmf

или

Melnikov15.wmf Melnikov16.wmf

pic_23.wmf

Рис. 2. Структурная схема нелинейной системы регулирования энергетической турбины в операторной форме

Последнее выражение представляет собой матричный оператор всей системы, который связывает спектральные характеристики входного и выходного процессов. Матричный оператор нелинейного элемента вычисляется по формуле [6]:

Melnikov17.wmf (1)

Для вычисления спектральных характеристик остальных процессов системы нужно воспользоваться структурными преобразованиями. В результате можно получить следующие зависимости:

Melnikov18.wmf Melnikov19.wmf Melnikov20.wmf Melnikov21.wmf Melnikov22.wmf

Алгоритм, определяющий процедуру решения задачи расчета параметров системы управления турбины, представлен на рис. 3.

pic_24.wmf

Рис. 3. Структурная схема алгоритма синтеза системы регулирования турбины в классе нелинейных систем

Вычислительный алгоритм синтеза системы регулирования первичного двигателя

Рассмотрим подробно алгоритм параметрического синтеза системы регулирования турбины ПТ-25/30-90, представленного на рис. 3.

Этап 1. При заданной структуре регулятора формируем вектор искомых параметров: p = [δ, Tz, Tμ]. Здесь же инициализируем вектор параметров начальными приближениями p0 = [1, 1, 1].

Этап 2. Накладываем ограничения на параметры, связанные с их физической реализуемостью. Эти ограничения будут пересланы в функцию оптимизации в качестве опции. Данная операция позволяет значительно сократить время поиска оптимальных параметров. Например, все полученные параметры регуляторов должны быть положительными.

Этап 3. Задаем эталонный переходный процесс при резком сбросе нагрузки (λ = –1):

Melnikov23.wmf

Стоит обратить особое внимание на то, чтобы эталонный переходный процесс соответствовал требованиям, предъявляемым в [12].

Этап 4. Выбор ортонормированного базиса на промежутке исследования [0, T]. В дальнейшем всегда будем использовать ортонормированный базис функций Уолша, упорядоченных по Адамару, и обозначать его H – квадратная симметричная матрица размера l×l. Вычисление спектральных характеристик входного и выходного воздействий.

Melnikov24.wmf Melnikov25.wmf

где λ, φЭ – вектор-столбцы отсчетов функций λ(t), φЭ(t) соответственно.

Этап 5. Определяем матричные операторы тех звеньев, параметры которых известны (как правило, это параметры объекта и часть параметров системы управления):

Melnikov26.wmf

Melnikov27.wmf

где AИ – матричный оператор интегрирования; I – единичная матрица; θ = 0,05 (коэффициент саморегулирования); для ПТ-25/30-90/10М: Tφ = 12 c, Tξ = 0,5 c. Эти матричные операторы, так же как и спектральные характеристики, рассчитываются один раз и хранятся в памяти ЭВМ.

Этап 6. Определяем матричные операторы звеньев, параметры которых подлежат определению (матричные операторы вычисляются при текущих значениях вектора искомых параметров p = [δ, Tz, Tμ]).

Melnikov28.wmf Melnikov29.wmf

Melnikov30.wmf

Здесь же можно задать начальное приближение для выходного сигнала (или его спектральной характеристики Melnikov31.wmf).

Этап 7. Вычисление матричных операторов нелинейных звеньев. Для этого необходимо предварительно найти спектральные характеристики процессов на входе нелинейных элементов:

Melnikov32.wmf Melnikov33.wmf

Melnikov34.wmf Melnikov35.wmf z = HCz;

Melnikov36.wmf

здесь Melnikov37.wmf – матричный оператор умножения функции представленной в скобках; z – вектор-столбец отсчетов функции z(t). Непосредственное применение формулы (1) для вычисления матричного оператора умножения является далеко не эффективным с точки зрения быстродействия его расчета. Так как система Matlab ориентирована на операции с матрицами, то для расчета матричного оператора умножения эффективно использовать следующую функцию:

function Ay=m_ymn1(H,f)

%Матричный оператор умножения на некоторую функцию %f – вектор-столбец

%H – матрица-базис функций Уолша

Ay=H*(H.*(ones(size(f))*f’)’);

Этап 8. Вычисление матричного оператора всей системы и спектральной характеристики выходного сигнала

Melnikov38.wmf

Melnikov39.wmf

Этап 9. Уточнение матричных операторов нелинейных элементов. Для этого проверяем решения задачи анализа при текущих значениях p: степень близости спектральных характеристик выходного сигнала, полученных на текущей и предыдущей итерациях.

Melnikov40.wmf (2)

Если (2) не выполняется, то Melnikov41.wmf, и переходим к этапу 7, если выполняется, то переходим к этапу 10.

Этап 10. Вычисление искомых параметров:

Melnikov42.wmf (3)

Реализация (3) представляет собой итерационный процесс, в котором параметры меняются в соответствии с выбранным алгоритмом нахождения минимума функции многих переменных. При очередном вычислении целевой функции (3) при конкретных значениях p, следует определять Cφ(p), т.е. возвращаться к этапу 6. Если минимум найден, то переходим к анализу системы при найденных параметрах регулятора.

Этап 11. Анализ системы регулирования турбины при найденных параметрах регулятора:

φ = HCφ; z = HCz; Melnikov43.wmf

μ = HCμ; ξ = HCξ,

где φ, Melnikov44.wmf z, μ , ξ  – вектор-столбцы отсчетов функций φ(t), Melnikov45.wmf z(t), μ(t), ξ(t) соответственно.

pic_25.wmf

Рис. 4. Эталонный и реальный переходные процессы относительного изменения частоты вращения ротора

Результаты синтеза системы регулирования турбины ПТ-25/30-90/10М представлены на рис. 4 при λ = –1, ε0 = 0,0167. Размер базиса составлял 512 функций Уолша. В качестве меры близости спектральных характеристик в выражениях (2) и (3) использовалась евклидова норма. Минимизация целевой функции (3) осуществлялась методом Гаусса – Ньютона, который реализован в среде Matlab с помощью функции lsqnonlin. Эта функция позволяет осуществлять условную минимизацию целевых функций, представляющих собой свертку вектора.

Для системы регулирования турбины ПТ-25/30-90/10М были найдены следующие параметры (ограничения на параметры не накладывались):

δ = 0,03; Tz = 0,0630 с; Tμ = 0,1014 с; (4)

Melnikov46.wmf

Для получения параметров (4) понадобилось 8 итераций. Степень нечувствительности сервомотора составляла ε0 = 0,0167, что соответствует 0,05 % нечувствительности системы регулирования относительно регулируемого процесса φ(t). В работах [1, 3, 4, 5, 14, 15] отражены детальные алгоритмы построения законов регулирования энергетических систем, в том числе и в робастной постановке, а также алгоритмы решения задач идентификации.

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и Правительства Калужской области (гранты № 14-48-03013, № 16-41-400701).