CC BY
17
УДК 621.311.001
Построение математической модели допустимого установившегося режима электроэнергетической системы
Н.П. Бадалян1, Е.А. Чащин2 1 ФГБОУВПО «Владимирский государственный университет имени А.Г. и Н.Г. Столетовых»,
Владимир, Российская Федерация 2 ФГБОУВПО «Ковровская государственная технологическая академия имени В.А. Дегтярёва»,
Ковров, Российская Федерация E-mail: norayrbadalyan@mail.ru; kanircha@list.ru
Авторское резюме
Состояние вопроса: Существующие методы расчета установившегося режима, основанные на решении системы нелинейных уравнений итерационным методом, в случае выполнения анализа электроэнергетической системы отличаются большим объемом производимых вычислений и существенными затратами машинного времени. Материалы и методы: Создание математической модели для расчета допустимого установившегося режима основано на Z-форме задания соответствия сети и введении ограничений типа неравенств, налагаемых на модули напряжений и реактивных мощностей независимых станционных узлов.
Результаты: Предложена математическая модель для анализа допустимого установившегося режима в большой электроэнергетической системе. Разработаны рекомендации по взаимному переходу от станционных узлов типа P-U к станционным узлам P-Q и обратно, обеспечивающие соответствие модулей напряжений и реактивных мощностей допустимым пределам, направленные на достижение допустимого установившегося режима электроэнергетической системы.
Выводы: Использование предложенной модели обеспечивает уменьшение числа итераций и повышение сходимости при решении системы уравнений, описывающей установившийся режим в электроэнергетической системе.
Ключевые слова: модель, мощность, параметр, станция, узел, матрица, нагрузка.
Mathematical Model Designing of Steady State Mode of Power Engineering System
N.P. Badalyan1, E.A. Chashchin2 Vladimir State University, Vladimir, Russian Federation 2Kovrov State Technological Academy, Kovrov, Russian Federation E-mail: norayrbadalyan@mail.ru; kanircha@list.ru
Abstract
Background: The existing calculation methods of steady-state mode, based on solving the system of nonlinear equations by using the iteration method in case of analysis of an electric power system, differ by very big volume of produced computing and significant costs of computer time.
Materials and methods: Designing the mathematical model for the calculation of the possible set mode is based on the Z form of accordance task of a network and reservation restriction type of inequalities, imposed on the modules of voltage and reactive powers of the independent station knots.
Results: The mathematical model for the analysis of the possible set mode in the large electro-power system is offered. Recommendations on the mutual transition from the station knots of P-U type to the P-Q station knots and back, providing accordance of the voltage modules and reactive powers of the possible limits, turned to the proper achievement of the possible set mode of the electro-power system.
Conclusions: The usage of the offered model helps to diminish the number of iterations and to increase the convergence when solving the equalizations system, describing the set mode in the electric power system.
Key words: мodel, capacity, parameter, station, knot, matrix, loading.
Основой для построения математической модели допустимого режима является предложенная ранее математическая модель расчета установившегося режима [1-9]. Однако решение задач, связанных с расчетом допустимого установившегося режима применительно к электроэнергетической системе (ЭЭС), традиционным путем, заключающимся в итерационном решении системы нелинейных алгебраических уравнений, связано с большими затратами машинного времени и рядом
других затруднений, вызванных большим объемом производимых вычислений.
Ниже для построения уравнения установившихся режимов принимается, что рассматриваемая ЭЭС состоит из источников (станционных) и потребителей (нагрузочных узлов), число которых составляет (М+1). После выбора первого источника в качестве базисного (балансирующего) ЭЭС будет состоять из М независимых узловых точек. Необходимо отметить, что в случае, когда к данному узлу
подключены как источник, так и потребитель, всегда можно их объединить и получить случай, когда к узлу подключен либо только источник, либо только потребитель.
При построении системы уравнений для расчета допустимого установившегося режима ЭЭС предполагается известной конфигурация исследуемой ЭЭС со всеми пассивными параметрами. Под допустимым установившимся режимом ЭЭС будем понимать такое состояние режимных параметров ЭЭС, при котором они удовлетворяют условиям типа неравенств, налагаемых на них. При этом условия типа неравенств формируются из допустимых предельных значений (нижних и верхних) отдельных режимных параметров. В этом случае, как обычно, ограничения налагаются на модули напряжений и реактивных мощностей станционных узлов.
Для изложения материала выбираем следующую систему индексов (рис. 1):
1. Для независимых станционных узлов т(п) = 1, 2, ..., Г, где Г - число независимых станционных узлов.
2. Для нагрузочных узлов к(1)=Г+1, Г+2, Г+3, ..., Г+Н, где Н- число нагрузочных узлов.
3. Для произвольных узлов /®=1, 2, ..., Г, Г+1, Г+2, ., Г+Н = М, где М - число независимых узлов в рассматриваемой ЭЭС.
Ф ■ = Р -^р/ Г1
Г
Г ( соэ ¥и, + О, эт Уи1 )
м
+ Х 1В 1 )
1=1
= о,
Фф = О/ -
г
и ( ( ¥и/ - О/ с°э ¥и/)-
( - я/1в/1)
1=1
= о,
(1)
где
А = ГГ [( + О/О1) ( - ) +
+(Р - 0Р ) ( -¥и)]; (2)
в=гг К+ОО1 )'п (-¥и )-
-(Р - О/р ) соэ( - ¥щ )].
В этом приближении математическая модель допустимого установившегося режима ЭЭС может быть получена на основании известной системы уравнений [1] при 2-форме задания соответствия сети.
Представим систему нелинейных алгебраических уравнений (1) в неявно выраженной форме:
ФР1 (Р,ОГ, ¥и ) = 0,1
Фф- (ОГ, ¥и ) = 0.}
Согласно выбранной системе индексов для независимых станционных т(п) и нагрузочных узлов к(1) , систему нелинейных алгебраических уравнений (3) представим как совокупность систем нелинейных алгебраических уравнений для станционных и нагрузочных узлов, в данном случае:
Фрт (Р,ОГ, ¥и ) = 0,1
ФЯт (,ОГ, ¥и) = 0,_[
Фрк (Р,ОГ, ¥ и ) = 0,|
Фяк (ОГ,¥и) = 0.[
Система неявных нелинейных алгебраических уравнений (4) относится к станционным узлам, а (5) - к нагрузочным узлам.
(3)
(4)
(5)
Рис. 1. Схематическое представление ЭЭС
Введем следующие обозначения:
Фрт (,О,и,Т) = 0,
рт
ФЯт (,О,и, Т ) = 0;
(6)
Фк =
(7)
Фрк (,О,и,¥и) = 0,
Фяк (Р,0,и, ¥и ) = 0.
При этом математическую модель допустимого установившегося режима ЭЭС удобно представить в следующем виде:
Фт (о,и,¥и) = 0;
Фк (,О,и,^) = 0;
(8)
(9)
(10)
(11)
ит,т'т — ит — ит,тах ;
°т,т1п — От — От,тах ,
где Um,min, От.т/п и ит,тэх, От.тэх — предельные минимальные и максимальные значения модулей напряжений и реактивных мощностей независимых станционных узлов.
Предположим, что все узлы являются узлами типа Р-О, тогда уравнения (8) и (9) изображают полную систему нелинейных алгебраических уравнений установившегося режима ЭЭС. Необходимо отметить, что при станционных узлах типа Р-и структура системы нелинейных алгебраических уравнений (8) и (9) не изменяется.
Если узлы являются типа Р-О, то это означает, что для каждого узла заданы активные и реактивные мощности и поэтому в конечном итоге уравнения (8) и (9) являются функцией от модулей и аргументов комплексных напряжений. Учитывая, что (8) и (9) получены из системы уравнений (1) и введены обозначения (6) и (7), можем написать следующее рекуррентное выражение на основании метода Ньютона-Рафсона:
"т " 1 ит И+1 "Тит ] И Г АТ ит ]
т с Тик АТик
ик ик Аик
ит ит Аит
(12)
где И - номер итерации.
Вторая столбцевая матрица правой части определяется на основании следующего матричного уравнения:
т 0°" <12 т В* д дФ рт дФрт
дТ ип дТи ди дип
АРт' дФрк дФрк дФрк дФрк
АРк дТ ип дТ и; ди і дип
АОк дФЯк дФЯк дФЯк дФЯк
АОт дТ ип дТ иі ди дип
дФЯт дФЯт дФЯт дФЯт
дТ ип дТи ди дип
АТ
\ л. ,,
АТ
ик
Аик
Аи„
(13)
где ДРт, ДОт - приращения активных и реактивных мощностей независимых станционных узлов; ДРк, ДОк - приращения активных и реак-
тивных мощностей нагрузочных узлов; ДТит, Дит - приращения аргументов и модулей комплексных напряжений независимых станционных узлов; ДТик, ДЦк - приращения аргументов и модулей комплексных напряжений нагрузочных узлов.
Введем следующие обозначения:
дФ,
рт
дФ
рт
дФ
рт
дФ
рт
тп рш
ит ип
а 9ФРк
кП д¥ ип
дТ,,
а 8ФРк
к1 дТ ■ ’ ит иі
т ди
дФРк
к1 ди; ’
дип
с 8ФРк
Скп дип '
ькп =
дФ
Як
Ьтп
дТ
ип
дФ
_ ич^ят
ькі =
дФ,
Як
дТ ■ иі иі
дФ
dki =■
дФ
Як
диі
дФ
dkn =
дФ
Як
дип
дТ
и 1
ь . = Ят d .=—ят
ті ’ ті я/1
дТиі диі
дФ
, =и ^ят
тп , ■
дип
При этом матричное уравнение (13) можно представить в следующем виде:
АРт
АРк
АОк
АОт
атп аті сті стп
акп акі скі скп
Ькп Ькі dki dkn
Ьтп Ьті <^ті <^тп
АТ
I—\ л. и
АТ
ик
А
Аи„
(14)
После обращения матрицы Якоби выражения (14) получим
АТ
і—\ л. її
АТ
ик
Аик
Аи„
атп аті Рт/ Ртп
акп акі Ркі Ркп
Укп Укі 8кі 8кп
У тп У ті 8ті 8тп
АРт
АРк
АОк
АОт
(15)
Выражения частных производных, входящих в вышеприведенные матрицы Якоби, определяются на основании аналитических выражений (4) и (5).
Предположим, что станционные узлы являются узлами типа Р-и, т.е. относительно этих узлов зафиксированы активные мощности и модули комплексных напряжений, тогда имеет место следующее соотношение:
Ди1 =Ди2 =... = Диг = 0,
или в матричной форме:
[Дит ] = 0. (16)
Тогда матричное уравнение (15) можно представить в виде двух подматричных уравнений:
АТ
1Л 1 и
АТ
ик
Аи„
атп аті Рті Ртп
акп акі Ркі' Ркп
Укп Укі 8кі 8 кп
АРт
АРк
АОк
АОт
(17)
[0] = [тп і У ті і 8ті і 8т
АРт
АРк
АОк
АОт
(18)
Из математического выражения (18) можно определить вектор реактивных мощностей независимых станционных узлов:
ат =
стп =
[АОт ] = -[8тп ]-1 х[Ат/ ], (19)
где
[Ат/] = [Ттп]х [АРт] + [тт!]х [АРк] + [8т/] х [А0к]■
На основании (19) можно установить численные значения приращений реактивных мощностей независимых станционных узлов:
[АОт] = [ [ А02 | ■■■ | АОг], (20)
где t - знак транспонирования.
Определяя численные значения приращений реактивных мощностей независимых станционных узлов, можно определить также их действительные значения, пользуясь вторым векторным уравнением (4). Определяя численные значения реактивных мощностей независимых станционных узлов, проверяется условие их допустимости согласно (11).
При проверке условия допустимости реактивных мощностей могут возникать следующие случаи:
1. Условие (11) полностью обеспечивается.
При этом непосредственно необходимо перейти к расчёту установившегося режима. На основании (17) вычисляются численные значения приращений АТит, А%к и Аик, затем, пользуясь рекуррентным выражением
т 1 um И+1 Tum" И AT um
т с Sr = Tuk - k u is
| k A k Sr 1
(21)
вычисляются их действительные значения.
Следует учитывать, что перед каждой итерацией необходимо проверять условие (11), и итерационный процесс считается завершенным, если искомые режимные параметры принимают желаемые численные значения с точки зрения их точности.
2. Условие (11) обеспечивается не полностью.
Заметим, что нарушение может быть двояко: реактивная мощность данного независимого станционного узла с индексом t окажется больше, чем верхний предел, или меньше, чем нижний предел. В обоих случаях данный станционный узел типа Рги{ заменяется на узел типа Р -ф, причем в первом случае принимается 01=0,тах, тогда как во втором случае принимается 01=01,тп.
Для таких станционных узлов искомыми переменными являются модули и аргументы комплексных напряжений. При этом перед каждой итерацией проверяется условие допустимости модулей комплексных напряжений согласно (10). Если нарушается условие допустимости модулей напряжений, то аналогичным образом из станционного узла типа Р-0 переходим к узлу типа Р-и.
Таким образом, путем перехода от станционных узлов типа P-U к станционным узлам P-Q и наоборот можно добиться получения допустимого установившегося режима ЭЭС, при котором модули напряжений и реактивные мощности находятся в пределах допустимости.
Список литературы
1. Хачатрян В.С., Аль-Исса Абрахим. Оптимизация режима электроэнергетической системы по реактивным мощностям прямым методом классического программирования при Z-форме задания состояния сети // Изв. НАН и ГИУ Армении. - 1977. - № 2. - С. 127-135.
2. Хачатрян В.С., Этмекчян Э.А. Метод расчета установившегося режима электрической системы // Изв. вузов СССР. Энергетика. - 1989. - № 5. - С. 12-18.
3. Хачатрян В.С., Хачатрян С.Ц., Сафарян В.С. Расчет установившегося режима электрических систем с применением матрицы Гессе при Z-форме задания состояние сети // Изв. вузов СССР. Энергетика. - 1990. -№ 1. - С. 20-23.
4. Хачатрян В.С., Этмекчян Э.А. Развитие гибридного метода расчета установившегося режима электрической сиситемы // Электричество. - 1991. - № 1. - С. 6-13.
5. Хачатрян В.С., Этмекчян Э.А., Аракелян В.П. Упрощенный метод расчета установившегося режима электроэнергетической системы // Электричество. - 1992. -№ 2. - С. 9-14.
6. Аракелян В.П., Хачатрян К.В., Аль-Дарвиш М.Б. Об одном методе расчета установившегося режима электроэнергетической системы // Изв. НАН и ГИУ Армении. Серия технических наук. - 1996. - № 1. - С. 17-22.
7. Бадалян Г.А. Об одном методе определения Z-матрицы обобщенных параметров электрических систем // Изв. НАН и ГИУ Армении. Серия технических наук. - 1996. -№ 1. - С. 11-16.
8. Бадалян Н.П., Хачятрян К.В. Новый метод определения обобщенных параметров установившегося режима электроэнергетической системы: сб. докл. II Меж-дунар. энергетической конф. в Армении. - Ереван, 2001. -С. 400-408.
9. Бадалян Н.П. Метод построения обращённой формы уравнений установившегося режима электроэнергетической системы // Изв. вузов и энергетических объединений стран СНГ. Энергетика. - 2004. - № 6. - С. 45-50.
References
1. Khachatryan, V.S., АГ-^Ба, Abrakhim. Izvestiya NAN i GIU Armenii, 1977, no. 2, рр. 127-135.
2. Khachatryan, V.S., Etmekchyan, E.A. Izvestiya vuzov SSSR. Energetika, 1989, no. 5, рр. 12-18.
3. Khachatryan, V.S., Khachatryan, S.Ts., Safaryan, V.S. Izvestiya vuzov SSSR. Energetika, 1990, no. 1, рр. 20-23.
4. Khachatryan, V.S., Etmekchyan, Е.А. Elektrichestvo, 1991, no. 1, рр. 6-13.
5. Khachatryan, V.S., Etmekchyan, E.A., Arakelyan, V.P. Elektrichestvo, 1992, no. 2, рр. 9-14.
6. Arakelyan, V.P., Khachatryan, K.V., Al'-Darvish, M.B. Izvestiya NAN i GIU Armenii, seriya tekhnicheskikh nauk, 1996, no. 1, рр. 17-22.
7. Badalyan, GA. Izvestiya NAN i GIU Armenii, seriya tekhnicheskikh nauk, 1996, no. 1, рр. 11-16.
8. Badalyan, N.P., Khachyatryan, K.V. Novyy metod opredeleniya obobshchennykh parametrov ustanovivshegosya rezhima elektroenergeticheskoy sistemy [The New Method of Defining Generalized Parameters of Steady-State Mode of Power Engineering System]. Sbornik dokladov II Mezhdu-narodnoy energeticheskoy konferentsii v Armenii [Reports Collection of the II International Engineering Conference in Armenia]. Erevan, 2001, pp. 400-408.
9. Badalyan, N.P. Izvestiya vuzov i energeticheskikh ob»edineniy stran SNG. Energetika, 2004, no. 6, pp. 45-50.
Бадалян Норайр Петикович,
ФГБОУВПО «Владимирский государственный университет им. Александра Григорьевича и Николая Григорьевича Столетовых»,
доктор технических наук, профессор кафедры электротехники и электроэнергетики, e-mail: norayrbadalyan@mail.ru
Чащин Евгений Анатольевич,
ФГБОУВПО «Ковровская государственная технологическая академия им. В.А. Дегтярёва», кандидат технических наук, доцент, зав. кафедрой электротехники, e-mail: kanircha@list.ru телефон (49232) 3-20-62.
