CC BY
245
УДК 621.311 Дата подачи статьи: 11.12.15
Б01: 10.15827/0236-235Х.113.126-133
ОБ АДРЕСНОМ РАСПРЕДЕЛЕНИИ МОЩНОСТЕЙ В ЭЛЕКТРИЧЕСКИХ СЕТЯХ
А.Б. Баламетов, д.т.н., профессор, руководитель научного направления, действительный член Академии электротехнических наук РФ, balametov.azniie@gmail.com; Э.Д. Халилов, к.т.н.., доцент, ведущий научный сотрудник, elmanxalilov2010@mail.ru; Т.М. Исаева, к.т.н., доцент, taranaiscteva@r~amЫer; Ф.Г. Исгендеров, аспирант, texnikishert@mcil.ru (Азербайджанский научно-исследовательский и проектно-изыскательский институт энергетики, ул. Зардаби Гасанбека, 94, г. Баку, Аз1012, Азербайджан; Азербайджанский государственный университет нефти и промышленности, просп. Азадлыг, 20, г. Баку, А21010, Азербайджан)
Реформы отрасли электроснабжения, возможность широкого использования возобновляемых источников электроэнергии, разделение услуг по передаче в электроснабжении, определение источников искажений синусоидальности тока и напряжений делают необходимым исследование задачи распределения потоков мощности и энергии между участниками рынка. В рыночных условиях стоимость электроэнергии у генерирующих станций и других поставщиков может быть разная, поэтому актуально определение узловых цен, в которые включается составляющая потерь в сети. Процесс отслеживания потоков электроэнергии позволяет устанавливать маршруты, связывающие конкретные генераторы с конкретными потребителями, и определять долевые вклады генераторов в покрытие нагрузки с учетом адресного распределения. Метод адресности позволяет также определять цену транзита электроэнергии. Правила присоединения новых потребителей к электрической сети обусловливают необходимость структурного анализа электроэнергетической системы на высших гармониках. В статье исследовано адресное разделение потоков и потерь активной мощности между участниками рынка электроэнергии. Разработаны алгоритм и программа для решения задачи адресности потоков и потерь мощности. В разработанной программе Z-матрица используется как для решения уравнений установившихся режимов, так и для расчета узловых потерь. Для исследования адресной передачи активных мощностей были использованы 6-узловая тестовая схема IEEE, 8-узловая и 26-узловая эквивалентные схемы Азербайджанской электроэнергетической системы. Результаты расчетов адресного распределения можно использовать при заключении двусторонних контрактов, для расчетов узловых цен, определения стоимости передачи и т.д.
Ключевые слова: потоки мощности, потери мощности, адресное распределение, моделирование, матрица, маргинальный метод, Z-узловой метод, программное обеспечение.
Внедрение конкурентных отношений на рынке электроэнергии требует решения ряда новых задач, например, о доле участия каждой станции в снабжении конкретной нагрузки и определении потоков мощности, текущих от каждого генератора по ветвям схемы замещения электрической сети, то есть какие генераторы снабжают конкретного потребителя, какая величина мощности каждого генератора передается по конкретным линиям электропередачи и какие потери мощности обусловлены передачей электроэнергии от каждого генератора к конкретной нагрузке.
В концепции конкурентного оптового рынка электроэнергии принято считать, что в интегрированной электроэнергетической системе (ЭЭС) невозможно проследить маршруты поставки электроэнергии от конкретного поставщика к конкретному потребителю. Однако в конце прошлого века в научных публикациях были предложены алгоритмы, основанные на использовании направлений потоков мощности и принципа пропорциональности [1-6].
Положенный в основу этих алгоритмов принцип пропорциональности основан на предположении, что в любом узле электрической сети каждая единица мощности выходит из узла в такой же пропорции [1], в какой мощность входит в этот узел.
Определение адресности является важным также для оценки
- показателей использования электрической сети различными участниками рынка;
- абонентской платы за пользование сетью вдоль пути протекания мощности;
- ответственности источников и потребителей электроэнергии за соблюдение режима (нахождение переменных режима в заданных технологических пределах);
- изменения топологии схемы сети таким образом, чтобы вся или большая часть производимой на определенной станции мощности поступала определенному потребителю: алгоритм адресности в этом случае будет определять эффективность принимаемых решений (при заключении прямых контрактов между производителем и потребителем на передачу электроэнергии);
- локальных балансов реактивной мощности.
В последние годы появилось большое число
публикаций, использующих как матричные, так и графовые алгоритмы решения задачи адресности. Решение проблемы адресности, использующее матричный подход, получено в [1].
Решение проблемы адресности реактивной мощности позволяет проследить пути передачи ре-
Рис. 1. Блок-схема решения задачи адресности потоков мощности Fig. 1. A flow chart for solving a power flow adressness problem
активной мощности от генераторов: проанализировать эффективность использования имеющихся в системе и установки дополнительных компенсирующих устройств.
На рисунке 1 приведена блок схема решения задачи адресности потоков мощности.
Прослеживание путей передачи мощности в сети может быть выполнено против их потоков, то есть от нагрузок к генераторам (прямое отслеживание), или вдоль потоков, то есть от генераторов к нагрузкам (обратное отслеживание).
Распределение активных мощностей без учета потерь
Задача определения коэффициентов адресности по значениям перетоков и инъекции активной или реактивной мощности является линейной задачей, использующей 1-й закон Кирхгофа. Прослеживание путей передачи мощности в сети может быть выполнено против их потоков, то есть от нагрузок к генераторам, или вдоль потоков, то есть от генераторов к нагрузкам.
Первый закон Кирхгофа можно представить в матричной форме: М х /в+ 3 = 0, где /в - матрица токов ветвей; 3 - матрица узловых (задающих) токов.
При решении задачи адресности электрическая сеть представляется в виде ориентированного графа, являющегося совокупностью узлов и ветвей, связанных между собой.
Рассмотрим сеть, состоящую из п узлов и т ветвей, и определим РУ как (пх1)-вектор узловых потоков мощностей (то есть сумма из узловых притоков или оттоков, включая местную генерацию и нагрузку соответственно), РОГ как (пх1)-вектор узловых генераций, РН как (пх1)-вектор узловых нагрузок, РВ как (тх1)-вектор потоков ветвей.
Для решения задачи адресности необходимо записать исходные данные, которые будут использоваться в расчете, а именно: столбцовые матрицы задающих генерирующих мощностей Рг, задающих нагрузочных мощностей РН, потоков мощностей в каждой ветви электрической сети РВ. Если узел генераторный, то записывается значение генерации, в других случаях значения равны нулю.
Первая матрица инциденции М, составленная для всех узлов схемы, включая и балансирующий узел, задает конфигурацию электрической сети. Число строк матрицы М равно количеству узлов, включая и балансирующий, а число столбцов соответствует количеству ветвей в сети. Таким образом, первая матрица инциденций М устанавливает взаимосвязь между узлами и ветвями направлен-
ного графа сети. Каждый элемент матрицы Mis, находящийся на пересечении i-й строки, соответствующей i-му узлу, и s-го столбца, соответствующего (—¿)-й ветви, равен:
1, если узел i - начало ветви s, Mis = <! 0, если узел i не относится к ветви s, -1, если узел i - конец ветви s.
Первая матрица инциденций M может быть представлена в виде суммы двух матриц [1]:
М=Мпр+Моб, (1)
где Мпр - матрица соединения ветвей во всех узлах схемы, элементами которой служат нули или единицы с положительными знаками, соответствует началу линии (обозначает направление потока); Моб - матрица, содержащая в качестве своих элементов нули или единицы с отрицательными знаками, соответствует концу линии.
Матрица смежности ориентированного графа С (nxn) определяется как матрица [С]¿=1, если поток направлен от узла i к узлу j, и может быть вычислена как
с = -МТр • Моб. (2)
Матрица Рпр может быть вычислена как
Рр =-МТр • diag(PB) • Моб, (3)
где diag(PB) - диагональная матрица потоков мощностей ветвей.
Матрица суммарных мощностей в узлах Ру может быть определена как сумма узловых оттоков или притоков:
Ру=Рн+Рпр • 1 или Ру=Рг+Рпр-1, (4)
где 1 - единичная матрица nx1, состоящая из единиц, число строк которой равно количеству всех узлов рассматриваемой сети.
Подставляя (3) в (4), после преобразований можно получить
Рн = Ру - Рпр • 1 = Ру + МТр • diag (Рв ) • Мо6 • 1 = = [l + Мтщ • diag (Рв ) • М0б • diag (Ру )-1 ]• Ру = (5)
= Ар • Ру,
Рг = Ру - РТ; ^ = Ру + Моб • diag (Рв ) • Мпр 4 = = [/ +Мто6 • diag (Рв ) • Мпр • diag (Ру )-1 ]• Ру = (6)
= Аоб • Ру,
предполагая, что узлы без узловых потоков, то есть те, для которых P,=0, были удалены из ориентированного графа так, чтобы (diagP)-1 существовала. Здесь I - единичная матрица nxn, состоящая из единиц в главной диагонали.
Матрицы Лпр и Лоб прослеживания путей вдоль и против потоков равны
Ар = I + МТр • diag (Рв ) • Моб • diag (Ру )-1,
Ло6 = I + Mis • diag(PB ) • Мпр • diag(Py)"'.
(7)
Матрицы ^4пр и Аоб могут быть сформированы непосредственно из потоков ветвей без использования матричных преобразований, показанных в (7).
Уравнения (5) и (6) позволяют проследить, как потоки мощности протекают в сети от отдельных генераторов или нагрузок:
Ру = Лр • Рн или Ру = ■ Рг . (8)
Доля узлового потока от местной генерации Рп равна Рп/Ри Следовательно, используя первое уравнение (8), локальную генерацию можно выразить как линейную комбинацию компонентов, доставленных отдельным нагрузкам:
P р Х[Ар]ik 'PHt • pi k=i
(9)
Моделирование задачи адресности
Для моделирования задачи адресности разработаны методика и алгоритм, реализованные в виде MATLAB-программы, работоспособность которой проверена на различных тестовых схемах IEEE.
Алгоритм адресного распределения потоков активной мощности включает семь этапов.
1. Составление схемы замещения электроэнергии.
2. Расчет установившегося режима ЭС.
3. Информация об узлах (генерация, нагрузка) и ветвях (нумерация и потоки мощности).
4. Формирование матриц узловых инъекций мощности, потоков мощности.
5. Формирование матриц инциденций.
6. Расчет коэффициентов адресности.
7. Определение вкладов генерации в потребление.
В отличие от адресного распределения потоков активной мощности моделирование адресности передачи реактивной мощности имеет следующие особенности:
- необходимость учета емкостных проводи-мостей ветвей на землю;
- необходимость учета индуктивных прово-димостей трансформаторов;
- увеличение числа узлов по сравнению с числом генераторов в схеме сети;
- увеличение числа нагрузочных узлов по сравнению с нагрузочными узлами.
Для учета особенностей адресности передачи реактивной мощности были разработаны дополнительные алгоритмы формирования фиктивных источников реактивной мощности, стоков реактивной мощности, фиктивных ветвей.
Для исследования адресной передачи активных мощностей были использованы 6-узловая тестовая схема IEEE, 8-узловая и 26-узловая эквивалентные схемы Азербайджанской ЭЭС.
Результаты моделирования
На рисунке 2 приведена 6-узловая тестовая схема IEEE. Данные по узлам и ветвям, а также результаты расчета установившегося режима тестовой схемы приведены в таблицах 1 и 2.
На рисунке 3 приведены результаты расчета адресности и орграф 6-узловой схемы, полученный по разработанной программе.
Таблица 1
Данные по узлам тестовой 6-узловой схемы
Table 1
Data for a 6bus test scheme
Таблица 2
Данные по ветвям тестовой 6-узловой схемы
Table 2
Data on the branches of a 6-noded test scheme
Результаты расчетов адресности активной мощности тестовой 6-узловой схемы по разработанной в среде МЛТЬЛБ программе приведены в таблице 3.
Таблица 3
Результаты расчетов адресности
Table 3
Adressness calculation results
Станция
1 2 3 4 5 6
1 100 0 0 0 0 0
2 25.96 50 0 0 0 0
5 F 3 0.76212 1.4679 60 0 0 0
4 52.365 21.505 0 0 0
5 41.35 12.242 16.468 0 0 0
6 9.0604 17.393 43.537 0 0 0
Узел
Узел 1 2 3 4 5 6
1 100 0 0 0 0 0
2 25.96 75.96 0 0 0 0
3 0.76212 2.23 62.23 0 0 0
4 52.365 32.67 0 73.87 0 0
5 41.35 18.598 17.08 3.88 70.06 0
6 9.0604 26.424 45.155 0.0033229 0.06 69.99
Станция
Нагрузка 1 2 3 4 5 6
1 0 0 0 0 0 0
2 0 0 0 0 0 0
3 0 0 0 0 0 0
4 49.622 20.378 0 0 0 0
5 41.315 12.232 1 6 .454 0 0 0
6 9.0617 17.396 43 . 543 0 0 0
Адресное распределение потерь между участниками рынка электроэнергии
Оценка потерь в электрических сетях различного класса напряжения и от транзита мощности, адресное разнесение потерь между продавцами и покупателями, нормирование потерь и т.д. являются приоритетными задачами [5, 7-10].
Для обоснования взаиморасчетов между получателями электроэнергии и владельцами транзитных сетей необходимо определить общий объем транзитных потерь в сети и распределить его между участниками энергообмена.
Таким образом, при распределении потерь активной мощности между участниками рынка электроэнергии возникают вопросы о доле участия электрических станций и потребителей электрической сети. Эти задачи возникают как на стадии заключения договоров, так и при оперативном планировании и на оптовом рынке электроэнергии и мощности.
В мировой практике реализации конкурентных отношений в электроэнергетике пока не существует единого подхода к этому вопросу.
Номер Напряжение, о.е. Рген, МВт Q^H, МВАр Рнаг, МВт Qm!', МВАр
1 1.050 100 33.84 0 0
2 1.050 50.00 59.16 0 0
3 1.070 60.00 83.76 0 0
4 1.007 0 0 70.00 70.00
5 1.007 0 0 70.00 70.00
6 1.018 0 0 70.00 70.00
Номер От узла К узлу От узла P(MW) От узла Q(MVAr) К узлу P(MW) К узлу Q(MVAr)
1 1 2 25.96 -1.59 -25.96 -1.59
2 1 4 41.20 22.03 -41.20 -22.12
3 1 5 32.83 13.40 -32.83 -16.05
4 2 3 2.23 -11.70 -2.23 5.13
5 2 4 32.67 44.67 -32.67 -43.91
6 2 5 16.27 13.34 -16.27 -16.20
Рис. 3. Орграф 6-узловой схемы Fig. 3. 6bus scheme orgraph
Суть пропорционального метода состоит в том, что суммарные потери между нагрузками и генераторами делятся поровну. Недостатками данного метода являются игнорирование сети и несправедливое распределение ответственности за потери между генераторами, далекими от центра нагрузок, и генераторами, расположенными в центре нагрузок.
Маргинальный метод заключается в использовании для разделения потерь коэффициентов чувствительности суммарных потерь к изменению генераций и нагрузок в узлах. Недостатком маргинального метода является то, что он не учитывает базисный узел, в нем велик дисбаланс между потерями, отнесенными к генераторным и нагрузочным узлам.
К числу специфических особенностей распределения потерь 2у-узловым методом относятся отрицательные доли потерь мощности. Существуют специальные приемы для исключения отрицательных потерь.
В [6] рассматривается метод определения коэффициентов адресности, основанный на поиске путей на ориентированном графе и позволяющий за число циклов, равное числу генераторных узлов, определить долю активной или реактивной мощности, текущей от каждого генераторного узла по ветвям схемы сети и в нагрузочные узлы.
Моделирование распределения потерь мощности
Матрица узловых сопротивлений 2у является обобщенным параметром электрической сети.
Обобщенные параметры находятся из решения системы уравнений, описывающей установившийся режим электрической сети в целом. Эти параметры могут быть получены на основе матричного и топологического методов [5].
Как известно, задача расчета установившегося режима заключается в определении значений узловых напряжений для заданных параметров электрической сети и значений узловых мощностей.
Для разделения потерь может быть использована матрица узловых сопротивлений:
"UБУ "
U' _ = + 1 Zy •
UБУ _ L y J _ I"
ДРЕ= Яе
V '=1 J=1
В целях моделирования распределения потерь 2у узловым методом разработан алгоритм, реализованный в виде программы, который опробован на различных тестовых и реальных электроэнергетических системах для локализации потерь при передаче. Алгоритм моделирования распределения потерь 2у узловым методом состоит из следующих этапов [10]:
- ввод информации;
- обращение 2у =(7у)-1-матрицы;
- решение системы линейных уравнений;
- вычисление потерь мощности, относимых в узлы.
Особенность методов матрицы 2у состоит в том, что надо получить, запомнить и применить в итерационных расчетах заполненную матрицу 2у. В связи с прямым решением уравнений сети напря-
жение каждого узла связано с токами всех узлов, что обусловливает быструю и надежную сходимость по сравнению с итеративными методами матрицы Уу. Напряжения узлов пересчитываются последовательными приближениями. 2у-матрица определяется в явном виде. В разработанной программе 2-матрица используется как для решения уравнений установившихся режимов, так и для расчета узловых потерь (рис. 4).
Рис. 4. Блок-схема разработанного ПО Fig. 4. A flow chart of developed software
Для проверки разработанного программного обеспечения были использованы тестовые схемы IEEE и эквивалентные 8- и 26-узловые схемы ЭЭС Азербайджана (рис. 5).
Проведены расчеты для эквивалентной 8-узло-вой схемы ЭЭС Азербайджана. Результаты этих расчетов по разработанной программе приведены в таблице 4.
Таблица 4
Потери мощности, отнесенные в узлы Zy-узловым методом
Table 4
Power losses that are distributed into nodes by Zybus method
Номер узла Потери мощности, МВт
2 -0,2120
3 44,6629
4 0,0225
5 15,6515
6 4,7038
7 -5,7621
8 7.3527
Суммарные потери мощности в ЭЭС составили 66.4192 МВт.
Таким образом, в нагрузочных узлах 3, 5, 6 и 8 имеются наибольшие доли потерь.
Разделение потерь мощности маргинальным методом
Маргинальные коэффициенты могут быть получены либо численным методом, либо на основе прямого вычисления частных производных. В численном методе с помощью программы расчета установившегося режима определяется изменение суммарных потерь при изменении нагрузки или генерации в узле на небольшую величину:
п п
АР = У а Р , АР = У а Р ,
/ I ген I ген I? / < наг I наг I >
1=1 1=1
~5АР дАР ' 'дАР дАР '
5Р dQ _ _ 56 5U _
Разработана программа распределения потерь активной мощности по частным производным потерь от активных мощностей. В таблице 5 приведены результаты распределения потерь активной мощности между генераторными и нагрузочными узлами для прямого определения частных производных от потерь на примере тестовой 6-узловой схемы IEEE.
Таблица 5
Распределение потерь активной мощности между генераторными и нагрузочными узлами
Table 5
Active power losses distribution between power generation buses and PQ buses
Производные от потерь (табл. 5) показывают, что генераторные узлы 2 и 3 получают компенсацию за потери. Наибольшее участие в потерях принимает 5-й узел. Участие в потерях нагрузочного узла 4 минимально.
На основании изложенного сделаем следующие выводы. Децентрализация процесса производства, передачи и распределения электроэнергии в электроэнергетике привела к необходимости адресного разделения потоков мощности и энергии между участниками рынка электроэнергии.
Но- Значе- Узло- Потеря мощно- Скоррек- Потеря, вы-
мер ние вая сти, вычислен- тирован- численная
узла част- мощ- ная по марги- ныи мар- по нормали-
ного ность, нальным коэф- гиналь- зованным
произ- МВт фициентам ный ко- коэффициен-
вод- эффици- там
ного ент по-
потерь терь
2 0,0451 50,0 2,2530 -0,03463 -1,7314
3 0,0542 60,0 3,2491 -0,04161 -2,4968
4 0,0671 70,0 -4,6997 -0,05159 3,6116
5 0,0804 70,0 -5,6267 -0,06177 4,3239
6 0,0774 70,0 -5,4181 -0,05948 4,1636
ЛРЕ= -10,2424 ЛРЕ= 7,8709
(-ъ) Azsrbaycan
Чг^ 2 lES
330\220 kV Parsabad YST
J
§irvan lES
330 kV
220\110 кВ Ssngsçal YST
Azsrbaycan lES
Ô
Sangaçal lES
£
Abçeron 4
YST330 kV _
6
T
6
220/110 кВ Mü§viq YST
500 kV
330 kV
500\330\220 kV Abçeron YST
220 kV
220\110 кВ Xirdalan YST
Рис. 5. Эквивалентная 8-узловая схема ЭЭС Азербайджана Fig. 5. An equivalent 8bus scheme of Azerbaijan electric-power system
8
3
7
5
9
В работе проведен анализ методов адресности потоков и потерь активных мощностей.
Для моделирования задачи адресности разработан алгоритм, реализованный в виде MATLAB-программы. Кроме того, проведена оценка адресного разделения потоков и потерь активной мощности на тестовых схемах IEEE и Азербайджанской ЭЭС с использованием разработанного программного обеспечения.
Результаты моделирования адресного разделения потоков мощности и энергии позволяют определить долю участия каждой станции в снабжении конкретной нагрузки и узла.
Литература
1. Bialek J. Tracing the flow of electricity. IEEE Proc. Generation, Transmission and Distribution, 1996, vol. 143, pp. 313-320.
2. Kirschen D., Allan R., and Strba G. Contributions of individual generators to loads and flows. IEEE Trans. Power Systems, 1997, 12, February pp. 52-60.
3. Bialek J.W., Kattuman, P.A. Proportional sharing
assumption in tracing methodology, Generation, Transmission and Distribution. IEEE Proc., 2004, vol. 151, iss. 4, pp. 526-532.
4. Conejo J.J., Arroyo M., Alguacil N. and Guijarro A.L. Transmission Loss Allocation: A Comparison of Different Practical Algorithms. IEEE Trans. Power Systems, 2002, vol. 17, no. 3, pp. 571-576.
5. Conejo A.J., Galiana F.D., Kockar I. Z-bus loss allocation. IEEE Trans, pp. 105-110.
6. Гамм А.З., Голуб И.И. Адресность передачи активных и реактивных мощностей в электроэнергетической системе // Электричество. 2003. № 3. С. 16-19.
7. Бартоломей П.И., Паниковская Т.Ю., Тихонов С.А. Анализ и разработка алгоритмов оценки ответственности субъектов рынка за потери электроэнергии // Энергосистема: управление, конкуренция, образование: сб. докл. III междунар. науч.-практич. конф. Екатеринбург. 2008. Т. 1. С. 623-627.
8. Баламетов А.Б., Халилов Э.Д. О прослеживании путей передачи мощностей в энергосистемах // Проблемы энергетики. 2011. № 4. С. 10-17.
9. Баламетов А.Б., Халилов Э.Д. О распределении потерь Z-узловым методом // Проблемы энергетики. 2011. N° 2. С. 30-38.
10. Баламетов А.Б., Халилов Э.Д. Методы прогнозирования потерь мощности электроэнергетических систем для анализа и управления // Электричество. 2013. № 7. С. 19-29.
DOI: 10.15827/0236-235X.113.126-133 Received 11.12.15
ON POWER TRACING DISTRIBUTION IN ELECTRIC NETWORKS Balametov A.B., Dr.Sc. (Engineering), Professor, Head of a Scientific Direction, balametov.azniie@gmail.com; Khalilov E.D., Ph.D. (Engineering), Associate Professor, Leading Researcher, elmanxalilov2010@mail.ru; Isaeva T.M., Ph.D. (Engineering), Associate Professor, taranaisaeva@rambler; Isgenderov F.G., Postgraduate Student, texnikishert@mail.ru;
(Azerbaijan Scientific-Research and Design-Prospecting Power Engineering Institute, Zardabi St. 94, Baku, Аз1012, Azerbaijan;
Azerbaijan State Oil and Industry University, Azadlyg Ave. 20, Baku, AZ1010, Azerbaijan)
Abstract. A research on the problem of power and energy flaws distribution between market participants is nesses-sary due to reforms of power supply sector, the possibility of wide use of renewed power sources, division of power and energy transfer services, definition of current and voltage harmonicity distortions sources. At market conditions electric power cost can be different in generating stations and other suppliers. Therefore, nodal pricing is an important problem. Nodal prices include losses in networks. The process of electric power flow tracing allows establishing the routes that connect particular generators with particular consumers, as well as defining shared contributions of generators to cover loading with regard to address distribution. The addressing method also allows defining a power transit price. The rules of joining new consumers to an electric network make it nessessary to carry out a structural analysis of a power system on ultraharmonics. The article investigates address distribution of active power flows and losses between the participants of the electric power market. The authors developed the algorithm and a program to solve the problem of power tracing and power losses. In the developed program Zbus matrix is used to solve the equations of the established modes, as well as to calculate nodal losses. For research of active power address transfer the authors used the IEEE 6bus test scheme, 8bus and 26bus equivalent schemes of Azerbaijan electric-power system. Results of address distribution calculations can be used when awarding two-sided contracts, calculating nodal prices, defining transfer costs etc.
Keywords: power flows, power losses, address distribution, modeling, matrix, marginal method, Zbus method.
References
1. Bialek J. Tracing the flow of electricity. IEEE Proc. Generation, Transmission and Distribution. 1996, vol. 143, pp. 313-320.
2. Kirschen D., Allan R., Strba G. Contributions of individual generators to loads and flows. IEEE Trans. Power Systems. 1997, no. 12, pp. 52-60.
3. Bialek J.W., Kattuman P.A. Proportional sharing assumption in tracing methodology. IEEE Proc. Generation, Transmission and Distribution. 2004, vol. 151, iss. 4, pp. 526-532.
4. Conejo J.J., Arroyo M., Alguacil N., Guijarro A.L. Transmission Loss Allocation: A Comparison of Different Practical Algorithms. IEEE Trans. Power Systems. 2002, vol. 17, no. 3, pp. 571-576.
5. Conejo A.J., Galiana F.D., Kockar I. Z-bus loss allocation. IEEE Trans. pp. 105-110.
6. Gamm A.Z., Golub I.I. Addressness of active and reactive power transfer in an electric-power system. Elektri-chestvo [Electricity]. 2003, no. 3, pp. 16-19 (in Russ.).
7. Bartolomey P.I., Panikovskaya T.Yu., Tikhonov S.A. Analysis and development of power loss responsibility evaluation algorithms of market participants. Energosistema: upravlenie, konkurentsiya, obrazovanie: sbornik dokladov III mezhdunar. nauch.-praktich. konf. [Proc. of the 3rd Int. Science and Practice Conf. "Power System: management, competition, generation"]. Ekaterinburg, 2008, vol. 1, pp. 623-627 (in Russ.).
8. Balametov A.B., Khalilov E.D. On tracking power transfer in power systems. Problemy energetiki [The Problems of Electrical Industry]. 2011, no. 4, pp. 10-17 (in Russ.).
9. Balametov A.B., Khalilov E.D. On losses distribution using Zbus method. Problemy energetiki [The Problems of Electrical Industry]. 2011, no. 2, pp. 30-38 (in Russ.).
10. Balametov A.B., Khalilov E.D. Methods of power loss forecasting in electric-power systems to analyse and control. Elektrichestvo [Electricity]. 2013, no. 7, pp. 19-29 (in Russ.).
Примеры оформления статьи в списке литературы
1. Баламетов А.Б., Халилов Э.Д., Исаева Т.М., Исгендеров Ф.Г. Об адресном распределении мощностей в электрических сетях // Программные продукты и системы. 2016. № 1 (113). С. 126-133.
2. Баламетов А.Б., Халилов Э.Д., Исаева Т.М., Исгендеров Ф.Г. Об адресном распределении мощностей в электрических сетях // Программные продукты и системы. 2016; DOI: 10.15827/0236-235X. 113.126-133.
3. Balametov A.B., Khalilov E.D., Isaeva T.M., Isgenderov F.G. On power tracing distribution in electric networks. Programmnye produkty i sistemy [Software & Systems]. 2016, no. 1, pp. 126-133 (in Russ.); DOI: 10.15827/0236-235X.113.126-133.
