CC BY
48
3
г л 2
1 i" Ч 1 Ч '
J 1 1
1 п 1 1
1 1 * ■ 1 1 1
1 I 1 1 1
1 1
if 1 1 1 1
\ , 1 \ \
.....-х
0 0.2 0.4 0.6 0.0 1
.0. tt
Рис. 2. Относительные временные зависимости изменения: 1) индуктивности обмотки, 2) генерируемого тока и 3) индукции в статоре генератора
СПИСОК ЛИТЕРАТУРЫ
1. Глебов И.А., Кашарский Э.Г., Рутберг Ф.Г. Синхронные генераторы кратковременного и ударного действия. - Л.: Наука, 1985. - 224 с.
Таблица. Расчетные параметры генераторов
Характеристики генераторов Масса генератора, т
0,1 1 10 100
Обороты ротора, об/мин 12000 6000 3000 1500
Период изменения индуктивности Т, мс 1,25 2,5 5 10
Кратность изменения индуктивности 68 110 145 210
Амплитуда тока 1т, кА 40 140 430 1350
Максимум мощности, МВт 6 60 460 3840
Длительность импульса тока, мс 0,25 0,5 1 2
Энергия импульса, кДж 1 20 300 4860
Энергия конденсатора, кДж 0,09 0,9 9 90
Повышение температуры обмотки, °С 1 3 7 16
Таким образом, электромашинные генераторы импульсов тока с периодически изменяющейся индуктивностью можно рассматривать как перспективные источники электромагнитной энергии для питания автономных электрофизических установок сериями импульсов, хотя их разработка, изготовление и применение требует решения многих научных и технологических проблем.
2. Физика и техника мощных импульсных систем / Сб. ст. под ред. акад. Е.П. Велихова. - М.: Энергоатомиздат, 1987. - 352 с.
3. Сипайлов Г.А., Лоос А.В., Чучалин А.И. Электромашинное генерирование импульсных мощностей в автономных режимах. - М.: Энергоатомиздат, 1990. - 168 с.
Установившийся режим электрической сети является вполне оправданной идеализацией, которая широко используется при моделировании систем любой природы. Практически все виды задач анализа, оптимизации, синтеза состояний и структуры электрических сетей связаны с расчетами их установившихся режимов [1].
В детерминированной постановке расчет установившегося режима, в зависимости от постановки задачи, сводится к нахождению значений зависимых параметров X для однозначно заданных независимых параметров У с последующим или одновременным использованием степеней свободы системы для получения физически реализуемого результата.
Таким образом, каждой ветви расчетной схемы ставится в соответствие свое разделение параметров на зависимые и независимые.
Рассмотрим в качестве обобщенной ветви расчетной схемы сети ветвь ¡- с комплексным сопротивлением —, где Д=Ке(—), Х=1т(—¡) и напряжениями Ц, Ц с фазовыми углами 5; и 81 (рис. 1).
ц 4 д-
Н-й—а < "I
Рис. 1. Схема обобщенной ветви расчетной схемы электрической сети
УДК 621.311.016.001.24
УЧЕТ ФИЗИЧЕСКИХ ОСОБЕННОСТЕЙ ПЕРЕДАЧИ ЭНЕРГИИ ПО ЭЛЕМЕНТУ ЭЛЕКТРИЧЕСКОЙ СЕТИ В АЛГОРИТМАХ РАСЧЕТА УСТАНОВИВШИХСЯ РЕЖИМОВ
С.Г. Слюсаренко
Томский политехнический университет E-mail: [email protected]
Исследован предельный по существованию режим передачи энергии по элементу обобщенной ветви расчетной схемы электрической сети. Проведен анализ неоднозначности математической модели ее режимного состояния. Предложены способ проверки существования решения задачи расчета установившегося режима и алгоритм управления итерационным процессом с целью получения физически приемлемых результатов.
Активная P¡j и реактивная Q¡j составляющие мощности — определяются выражениями:
и2 U.U.
P.j =cosT7 cos(ö — ö. + T7 ), (1)
'J у 7 - 7 . J 7 '
U2 U.U.
Q'j = TT sin sin(ö — öj + ),
(2)
X
где T . = arctg R-.
В табл. 1 приведены возможные варианты разделения параметров ветви ¡- на зависимые и независимые. В качестве одного из параметров выступает сдвиг фазовых углов ^.¡¡=8—8.
Таблица 1. Разделение параметров на зависимые и независимые
№ Y X № Y X
1 Uj¡j Pj'Qj 6 UiPijQij Ui Ai
2 UUiP A¡pQj 7 UiAjAj UQ
3 U¡,Uj,Qij ii 8 UjAij,Qj UP
4 UiAj,Pj UAi 9 UiPj Q UiAj
5 UAjQj UP
Фазовый угол напряжения 8¡ определяется на основе векторной диаграммы напряжений:
Q,JR,J - рх„
ö = ö + arctg
U2 — PR — Q.X.. ''<0
' 'J 'J J 'J
Если исключить из математических решений (1, 2) относительно X, приведенных в табл. 1, те, у которых модули напряжений отрицательны, то для первых восьми вариантов разделения параметров на зависимые и независимые решения будут однозначными, а для девятого варианта состава X будет иметь место неоднозначность. Для большинства постановок задач расчета установившегося режима значительная часть ветвей оказывается в условиях, когда У={и,РфО}, поэтому далее рассмотрим эту ситуацию более подробно.
где 1=Щ-Р$гО£цг
Знак минус перед внешним корнем (4) физического смысла не имеет, но математические решения с отрицательным значением и возможны.
Обозначим через к выражение под внутренним корнем (4). Обратимся к двум решениям (4), которые имеют и>0. Назовем первым решение ил, соответствующее знаку плюс перед внутренним корнем (4), тогда второе решение и2 будет соответствовать знаку минус. Эти решения возможны, если к>0. Значения и и и2 становятся кратными при к=0 и решение системы (1, 2) отсутствует, если к<0. Поскольку ветвь ¡-с У={Ц,Рр2^ рассматривается как элемент расчетной схемы электрической сети, то от варианта удовлетворения ее параметрами выражению (4) будет зависеть исход процедуры расчета установившегося режима сети.
Физическая однозначность режимного состояния ветви —
Результатам решения уравнений (1, 2) соответствуют различные входные сопротивления —вх.я и —'ш.й, вещественные и мнимые составляющие которых определяются выражениями:
-и 212
вх.,(1>2)) = '(1,2) сое %, -и 212
1т(Квх. ,(1,2) ) = (1,2) ЯШ % ,
Математическая неоднозначность модели
режимного состояния элемента электрической сети
Проблеме неоднозначности результатов моделирования установившихся режимов посвящено большое количество работ. Основные принципы и результаты исследований отражены в [2-5]. Но эффективных алгоритмов получения правильного решения задачи предложено не было. Рассмотрим причины появления неоднозначности при моделировании режима на примере двухузловой сети.
Рассмотрим выражения (1, 2) как систему уравнений и решим ее относительно и и 8. Исключив 8, получим уравнение относительно и:
и,4 - 2и (2 и2+р^ + )+(Р + Q2 К = °(3)
Решение (3) имеет вид:
U . = ±
U2
—- + P..R.. + Q..X.. ±
2 .- .- ^- .-
U
+ P.R.. 2 . .
(4)
-Q..X..
— (P-+Q-)Z-
где %ц= гтЛдр + я\р<°.
Режим ветви ¡- определяется значениями и и входного сопротивления -вх.; (рис. 2). Значения и, Р., <2. являются следствиями. Поскольку —ш.я и —ш.й физически не могут быть альтернативными, то и1 и и 2 тоже не являются альтернативными.
Традиционно электрическая нагрузка характеризуется значениями потребляемых активной и реактивной мощностей (Рн и 2н) и желаемым значением напряжения нагрузки Ц^.
Рис. 2
Схема замещения — входным сопротивлением
Физически, приемлемое (желаемое) значение напряжения на шинах нагрузки обусловлено техническими условиями эксплуатации токоприемников, а значения Рн и 2н - их результирующим сопротивлением —н, которое является однозначным:
и
р =_нж С08 у ^
и2
а=и^ш
Следовательно,
и
•] п¡р7^]
(5)
и2
^ = 7.' = 7 , .
1, вх. 1 пр.
и2
—^ + Р.' 7.. 2 пр.
С08(у- ) . Р, 7
-= 0,
где х¥2н - аргумент —н.
Поэтому результат расчета режима сети нельзя считать приемлемым, если напряжения на шинах нагрузок существенно отличаются от желаемых для них значений, поскольку он будет соответствовать тогда другим значениям их эквивалентных сопротивлений и, следовательно, другому составу токоприемников.
Следует отметить, что математические модели режимов, в которых электрические нагрузки представлены -н=сош1;, являются линейными и результат расчета режима можно считать полученным, если напряжения на шинах всех нагрузок близки к желаемым значениям. В противном случае будут наблюдаться отклонения мощностей нагрузок модели от заданных значений. В таких моделях остаются актуальными проверка существования решения и введение напряжений на шинах электрических нагрузок в допустимую область.
Параметры предельного по существованию режима
Обратимся вновь к выражению (4) с тем, чтобы выяснить, при каких условиях достигается предельный по существованию режим для ветви I-}.
Значения Р} и для которых выражение под внутренним корнем (4) обращается в нуль, обозначим пр. и <2} пр. и назовем предельными по существованию, а соответствующие им значения Ц обозначим Ц кр. и назовем критическими.
Из (4) следует, что для заданных значений Ц и 2 предельному по существованию режиму ветви —
соответствует равенство:
и' =
— + Р1] прЛ + а1] пр Д] = >% прЛ •
из которого
р]
и 2 С0в ф,
^4=^ 27//(81впР;/ -С08(Тг -ф/))
. (6)
Для р<0 выражение (6) идентично формуле предельной мощности [6].
Из (6) и (5) следуют выражения для пр. и Ц пр.:
и 2
1] пр. ф,, = С0ПЙ
\22/.(81ёпРР/ -С08(У2., -ф/))
и1
У
и 2 • Р
2^Р/, - С08(Тг -ф/,))
(7)
Из (6) и (7) следует, что р пр^/р пр.!!р<0 и Ц кр.!р}>0>! Ц кр.!!р}<0, а при ЧггФц и Р}>0 Рпр =да и Ц кр.'= да.
Режимные характеристики обобщенной ветви
Как отмечено ранее, режим ветви — для заданных Ц, предопределяется значением входного сопротивления 2вх. ; (рис. 2). Задавшись неизменным значением ф} и изменяя 2,х. ; от 0 до да можно рассмотреть множество возможных режимных состояний как для нагрузочного, так и для генераторного характера и получить изображения зависимостей р.(Ц), <(Ц), р(Ц), еда ,(Ц) и т.п.
В данном случае представляют интерес зависимости Р}( Ц), 2Ж. ( Ц). Эти зависимости можно также получить, используя значения Ц, вычисляемые по [4] для ф}=сош1 и р, изменяющегося в диапазоне от предельного значения нагрузочного, до предельного значения генераторного режимов. На рис. 3 представлены зависимости Р}( Ц) и 2вх.;(Ц) при 1ёф}=0,18. Кривые 1, 1' и 2, 2' показывают характер изменения входных сопротивлений 2,х. ; и соответствующих им потоков мощности Р} в генераторном режиме соответственно в областях 2,х. ;<2 и 2,х. >2}.
Аналогичные зависимости представлены кривыми 3, 3' и 4, 4' для нагрузочного характера Р}.
Поскольку для реальных сетей имеет достаточно узкий диапазон значений, то аналитическую оценку предельного режима лучше всего получить для условий Ц=сош1:, ф}=сош11.
После замены на Р^ёф} выражение (4) принимает вид:
и,
и2 С08(Т + Р 7 -
2 11 11
-ф„)
ф,,
( и2 С08(Т7 -±. 1-2- + РЛ-^
II 2 С0ф
ф,)
Р111
. 1] 1] С0^2 ф,
Тогда предельному по существованию режиму будет соответствовать равенство:
Р,
пр. 1
%
пр. 2
Р, г \ 2
/ ' / 1' 3 , / ' 1 4 У 1 2'
^ —------ _ V- - — ' / ( и
и ^ икр. 2 и = и
3' 4'
Рис. 3. Режимные характеристики обобщенной ветви
и
Таблица 2. Зависимости Рг/пр. и Ц от
Р/<0 Р/>0
Рпр., МВт Ц,кр., кВ Р/ пр., МВт Цкр., кВ
0,5 -108,6 63,3 513,5 137,6
0,2 -137,7 68,0 343,7 107,5
0 -161,3 73,0 264,7 93,4
-0,2 -187,7 79,4 206,4 83,3
-0,5 -229,0 91,9 147,3 73,7
Наиболее простой проверкой существования решения является выполнение для каждого узла г расчетной схемы условия (рис. 4):
Напряжения Цкр1 и Цкр.2 и мощности Рпр1 и Рпр.2 соответствуют предельным по существованию режима значениям Р/. Вид режимных характеристик, приведенных на рис. 3, существенно зависит от \%ф/. В табл. 2 для демонстрации этого факта приведены предельные значения Рг/ и Ц кр. для ветви с —=8+/32, Ом, Ц=115 кВ для некоторых значений ^щ.
Контроль отсутствия результата расчета
установившегося режима
Нередко расчетчик режимов сталкивается с тем, что программа выдает сообщение о том, что решение задачи не получено из-за отсутствия сходимости итерационного процесса. При этом ему не указывается причина такого результата.
Характер итерационного процесса расчета установившегося режима определяется многими факторами. Но самым важным из них является некорректность исходных данных, в результате которой для одной или нескольких ветвей схемы сети формируется требование передать энергию мощностью, превышающей предельное значение. Поскольку решение задачи в такой ситуации отсутствует, то расчетчику желательно указать те элементы, из-за параметров которых режим не может существовать. Проанализировав сложившуюся ситуацию, он может выйти из нее, имитируя работу регуляторов напряжения в сторону увеличения предельных значений передаваемых мощностей или исправляя некорректные данные. Достаточно часто некорректные ситуации возникают при имитации вариантов резервирования питания или развития электрических сетей.
Для выявления отсутствия решения задачи удобно использовать узловые эквиваленты (рис. 4) и, при необходимости, редактировать параметры в нужном направлении. Параметры узловых эквивалентов на итерациях должны пересчитываться.
Л
а
Рис. 4.
У
> У
(8)
Схемы обобщенного узла сети и его эквивалента: —, - нагрузочная или генераторная мощность узла г, —0 - результирующая проводимость узла г на шину нулевого потенциала, Ц и — - эквивалентные параметры ветвей, примыкающих к узлу г
Удерживание итераций в области притяжения к правильному математическому решению задачи расчета установившегося режима
Расчет установившегося режима является базовой задачей для решения всего комплекса других не менее важных задач, связанных с оценкой экономичности и режимной надежности функционирования электрических систем. При использовании нелинейных моделей установившихся режимов существует опасность получения неправильных начальных условий для их решения и тогда полученные результаты можно рассматривать как опасные для принятия решений в управлении режимами. В связи с этим получению правильных или физически реализуемых результатов расчета установившихся режима электрических систем уделяется большое внимание [2-5].
Наиболее эффективное решение этой проблемы, по мнению автора, является управления итерационным процессом в направлении получения правильного результата. В алгоритме решения этой задачи предлагается использовать следующие положения:
1. Если для всех элементов электрической сети значения параметров физически приемлемые, то они физически приемлемы и для сети в целом.
2. Результат моделирования установившегося режима электрической системы является приемлемым, если выполняется первое условие и напряжения на шинах нагрузок допустимо близки к целесообразным (желаемым) значениям. При этом, токовые перегрузки элементов сетей контролируются, но не фигурируют в качестве основных нарушений ограничений.
3. Получение правильного результата расчета установившихся режимов должно обеспечиваться начальной ориентацией и удерживанием параметров модели на итерациях с использованием степеней свободы в области притяжения к физически приемлемому решению.
Важно отметить, что многочисленные проверки неоднозначных решений показали, что только в одном, физически приемлемом, решении параметры всех ветвей удовлетворяют выражению (4) со знаком плюс, перед внутренним корнем, а у остальных имеются одна или несколько ветвей, параметры которых удовлетворяют тому же выражению со знаком минус.
Независимо от для режимных состояний с 2,х. ¡<2!/ максимальное значение Ц равно Ц кр. Это обстоятельство позволяет использовать следующий формальный критерий для проверки физической реализуемости параметров ветвей расчетной схемы электрических сетей на итерациях расчета
V1
установившихся режимов: если гвх 1 = —— > , то
и
параметры модели режимного состояния ветви должны удовлетворять выражению (4) со знаком плюс перед внутренним корнем.
Многократные экспериментальные исследования подтвердили состоятельность указанной выше концепции. Для контроля физической реализуемости вектора зависимых переменных на итерациях оказалось достаточным обеспечение физической реализуемости параметров для эквивалентов (рис. 4) всех расчетных узлов.
В нормальных условиях функционирования электрических сетей все элементы, кроме шунтов на шину нулевого потенциала, имеют режимные состояния, близкие к холостому ходу и, следовательно, их параметры находятся в области 2,х. >2. Реальное попадание параметров одной или нескольких ветвей в область 2*. ¡<2 связано с существенными изменениями режима. Такое может быть с протяженной линией электропередачи, попавшей в условие устойчивого короткого замыкания в питаемой сети, в результате чего ее эквивалентное входное сопротивление может стать по модулю меньше продольного сопротивления. При моделировании режимов схем с предварительным преобразованием сильно загруженных связей, соединенных в звезду, в эквивалентный треугольник элементы треугольника также могут оказаться в условиях значительного уменьшения входного сопротивления.
При переходе из области режимных состояний с 2,х. ¡<2 в область 2,х. >2/ вблизи предела по существованию происходит смена знака регули-ВРГ
рующего эффекта -дЦ-■ С уменьшением Р9 соотношения регулирующих эффектов зависят от
Изложенные обоснования необходимости контроля условий существования решений уравнений установившегося режима программно реализованы и протестированы с использованием программного комплекса «Мустанг», применяемого в энергосистемах [4]. В указанном комплексе используется модель установившегося режима в форме уравнений баланса мощностей для узлов сети. Система уравнений решается методом Ньютона-Рафсона и по параметру. Тестирование позволило выявить
следующее. При задании исходных данных предельного и за пределом существования режима, программа «Мустанг» выдает сообщение об отсутствии сходимости итерационного процесса. Тестируемая программа указывает, в каких узлах сложилась ситуация противоположная (8), и предлагает результат расчета с измененными значениями —¡.
При вариациях значений вектора начальных приближений зависимых переменных X методы Ньютона-Рафсона и по параметру демонстрировали неоднозначность решения уравнений установившегося режима. Тестируемая программа в этих же условиях обеспечивала один результат, параметры которого удовлетворяли условию 2х. >2 для всех ветвей расчетной схемы. Тестированная программа 1пёогС1гсш\ реализует модификацию двухэ-тапного метода [5] и внедрена в филиалах «Северные электрические сети», «Центральных электрические сети» ОАО «КузбассЭнерго», в электрических сетях гг. Томск, Северск, Абакан и ряда крупных промышленных предприятий.
Заключение
Управление итерационным процессом расчета установившихся режимов электрических сетей с целью получения физически приемлемых решений должно осуществляться с использованием степеней свободы системы (независимых переменных). Также можно использовать дополнительную процедуру введения параметров модели в допустимую область.
Моделирование установившихся режимов электрических сетей должно опираться на причинную основу физических явлений, происходящих в системе-оригинале. Контроль отсутствия результата моделирования установившегося режима должен быть обязательной процедурой, поскольку она существенно упрощает поиск причины отсутствия сходимости итераций.
Изложенные способ диагностики выхода параметров модели электрических сетей из области существования режима и алгоритм управления итерационным процессом расчета установившегося режима с целью получения правильного результата моделирования режима применимы также и для трубопроводных инженерных сетей, т.к. модели их режимов базируются на принципе подобия электрических сетей постоянного тока.
СПИСОК ЛИТЕРАТУРЫ
1. Федин В.Т., Латышева И.Г. Анализ эксплуатационных программ расчета установившихся режимов энергосистем // Известия вузов. Энергетика. - 1973. - № 9. - С. 3-8.
2. Идельчик В.И., Лазебник А.И. Аналитическое исследование сходимости решения уравнений установившегося режима электрических систем. // Известия АН СССР. Энергетика и транспорт. - 1972. - № 3. - С. 47-56.
3. Борисов Р.И., Готман В.И., Слюсаренко С.Г. Экспериментальное исследование сходимости и неоднозначности решения
уравнений установившегося режима электрических систем // Известия вузов. Энергетика. - 1978. - № 9. - С. 14-20.
4. Идельчик В.И. Расчеты и оптимизация режимов электрических сетей и систем. - М.: Энергоатомиздат, 1988. - 288 с.
5. Слюсаренко С.Г. Модификация алгоритма расчета установившихся режимов электрических систем на основе поочередного уточнения потокораспределения и напряжений // Процессы и режимы электрических систем. Межвузовский научно-технический сборник. - Томск: Изд-во ТПИ, 1990. - С. 37-43.
6. Зевеке Г.В., Ионкин П.А., Нетушил А.В., Страхов С.В. Основы теории цепей. - М.: Энергоатомиздат, 1989. - 528 с.
