Лыкин А.В. Электрические системы и сети - файл Part_3.doc

приобрести
Лыкин А.В. Электрические системы и сети
скачать (2874.2 kb.)
Доступные файлы (11):
Part_3.doc3312kb.08.06.2006 21:19скачать
Part_4.docскачать
Part_5.doc2076kb.16.12.2001 18:28скачать
Part_6.doc256kb.16.12.2001 18:29скачать
Part_7.doc1427kb.16.12.2001 18:33скачать
n6.doc43kb.28.01.2002 12:33скачать
Ps-8_Literatura.doc32kb.16.12.2001 19:05скачать
n8.doc743kb.08.06.2006 21:19скачать
n9.doc752kb.08.06.2006 21:19скачать
n10.doc1115kb.08.06.2006 21:19скачать
n11.doc850kb.08.06.2006 21:19скачать

Part_3.doc

  1. моделирование режимов сложных схем электрических сетей

    1. применение теории графов для моделирования схем электрических сетей


Электрические сети современных ЭЭС насчитывают сотни и даже тысячи ЛЭП и трансформаторов. Расчеты режимов сложных схем электрических сетей требуют специальных моделей представления схем и компактной записи уравнений. Такими моделями являются графы и матрицы.

Линии, трансформаторы и другие элементы электрической сети представляются в расчетах своими схемами замещения, состоящими из ветвей с сопротивлениями и проводимостями. Все шины электрических станций и подстанций являются узловыми точками сети. Количество этих узловых точек или узлов схемы сети обозначим буквой n, а количество ветвей, соединяющих эти шины,  m. Если сеть не содержит замкнутых контуров, то количество узлов и ветвей различаются на 1: n = m + 1. При наличии контуров n = m + 1 – k, где k – количество независимых контуров.

Графы являются топологическими моделями схем электрических цепей.

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

Элементами ЭЭС, которые моделируются ребрами графа, являются ЛЭП, трансформаторы, реакторы, батареи конденсаторов и др. Как правило, все они представляются П-образными схемами замещения и поэтому имеют элемент связи между двумя граничными узлами – продольная ветвь, и элементы, связывающие узлы с нейтральной точкой системы N  поперечные ветви, рис. 3.1.



Рис. 3.1. П-образная схема замещения – а и ее граф  б
Для ЛЭП (см. раздел 1, (1.13)):




,

(3.1)







(3.2)

Обычно Z = (r0 + jx0)l и .

Для трансформатора (см. раздел 1, (1.42))







(3.3)

при kт > 1. Если kт = 1, то из (3.3) получается Г-образная схема замещения трансформатора.

Для реакторов и батарей конденсаторов, включенных в виде продольных элементов сети, параметры схемы замещения: Z = jXр и Z = jXc. Y1 = Y2 = 0 (Y1 или Y2 может быть отлично от нуля и моделировать потери активной мощности в реакторе или батарее конденсаторов). В случае их включения в виде поперечных ветвей: Z = 0, а Y1 и Y2 представляются одной поперечной ветвью – Y шунта: . Аналогично могут представляться своими схемами замещения электрические нагрузки (см. раздел 1, (1.46)).

Рассмотрим пример схемы электрической сети, состоящей из ЛЭП и трансформатора, рис. 3.2,а. Ее схема замещения есть две соединенные между собой П-образные схемы замещения ЛЭП и трансформатора, а граф будет состоять из двух графов П-образных схем, рис. 3.2,б.



Рис. 3.2. Схема простой электрической сети – а и ее граф – б

Для более сложных схем, например схемы на рис. 3.3,а, удобно ввести в рассмотрение нейтральную плоскость в сети и рассматривать узлы графа сети «висящими» над нейтральной плоскостью N и соединенными с ней поперечными ветвями, рис. 3.3,б.



Рис. 3.3. Схема электрической сети из восьми узлов и десяти ветвей – а и ее граф – б

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



Рис. 3.4. Граф сети без изображения ребер, связанных с нейтральной плоскостью

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

Первая матрица инциденций M определяет связи узлов и ветвей и составляется для направленного графа. Она состоит из нулей и единиц. Узлам графа ставятся в соответствие строки матрицы, а ветвям – столбцы. В каждой строке связанного графа есть хотя бы одна единица. Она указывает на наличие связи между соответствующим узлом и ветвями графа. Если такая ветвь имеется, то в соответствующем столбце ставится единица. Знак плюс означает, что направление на ветви задано «от узла», и знак минус, что направление задано «к узлу». Так для графа на рис 3.4 первая матрица соединений выглядит следующим образом




.

(3.4)

Количество единиц в любой строке матрицы указывает на количество присоединенных к данному узлу ветвей. В каждом столбце матрицы M находится ровно две единицы, которые показывают, из какого узла выходит (по направлению) и в какой узел входит ветвь, описанная в данном столбце. Соответственно говорят о начале ветви (узел, отмеченный положительной единицей) и конце ветви (отрицательная единица).

Использование матрицы M в практических расчетах очень неэффективно, так как требует хранение в памяти компьютера большого числа нулей (нулевых элементов) и занимает много времени в процессе расчета для поиска ненулевых элементов. Обычно используется компактная форма этой матрицы – так называемая ленточная форма L




.

(3.5)

В первой строке матрицы L указывается номер (имя) начального узла, а во второй, того же столбца, – номер (имя) конечного узла. Пара номеров узлов в столбце образует имя ветви, например, для ветви b это 1 – 2.

Ленточная форма графа является самой компактной и используется при описании топологии схемы электрической сети.

Кроме матрицы инциденций M для моделирования режимов электрических сетей используется матрица инциденций N (вторая матрица соединений). В этой матрице устанавливается связь между независимыми контурами и ветвями графа. Строкам матрицы ставятся в соответствие независимые контуры, столбцам – ветви графа. Для составления матрицы N необходимо задать направления обхода контуров. При поочередном обходе каждого контура проверяется на совпадение направления обхода контура и направления ветви. Если они совпадают, то в соответствующем столбце матрицы N ставится положительная единица, в противном случае – отрицательная единица. Если ветвь не входит в рассматриваемый контур, то в этом столбце ставиться ноль. Так для графа сети, изображенного на рис. 3.4, можно выбрать в качестве независимых контуры: I  {a, e, h}, II – {a, b, c, f, i, e}, III – {f, g, h}. При обходе указанных контуров по часовой стрелке получается следующая матрица N:




.

(3.6)




Рис. 3.5. Многослойный граф
В некоторых случаях можно использовать многослойные графы, в которых сеть каждого напряжения располагается в отдельном слое. Получается, что в горизонтальных слоях находятся ветви, моделирующие линии электропередачи, а между ними вертикально изображаются трансформаторные связи, рис. 3.5. Таких слоев может быть столько, сколько ступеней номинального напряжения имеется в сети.

Пример. Восстановить матрицу соединений M по ее ленточной форме L для графа сети на рис. 3.4.

Расчет выполним в системе Mathcad.











    1. матричные формы моделей электрических сетей и их режимов


Каждая продольная ветвь в графе электрической сети характеризуется сопротивлением Zj = Rj + jXj, а поперечная ветвь – проводимостью Yi = Gi + jBi (j = 1,2,…,m; i = 1,2,…,n), которые образуют матрицы параметров электрической сети – матрицу сопротивлений продольных ветвей и матрицу-столбец проводимостей поперечных ветвей – шунтов:




.

(3.7)

Здесь Zвjj = Zj, а YNi = Yi . Недиагональные элементы матрицы Zв обычно равны нулю, хотя в некоторых случаях учитывают взаимные сопротивления ветвей, которые могут быть отличны от нуля, например для близко расположенных ЛЭП возможно наличие взаимной индукции.

Кроме пассивных ветвей в сети существуют активные ветви, включающие источники ЭДС и тока. Эти ветви, как правило, являются поперечными и моделируют генераторы электрических станций (ЭДС) и потребителей электрической энергии (источники тока), рис. 3.6,а.

В

.

(3.8)



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




Рис. 3.6. Изображения поперечных ветвей
Принято не изображать на графе сети не только шунтирующие проводимости, но и активные поперечные ветви с ЭДС и источником тока, однако источник тока все же задают упрощенным изображением в виде стрелочки, направленной в узел, рис. 3.6,б. Это показывает, что в сеть «вливается» извне ток генерации или нагрузки (с обратным знаком). Такие токи называются токами инъекции (injection current) или задающими токами.

Матрицы E и J задают режим работы электрической сети и являются векторами независимых переменных. Они относятся к режимным параметрам электрической сети. Другие режимные параметры называются зависимыми переменными. К ним относятся напряжения в узлах, токи и напряжения в продольных ветвях и ряд других параметров режима:

U – матрица напряжений в узлах (узловые напряжения);

I – матрица токов ветвей;

U – матрица напряжений в ветвях (падения напряжения на сопротивлениях ветвей);

Sв(н) – матрица потоков мощности в начале ветвей;

Sв(к) – матрица потоков мощности в конце ветвей;

Sв – матрица потерь мощности в ветвях.

    1. Узловые уравнения установившегося режима





Рис. 3.7. Пример графа электрической сети
Рассмотрим пример направленного графа электрической сети, изображенного на рис. 3.7. Для удобства записи в матричной форме параметров ветвей присвоим каждой ветви ее порядковый номер (на рис. 3.7 курсив). Составим матрицу соединений M для этого графа

У

.

(3.9)




множим эту матрицу на матрицу токов ветвей, будем иметь:





.

(3.10)

Полученное соотношение является первым законом Кирхгофа в матричной форме записи




.

(3.11)

Так как к узлам графа электрической сети еще присоединены другие поперечные ветви с ЭДС и проводимостью шунта, то задающий ток в (3.11) включает в себя также токи данных ветвей




.

(3.12)

Здесь: JГ – матрица токов генерации (ветви с ЭДС), которые определяются через мощности генерации; JН – матрица токов нагрузки, которые определяются через мощности нагрузки (имеет обратное направление – от узла); JY – матрица токов в проводимости шунтов, которые зависят от проводимости шунта из матрицы YN и напряжения в узле из матрицы U (также имеет обратное направление – от узла, так как моделирует потребление мощности).

Умножим транспонированную матрицу соединений МT на матрицу узловых напряжений, получим:






(3.13)

или




.

(3.14)

По закону Ома в матричной форме записи имеем




.

(3.15)

или




.

(3.16)

Подставим в (3.11) выражение для матрицы токов ветвей (3.12) и затем (3.14), получим




.

(3.17)

Введем обозначение




,

(3.18)

тогда (3.17) приобретет вид




.

(3.19)

Полученное соотношение является уравнением узловых напряжений (потенциалов) в матричной форме записи. Матрицу Y называют матрицей узловых проводимостей электрической сети. Рассмотрим структуру этой матрицы, для чего выполним матричные перемножения в (3.18). Заметим, что обратная матрица сопротивлений ветвей легко получается в силу своего диагонального вида – ее элементы суть обратные величины к сопротивлениям ветвей.

В начале перемножим первые две матрицы матричного произведения:









.

(3.20)

Полученную матрицу умножим справа на матрицу MT. В результате получим:

.

(3.21)

Из полученной матрицы можно сделать следующие выводы о вычислении ее элементов:

1) Элементы, расположенные на диагонали матрицы, вычисляются как сумма проводимостей ветвей, подходящих к соответствующему узлу:




,

(3.22)

где Yij – диагональный элемент матрицы Y;

Zj – сопротивление j-й ветви;

i – множество номеров узлов, связанных с i-м узлом.

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




.

(3.23)

Запишем узловое уравнение для узла с номером i:




.

(3.24)

Объединив подобные члены, получим, что в диагональные элементы матрицы Y войдут дополнительные слагаемые YNi:




,

(3.25)

т. е. диагональный элемент будет равен сумме проводимостей всех подходящих к i-у узлу ветвей, включая поперечную ветвь – шунт YNi.

Задающие токи узлов в (3.19) будут состоять только из токов генерации и токов нагрузки.

В случае отсутствия связей с нейтральной плоскостью N система уравнений (3.19) не имеет единственного решения, так как в этом случае определитель матрицы Y равен нулю. Сумма всех задающих токов в такой сети равна нулю:




.

(3.26)

Следовательно, среди всех n узлов можно выделить узел, например, с номером n, ток в котором равен:




.

(3.27)

Для уравнений узловых напряжений это означает, что одно уравнение лишнее, т. е. зависит от остальных уравнений и может быть получено через сумму всех остальных уравнений. Так как ток в этом узле может быть получен из баланса токов в сети (3.27), то его называют балансирующим. Обычно это шины мощной электростанции или системы.

Таким образом, из системы (3.19) исключается одно уравнение и тогда получается система независимых линейных уравнений порядка n – 1. Однако, так как число неизвестных напряжений по-прежнему равно n, то в одном из узлов следует задать напряжение по величине и фазе, так, чтобы все напряжения вычислялись относительно этого известного напряжения. Такой узел в сети называется базисным. Обычно фазу напряжения базисного узла принимают равной нулю, т. е. вектор напряжения базисного узла совмещают с действительной осью. Остальные узлы называют независимыми узлами.

Во многих случаях балансирующий узел и базисный узел совмещают, и в дальнейшем будем считать, что это один и тот же узел.

Таким образом, с исключением уравнения для базисного балансирующего узла с номером n, будем иметь систему уравнений (3.19) с числом уравнений n – 1, однако в эти уравнения будет входить слагаемое с заданным напряжением базисного узла.

Изменим номер базисного балансирующего узла. Пусть его номер есть 0 (ноль). Тогда уравнение (3.19) приобретет следующий вид:




.

(3.28)

где Y0 – матрица проводимостей ветвей, связывающих независимые узлы с базисным балансирующим узлом;

U0 – напряжение базисного узла (скаляр).

Матрица узловых проводимостей в (3.28) имеет порядок n – 1 и определятся через матрицу инциденций M, в которой нет одной строки, соответствующей балансирующему узлу.

Необходимо заметить, что во всех уравнениях, где одновременно присутствуют токи и напряжения: (3.15), (3.16), (3.17), (3.19), (3.24) и (3.28) напряжения даны в фазных значениях, хотя индекс (буква «ф») для простоты не записывался. Эти же уравнения можно считать записанными и для линейных напряжений, однако токи будут увеличенными в раз и для вычисления истинных токов их следует уменьшать в .

    1. формы линейных уравнений установившегося режима и их решение


В качестве известных независимых переменных в уравнениях установившегося режима могут выступать задающие токи узлов и напряжение базисного узла. В этом случае решение уравнения (3.28) может быть записано в виде




.

(3.29)

Здесь Z – матрица узловых сопротивлений.

Численное решение системы уравнений (3.28) выполняется методом Гаусса или другим методом решения системы линейных алгебраических уравнений.

В случае, когда если известны мощности в узлах сети – задающие мощности Si, то токи можно вычислить приближенно , (i = 1,…,n – 1). Задающие мощности также как и токи складываются из мощности генерации и мощности нагрузки:




.

(3.30)

Другой приближенный подход связан с представлением задающих токов через напряжения и проводимости , где YSi – проводимость нагрузки или/и генерации (схема замещения). Для i-о узла имеем:






(3.31)

Объединив подобные члены, получим






(3.32)

где в элемент Yii входит проводимость YSi. Знак перед этой проводимостью зависит от того, какая мощность преобладает в узле: плюс, если нагрузка и минус, если генерация. В матричной форме записи:




.

(3.33)

Решение матричного уравнения (3.29) запишется в виде:




.

(3.31)

Комплексную матрицу узловых проводимостей Y иногда представляют в блочной форме через ее вещественную G и мнимую B составляющие и тогда система уравнений (3.33) становится системой с вещественными величинами:




.

(3.35)

После перемножения двучленов в (3.31), будем иметь:




.

(3.36)

Приравняем отдельно вещественные и мнимые части полученного уравнения и получим два матричных уравнения с вещественными величинами:






(3.37)

или в компактной форме:




.

(3.38)

Решение (3.38) запишется в виде:




.

(3.39)

Пример. Рассчитать напряжения в узлах и токи в ветвях схемы электрической сети, граф которой изображен на рис. 3.7. Исходные данные для расчета и расчет представлен в системе Mathcad.















Модель электрической сети





























    1. нелинейные уравнения установившегося режима


Так как во многих случаях расчеты ведутся при заданных мощностях нагрузок и генерации, то их следует ввести в уравнения установившегося режима.

Мощность в трехфазной сети в симметричных режимах выражается суммарной мощностью всех трех фаз:




.

(3.40)

В матричной форме это выражение можно записать, используя операцию диагонализации матрицы U. Матрица diag{U} есть квадратная матрица, в которой элементы матрицы U расположены по главной диагонали, а все остальные элементы равны нулю. Тогда




.

(3.41)

Уравнение установившегося режима









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






(3.42)

В левой части этого уравнения после умножения на напряжения стали линейными.

Умножим левую и правую части уравнения (3.42) слева на матрицу diag{U}, получим




.

(3.43)

Система уравнений (3.43) является системой нелинейных уравнений установившегося режима. В зависимости от формы представления комплексных величин применяют две основные формы этой системы уравнений.

В начале рассмотрим алгебраическую форму записи. Для i-о узла имеем:




.

(3.44)

После перемножения двучленов и разделения уравнения на два уравнения с вещественными величинами, получим систему 2(n – 1) алгебраических уравнений.






(3.45)

Здесь i = 1,…,n – 1.

Тригонометрическая форма нелинейных уравнений установившегося режима может быть получена, если комплексные величины в уравнении (3.39) записать в виде:




,

(3.46)

тогда




.

(3.47)

Уравнение (3.47) в тригонометрической форме запишется как




.

(3.48)

и после разделения на два вещественных уравнения






(3.49)

Обычно вместо угла ij используют дополняющий до 90 угол ij. ij = 90  ij, ij = 90  ij.

Тогда cos(i – j – ij) = cos(i – j – 90 + ij), а с учетом четности функции косинус cos(i – j – 90+ ij) = cos(90 – i + j – ij). Имея в виду, что cos(90 – ) = sin(, получим: cos(90 – i + j – ij) = sin(i – j + ij).

Аналогично sin(i – j – ij) = sin(i – j – 90 + ij) = –sin(90 – i + j – ij), в силу нечетности функции синус. Так как sin(90 – ) = cos(), получим:

–sin(90 – i + j – ij) = –cos(i – j + ij). Подставляя полученные соотношения в (3.49), будем иметь:






(3.50)

В полученной системе нелинейных уравнений установившегося режима искомыми переменными являются модули и фазовые углы, напряжений, в то время как в уравнениях (3.45) неизвестными являются вещественная и мнимая составляющие напряжений.

Пример. Рассчитать напряжения в узлах и потоки мощности в ветвях схемы сети, граф которой изображен на рис. 3.7. Исходные данные для расчета и расчет представлен в системе Mathcad.

















































































моделирование режимов сложных схем электрических сетей
Учебный материал
© nashaucheba.ru
При копировании укажите ссылку.
обратиться к администрации