Технические науки — от теории к практике № 2 (50), 2016г_'_
СЕКЦИЯ
«ЭНЕРГЕТИКА И ЭНЕРГЕТИЧЕСКИЕ ТЕХНИКА И ТЕХНОЛОГИИ»
ПРАКТИЧЕСКИЕ МЕТОДЫ РАСЧЕТА ТОКОВ КОРОТКОГО ЗАМЫКАНИЯ В ЭЭС ПРИ СЛОЖНЫХ ВИДАХ ЗАМЫКАНИЙ С ИСПОЛЬЗОВАНИЕМ ФАЗНЫХ КООРДИНАТ
Осак Алексей Борисович
научный сотрудник, Институт систем энергетики им. Л.А. Мелентьева (ИСЭМ) СО РАН, РФ, г. Иркутск E-mail: osakalexey@mail. ru
PRACTICAL METHOD OF CALCULATION OF SHORT CIRCUIT IN THE POWER SYSTEM UNDER DIFFICULT TYPES OF FAULT WITH PHASE COORDINATES
Alexey Osak
researcher, Melentiev Energy Systems Institute of Siberian Branch of the Russian Academy of Science (ESI SB RAS),
Russia, Irkutsk
АННОТАЦИЯ
Для выполнения расчетов токов короткого замыкания (ТКЗ) при сложных видах повреждений и несимметрии элементов сети целесообразно использовать метод фазных координат. В статье представлена методика расчета ТКЗ в фазных координатах в матричном виде и даны рекомендации по определению параметров схемы замещения для практических задач. Данная методика с 1998 года используется в лабораторных и научно исследовательских работах в Иркутском государственном техническом университете [2; 3], которая здесь уточнена в части задания трансфор-
Технические науки — от теории к практике _№ 2 (50), 2016г
маторных ветвей в первой матрице инциденций. В ИСЭМ СО РАН [11] данная методика была дополнена учетом предшествующего аварии режима. Появление современных математических программных пакетов позволяет широко использовать данную методику для расчета ТКЗ в фазных координатах без использования специализированного программного обеспечения.
ABSTRACT
To perform the calculation of short circuit current in complex forms of fault, and the asymmetry of grid elements it is expedient to use the method of phase coordinates. The method of calculating the short-circuit current in phase coordinates in the matrix form. The recommendations for the definition of equivalent circuit parameters for practical works. This method is used since 1998 in laboratory and scientific research works in the Irkutsk State Technical University [2; 3], which has been added of the calculation of transformer branch in a first incidence matrix. In the ESI SB RAS [11], this method has been added taking into account the previous regime of the accident. The advent of modern mathematical software packages allows extensive use of this method to calculate the short-circuit current in phase coordinates without the use of specialized software.
Ключевые слова: токи короткого замыкания; фазные координаты; ЭЭС.
Keywords: current of short circuit; phase coordinates; power systems.
Задача расчета токов короткого замыкания (ТКЗ) является одной из типовых задач расчета режимов электроэнергетических систем (ЭЭС), которая необходима как для проектирования, так и для эксплуатации. Задача расчета ТКЗ применяется для выбора характеристик первичного и вторичного оборудования, расчета уставок устройств релейной защиты и автоматики (РЗА), проверки селективности РЗА, разбора и анализа фактических коротких замыканий (КЗ) и определения места повреждения (ОМП).
Для решения задачи расчета ТКЗ широко применяются промышленные программные комплексы отечественной и зарубежной разработки, такие как: АРМ СРЗА (г. Новосибирск), блок RastrKZ в составе RastrWin (г. Екатеринбург), EnergyCS ТКЗ (г. Москва), блок ТКЗ в составе ПВК АНАРЭС (г. Иркутск) [11].
В тоже время имеется ряд случаев, когда расчет ТКЗ на промышленных программных комплексах приводит к большим трудностям. К таким случаям можно отнести [9, с. 110; 12, с. 29]: сложные виды повреждений (многократные КЗ, КЗ между провод-
Технические науки — от теории к практике № 2 (50), 2016г_
никами разных классов напряжения, сочетания КЗ и обрывов проводов и др.) или несимметрию нескольких элементов ЭЭС (неучет которой для конкретной практической задачи является неприемлемым). В этих случаях расчет ТКЗ проще всего выполнять с использованием фазных координат (ФК).
Метод ФК [9, с. 182] - это метод представления схемы замещения ЭЭС в пофазной постановке, а не в виде параметров схемы замещения прямой, обратной и нулевой последовательности как это осуществляется в широко применяемом методе симметричных составляющих (СС). При этом в ФК количество рассматриваемых фаз и схемы пофазного соединения всех элементов сети соответствуют фактическим схемам электроустановок (1 или 3 фазы, наличие или отсутствие нулевого провода и заземления нейтрали, наличие грозотроссов и их заземление). Продольные и поперечные сопротивления каждого элемента сети, а также взаимоиндукция между фазами и между элементами задается для каждой фазы отдельно [12, с. 248]. Вышеуказанные свойства ФК как раз и позволяют легко решать задачи расчета ТКЗ при сложных видах повреждений и несимметрии элементов сети.
Метод расчета режимов ЭЭС в ФК напрямую вытекает из расчета электрических цепей по законам Ома и Кирхгофа. Поэтому метод ФК известен уже более века. Но по сравнению с методом расчета в СС (при общепринятом допущении о симметричности элементов электрической сети матрицы сопротивлений в СС диагональные) метод ФК существенно более требователен к вычислительным ресурсам, потому что матрицы сопротивлений электрической сети в ФК не являются диагональными. Поэтому до появления высокопроизводительной вычислительной техники метод ФК не имел широкого применения для расчета режимов ЭЭС, в т. ч. для расчета ТКЗ.
С середины 1990-х годов уровень развития персональных компьютеров достиг уровня (по быстродействию и объему оперативной памяти), при котором применение метода ФК стало целесообразным. Имеется большое число публикаций на тему применения ФК для расчета различных режимов ЭЭС [1-9; 11-12], в т. ч. для расчета ТКЗ [9; 11-12].
В настоящее время широкое распространение получили универсальные математические программные комплексы (программные пакеты), позволяющие решать математические задачи различной сложности, в том числе задачи расчета режимов ЭЭС. Например:
• коммерческие: MATLAB, Mathcad, Mathematica;
• свободные: Scilab, GNU Octave, FreeMat.
ССибАК
Технические науки — от теории к практике и'и'и'.\//)^с. ш/о_№2 (50), 2016г
Появление вышеуказанных универсальных математических программных комплексов позволяет качественно повысить доступность метода ФК для практических расчетов.
В настоящей статье представлена математическая постановка задачи расчета ТКЗ в ФК, которая может быть реализована в любом из универсальных математических пакетов или на языке программирования высокого уровня.
В ряде работ для выполнения расчетов режимов в ФК предлагаются достаточно сложные и индивидуальные алгоритмы формирования матрицы узловых потенциалов для разных элементов сети [12, с. 248; 9, с. 183; 7, с. 28]. Такой подход справедлив и работоспособен, но существенно усложняет алгоритм формирования схемы большой размерности. Используя современные математические пакеты или современные математические библиотеки для работы с матрицами и комплексными числами можно многие практические расчеты ТКЗ выполнять без каких-либо специальных и сложных алгоритмов, если не требуется учет индивидуальных особенностей отдельных элементов сети. Ниже излагается такая методика расчета ТКЗ в ФК в матричном виде.
Схема замещения для расчета ТКЗ в ФК.
Для расчета режима в ФК составляется схема замещения ЭЭС в пофазной (в 3-х фазной) постановке в виде графа сети (узлов и ветвей), топология которого может быть описана в матричном виде первой матрицей инциденций М [10, с. 37].
М представляет собой матрицу, каждая строка которой соответствует одному из узлов, а каждый столбец одной из ветвей. В ячейках матрицы записывается «0» - если ветвь не связана с узлом, «+1» - если узел является началом обычной ветви, «-1» - если узел является концом ветви. Для трансформаторных ветвей вместо «+1» со стороны узла низкого напряжения нужно записывать Ктр
(где К = ив/ин ), знак «+» или «-» указывается как и для обычных ветвей. Каждый столбец матрицы может содержать только одно положительное число («+1» или «+ К») и одно отрицательное
(«-1» или « -Ктр »).
В расчетах в ФК базисным узлом предпочтительно считать узел «земля», потенциал которого равен нулю, соответственно напряжение узла будет равно фазному напряжению узла относительно земли. Строку матрицы М, соответствующую базисному узлу нужно исключить (чтобы размерность матрицы соответствовала числу
С1
Сибдк
Технические пауки — от теории к практике №2 (50), 2016г_www.sibac.info
независимых узлов), а при необходимости эту строку всегда можно восстановить.
Первую матрицу инциденций М можно представить в виде блочной матрицы М =|М. М р||, Где М а - первая матрица инциденций для дерева графа сети (квадратная матрица). М р - первая матрица инциденций для хорд.
Если требуется выполнять расчеты методом контурных токов, то необходимо составить вторую матрицу инциденций N. Строки матрицы соответствуют независимым контурам, а столбцы ветвям. В ячейках матрицы записывается «0», если ветвь не входит в контур, «+1», если обход контура совпадает с направлением тока ветви, «-1», если не совпадает. Вторую матрицу инциденций N [10, с. 43] также
можно представить в виде блочной матрицы N =||N. N ||, где Nа -соответствует ветвям дерева графа сети, а Nр - соответствует ветвям-хордам (квадратная матрица). Если выбирать независимые контуры таким образом, чтобы каждая хорда входила только в один контур, и направление тока в хорде совпадало с направлением обхода контура, то можно получить матрицу Nр - единичной, тогда матрицу N а
можно получить алгебраически N а =
-м р\м а.
Для расчета схемы в матричной форме необходимо вначале составить матрицы параметров схемы замещения: 2_в, И . J .
Матрица ЭДС ветвей Ев является матрицей вектор-столбцом, содержащей в каждой ячейке значение ЭДС в комплексном виде соответствующей ветви, причем знак перед ЭДС ставится в соответствии со следующим правилом. Если направление ЭДС совпадает с направлением тока, то ставится знак «+», иначе «-» (принятое направление тока из начала в конец ветви). Матрица
задающих токов узлов J является матрицей вектор-столбцом, элементами которой являются задающие токи в комплексном виде соответствующих узлов, причем знак «+» ставится перед током, если он входит в узел, «-» - если выходит из узла. Матрица сопротивлений ветвей является квадратной матрицей. В главной диагонали ^ находятся значения собственных сопротивлений ветвей в комплексном виде со знаком «+». В недиагональных элементах находятся значения взаимных сопротивлений между ветвями. Причем знак «+» ставится при согласном включении, «-» при встречном. Можно использовать и другие варианты направления векторных величин, главное, чтобы была единая система направлений для всей схемы, с учетом знаков
Технические науки — от теории к практике и'и'и'.\//)^с. ш/о_№2 (50), 2016г
в системе уравнений. В нижеизложенной методике используется следующие матричные уравнения.
Первый закон Кирхгофа [10, с. 58]:
М 1Б = }
Второй закон Кирхгофа:
N гв1Б = Ы -Ев
Закон Ома:
и В~ I в~ Е в
Расчет методом узловых потенциалов.
Матрица узловых проводимостей:
КУ=М КВ -м'=М ■ •М' Уравнение узловых напряжений [10, с. 76]:
-М -Ув-Ёв иА=1' 0 -М -1В-ЁВ)
Выражение токов ветвей через узловые напряжения:
1в=1Лм'-иА+Ёв)
Падение напряжения в ветвях схемы находится по выражению:
йв=М'йА
Потери мощности в только одной ьой ветви находятся по выражению:
=1в1-иВ1
*
Технические науки — от теории к практике № 2 (50), 2016г_
Поток мощности в i-ой ветви со стороны j-ого узла (начало или
конец ветви) определяется по выражению:
*
где: ид. - напряжение]-ого узла относительно базисного.
Расчет методом контурных токов.
Матрица контурных ЭДС:
Ек = N -Бв
Матрица контурных сопротивлений:
г к=N • г в • N
Уравнение контурных токов [10, с.74]:
г к • 1к=Ек-N • г в•
I к г к
Ек - N • гв •
М а
0
М -
0
• J
• J
Выражение токов ветвей через контурные токи:
1в = N I К + |М а\и
Падения напряжения в ветвях, входящих в дерево схемы находится по выражению:
и а аа -1 а а
Напряжения в узлах схемы относительно базисного определяется по выражению:
и Л = (М I)-1 и а
ССибАК
Технические науки — от теории к практике www.sibacinfo_№2 (50), 2016г
Практические рекомендации по расчету ТКЗ.
Расчет ТКЗ в ФК можно выполнять как методом узловых потенциалов, так и методом контурных токов. Если схема замещения содержит ветви с нулевым сопротивлением (например, соответствующие выключателю), то при расчете методом узловых потенциалов
матрицы проводимостей ветвей по выражению YВ = ZВ произойдет деление на ноль. Поэтому у такой ветви необходимо нулевое сопротивление заменять некоторой малой величиной, что может несколько снижать точность расчета. В методе контурных токов такая проблема отсутствует, если только не весь контур содержит ветви с нулевым сопротивлением. Если схема замещения содержит ветви с очень большим сопротивлением (соответствующие отключенной ветви или обрыву, когда считается режим нескольких несвязанных подсистем), то в методе контурных токов возможно снижение точности за счет погрешности математических операций, тогда лучше покажет себя метод узловых потенциалов.
Составление матриц с параметрами схемы замещения для практических расчетов ТКЗ.
При наличии готовой расчетной схемы для расчета ТКЗ в любом из промышленных программных комплексов можно для решения практических задач получить параметры схемы замещения в ФК следующим образом. Вначале поэлементно выполняем преобразования матрицы с параметрами схемы замещения одного элемента сети в СС в матрицу с параметрами схемы замещения этого же элемента сети в ФК по следующим выражениям (обращая внимание на порядок расположения в матрицах параметров соответствующих последовательностей и фаз):
Е ЛЕС ~ Eoil
J Л ВС — J 012
¿L
[111 1 .2*
где: S = i a2 a , a = e 3 , a2 = e 3 - операторы поворота вектора.
i 2 1 a a
Если имеются данные о взаимоиндукции нулевой последовательности между двумя или несколькими ЛЭП, то формируется матрица
Z012 не для одной ЛЭП (размерностью 3х3), а сразу для нескольких ЛЭП (размерностью 6х6, 9х9 и т. д.), в соответствующих недиагональных элементах указывается сопротивление взаимоиндукции нулевой последовательности между соответствующими ЛЭП.
С1
Сибдк
Технические пауки — от теории к практике №2 (50), 2016г_www.sibac.info
Далее формируется матрица параметров схемы замещения:
Е . ./ , как блочные матрицы отдельных элементов сети.
После выполнения расчета в ФК можно параметры режима в ФК конкретного элемента преобразовать в параметры режима с СС:
Г * Г / И
1 «
АБС
-1 ■
АБС
Определение параметров схемы замещения генераторов для учета предшествующего режима при расчете ТКЗ.
При разборе аварий, анализе фактических коротких замыкания (КЗ), при проверке селективности РЗА, решении задачи определения места повреждения (ОМП) для повышения точности расчета ТКЗ требуется учет предшествующего аварии установившегося режима. Это позволяет уточнить модули и фазы ЭДС для генераторов. Если предшествующий режим был тяжелым, т. е. перетоки по отдельным электропередачам приближались к допустимым пределам по устойчивости, то влияние учета фаз ЭДС может существенно повлиять на результаты расчета ТКЗ.
Особенностью математического моделирования вращающихся электрических машин (в нашем случае синхронных генераторов) является то, что общепринятая система координат q" повернута относительно общесхемной системы координат на угол 5, который является углом поворота ЭДС генератора и определяется предшествующим режимом.
Для неявнополюсных машин (турбогенераторов) все особенности учета предшествующего режима ограничиваются только определением параметров ЭДС генератора (прямой последовательности) по выражению:
Р =Еп=и
СГ( 1) О пр
и 2
пр
Активным сопротивлением генератора можно пренебречь, либо определить его через постоянную времени обмотки статора (справочный параметр):
X''
Ксг -
1 а
®0 ■ Та
ССибАК
Технические науки — от теории к практике и'и'и'.\//)^с. ш/о_№2 (50), 2016г
Соответственно комплексное сопротивление прямой последовательности генератора для расчета ТКЗ:
7 СГ (1) = КСГ + jXd
Для обратной последовательности:
7СГ (2) = jXd
Обмотка статора синхронных генераторов обычно включается в звезду, причем нулевая точка в малых машинах изолирована, а в крупных машинах с целью выполнения релейной защиты от замыканий на землю заземляется через большое сопротивление. Поэтому токи нулевой последовательности либо отсутствуют, либо весьма невелики. Соответственно если в справочных данных отсутствуют параметры нулевой последовательности, то для перевода сопротивлений генератора в ФК можно сопротивление нулевой последовательности приближенно принять равным сопротивлению обратной последовательности.
Для явнополюсных машин (гидрогенераторов) при учете предшествующего режима необходимо учесть неравенство сверхпереходных сопротивлений по осям и
Активное сопротивление гидрогенератора определяется через постоянную времени обмотки статора (справочный параметр):
2 • X"-X"
®о • (К + X" ) • т„
Сопротивление гидрогенератора прямой последовательности в комплексном виде представить невозможно (из-за неравенства сверхпереходных сопротивлений по осям и "q"), поэтому его нужно представлять в виде матрицы 2х2. Матрица сверхпереходных сопротивлений генератора в собственных осях генератора:
7"
7 с
КСГ ~Х1
X" япг