CC BY
22
УДК 621.311.018
DOI: 10.14529/ро\ег180104
ОСОБЕННОСТИ ПРИМЕНЕНИЯ МАТРИЦ
ПРЕОБРАЗОВАНИЯ И УРАВНЕНИЙ СВЯЗИ
ПРИ АНАЛИЗЕ НЕСИММЕТРИЧНЫХ ПОВРЕЖДЕНИЙ
П.Л. Воронов
Чувашский государственный университет им. И.Н. Ульянова, г. Чебоксары, Россия, ООО НПП «ЭКРА», г. Чебоксары, Россия
Тензорно-топологический метод Крона, развиваемый многими отечественными и зарубежными исследователями, открывает широкие возможности для решения прикладных задач в сложных электроэнергетических системах с помощью тензорного анализа и диакоптики. Его применение к анализу и расчету несимметричных и сложных видов повреждений, возникающих в системах электроснабжения (СЭС), имеет свои особенности, обусловленные необходимостью преобразований переменных (систем координат). Помимо задач определения матриц преобразования и уравнений связи появляются проблемы объединения отдельных частей системы, которые представлены в различных координатах. Затруднения возникают вследствие необходимости преобразования этих матриц при построении схем замещения.
В предлагаемой статье на простейших примерах рассмотрены некоторые особенности совместного применения матриц преобразования координат, матриц соединения и уравнений связи, а также их различения при моделировании и анализе несимметричных повреждений в СЭС.
Ключевые слова: несимметричные повреждения, матрицы преобразования, симметричные составляющие, схемы замещения последовательностей, уравнения связи.
Введение
Физическая структура реальных сложно-разветвленных электрических сетей и взаимосвязанных машин систем электроснабжения (СЭС) является исключительно сложной системой как с точки зрения описания ее динамических свойств посредством уравнений электромагнитного поля, так и с помощью сетевых моделей и эквивалентных электрических схем замещения. Поэтому все более настоятельным является привлечение для анализа установившихся и переходных режимов СЭС тензорно-топологического метода [1-5], который в настоящее время с развитием и совершенствованием программно-вычислительных комплексов становится по преимуществу средством инженерных расчетов.
В зависимости от целей анализа при расчетах на ЭВМ, как правило, используется исходная информация в виде заданных параметров системы и режима, а также топологической структуры или схемы замещения. Часто параметры разных элементов СЭС задаются в различных системах координат (А, В, С; <1, q, 0; 1, 2, 0; / Ь, 0 и т. д.) [6-8], поэтому при конкретных расчетах нередко приходится проводить предварительные преобразования параметров, например, к одной системе координат. Между тем даже при простом переходе от фазных переменных к симметричным составляющим и, наоборот, необходимо преобразовывать каждый элемент трехфазной цепи в отдельности. Применение метода преобразования координат к отдельным элементам СЭС и к анализу электрических сетей с распределенными источниками энергии подробно проиллюстрировано в статье [9].
Общая теория расчетов трехфазных систем при несимметричных и сложных видах повреждений дана в [10-12]. Она базируется на применении метода симметричных составляющих и построении комплексных схем замещения с промежуточными трансформаторами. В [12-16] представлены различные модификации методов анализа несимметричных и неполнофазных режимов в базисах, симметричных составляющих и в фазных координатах. В [17, 18] развит метод расчета несимметричных повреждений на основе тензорной методологии Крона, сочетающий в себе два преобразования: соединение схем последовательностей и учета следствий из граничных условий, соответствующих конкретным видам несимметрии.
Тензоры и матрицы преобразования
электрических систем
Первым шагом в развитии тензорно-топологи-ческого метода анализа электрических систем стало применение матриц и многомерных векторов. Тензоры всегда рассматривались как инвариантные объекты относительно преобразования координат, выражаемые через свои компоненты в форме матриц и векторы базиса, соответствовавшие координатным системам, в которых эти компоненты определялись. Системы координат, вводимые произвольно исследователями, устанавливают соответствие между числами и топологическими (геометрическими) пространствами. В каждой точке такого пространства в зависимости от его размерности имеют место координатные линии. Эти линии выбираются по усмотрению, но не связаны непосредственно с характеристиками изучаемого физического явления. Однако важно, что
законы изменения физических процессов, рассматриваемые исследователем относительно этого произвольного выбора систем координат, всегда инвариантны. Данное обстоятельство накладывает существенные ограничения на формулы преобразования тензорных величин при изменении систем координат. Матричный анализ и формулы (законы) преобразования переменных (тензоров) нашли в настоящее время широкое применение в прикладной электротехнике, но все же до сих пор наблюдается оппозиция слову «тензор». Это объясняется не только тем, что при анализе электрических сетей оперируют лишь тензорами 0, 1, 2 рангов, компоненты которых выражаются числами и матрицами первого и второго порядка, но и представлениями, что топологическая структура возбужденной электрической цепи должна быть непременно идентична ее ориентированному графу (структуре алгебраической топологии). На самом же деле она таковой не является, поскольку матрицы соединения и матрицы преобразования координат кардинально отличаются от матриц инци-денций графа цепи. Матрицы инциденций связывают разно-мерные пространства и автоматически удовлетворяют теореме Стокса, а не законам Кирхгофа. Конечно, сами по себе матрицы не имеют физического смысла в отличие от векторов (тензоров), представляющих собой объективную, предметную и измеряемую сущность, обладающую множеством групп, составляющих в бессчетном числе систем координат. Каждый элемент из этих групп может быть записан в форме матрицы, но она лишь столбец, таблица, многомерный куб чисел или составляющих тензора в определенной системе координат.
Любой процесс или движение определяется посредством выделяемой системы координат, но движение среды подчиняется определенным законам, причем для практики вполне достаточно знать только некоторые средние, суммарные или глобальные характеристики его, описываемые тензорными инвариантными уравнениями с переменными, принимающими обобщенные значения, например, переменными Лагранжа. Эти переменные представляют сопутствующую систему координат. Напомним, что в уравнениях Лагранжа - Максвелла электрические токи рассматриваются как обобщенные скорости, а в качестве соответствующих переменных или координат принимаются количества электрических зарядов в данный момент времени, причем эти обобщенные координаты не входят в уравнения Максвелла в явной форме. Между тем определение вектора только его компонентами является недостаточным, даже если подчеркивается, что они преобразуются при замене системы координат по определенным правилам. Приводя их значения в форме матриц, надо указывать и базис, в котором они заданы. Векторы базиса в геометрическом смысле направлены по касатель-
ным к координатным линиям. Они могут тоже изменяться при замене координат посредством соответствующих матриц преобразования. Например,
1 2
если заданы две системы координат X , X и
12
У ,У на плоскости, причем существует соответствие X' = X' (У\У2), то выполняется равенство дХ1 / дУ1 = С' ■, причем индекс 1 соответствует строке, а 1 - столбцу матрицы С'^. Взаимно обратной к ней будет матрица дУ-1 / дХк = , а произведение этих матриц является единичной матрицей Ъ\к, называемой символом Кронекера. Заметим, что обозначения с помощью индексов и точек, показывающих порядок следования индексов, существенно облегчают выполнение операций с тензорами, позволяют избежать ошибок при преобразованиях и решении тензорных уравнений. Они также указывают ранг тензора и степень его ковариантности или контравариантности. В прямой нотации используется обычная запись [С], [А]. Из определения вектора (тензора первого ранга) следует, что скорость или электрический ток - это контравариантные тензоры (у них верхние индексы), а вектор-градиент скалярной функции или электрическое напряжение - это ковариантные тензоры (у них нижние индексы). Скаляр - это тензор нулевого ранга, характеризующийся числом. Скаляром является мощность электрической сети. Помимо этих величин, связанных с возбужденной электрической сетью, инвариантными объектами относительно преобразования систем координат являются также тензоры второго ранга: полное (кажущееся) сопротивление или импеданс ] и проводимость или адмиттанс [У1' ]
сети, представляемые в различных системах координат матрицами второго порядка. Для определения их компонент при смене системы координат требуется уже две матрицы преобразования: \1'] = [С][1 ][С], [У'] = [А][У][А]. Количество матриц, необходимых для преобразования тензора, равно его рангу или валентности.
В общем случае вся совокупность матриц преобразования, являющихся ключом к тензорному анализу, обладает свойством группы [С ] = [С1][С2][С3]^, но каждая отдельная из них
представляет собой смешанный тензор второго ранга, у которого один индекс ковариантный, а второй -контравариантный. Он выражает отношение между новыми и старыми переменными и имеет особый закон (формулу) преобразования. Его преобразование осуществляется с помощью трех матриц. Например, если для токов заданы соотношения:
[12] = [С.2Ш1], [12] = [С2з][13], [1'1] = [С^р4], то из них можно найти путем последовательных под-
становок и преобразований матрицу [С.4], преобразующую токи [13] в [14]. Она имеет вид: [С.34 ] = [С.23 ]-1 [С.2! ] [С4 ]. Следовательно, полную совокупность всех матриц преобразования можно обоснованно называть тензором, поскольку каждый тензорный или геометрический объект объединяется с целой группой матриц преобразования, причем каждая новая система координат, в которой исследуется объект, характеризуется своей собственной матрицей преобразования.
Таким образом, если известно уравнение состояния физической системы в матричной форме, то оно будет таким же и для других систем той же физической природы, но только при условии существования группы матриц преобразования и замене всех величин в этом уравнении тензорами. В тензорных уравнениях, являющихся инвариантом относительно определенной группы матриц преобразования, важно всегда одновременно рассматривать три понятия: преобразование, инвариантность и группа преобразований. Исходя из этих представлений, можно ассоциировать все стационарные электрические сети как производные от некоторой исходной или «элементарной» (примитивной) сети и реализовать на основе простейших алгоритмов процесс составления их уравнений и решения по частям с помощью соответствующих матриц преобразования, элементами которых являются только значения 0, 1, -1.
Особенности матриц соединения
и преобразования координат
Во многих случаях аналитического и численного решения на ЭВМ сложно-замкнутых электрических сетей приходится сталкиваться с рядом специфических особенностей при замене координат. Дело в том, что матрица преобразования дХ1 / дY1 — С' ■, являющаяся матрицей соединения
отдельных ветвей и частей (подсистем) в результирующую систему, часто оказывается сингулярной (вырожденной и не имеющей обратной матрицы). Но поскольку переход от «элементарной» все-контурной сети, представляющей собой в простейшем случае отдельные ветви, замкнутые на себя, не требует применения обратных матриц, то возникает вопрос о правомерности отнесения таких преобразований к тензорным. Суть вопроса состоит в том, что в тензорном анализе используются несингулярные матрицы, а операции замыкания, размыкания и соединения, например, ветвей для получения множества сетей, состоящих из одной и той же их совокупности, но связанных различным образом, до работ Крона не рассматривались как преобразования координат. Однако несингулярная матрица для любой из образованного таким путем счетного множества сетей на самом деле всегда существует, если принять во внима-
ние, что все электрические сети являются ортогональными сетями [19]. В принципе любая из этого множества сеть может быть принята за «элементарную», а переход от ее токов, напряжений и сопротивлений к другой сети может быть выполнен по формулам:
[1Э] = [С.эк][/к],
К]=[с;э][иэ], (1)
^=[с;э]^ээ][с.эк].
С практической точки зрения и удобства составления матриц преобразования и вычислений рекомендуется применять, как правило, в качестве исходной «элементарной» сети все-контурную цепь. Тогда матрица соединения [С.эк] при переходе от параметров «элементарной» сети к новой сети с меньшим числом контуров выполняет помимо функции соединения ветвей еще и функцию наложения ограничений на новую сеть, изменяя ее степень свободы. В таком преобразовании она представляет собой лишь часть, относящуюся к замкнутым путям (контурам) токов, несингулярной матрицы ортогональной сети.
Одновременно с матрицами соединения при анализе повреждений в СЭС применяются и матрицы, с помощью которых непосредственно осуществляется замена переменных (координат). Чаще всего для трехфазных сетей используются матрицы преобразования от фазных координат к симметричным составляющим 1, 2, 0 и от вращающейся системы действительных координат d, q, 0 к комплексным составляющим (прямым и инверсным координатам) вращающегося поля /, Ь, 0. Эти матрицы при соблюдении условия инвариантности мощности имеют вид:
[С0,1,2 ]
[С 0, :Г,ь ] =
0 1 2
Л 1 1 1 1
в 1 а 2 а '73
С 1 а а 2
0 / Ь
0 л/2 0 0 1
d 0 1 ■ 72
q 0 1 - ■
Они, преобразуя параметры сети, принципиально отличаются от матриц соединения, содержащих в себе скрытую и полезную информацию о топологических свойствах и поведении сети. Эта информация в зашифрованном виде выражает собой законы Кирхгофа. И хотя сами матрицы соединения не входят непосредственно в уравнения состояния сети, они в отличие от параметров системы и параметров режима (тока, напряжения мощности, частоты и т. д.) характеризуют топологию сети и определяют координатные линии,
вдоль которых изменяются переменные. Матрицы же преобразования координат, не затрагивая конфигурацию (соединение ветвей) сетей, преобразуют только параметры режима, например, токи одной совокупности переменных в токи другой совокупности переменных или координат. Заметим, что их элементы могут быть комплексными числами, поэтому при выполнении преобразований
г--,0 1 2 -1*
транспонированная матрица [СА ^С ] принимается комплексно-сопряженной, причем она равна обратной матрице [С°/вС ] 1. Группы матриц преобразования, если они имеют в качестве своих компонент только целые числа в алгебраической топологии, называются группами Бетти. Но помимо таких групп преобразований в расчетах электрических сетей и особенно машин используются группы, в которых [С] являются функциями переменных (координат). Примером такого преобразования является переход от фазных переменных к переменным в координатных осях d, q, 0, жестко связанных с ротором машины, и к синхронно вращающимся осям ds, qs, 0. Процедуры такого вида преобразований для электромеханических систем разработаны на основе тензорной теории анализа обобщенных («элементарных») машин разных типов. Уравнения идеализированных обобщенных машин в различных системах координат подробно представлены и проанализированы в [20]. Заметим, что тензорный метод не уменьшает количество вычислений, но экономит время благодаря применению специальных алгоритмов и приемов, стандартизирующих и автоматизирующих процесс анализа режимов СЭС на ЭВМ, исключает необходимость изобретения новых методов для решения появляющихся сложных задач. Преимущества его особенно наглядно проявляются при анализе сложно-замкнутых сетей, содержащих электрические машины. Поскольку физические процессы в обобщенных машинах описываются с помощью эквивалентных параметров относительно их электрических и механических зажимов на основе фундаментальных предположений, то они могут быть представлены совокупностью движущихся линейных электрических сетей с сосредоточенными параметрами. В таком виде обобщенная машина отличается от «элементарной» неподвижной сети только в отношении расположения ее обмоток (токовых слоев) по двум взаимно перпендикулярным осям и в отношении их пространственного перемещения относительно друг друга. Но поскольку при таком расположении обмоток взаимный импеданс имеет место только между обмотками по одной оси, то модель такой машины можно представить в виде двух отдельных многообмоточных трансформаторов, расположенных ортогонально. Данная модель вполне допустима при определении токов короткого замыкания (КЗ) и рас-
чета других видов повреждений для начального и произвольного момента времени при условии введения расчетных ЭДС Е"ч, Е'ч, Е1, определяемых
из векторных диаграмм для соответствующих режимов работы машин. Для исследования режима качаний и переходных процессов при несимметричных повреждениях целесообразнее всего использовать координаты/, Ь, 0 [18].
Примеры использования различных матриц преобразования и соотношения между ними Проиллюстрируем на простейших примерах особенности использования различных матриц преобразования и соотношений между ними при расчетах несимметричных повреждений в СЭС.
На рис. 1 приведена трехфазная схема замещения сети с несимметричной нагрузкой (в общем случае пусть 2во ф 1рр ф ) или сети с компенсированной нейтралью, широко используемой в качестве модели для исследования режима замыкания одной фазы на землю, когда под этими сопротивлениями можно понимать емкости фаз относительно земли 2С = 1/ юС .
Рис. 1. Трехфазная схема замещения сети с однофазным замыканием
При известных параметрах 7 ветвей, обозначенных латинскими буквами, и ЭДС нетрудно в отсутствие замыкания К(1) (обозначенного пунктирной линией) составить ее уравнения для трех независимых замкнутых контуров а, р, у, выполнив преобразование от «элементарной» сети с помощью матрицы соединения
[CL -
A B C D F G N
а 1 1 1
Р 1 1 1
У 1 1 1
(2)
и формул преобразования (1). Представим эти уравнения в базисе фазных переменных в матричной форме:
еа
еР =
еу
% АЛ + %ВВ + %NN
%NN + ^ИИ +
%СС + +
'у
Решив их относительно контурных токов, можно затем определить и токи во всех ветвях. Конечно, для такой простой сети эти уравнения можно было бы записать и непосредственно, используя правила составления контурных матриц. Однако для сложно-разветвленных сетей с взаимоиндукцией между ветвями справиться с такой задачей без ошибок крайне затруднительно. Соединение ветвей этой схемы в многоконтурную сеть можно осуществить и с помощью уравнений связи [М ] •['] = 0 . Под уравнениями связи понимают соотношения между переменными, например, между токами в узле. Тогда [М] - матрица алгебраического отображения сети или узловая матрица соединений (на рис. 1 узлы связи обозначены цифрами). Для рассматриваемой сети эта матрица имеет вид:
[М ] = 2
А В С В И о N
1 1 -1
2 1 -1
3 1 -1
4 1 1 1 -1
(3)
Поскольку матрицы [С.эк ] и [М] выполняют одну и ту же функцию соединения ветвей, то должно существовать определенное соотношение между ними. В [21] такое соотношение было получено в форме двух уравнений: [М] = [МдМс ], [С;3] = [-Мй (Мд1)^ |1]. В этих уравнениях подматрица [Мд] - неособенная матрица некоторого дерева графа исследуемой схемы замещения, а [Мс ] - подматрица, соответствующая хордам этого же графа. Между тем искомое соотношение между матрицами [С.эк] и [М] можно найти существенно проще, не прибегая к теории графов. Для этого достаточно приравнять независимые контурные токи токам фаз 'а = 'А, № = 'В, iу = Iе и дополнить эти три независимых уравнения еще тремя зависимыми: 'а = '°, ^ = '¥, 'у = ^ - 'а - ^. Из полученных шести уравнений и определяется матрица (2), а также соотношение между матрицами (2) и (3).
Однако в практических расчетах, как правило, требуется обратный переход от уравнений [/э] = [С.эк][/'к] к уравнениям связи [М] [г] = 0 . В этом случае проще всего взять транспонированную матрицу (2) и представить ее посредствам перестановки строк в виде сложной матрицы, разделив её координатные оси на относящиеся к независимым и зависимым контурам (k' ,k") и образующим их ветвям (Ь', Ь"). Вычитая из нее затем единичную квадратную матрицу, находим
с;3] - [1] = ь
V г V V" V V"
V 1 0 - V 1 0 = Ь ' 0 0
Ь" С 0 V 0 1 Ь" С -1
(4)
После исключения нулевых строк матрицы (4) получаем соотношение [С ' 11] = [М ], которое нетрудно проверить. Его удобно использовать при анализе несимметричных повреждений.
Обратимся вновь к схеме (см. рис. 1) и при замыкании фазы на землю выполним ее расчет, воспользовавшись методом симметричных составляющих. Предположим, что ЭДС и сопротивления генератора уже заданы в системе симметричных составляющих, а сопротивления емкостной нагрузки 2С и сопротивление 2Ш в нейтрали даны в фазных координатах. При переходе к симметричным составляющим помимо исходной (индекс «к») и «элементарной» (индекс «э») сетей будут вводиться еще две: «элементарная» сеть последовательностей (индекс «п») для отдельных элементов исходной трехфазной сети, преобразованных к симметричным составляющим посредством матрицы [С^ ВС ], и комплексная схема замещения (индекс
«5»), определяемая обычно следствиями из граничных условий для различных повреждений. Между токами любой пары данных сетей имеют место соотношения вида:
[''э] = [С.эк][/'к], [г3] = [Сп][/'п], [г'к] = [С5]['5], [Vп] = С ]['5].
Следовательно, матрица преобразования [С." ] с учетом формулы преобразования для [С.34], приведенной выше, тогда будет [С.5 ] = [С.эп]-1[С.эк][С.к5 ]. Это преобразование выполнимо, поскольку матрица
[С.эп] несингулярная и имеет обратную. Она состоит из стольких матриц [СА'вС ], сколько трехфазных
элементов содержится в исходной сети. Зная эту матрицу, можно без каких-либо затруднений определить по формулам (1) уравнения комплексной схемы замещения в симметричных составляющих. Для этого необходимо сначала найти преобразованные к симметричным составляющим соответствующие матрицы каждого из элементов «элементарной» сети последовательностей. Тогда полная матрица сопротивлений этой сети в симметричных составляющих
1 2 0 0'' 1 2 0
1' 3Z1G
2' 3Z2G
0' 3Z0G
0'' 2Zq -Zc -Zc
1'' -Zc 2Zc -Zc
2 -Zc -Zc 2Zc
0 9Zn
Следовательно, для рассматриваемого примера матрица преобразования имеет вид
[Сп] = f
0' 1' 2' 0'' 1 2 0
A 1 1 1
B 1 a 2 a
С 1 a a2
D 1 1 1
F 1 a2 a
G 1 a a2
N 3
(5)
Чтобы совершить теперь переход к матрице [С" ], нужно найти матрицу [С.* ]. Однако ее можно определить только с помощью пока неизвестной комплексной схемы замещения. Поэтому на этом шаге целесообразно прибегнуть к уравнениям связи. Поскольку при этом осуществляется переход к токам комплексной схемы замещения в симметричных составляющих, то уравнение связи видоизменяется на уравнение
М ][С.эп][/' ] = [М'][/' ] = 0.
Оно в развернутой форме с учетом (4) и разделения матрицы (5) на четыре (отмеченные в (5) жирными линиями) подматрицы [£>,■], которые соответствуют независимым и зависимым токам ¿1 и ¿2 , приобретает вид:
[M' ][/' ] =
0 0
С -1
A D2
D 3 d4
= 0.
Решая это уравнение относительно зависимых токов, находим
¿2 =-(СЪ2 -£4)-1(С'Ц -Б3>1 = С"¿1. (6)
Таким образом, с учетом формы записи транспонированной матрицы (2) представим и получаемую из (6) матрицу [С." ] тоже в виде транспонированной
[С ] =
С
1 2 0 0'' 1 2 0
1 1 -1
2' 1 -1
0' 1 -1 -1
Тогда матрица сопротивлений, вычисленная по формуле преобразования ] = [ССП][7пп][С° ], и контурные уравнения искомой комплексной схемы, представятся в виде
1
1 2 0
1 7\о +17С 7с 7с
1 Е\о 3 3 1 I'
2 0 " 2 7с + 3 7с 7с • 2 12
0 0 3 3 0 10
0 _ 3 _ 3 700 + 3 7 С + 37N
(7)
Данной системе уравнений соответствует схема замещения, показанная на рис. 2.
Z1G 7С
Е,
Рис. 2. Схема замещения, соответствующая уравнениям (7)
Заметим, что она совпадает с аналогичной схемой, построенной для случая однократной продольной несимметрии, когда сопротивления ZC имеют место лишь в двух фазах. Конечно, можно было бы получить и другую комплексную схему (рис. 3), представив матрицу Zпп в измененном виде.
Рис. 3. Комплексная схема при замыкании одной фазы на землю
В качестве второго примера, наглядно иллюстрирующего суть совместного использования матриц преобразования координат, матриц соединения и уравнений связи, рассмотрим расчет сложного вида повреждения в трехфазной сети, представленной на рис. 4.
ел
ZA
-о о
7
ев
7в
1л
7
нг
7с
и
К2)
7Н.
1с
/иг
Рис. 4. Схема сети с двумя повреждениями
е
с
Пусть генератор, питающий нагрузку, задан параметрами прямой, обратной и нулевой последовательностей 2Ь 22, 20, нагрузка - действительными сопротивлениями фаз 2НГ, а несимметрия вызывается одновременными повреждениями: разрывом фазы А (В1-1) и двухфазным КЗ между фазами В и С (К(2)). Требуется в аналитической форме определить матрицу решения, с помощью которой можно затем вычислить токи при любых действительных параметрах системы и значениях ЭДС генератора.
Обратимся вновь к тензорному методу, в котором ключевую роль играют матрицы преобразования. Обозначим матрицу перехода к симметричным составляющим через [Сх], а матрицу соединения, характеризующую связи между элементами, через [Сс]. Поскольку параметры генератора заданы в симметричных составляющих, а решение надо найти в базисе действительных фазных координат, то сначала преобразуем матрицы сопротивлений и ЭДС генератора к фазным переменным с помощью матрицы [С]:
[Zgs ] ="
0 1 2
0 Zo
1 Zi
2 Z 2
[eGS ] =
0 1 2
0 1 2 [C ] 1 ^ 1 1 1
0 Ei 0 , [Cs ]= 73 5 1 a 2 a
C 1 a a2
к фазным переменным с помощью матрицы [С] по формулам:
[2G ] = [С5 t[2GS ]С ^ [eG ] = [С ]* [eGS ] .
В результате получаем:
[Zg ] —•
A В C
A Z0 + Z1 + Z 2 Z0 + a2 Z1 + aZ2 Z0 + aZ1 + a2 Z 2
В Z0 + aZ1 + a2Z 2 Z0 + Z1 + Z2 Z0 + a2 Z1 + aZ 2
C Z0 + a2 Z1 + aZ 2 Z0 + a2 Z1 + aZ2 Z0 + Z1 + Z 2
1 A В C
73 • E aE1 a2 E
Чтобы найти теперь матрицы сопротивлений и ЭДС в связанной (исходной) системе с повреждениями в действительных фазных координатах, необходимо воспользоваться уравнениями связи для токов и матрицей преобразования [С]. Выбрав в качестве независимых токов токи фазы В и нагрузки: = iB ; Н = /НГ , определяем через них зависимые токи: i'A = 0, С = —В - Н и находим матрицу соединения в виде
A В C НГ
В' 1 -1
НГ' -1 1
Сложив матрицы сопротивлений генератора и нагрузки и выполнив операции умножения с помощью матрицы соединения согласно формулам преобразования, приведенным выше, после ряда упрощений находим искомые матрицы сопротивлений и ЭДС всей цепи в фазной системе координат относительно независимых токов в виде:
[e]
V3
В' НГ '
В' 3(Z1 + Z2) Z1(1 - a2) + Z2(1 - a)
НГ ' Z1(1 - a) + Z2(1 - a2) Z0 + Z1 + Z2 + 3Z НГ
В' НГ '
a(1 - a)E1 -a2 E
G
Выводы
1. Совместное использование уравнений связи матриц соединения и преобразования координат расширяет возможности метода симметричных составляющих для расчета несимметричных повреждений.
2. Изложенная методика расчетов позволяет находить комплексные схемы замещения для различных видов несимметрии в электрических системах.
3. Она эффективно сочетается с методом расчета сложных систем и схем отдельных последовательностей по частям.
Литература
1. Крон, Г. Тензорный анализ сетей /Г. Крон. -М.: Сов. радио, 1978. - 720 с.
2. Happ, H.H. Foundations of tensor network theory / H.H. Happ // Journal of the Franklin Institute. - 1968. - Vol. 286, no. 6. - P. 561-564. DOI: 10.1016/0016-0032(68)90392-X
3. Петров, А.Е. Тензорный метод двойственных сетей / А.Е. Петров. - М.: ООО ЦИТиП, 2007. - 602 с.
4. Сохор, Ю.Н. Тензорный анализ сетей и диа-коптика в инженерных расчетах / Ю.Н. Сохор. -М.: LAP Lambert Academic Publishing, 2012. - 200 с.
5. Amari, S. Topological foundations of Kron's Tearing of electrical networks / S Amari // RAAG Memoirs. - 1962. - Vol. 3. - P. 322-350.
6. Park, R.H. Two-reaction theory of synchronous machines. Part II / R.H. Park // Trans. AIEE. - 1933. - Vol. 52. - P. 352-354. DOI: 10.1109/T-AIEE.1933.5056309
7. Clarke, E. Circuit analysis of A-C power systems /E. Clarke, J. Wiley. - New York, 1943. - Vol. 1. -540 p.
8. Ky, Y.H. Transient analysis of rotating machines and stationary networks by means of rotating reference frames. Part I / Y.H. Ky // Trans. AIEE. - 1951. - Vol. 70. - P. 943-957. DOI: 10.1109/T-AIEE.1951.5060505
9. Воронов, П.Л. Применение метода преобразования координат к анализу электрических сетей с распределенными источниками энергии / П.Л. Воронов, В.А. Щедрин // Региональная энергетика и электротехника: проблемы и решения: сб. науч. тр. - 2015. - Вып. XI. - С. 42-65.
10. Ульянов, С.А. Электромагнитные переходные процессы в электрических системах / С.А. Ульянов. -М.: Энергия, 1964. - 704 с.
11. Han, Z.X. Generalized Method of Analysis of Simultaneous Faults in Electric Power Systems / Z.X. Han // IEEE Transactions on Power Apparatus
and Systems. - 1982. - Vol. 101, no. 10. - P. 39333942. DOI: 10.1109/TPAS.1982.317045
12. Лосев, С.Б. Вычисление электрических величин в несимметричных режимах электрических систем / С.Б. Лосев, А.Б. Чернин. - М.: Энерго-атомиздат, 1983. - 528 с.
13. Авербух, А.М. Примеры расчетов непол-нофазных режимов и коротких замыканий / А.М. Авербух. - Л.: Энергия, 1979. - 184 с.
14. Лямец, Ю.Я. Эквивалентирование многопроводных систем при замыканиях и обрывах части проводов / Ю.Я. Лямец, Д.Г. Еремеев, Г.С. Нудельман //Электричество. - 2003. - № 11. -С. 17-27.
15. Anderson, P.M. Analysis of Faulted Power Systems / P.M. Anderson. - Ames, Iowa: The Iowa State University Press, 1973. - 513 p.
16. Anderson, P.M. Analysis of Simultaneous Faults by Two-Port Network Theory / P.M. Anderson // IEEE Transactions on Power Apparatus and Systems. - 1971. - Vol. 90, no. 5. - pp. 2199-2205. DOI: 10.1109/TPAS.1971.293041
17. Lewis, W.E. Matrix Methods for the Evaluation of Simultaneous Faults in Three-phase Systems / W.E. Lewis, J.H. Banks // Institution Monograph. -1955. - pp. 231-243.
18. Воронов, П.Л. Расчет по частям трехфазных сетей при несимметричных и сложных видах повреждений / П.Л. Воронов // Вестник Чувашского университета. - 2017. - № 1. -С. 76-87.
19. Kron, G. Diacoptics: The Piecewise Solution of Large-Scale Systems / G. Kron. - London, MacDonald, 1965. - 166p.
20. Уайт, Д. Электромеханическое преобразование энергии / Д. Уайт, Г. Вудсон. - М.: Энергия, 1964. - 537 с.
21. Ионкина, П.А. Теоретические основы электротехники / П.А. Ионкина. - М.: Высшая школа, 1976. - Т. 1. - 544 с.
Воронов Павел Леонидович, аспирант кафедры электроснабжения промышленных предприятий им. А.А. Федорова, Чувашский государственный университет им. И.Н. Ульянова; инженер отдела разработки станционного оборудования, ООО НПП «ЭКРА», г. Чебоксары; plv911@mail.ru.
Поступила в редакцию 26 ноября 2017 г.
DOI: 10.14529/power180104
ON USING TRASFORMATION MATRICES AND CONSTRAINT EQUATIONS FOR ANALYSIS OF NON-SYMMETRIC FAULTS
P.L. Voronov, plv911@mail.ru
I.N. Ulyanov Chuvash State University, Cheboksary, Russian Federation,
EKRA Research and Production Enterprise Ltd, Cheboksary, Russian Federation
Kron's tensor-topological method being advanced by many Russian and foreign researchers opens wide opportunities for solving applied problems in complex electric power systems by means of tensor analysis and diakoptics. Its application to the analysis and calculation of asymmetric and complex types of faults occurring in power supply systems has its own peculiarities which are due to the need for transformations of variables (coordinate systems). In addition to the problems of defining transformation matrices and constraint equations, there are problems of combining individual parts of the system represented in different coordinates. Difficulties arise from the need to transform these matrices when constructing equivalent circuits.
The paper discusses simplest examples to show some features of how coordinate transformation matrices, coupling matrices and constraint equations can be jointly applied, as well as how they can be differentiated when modeling and analyzing non-symmetric faults in power supply systems.
Keywords: non-symmetric faults, transformation matrices, symmetric components, sequence equivalent circuits, coupling equations.
References
1. Kron G. Tenzornyy analiz setey [Tensor Analysis of Networks]. Moscow, Sov. Radio Publ., 1978. 720 p.
2. Happ H.H. Foundations of Tensor Network Theory. Journal of the Franklin Institute, 1968, vol. 286, no. 6, pp. 561-564. DOI: 10.1016/0016-0032(68)90392-X
3. Petrov A.E. Tenzornyy metod dvoystvennykh setey [Tensor Method of Dual Networks]. Moscow, TsITiP Publ., 2007. 602 p.
4. Sohor Ju.N. Tenzornyy analiz setey i diakoptika v inzhenernykh raschetakh [Tensor Analysis of Networks and Diakoptics in Engineering Calculations]. Moscow, Lambert Academic Publishing, 2012. 200 p.
5. Amari S. Topological Foundations of Kron's Tearing of Electrical Networks. RAAG Memoirs, 1962, vol. 3, pp. 322-350.
6. Park R.H. Two-Reaction Theory of Synchronous Machines. Part II. Trans. AIEE, 1933, vol. 52, pp. 352354. DOI: 10.1109/T-AIEE.1933.5056309
7. Clarke E., Wiley J. Circuit Analysis of A-C Power Systems. New York, 1943, vol. 1. 540 p.
8. Ky Y.H. Transient Analysis of Rotating Machines and Stationary Networks by Means of Rotating Reference Frames. Part I. Trans. AIEE, 1951, vol. 70, pp. 943-957. DOI: 10.1109/T-AIEE.1951.5060505
9. Voronov P.L., Shchedrin V.A. [Application of the Coordinate Transformation Method to the Analysis of Electrical Networks with Distributed Sources of Energy]. Regional'naya energetika i elektrotekhnika: problemy i resheniya: sb. nauch. tr. [Regional Power Engineering and Electrical Engineering: Problems and Solutions: a Collection of Sci. Papers], 2015, iss. XI, pp. 42-65. (in Russ.)
10. Ul'yanov S.A. Elektromagnitnye perekhodnye protsessy v elektricheskikh sistemakh [Electromagnetic Transients in Electrical Systems]. Moscow, Energiya Publ., 1964. 704 p.
11. Han Z.X. Generalized Method of Analysis of Simultaneous Faults in Electric Power Systems. IEEE Transactions on Power Apparatus and Systems, vol. 101, no. 10, pp. 3933-3942. DOI: 10.1109/TPAS.1982.317045
12. Losev S.B., Chernin A.B. Vychislenie elektricheskikh velichin v nesimmetrichnykh rezhimakh elektri-cheskikh sistem [Calculation of Electrical Values in Nonsymmetric States of Electrical Systems]. Moscow, Energoatomizdat Publ., 1983. 528 p.
13. Averbuh A.M. Primery raschetov nepolnofaznykh rezhimov i korotkikh zamykaniy [Examples of Calculating Incomplete Phase States and Short Circuits]. Leningrad, Energiya Publ., 1979. 184 p.
14. Lyamets, Yu.Ya., Eremeev D.G., Nudel'man G.S. [Equivalence of Multi-Wire Systems in the Case of Short-Circuits and Partial Wire Breakage]. Elektrichestvo [Electricity], 2003, no. 11, pp. 17-27. (in Russ.)
15. Anderson P.M. Analysis of Faulted Power Systems. Ames, Iowa: The Iowa State University Press, 1973. 513 p.
16. Anderson P.M. Analysis of Simultaneous Faults by Two-Port Network Theory. IEEE Transactions on Power Apparatus and Systems, 1971, vol. 90, no. 5, pp. 2199-2205. DOI: 10.1109/TPAS.1971.293041
17. Lewis W.E., Banks J.H. Matrix Methods for the Evaluation of Simultaneous Faults in Three-Phase Systems. Institution Monograph, 1955. pp. 231-243.
18. Voronov P.L. [Piecewise Calculation of Three-Phase Networks with Asymmetric and Complex Types of Damage]. Bulletin of the Chuvash University, 2017, iss. 1, pp. 76-87. (in Russ.)
19. Kron G. Diacoptics: The Piecewise Solution of Large-Scale Systems. London, MacDonald, 1965. 166 p.
20. White D., Woodson G. Elektromekhanicheskoe preobrazovanie energii [Electromechanical Energy Conversion]. Moscow, Energiya Publ., 1964. 537 p.
21. Ionkina P.A. Teoreticheskie osnovy elektrotekhniki [Theoretical Foundations of Electrical Engineering]. Moscow, Vysshaya Shkola Publ., 1976, vol. 1. 544 p.
Received 26 November 2017
ОБРАЗЕЦ ЦИТИРОВАНИЯ
FOR CITATION
Воронов, П.Л. Особенности применения матриц преобразования и уравнений связи при анализе несимметричных повреждений / П.Л. Воронов // Вестник ЮУрГУ. Серия «Энергетика». - 2018. - Т. 18, № 1. - С. 27-37. DOI: 10.14529/power180104
Voronov P.L. On Using Trasformation Matrices and Constraint Equations for Analysis of Non-Symmetric Faults. Bulletin of the South Ural State University. Ser. Power Engineering, 2018, vol. 18, no. 1, pp. 27-37. (in Russ.) DOI: 10.14529/power180104
