1 Расчет установившихся режимов Математические модели Методы Санкт-Петербург

Описание презентации 1 Расчет установившихся режимов Математические модели Методы Санкт-Петербург по слайдам

1 Расчет установившихся режимов Математические модели Методы Санкт-Петербург НОУ ЦПКЭ, СПб. ГПУ доц. ,1 Расчет установившихся режимов Математические модели Методы Санкт-Петербург НОУ ЦПКЭ, СПб. ГПУ доц. , к. т. н. Шаргин Ю. М.

2 Линии коэффициент распространения электромагнитной волны волновое сопротивление Параметры П-схемы замещения L до ~3002 Линии коэффициент распространения электромагнитной волны волновое сопротивление Параметры П-схемы замещения L до ~300 км: L >300 км: Yi j Yj Zij i ly 2 1 YY 0 ji lz. Z 0 ij )lsh(ZZ 0 вij 2 l th Z 1 YY 0 в ji 0 0 в y z Z 000 jbgy 000 jxrz 000 yz

3 Трансформаторы  j i  Ui Uj Uj Y Zij  kij =3 Трансформаторы j i Ui Uj Uj Y Zij kij = Uj Ui Ii Ij j Ui Uj Yj Zij Ij i Yi Y Iii Ijj jbg. Y 3 2 ][. . 10 к. Вном. ВН к. Втхх U P g 2 ][. 100 % к. Вном. ВН АМВномхх m U SI b Ом S UP r АМВном к. Вном. ВНк. Втзк T 3 2 2. . 10 ][22 Омrzx. ТТТ ][ 100 %2 Ом S Uu z ном отвk Т 100 1. %U n. UUотв отвном. ВНотв ном. НН отв ij U U k.

4 Параметры пассивной П-схемы замещения Г-схема и П-схема эквивалентны,  если. Параметры П-схемы определяются4 Параметры пассивной П-схемы замещения Г-схема и П-схема эквивалентны, если. Параметры П-схемы определяются через параметры Г-схемы. Д ействительный k ij , тогда: П-схема позволяет заменить две или более параллельно включенные трансформаторные ветви с различными коэф-фициентами трансформации одной эквивалентной ветвью. ij jiji j 2 ij ij ij j i i Z Uk. U I, k Z k U U I jj ij ji j ii ij ji i YU Z UU I ‘ , ‘ i j ij ij k Z Z)1( 1 )1( ij ij ij ik Z k Y )k 1( Z 1 Z 1 Yij ijijij j Пj. Гj Пi. Гi II II. .

5 Активная П-схема замещения трансформаторной ветви Если коэффициент трансформации комплексный (при попе-речном регулировании напряжения5 Активная П-схема замещения трансформаторной ветви Если коэффициент трансформации комплексный (при попе-речном регулировании напряжения автотрансформаторов), вращение векторов напряжений осуществляется за счет воздействия задающих токов: Трехобмоточные трансформаторы и автотрансформаторы имеют трехлучевую схему замещения. Ветви среднего и низкого напряжений представляют Г-схемой с идеальными трансформаторами или активной П-схемой замещения. То же, трансформаторы с расщепленной обмоткой НН. ijij ij Z k Y )1 k(YYijiji ij ij j Z )k 1( Y 1 kk YUJ ijij ijji 1 kk YUJ ijij j Ui Uj Yj Yij Ij i Yi Y Ii Ji Jj Iii Ijj

6 Расчет электрического режима ветвей, моделируемых П-схемой замещения  j  Ui Uj Yj6 Расчет электрического режима ветвей, моделируемых П-схемой замещения j Ui Uj Yj Yij Ij i Yi Ji Jj Ii Ijj Iii ijjiij. YUUI)( iiiii UYJI jjjjj UYJI ijiii. III ijjjj. III iii. IU 3 S jjj. IU 3 S jiij. SSS

7 Нагрузка СХН: естественные  с учетом регулирования напряжения  для длительных режимов f7 Нагрузка СХН: естественные с учетом регулирования напряжения для длительных режимов f = f ном U = U крf. Up. Pнг, f. Uq. Qнг, 2. . 1. номнг 2 номнг 0 номнгнг U U aa. PP 2. . 1. номнг 2 номнг 0 номнгнг U U bb. QQ 0, 1 210 aaa 0, 1 210 bbb

8 Типовые СХН 8 Типовые СХН

9 P нг =P нг (U)    н ом. н г 129 P нг =P нг (U) н ом. н г 12 ln 2 e xp~ U U aa. Up

10 Q нг =Q нг (U)   номнг. U U bb. Uq. 1210 Q нг =Q нг (U) номнг. U U bb. Uq. 12 ln 2 e xp~

11 СХН с учетом регулирования напряжения отклонения напряжения от номинального значения в сторону уменьшения11 СХН с учетом регулирования напряжения отклонения напряжения от номинального значения в сторону уменьшения и увеличения, компенсируемые АРНТ и/или действиями оперативного персонала , , , . 2. . 1. UUUесли U UU bb. Q UUUUUесли. Q UUUесли U UU bb. Q Q номнг 2 номнг 0 номнгномнг 2 номнг 0 номнг нг UU ,

12 СХН с учетом регулирующего эффекта нагрузки где К нг - полный регулирующий эффект12 СХН с учетом регулирующего эффекта нагрузки где К нг — полный регулирующий эффект нагрузки где k нг — регулирующий эффект нагрузки при U=const Зависимостью реактивной мощности нагрузки от частоты либо пренебрегают, считая мощность постоянной, либо пересчитывают по формуле (старые руководящие указания по устойчивости энергосистем). а) б) в) номномнг 2 номнг 0 номнгнг f f a U U aa. PP 3 2. . 1. ном нгномнгнг f f KPP 1. f U U P f P Kнгнг нг ном нгнгнг f f k. UPP 1 номном нг. номнгнг нг fff PPP f P k нгнгнг tg. PQ const номнгнг .

13 Математические модели нагрузки а – шунт; б – идеальный источник тока; в –13 Математические модели нагрузки а – шунт; б – идеальный источник тока; в – реальный источник тока; г – реальный источник ЭДС

14 Математические модели нагрузки Шунт: Идеальный источник тока: Реальный источник ЭДС: Преобразование реального источника14 Математические модели нагрузки Шунт: Идеальный источник тока: Реальный источник ЭДС: Преобразование реального источника ЭДС в реальный источник тока и наоборот: нгнгнг нг S UU I U Z 233 23 U S U I Yнгнг нг. YU Z U S 2 2 3 3 U 3 S Jнг нг нгнг. IZUE UYIJнгнгнг UY U S Jнг нг нг 3 EY Z E J Z 1 Y

15 Генераторы PQ -генератор.  Нерегулируемый генератор с заданными значениями генерации активной и реактивной15 Генераторы PQ -генератор. Нерегулируемый генератор с заданными значениями генерации активной и реактивной мощности, не зависящими от напряжения в узле (от режима сети). Свободными параметрами режима генераторного узла, которые определяются расчетом, являются модуль и угол напряжения в узле, т. е. вектор напряжения. PU -генератор. Регулируемый генератор с заданными уставками активной мощности и модуля напряжения в узле. Активная мощность генератора постоянна и не зависит от режима сети. Для поддержания заданного модуля напряжения изменяется реактивная мощность генератора в пределах регулировочного диапазона от Qmin до Qmax. Свободными параметрами режима являются реактивная мощность генератора и угол напряжения узла. При выходе на ограничения по реактивной мощности фиксируется ее граничное значение (изменяется тип генератора на PQmin или PQmax ), свободными параметрами режима становятся модуль и угол напряжения. PU -генераторы участвуют в сведении баланса реактивной мощности в схеме. В зависимости от поведения генераторов в процессе балансирования режима они могут быть регулируемыми и изменять свою активную и/или реактивную мощность, участвуя в сведении/замыкании баланса мощности в схеме, или нерегулируемыми с фиксированными значениями активной и реактивной мощности. Различают следующие типы генераторов:

16 Типы генераторных узлов PEq -генератор.  Регулируемый генератор с заданными уставками активной мощности16 Типы генераторных узлов PEq -генератор. Регулируемый генератор с заданными уставками активной мощности и модуля напряжения в узле. Для этого типа генератора дополнительно должны быть заданы синхронные индуктивные сопротивления по осям d и q. При изменении режима сети изменяются возбуждение генератора (модуль синхронной ЭДС Eq ) в пределах регулировочного диапазона от Eqmin до Eqmax и угол между векторами ЭДС и напряжения узла таким образом, чтобы получить заданные значения активной мощности и модуля напряжения. Свободными параметрами режима генераторного узла, которые определяются расчетом, являются реактивная мощность генератора и угол напряжения узла. При выходе на ограничения Eqmin или Eqmax модуль Eq фиксируется, при этом свободным параметром режима становится модуль напряжения в узле, а реактивная мощность генератора и угол напряжения по-прежнему остаются свободными параметрами режима. Если ток якоря генератора превышает номинальный ток, то, если не достигнуты ограничения по возбуждению, изменяется Eq таким образом, чтобы уменьшить ток якоря до номинального, в противном случае уменьшается активная мощность генератора. PEq -генераторы участвуют в регулировании баланса реактивной мощности в схеме. В отличие от PU -генератора, при регулировании PEq -генератора учитывается взаимосвязь процессов регулирования активной и реактивной мощности, изменение модуля вектора синхронной ЭДС влияет на угол и наоборот.

Типы генераторных узлов 17 Типы генераторных узлов

18 Типы генераторных узлов Векторная диаграмма PEq-генератора)cos(UE q qd q d q. U x18 Типы генераторных узлов Векторная диаграмма PEq-генератора)cos(UE q qd q d q. U x xx E x x E maxmi nqqq.

19 а) шунт: б) идеальный источник тока: в) реальный источник тока: г) реальный источник19 а) шунт: б) идеальный источник тока: в) реальный источник тока: г) реальный источник ЭДС: Математические модели генераторовнгнгнг нг S UU I U Z 233 23 U S U I Yнгнг нг U 3 S JГ Г ГГIZUE UYIJГГГ UY US JГ Г Г

20 Практические методы расчета режимов энергосистем Пусть электрическая схема имеет  r  узлов20 Практические методы расчета режимов энергосистем Пусть электрическая схема имеет r узлов и s ветвей. Если заданы ЭДС всех источников и параметры всех ветвей, то неизвестными являются s токов ветвей. После нахождения токов, можно вычислить напряжения в узлах По первому закону Кирхгофа можно составить r уравнений , но только r-1 из них независимы друг от друга. При расчете электрических режимов используются первый и второй законы Кирхгофа n — число ветвей, присоединенных к узлу m — число ветвей, входящих в контур По второму закону Кирхгофа можно составить s-r+1 независимых уравнений. Законы Кирхгофа позволяют составить систему из (r-1)+(s-r+1)=s независимых уравнений, решив которую можно найти токи в s ветвях. Для уменьшения числа уравнений можно упростить схему, используя эквивалентные преобразования, или использовать специальные методы контурных токов , узловых напряжений. 0 1 n k ki m k keu 11 0 1 n k k. I m k k. EU

21 Метод узловых напряжений Минимальная исходная информация для проведения расчетов: схема замещения электрической сети21 Метод узловых напряжений Минимальная исходная информация для проведения расчетов: схема замещения электрической сети и параметры ее элементов, текущие значения коэффициентов трансформации, отключенные с одной или двух сторон ветви схемы, полная мощность нагрузки в нагрузочных узлах, полная мощность нерегулируемых генераторов, уставки активной мощности и напряжения регулируемых генераторов, располагаемый диапазон изменения реактивной мощности PU -генераторов, зависящий от числа включенных в узле генераторов и их загрузки по активной мощности, или пределы изменения модуля синхронной ЭДС PEq -генераторов. По этому методу составляется система уравнений для определения напряжений узлов относительно одного из узлов схемы, выбранного в качестве опорного узла. Напряжение опорного узла выбирается произвольно и может быть равно нулю. Тогда узловые напряжения являются фазными напряжениями. Определив искомые узловые напряжения, которые являются фазными напряжениями, можно найти напряжения на ветвях, токи ветвей, генерацию, потребление, перетоки мощности, потери в сети и другие величины, необходимые для анализа режима.

22 Схема замещения электрической сети Трансформаторы представлены П-схемой, тогда Схема имеет пять  узлов.22 Схема замещения электрической сети Трансформаторы представлены П-схемой, тогда Схема имеет пять узлов. Узел нулевого потенциала (“земля”) имеет номер 0. 2 4 1 3 Y 24 Y 34 Y 13 Y 23 Y 12 JГ 1 Y 20 Y 40 Y 10 Y 30 JГ 3 JГ 2 JГ 4 Jнг 1 Jнг 2 Jнг 3 Jнг 4 43344224322331132112 , , YYYYY

23 Формирование системы узловых уравнений где, 0)()(1110131312121 нг. ГJJUYUUYUUY , 0)()()(22220242423232121 нг. ГJJUYUUYUUYUUY ,23 Формирование системы узловых уравнений где, 0)()(1110131312121 нг. ГJJUYUUYUUY , 0)()()(22220242423232121 нг. ГJJUYUUYUUYUUY , 0)()()(33330343432323131 нг. ГJJUYUUYUUYUUY. 0)()(4444043434242 нг. ГJJUYUUYUUY 4444343242 3434333232131 2424323222121 1313212111 , , , JUYUYUYUY JUYUYUY , 10131211 YYYY , 2024232122 YYYYY , 3034323133 YYYYY 40434244 YYYY iнгi. Гi. JJJ

24 Постановка задачи расчета установившегося режима Система узловых уравнений является нелинейной системой уравнений, поскольку24 Постановка задачи расчета установившегося режима Система узловых уравнений является нелинейной системой уравнений, поскольку результирующие задающие токи зависят от напряжений, которые нужно найти (слайды 14 и 18). Если в узлах схемы включены нелинейные шунты, их параметры тоже зависят от режима сети. Нелинейная система может быть решена только методом последовательных приближений. На первый взгляд кажется, что для получения решения достаточно в качестве начальных значений напряжений узлов задать номинальные напряжения и ориентировать вектора напряжений по оси вещественных комплексной плоскости, задать активную и реактивную мощности нагрузки, активную и реактивную мощности генераторов задать равной мощности нагрузки с учетом примерной оценки потерь мощности в сети, и применить какой-либо метод решения системы нелинейных уравнений. Итерации прекращаются по достижении заданной точности сведения балансов токов. При такой постановке задачи итерационный процесс или расходится, или дает решение, как правило, не соответствующее условиям задачи, т. к. система нелинейных уравнений может иметь несколько решений. Для получения решения, имеющего смысл, необходимо учесть физические свойства объекта, для которого составлена система уравнений.

251.  Чтобы получить напряжения узлов, соответствующие условиям задачи,  необходимо в ходе итеративного251. Чтобы получить напряжения узлов, соответствующие условиям задачи, необходимо в ходе итеративного решения системы уравнений регулировать напряжение в узлах схемы. Уже давно установлено, что представление генерации в схеме только PQ -генераторами препятствует получению заведомо существующего решения и является причиной неоднозначного решения нелинейной системы узловых уравнений или нарушения условий существования решения. Изменение модуля вектора задающего тока PQ -генератора в пределах диапазона допустимых значений при изменении напряжения узла не контролируется и не ограничивается, а по углу этот вектор жестко привязан к вектору напряжения. Поведение PQ -генератора в итеративном процессе не соответствует поведению реального генератора при изменении режима сети. Для регулирования напряжения в одном или нескольких узлах схемы PQ -генераторы необходимо заменить регулируемыми PU — или PEq -генераторами, которые, поддерживая заданное значение модуля напряжение, позволяют векторам задающего тока и напряжения узла поворачиваться на различные углы в пределах, определяемых заданным регулировочным диапазоном реактивной мощности генератора (возбуждения генератора). Регулируемые генераторы обеспечивают уровень напряжений, соответствующий условиям задачи, и участвуют в сведении баланса реактивных составляющих токов (мощностей) в рассматриваемой схеме. Генераторные узлы с регулируемыми генераторами называют опорными узлами. Постановка задачи расчета установившегося режима i=1, 2, 3, 4 2. Чтобы исключить неконтролируемое вращение векторов в ходе итераций, препятствующее сведению баланса токов, и обеспечить однозначное расположение векторов задающих токов и напряжений на комплексной плоскости достаточно в одном из опорных узлов задать фазу напряжения. i. Uj ie. Ji. J

26 Постановка задачи расчета установившегося режима PU -генератор. PQ -генератор. I U +1 +j26 Постановка задачи расчета установившегося режима PU -генератор. PQ -генератор. I U +1 +j I U p. I а. I

Постановка задачи расчета установившегося режима Peq - генератор 27 Постановка задачи расчета установившегося режима Peq — генератор

28 Базисный, балансирующий и опорный узлы Базисный узел – узел с заданным комплексом напряжения28 Базисный, балансирующий и опорный узлы Базисный узел – узел с заданным комплексом напряжения ( исключаютс я уравнения балансов Р и Q ). Заданы модуль напряжения и его угол, отсчитываемый от оси вещественных комплексной плоскости против часовой стрелки, или реальная и мнимая части комплекса. Обычно комплекс напряжения в базисном узле ориентируют по оси вещественных комплексной плоскости, тогда угол или мнимая часть комплекса напряжения базисного узла равны нулю. Балансирующий узел – узел, замыкающий баланс активной мощности в схеме ( исключаетс я уравнение баланса Р ). Уравнение балансирующего узла исключается из системы узловых уравнений. Результирующий задающий ток в этом узле определяется после расчета напряжений в узлах непосредственно по его уравнению. (Если базисный узел одновременно является балансирующим, его целесообразно называть базисно-балансирующим, ББУ ). Опорный узел – узел с генератором, регулирующим напряжение в узле (тип PU ). Уравнрение баланса Q исключается. Опорные узлы необходимо выбирать таким образом, чтобы фазовые углы между напряжениями любого узла с PQ -генератором или нагрузочного узла и ближайшего по “дереву” схемы опорного узла не превышал ~30 -35 . Это предупреждает затруднения, связанные с неоднозначностью и нарушением условий существования решения. NB: В базисно-балансирующем узле в отличие от балансирующего узла вектор напряжения при итеративном решении системы уравнений не может вращаться. Это нужно учитывать при выборе опорных узлов, если на линиях к ББУ получаются углы более 30.

29 При нумерации узлов, как правило, базисный узел получает последний номер, а балансирующий —29 При нумерации узлов, как правило, базисный узел получает последний номер, а балансирующий — предпоследний, это позволяет избежать перенумерации узлов при формировании системы уравнений. Система уравнений Узел 4 – базисный, узел 3 – балансирующий: Узел 4 – базисно-балансирующий: Классическая постановка задачи расчета установившегося режима методом узловых напряжений не позволяет иметь несколько балансирующих узлов, т. к. число искомых напряжений должно быть равно числу уравнений, но можно назначить несколько базисно-балансирующих узлов. Однако в этом случае необходимо корректно задать модули и углы напряжений в этих узлах (последовательные приближения). Исходная система уравнений: . , , 4444343242 4242323222121 1313212111 UYJUYUYUY . , , 4343333232131 4242323222121 1313212111 UYJUYUYUY 4444343242 3434333232131 2424323222121 1313212111 , , , JUYUYUYUY JUYUYUY

30 Методы решения системы нелинейных узловых уравнений Основными характеристиками итерационного метода являются:  условия30 Методы решения системы нелинейных узловых уравнений Основными характеристиками итерационного метода являются: условия и надежность сходимости к решению, “чувствительность” к точности начального приближения, скорость сходимости, характеризуемая необходимым числом итераций для получения решения с заданной точностью, характер сходимости — апериодический или колебательный, сложность алгоритма и потребность в ресурсах компьютера при его реализации. При так называемом “плоском” старте в качестве начального приближения принимают номинальные напряжения в узлах, а вектора напряжений ориентируют на комплексной плоскости так же, как вектор напряжения в базисном узле. Зная напряжения в узлах, можно определить задающие токи узлов, что дает возможность решить систему уравнений и найти новые значения напряжений. В процессе итераций происходит поворот векторов задающих токов и напряжений на комплексной плоскости (уточняется ориентация) и корректируются их модули, контролируются ограничения режима регулируемых генераторов. Итерации прекращаются по достижении заданной точности расчета.

31 Метод Зейделя (поузловое определение напряжений) k – номер итерации    31 Метод Зейделя (поузловое определение напряжений) k – номер итерации . 1 , 1 3 434232131 33 3 2 424323121 22 2 1 313212 11 1 U SS UYUYUY Y U U SS UYUY Y U нг 3 Г 3 нг 2 Г 2 нг 1 Г 1 ii i iнгi. Гn ij k jiji j k jij k i Y U SS UYUY U 1 )(1 1 )1(

32 Метод Ньютона (одновременное определение напряжений узлов) Функции балансов содержат сопряженный комплекс напряжения, поэтому32 Метод Ньютона (одновременное определение напряжений узлов) Функции балансов содержат сопряженный комплекс напряжения, поэтому не являются аналитическими, не имеют комплексной производной. . 0 3 , , , 0 3 , , 3 4343332321313313 22 4243232221213312 11 3132121113311 3 нг 3 Г I 2 нг. Г I 1 нг. Г I U SS UYUYUYUYUUUw U SS UYUYUYUUUw . 0 3 , , , 0 3 , , 3 434 2 3332321313313 22 24242323 2 22221213312 11 13131212 2 1113311 нг 3 Г 333 S нг. Г S SS UUYUYUUYUUYUUUw SS UUYUUYUYUUUw UU jj. UUe. U Uj. UU U sincos

33 Запись уравнений в форме баланса токов с переменными U’ и U” функция баланса33 Запись уравнений в форме баланса токов с переменными U’ и U” функция баланса активных составляющих токов функция баланса реактивных составляющих токов , 3 3 2 2 4 , 1 i iiiнг i ii. Г ij j jijjijiiiiiii. I U UUQUUP U UQUP Ub. Ugw . UUU где U UUPUUQ U UPUQ Ug. Ubw iii i iiiнг i ii. Г ij j jijjijiiiiiii. I 222 2 2 4 , 1 ,

34 Запись уравнений в форме баланса мощности c переменными U,  В программах расчета34 Запись уравнений в форме баланса мощности c переменными U, В программах расчета установившихся режимов чаще используют систему уравнений в форме баланса мощности с переменными U, U , при этом функции балансов полной мощности разделяются на балансы активной и реактивной мощностей: угол, отсчитываемый от вектора напряжения i -го или j-го узла к вектору напряжения базисного узла против часовой стрелки Эта форма записи позволяет уменьшить порядок системы узловых уравнений, для узлов с PU -генераторами исключаются уравнения баланса реактивной мощности. При выходе на ограничения по Q и смене типа генераторных узлов требуется переформирование системы уравнений. δ i; j —

Запись уравнений в форме баланса мощности c переменными U,  35 Дальнейшего снижения трудоемкостиЗапись уравнений в форме баланса мощности c переменными U, 35 Дальнейшего снижения трудоемкости расчета можно достигнуть с помощью следующей замены переменных: Тем самым устраняется необходимость многократного обращения к подпрограммам вычисления тригонометрических функций, и далее задача решается в переменных U , . Окончательно фазы напряженийȠ в узлах определяются после сходимости к решению. Ƞ = tg (δ i; j / 2 ) и далее 1 — Ƞ 2 2 Ƞ sin δ i; j = —— ; cos δ i; j = ——. 1+ Ƞ 2 1 + Ƞ

36 матрица Якоби  на k-й итерации вектор поправок напряжений вектор небалансов на k-й36 матрица Якоби на k-й итерации вектор поправок напряжений вектор небалансов на k-й итерации k – номер итерации)( )( k k UWU U W U )()( kk II II UU UW UW UW U U U 3 2 1 )( 3 )1( 3 )( 2 )1( 2 )( 1 )1( 1 )()1( U UU UU UU kk kk. UUU 321 )( 3)1( 3 )( 2)1( 2 )( 1)1( 1 )()1( UU UU UU kk kk UUU )( )( )( kk k UW UW UW )( 3 )( 2 )( 1 )( k I k I k w w w UW )( 3 )( 2 )( 1 )( k I k I k ww w UW

37 Матрицы-клетки матрицы Якоби    3 3 2 3 1 3 337 Матрицы-клетки матрицы Якоби 3 3 2 3 1 3 3 2 2 2 1 2 3 1 2 1 1 1 U w U w U w U w U w III III III I U W 3 3 2 3 1 3 3 2 2 2 1 2 3 1 2 1 1 1 U w U w U w U w U w III III III I U W

38 Элементы матриц-клеток Диагональные: Недиагональные:   2433 2222 i iiiнг ii ii. Гiii.38 Элементы матриц-клеток Диагональные: Недиагональные: 2433 2222 i iiiнг ii ii. Гiii. Г ii i i. I U UUQUUP UU UUQUUP g U w 2433 22 22 i iiiнг ii iii. Г ii i i. I U UUQUUP UU UUPUUQ b U w 2433 22 2 2 i iiiнг ii ii. Гiii. Г ii i i. I U UUPUUQ UU UUPUUQ b U w 2433 22 22 i iiiнг ii ii. Гiii. Г ii i i. I U UUPUUQ UU UUQUUP g U w i j j i. Ig U w i j j i. Ib U w i j j i. Ig U w

39 Каждая итерация состоит из нескольких шагов:  1.  Вычисляется вектор небалансов. Если39 Каждая итерация состоит из нескольких шагов: 1. Вычисляется вектор небалансов. Если достигнута заданная точность сведения баланса, итерации прекращаются. 2. Вычисляются элементы матрицы Якоби с контролем ограничений режима регулируемых генераторов. 3. Решается система уравнений методом Гаусса и определяется вектор поправок. 4. Определяются напряжения узлов.

40 Расчет установившегося режима без балансирующего узла Обычно совмещают вектор напряжения в базисном узле40 Расчет установившегося режима без балансирующего узла Обычно совмещают вектор напряжения в базисном узле (узел 4) с осью вещественных, тогда угол или мнимая составляющая комплекса напряжения в этом узле равна нулю. Число переменных получается на единицу меньше числа уравнений. Вместо нее вводится новая переменная — относительное отклонение частоты от номинальной Система уравнений для определения поправок реальной и мнимой составляющих напряжений и отклонения частоты принимает вид: Используется полная исходная система узловых уравнений. Эта система уравнений преобразуется к системе из 8 уравнений балансов активных и реактивных составляющих токов или мощностей. В базисном узле фиксируется только положение вектора напряжения (угол или мнимая составляющая). 4444343242 3434333232131 2424323222121 1313212111 , , , JUYUYUYUY JUYUYUY ном f ff f *

41 Если активные мощности генераторов постоянны, режим балансируется за счет регулирующего эффекта нагрузки, поэтому41 Если активные мощности генераторов постоянны, режим балансируется за счет регулирующего эффекта нагрузки, поэтому СХН должны учитывать зависимость нагрузки от частоты, кроме того, целесообразно учесть зависимость реактивных параметров элементов схемы замещения от частоты: Для регулирования частоты можно выделить один или несколько генераторов и распределять возникающий при отклонениях частоты небаланс активной мощности между этими генераторами: — коэффициент крутизны статической характеристики генератора по частоте 4 3 1 * 4 3 4 1 4 * 3 3 3 1 3 * 1 3 1 1 1 I III III w w w f U U f w U w U w f w U w *f 1 b fb max. mi n. *, ГГГГ. номс. ГГPf. PKPf. P номном Г. номном. ГГ c fff PPP K .

Метод эквивалентных преобразований Если заданы параметры элементов схемы замещения,  причем генерация и нагрузкаМетод эквивалентных преобразований Если заданы параметры элементов схемы замещения, причем генерация и нагрузка представлены реальными источниками ЭДС (тока), напряжения узлов и токи продольных и поперечных ветвей схемы можно найти, используя эквивалентные преобразования схемы замещения. Преобразование является эквивалентным, если электрические режимы узлов и ветвей в непреобразованной части схемы до и после какого-либо эквивалентного преобразования одинаковы. Задачу сведения балансов токов (балансов активной и реактивной мощностей в узлах схемы) при итеративном расчете установившегося режима с использованием метода эквивалентных преобразований можно решить, если найти способ регулирования источников ЭДС на поддержание заданных уставок активной мощности и напряжения в процессе итераций.

43 Характеристики:  Безматричный метод, не требует формирования, хранения и обработки слабозаполненных матрицы узловых43 Характеристики: Безматричный метод, не требует формирования, хранения и обработки слабозаполненных матрицы узловых проводимостей и матрицы Якоби при решении системы узловых уравнений. Применим для расчета как установившихся режимов, так и мгновенных (на шаге численного интегрирования) электрических режимов при расчете электромеханических переходных процессов. Позволяет построить компактный алгоритм и вычислительное ядро программы в комплексном виде. Обеспечивает простую и эффективную обработку разделения схемы на произвольное число частей, вплоть до отключения отдельных узлов, с последующим поочередным или одновременным слиянием частей. Позволяет реализовать корректное эквивалентирование схемы исключением части узлов, поскольку в основе самого метода расчета установившегося режима лежит исключение узлов. Обладает надежной сходимостью при балансировании установившихся режимов и дает однозначное решение за счет применения адекватной PEq -модели генератора. Эта модель органично связана с динамической моделью генератора, обеспечивая корректный расчет начальных условий на основе векторной диаграммы генератора. Метод эквивалентных преобразований

44 Не критичен к начальному приближению, от начального приближения слабо зависит число итераций, 44 Не критичен к начальному приближению, от начального приближения слабо зависит число итераций, необходимых для поворота векторов ЭДС на требуемые углы при балансировании установившегося режима. Предоставляет полный и точный контроль ограничений режима генераторов при балансировании установившегося режима. Дает возможность реализовать различные законы управления регулируемыми компенсирующими устройствами в ходе итеративного балансирования режима с полным и точным контролем ограничений их режима. Недостатком метода эквивалентных преобразований является большое число итераций, необходимое для балансирования установившегося режима, по сравнению с решением системы узловых уравнений методом Ньютона. На расчет задачи может потребоваться несколько сотен итераций. Метод эквивалентных преобразований

45 Метод эквивалентных преобразований Расчет электрического режима методом эквивалентных преобразований состоит из следующих этапов:45 Метод эквивалентных преобразований Расчет электрического режима методом эквивалентных преобразований состоит из следующих этапов: Замещение рассматриваемой электрической схемы эквивалентным источником тока в узле свертки (свертка схемы к одному из узлов). Определение напряжения в узле свертки (при расчете мгновенного электрического режима в процессе интегрирования) из условия, что ток эквивалентного источника тока в узле свертки равен нулю. Или определение тока эквивалентного источника, включенного в базисный узел, выбранный за узел свертки, при расчете установившегося режима. Определение напряжений в узлах при развертывании схемы в порядке, обратном исключению узлов. Расчет электрического режима продольных и поперечных элементов схемы.

46 Эквивалентные преобразования Для свертки схемы необходимы следующие эквивалентные преобразования:  замена параллельных ветвей46 Эквивалентные преобразования Для свертки схемы необходимы следующие эквивалентные преобразования: замена параллельных ветвей одной эквивалентной ветвью исключение узла слиянием узлов ветви с нулевым сопротивлением (включенного выключателя) исключение узла при помощи преобразования активной многолучевой звезды в эквивалентный многоугольник

47 Слияние узлов включенного выключателя Хорошо работает и обычное замещение ветви выключателя очень небольшим47 Слияние узлов включенного выключателя Хорошо работает и обычное замещение ветви выключателя очень небольшим и очень большим сопротивлением во включенном и отключенном состояниях i j m Y im Y jm Jm 0 Ymn 0 Ym 0 0 n а) б) Yn 0 0 i j m Y im Y jm 0 Jmn 0 Jn 0 l k Y kn Y ln

48 Исключение одиночного узла Примечание. Это частный случай исключения узла при помощи более общего48 Исключение одиночного узла Примечание. Это частный случай исключения узла при помощи более общего преобразования активной многолучевой звезды в эквивалентный многоугольник i j m Y im Y jm Y mn Y n 0 0 J n 0 i j m n e Yim Y jm Z mn Z n 0 0 0 n E i j m Y im Y jm Y m 0 0 J m 0 n 0 Y m 0 J m 0 ( n ) I nm а) б) в )

49 Исключение узла преобразованием активной многолучевой звезды в эквивалентный многоугольник  1  249 Исключение узла преобразованием активной многолучевой звезды в эквивалентный многоугольник 1 2 i j m n 1 2 i j m e e Z 1 n Z 2 n Z in Z jn Z mn Z n 0 Z ij Z ie 0 0 0 n. E a ) б ) Z 1 e Z m e Z 12 Z jm 1 2 i j m n 1 2 i j m Y 1 n Y 2 n Yin Yjn Ymn Yn 0 Yij 0 J 20, Y 20 Ji 0, Yi 0 J 10, Y 10 Jm 0, Ym 0 Jj 0, Yj 0 Jn 0 в) г)

50 Исключение узла преобразованием активной многолучевой звезды в эквивалентный многоугольник Следует подчеркнуть, что преобразование50 Исключение узла преобразованием активной многолучевой звезды в эквивалентный многоугольник Следует подчеркнуть, что преобразование активной многолучевой звезды в эквивалентный многоугольник не существует, если генерация и/или нагрузка в исключаемом узле представлена идеальными источниками тока, поскольку для получения формул использовалось преобразование источника тока в источник ЭДС и обратно. Для идеальных источников тока и ЭДС взаимное преобразование не определено, их воздействие на режим схемы различно — в узле с идеальным источником ЭДС напряжение постоянно, а в узле с идеальным источником тока — нет. Если формально выполнить исключение узла с идеальным источником тока, искажается режим внешней сети, примыкающей к узлам, смежным исключаемому узлу. Это проявляется в том, что потери мощности в ветвях исходной и эквивалентной схем не совпадают. nn jnin ij Y YY YY)( nn in n n i Y Y YY 0 )( 0 nn in n n i Y Y JJ 0 )( 0 021. . . nmnjninnnnn. YYYYYYY

Регулирование P eq  -генераторов Разложив выражения для активной мощности и напряжения в рядРегулирование P eq -генераторов Разложив выражения для активной мощности и напряжения в ряд Тейлора, можно получить линеаризованную систему уравнений для определения поправок модуля и угла вектора ЭДС при регулировании генератора на поддержание заданных уставок активной мощности и напряжения 51, . YUYEUYUEUIUj. QPS 233)(33 g. UEUy. P YUE 23)cos(3 b. UEUy. QYUE 23)sin(3 Ej. Ee. EUj. Ue. U Yj q yejbg Z Y 122 bgy b b. Qy. EEy U UEUE 6 12 cos 9 cos

Регулирование P eq  -генераторов 52, .    , , 0 0Регулирование P eq -генераторов 52, . , , 0 0 UU U E E U PP P

53 Регулирование генераторов 53 Регулирование генераторов

54 Регулирование генераторов  Выполняется итерация. После итерации модуль вектора напряжения в узле уменьшился,54 Регулирование генераторов Выполняется итерация. После итерации модуль вектора напряжения в узле уменьшился, а по углу повернулся на угол по часовой стрелке (красный цвет). Определяются: При выходе на ограничения ЭДС Eq фиксируется, но реактивная мощность остается свободным параметром режима. При перегрузке генератора по току статора уменьшается Eq , если она не вышла на нижнее ограничение, в противном случае уменьшается активная мощность. Таким образом, выполняется полный и точный контроль по ограничениям режима генератора, т. е. балансирование режима схемы осуществляется с контролем возможности существования режима самих генераторов. NB: с расчетчика снимается обязанность контроля реактивной мощности генераторов. Процесс балансирования приобретает физический смысл – регулирование генераторов при безынерционном изменении возбуждения и положения ротора генератораотклонения поправки 1) 2) 3) 4) YUEU 3 IU 3 j. QP PP 0 UU 0 E E EEEE sin. EE, cos. EEComplex. E UE q qd q d qcos. U x xx E x x E maxqqminq.

Определение предельных перетоков по линиям и сечениям методом непрерывного утяжеления режимов Теория статической апериодическойОпределение предельных перетоков по линиям и сечениям методом непрерывного утяжеления режимов Теория статической апериодической устойчивости – самостоятельно по книге Г. В. Меркурьев, Ю. М. Шаргин Устойчивость энергосистем, изд.

56 Определение предельных перетоков по линиям и сечениям методом непрерывного утяжеления режимов  56 Определение предельных перетоков по линиям и сечениям методом непрерывного утяжеления режимов

57 Определение предельных перетоков по линиям и сечениям методом непрерывного утяжеления режимов Нужно найти57 Определение предельных перетоков по линиям и сечениям методом непрерывного утяжеления режимов Нужно найти предельный переток по линии 6370 -470 ПС Ленинградская – Калининская АЭС 750 к. В и в сечении 1 , 2 , в которые входят также линии 330 и 110 к. В. Генерация электростанций ОЭС Северо-Запада, задействованных в изменении режима схемы, в исходном режиме составляет: ЛАЭС — 2000 МВт: первая очередь (узел 1316) не работает; вторая очередь (узел 1376) работают четыре агрегата по 500 МВт, суммарная мощность 2000 МВт. Киришская ГРЭС — 900 МВт, работают три агрегата по 300 МВт. ОЭС Балтии принимает из ОЭС Северо-Запада 534 МВт по сечению 3. ОЭС Белоруссии принимает 1332 МВт по сечению 4 из ОЭС Центра. В подпорном узле 470 на Калининской АЭС работает один энергоблок на номинальной мощности 1000 МВт. Загрузка сечений осуществляется за счет разгрузки второй очереди ЛАЭС до 1000 МВт по линейному закону за 1000 шагов со скоростью 1 МВт/шаг, Киришская ГРЭС разгружается от 900 до 500 МВт за 2000 шагов со скоростью 0, 2 МВт/шаг. До середины процесса утяжеления режима скорость перераспределения генерации активной мощности между Северо-Западом и Центром составляет 1, 2 МВт/шаг. После завершения разгрузки ЛАЭС приближении к предельному режиму и при нарушении апериодической устойчивости скорость перераспределения баланса мощности составляет 0, 2 МВт/шаг, что обеспечивает плавное «сползание» режима с четким обозначением предельных перетоков по сечениям 1 и 2. Файл nw 272_ss

58 Изменение мощности генераторов ЛАЭС и Киришской ГРЭС 58 Изменение мощности генераторов ЛАЭС и Киришской ГРЭС

59 Изменение перетоков по сечениям 59 Изменение перетоков по сечениям

60 Изменение напряжений 60 Изменение напряжений