Упрощенная модель определения узловых цен на рынке электроэнергии
Автор: Булатов Борис Георгиевич, Каркунов Вадим Олегович
Журнал: Вестник Южно-Уральского государственного университета. Серия: Энергетика @vestnik-susu-power
Рубрика: Электроэнергетика
Статья в выпуске: 34 (167), 2009 года.
Бесплатный доступ
Разработана математическая модель, в которой учет режимных ограничений осуществлен по методу Лагранжа, а для учета потерь использованы R-схема и коэффициенты токораспределения. Проведены сравнительные расчеты для виртуальной энергосистемы и для первой ценовой зоны, подтвердившие возможность использования метода для целей прогнозирования узловых цен.
Рынок электроэнергии, узловые цены, метод лагранжа, коэффициенты токораспределения
Короткий адрес: https://sciup.org/147158077
IDR: 147158077
Текст краткого сообщения Упрощенная модель определения узловых цен на рынке электроэнергии
С 1 сентября 2006 года постановлением Правительства Российской Федерации введены новые правила функционирования оптового рынка электроэнергии (мощности), которые меняют всю систему взаимоотношений покупателей и поставщиков электрической энергии и мощности, усиливая роль свободных цен. В новой модели оптового рынка предусмотрены два способа торговли электроэнергией - это свободные двусторонние договоры и рынок «на сутки вперед» (РСВ). В рамках свободных двусторонних договоров участники рынка сами определяют контрагентов, цены и объемы поставки. Основой рынка «на сутки вперед» является проводимый администратором торговой системы (АТС) конкурентный отбор ценовых заявок поставщиков и покупателей за сутки до реальной поставки электроэнергии с определением цен и объемов поставки на каждый час суток. Для проведения отбора ценовых заявок системный оператор (СО) предоставляет АТС актуализированную расчетную модель и планируемое почасовое потребление.
Успешная экономическая деятельность участников оптового рынка электроэнергии зависит от объективного учета складывающейся на рынке ценовой ситуации при формировании заявок на покупку/продажу электроэнергии. В данных условиях особую актуальность приобретает задача оценки значений равновесных цен на электроэнер- гию в узлах расчетной модели, определяющих ситуацию на оптовом рынке.
Целью работы является создание практического инструмента для участников ОРЭ по определению узловых цен на рынке «на сутки вперед», который будет являться упрощенной альтернативой расчетной модели, применяемой АТС.
Методика АТС
Конкурентный отбор ценовых заявок на сутки вперед, заключающийся в определении планового почасового производства и потребления электро энергии участниками оптового рынка, проводится администратором торговой сети с использованием расчетной модели, предоставленной СО, в соответствии с Регламентом [1]. Математическая мо дель формируется как оптимизационная задача нелинейного программирования. Целевая функция задачи расчета узловых цен и объемов РСВ задает ся следующим образом:
где Pg(l) - объем в 1-й ступеньке ценовой заявки генератора g; Qg - реактивная мощность генератора g; Рс (1) - объем в 1-й ступеньке ценовой заявки потребителя с; Uj, dj — модуль и фаза напряжения в узле j соответственно.
Сложность задачи определяется большим числом ограничений в форме равенств и неравенств. Рассмотрим лишь основные из них.
Ограничения в форме равенств определяются условиями баланса активной и реактивной мощности в каждом узле]:
ZPu+ZPg-ZP==0 для каждого],
‘ 8 ° (2)
E^j+ZQg-ZQ™""0 Для каждого], i g с где Pg = ^Pg(l)-активная мощность, выдаваемая
1 поставщиком;
РС=^РС(1) - активная мощность,
1 .
принимаемая потребителем.
Учет зависимостей потоков активной ру = f(Ui, Uj, d1, dj) и реактивной мощности Qy = f(Uj, Uj, dj, dj) по ветвям сети от модулей и фаз напряжений требует решения задачи расчета режима энергосистемы на каждой итерации поиска экстремума.
Ограничения в форме неравенств более разнообразны и учитывают:
-
- ограничения по производству и потреблению электроэнергии, исходя из заявок участников рынка
О < Р (1) < Pb,d для всех g, 1;
О < Рс (1) < Pb,d для всех с, 1, где Pbld - активная мощность, соответствующая объему 1-й ступени ценовой заявки поставщика; Pbld - активная мощность, соответствующая объему 1-й ступени ценовой заявки потребителя;
-
- ограничения на производство реактивной мощности
Q7 < Qg < Q^ для всех g;(4)
-
- ограничения потоков активной мощности по контролируемым сечениям или ветвям
p™n ^ Xpu - р^ дая всех с >< ст где с - контролируемые сечения;
-
- ограничения на уровни напряжения
0,9ЦНОМ < Uj < IJU^ для всехj.(6)
Результатом решения оптимизационной задачи являются оптимальные значения Pg(l), Рс(1), определяющие объемы производства и потребления электрической энергии. В ходе решения определяются и двойственные переменные к ограничениям [2]. Множители Лагранжа Х^ к ограничениям (2) по балансу активной мощности в узлах с экономической точки зрения интерпретируются как равновесные узловые цены. Множители Лагранжа ца к активным ограничениям (5) по контролируемым сечениям с экономической точки зрения интерпретируются как цены последнего МВт пропускной способности сечения.
Используя математический пакет MATLAB [3,4], был проведен расчет узловых цен в тестовой сети, результаты которого приведены на рис. 1. В примере все нагрузки приняты ценопринимающими, ценовые заявки генераторов имеющими одну ступень, балансирующий узел б с U=115 кВ и d=0, контролируемые потоки в ветвях 6-5 и 3-4.
Упрощенная методика по R-схеме
Реализация метода АТС требует значительных объемов информации для актуализированной г2=520
S2 = 27,53+jl4,14
Л Г2
113,25
*2С2=520
^27^3+j 14,14
|27,32+jl3,48
▼ 21,51+j6,19
С, = 525
Л4
Л1
сз= 526,25 l3 = 111,51 3

S3=50+j20
8,6+j5^1
лз v i .<'6
Sj = 14,28+J8,75
^0 = 525
Л2
Предельно 93 МВт
П
Л5
s5= 6(^24
21,6+j6,49

^ 110+j42,9
*С6 = 500
5 IL-112,39
С5 = 521,47 ^108,71+j 38,69
предельно tj« ПО МВт
Л6
8^7+j5,42i
С4 = 527,75
S4= 40+16
Рис. 1. Результаты расчета по методике АТС

Гб =500
S6= 110+42,9
модели, поскольку узловые цены определяются на основе уравнений баланса, требующих расчета режима. Природа режимной поправки цены в узлах заключается в различном влиянии последнего МВт узловой мощности на потери в сети, количественной мерой которого является удельный прирост потерь [5, 6]. Для приближенного определения потоков мощности в ветвях, потерь в сети и приростов потерь можно использовать коэффициенты токораспределения, что не требует итерационных расчетов.
Таким образом, задача может решаться в два этапа. На первом определяются равновесные узловые цены и оптимальный план производст-ва/потребления электроэнергии, т.е. величин Pg(l), Рс(1), с учетом баланса в системе, ценовых заявок и ограничений по сечениям, но без учета потерь в сети. На втором этапе по найденным узловым мощностям проводится учет потерь, уточнение значений Pg(l), Рс(1) и корректировка узловых цен.
Итак, на первом этапе целевая функция задачи определяется следующим образом:
£№(1)-Z№,(l)^p™x ,. (7)
с 1 g 1
где Pg(l) - объем в 1-й ступени ценовой заявки генератора g;
Рс(1) - объем в 1-й ступени ценовой заявки потребителя с.
В составе ограничений задачи оптимизации учитываются:
-
- баланс активной мощности в системе без учета потерь
ZPg-ZPe=0> (8)
g с где Pg = ^Pg(l) - активная мощность, выдаваемая 1
поставщиком;
Рс = ^ Рс (1) - активная мощность, принимаемая потребителем;
-
- ограничения по производству и потреблению электроэнергии, исходя из заявок участников рынка:
О < Р (1) < Pb,d для всех g, 1, . (9)
О < Рс (1) < Pbld для всех с, 1.
где Pb,d, Pbld - активная мощность, соответствующая объему 1-й ступени ценовой заявки поставщика/потребителя соответственно;
-
- ограничения перетоков активной мощности по контролируемым сечениям
-
р ™" < ^ р.. < р™ для всех о, (10)
СТ где все р^ определяются с помощью предварительно найденной матрицы коэффициентов потокораспределения [а].
Результатом первого этапа решения являются двойственные переменные X, F, р, по которым и определяются узловые цены без учета потерь [2]:
С = X + р + [а]т т ц .
(И)
Найденные таким образом цены определяют приращение целевой функции при изменении узловой мощности на 1 МВт.
Задачей второго этапа является корректировка узловых цен с учетом потерь в сети. В соответствии с методикой оптимизации режимов поправку на потери можно осуществить по формуле
1-- 5Р,
Эл где - удельные приросты потерь.
Приближенно потери активной мощности в сети
^(Ps2+Qs)Rs OW)^^
Л=-----2
и2
U2
где UH - номинальное напряжение расчетной схемы;
tg(p - средний коэффициент мощности.
Потоки мощности по ветвям Ps по найденным на первом этапе узловым мощностям Р определяются по матричному уравнению
Ps= [а] Р.
Как известно из теории оптимизации режимов, удельные приросты потерь можно приближенно определить по формуле
2 " 22 Ps " Rs ’
UH s
в которой сумма вычисляется по любой группе ветвей, связывающих рассматриваемый узел с балансирующим узлом, и будет пропорциональна потере напряжения АТТ = Ц - ПБУ в группе. Потерю напряжения можно найти путем решения узлового уравнения для R-схемы замещения сети. Тогда выражение (15) примет вид
9Р
U2
где
обращенная матрица узловых проводимостей R-схемы;
-
Р - вектор узловых мощностей.
Для уточнения плана генерации/потребления необходимо потери по ветвям, определенные по выражению (13), учесть в уравнении баланса. При повторном расчете все потери мощности перенесутся в балансирующий узел, что в рассмотренном примере может привести к нарушению ограничения по ЛЭП. Причиной этого является формальное отсутствие связи между потоками мощности по ветвям расчетной схемы и мощностью балансирующего узла, поскольку коэффициенты токораспределения от балансирующего узла равны 0.
Для исключения несоответствия примем в качестве балансирующего фиктивный узел № 7 с ценопринимающей нагрузкой в 0 МВт, включенный в рассечку линии Л7. Фиктивная нагрузка должна быть нулевой, иначе в ходе оптимизации потери мощности будут переноситься на балансирующий узел, а не на все электростанции. Так как Л7 имеет ограничение по пропускной способности, то приблизим фиктивный узел № 7 к нагрузочному узлу № 5 на малое расстояние и с участка 5-7 снимем ограничение по контролируемому сечению Л7. Результаты расчета по упрощенной методике приведены на рис. 2.
Сравнение результатов расчета по рассмотренным методам дается в табл. 1.
Решение рыночной задачи для виртуальной энергосистемы
Проверим работоспособность предложенной методики определения узловых цен для небольшой виртуальной энергосистемы, в составе которой КЭС, ТЭЦ, 7 потребителей и п/ст Север связи с магистральной сетью, выступающая на аукционе в качестве поставщика энергии. Потребители, пи тающиеся от п/ст 1, п/ст 2, п/ст 3, п/ст 5, п/ст 6, рассматриваются как гарантирующие поставщики (энергосбытовые организации), следовательно, объемы покупки электроэнергии определяются ценопринимающими заявками (*). Потребители п/ст 4 и п/ст 7 помимо ценопринимающей ступени заявки имеют дополнительную, отражающую предельную стоимость требуемого объема электрической энергии. Объемы регулируемых договоров и технологический минимум выработки электроэнергии электростанций учитываются в ценопринимающей ступени (*). Ценовые заявки некоторых участников рынка представлены на рис. 3.
Экранная страница программы, реализующей исследуемые методы проведения аукциона АТС, показана на рис. 4. Здесь представлены результаты расчета узловых цен по методу АТС для случая ремонта двух ЛЭП. Следует отметить, что изменение оперативной схемы на экране вызывает соответствующую перестройку расчетной модели, расчет и отображение результатов.
Таблица 1
Методика |
АТС |
R-схема |
Узловые цены, руб./МВт ч |
1-525 2-520 3 - 522,5 4-527,75 5-521,41 6-500 |
1-525 2-520 3 - 526,5 4 - 528,02 5-521,53 6-500 |
План производства, МВт |
1 - 14,28 2-27,53 6-110 |
1 - 15,36 2-26,41 6-110 |

Сравнительные результаты расчета


Рис. 3. Ценовые заявки
Файл Просмотр Настройки Инструменты Справка
Напряжения в у? л ах, кВ
п/ст1 п/ст2 п/стЗ п/ст4 пЛт5 п/стб п/ст7
КЭС-1 ТЭЦ п/ст Север
НОРЭМ
522.36
110 кВ
НОРЭМ—*|
Транзит
220 кВ
Р, МВтч С, руб
Ценопр.

1101®
1101®
110 кВ
10 *в
110 кВ
520.92 «*В
110 кВ яйап» 7
10 кВ
Р, МВтч ' С, руб |
20 Ценопр. j
10 555^
564.6
Р, МВтч С, руб 20 Ценопр.
10 550
io е
*т2
10 кВ
518.51
111.2
111.7
113.9
105.6
105.7
111.7 115
188.1
112.9
229.8
110в
514.03
Р. МВтч ■ С, руб | 45 Ценопр. !
0?
Р, МВтч С5 руб 45 Ценопр.
0 С
р, МВтч с, руб ;
30 Ценопр. :
Расчет
Р, МВтч
о
о
о
о
е
КЭС»1
Р,МВтч; С. руб 45 Ценопр.
Р, МВтч] С, руб 45^Ценопр •
МО кВ

513.86
Ncm3
W€m4
552,33 -- но «в
10 кВ
Ценопр
Рис. 4. Схема виртуальной энергосистемы
Р, МВтч С, руб
45 Ценопр.
о
о
Результаты расчета узловых параметров по методике АТС представлены в табл. 2. Значение целевой функции благосостояния F= 291 320.
Для тех же условий был проведен расчет по методу R-схемы, общие результаты которого представлены в табл. 3. Значение целевой функции F= 291 690.
Анализируя результаты конкурентного отбора ценовых заявок по разным методикам, приходим к выводу, что при отключении ЛЭП изменяется конфигурация сети, в результате чего увеличиваются пути от дешевых источников электроэнергии, что приводит к загрузке более дорогого генерирующего оборудования и к росту значений узловых цен в узлах, находящихся на значительном расстоянии от поставщиков энергии с низким уровнем ценовых заявок. Отключение ЛЭП привело к существенному увеличению значения узловой
Таблица 2
Результаты расчета аукциона электроэнергии по методу АТС
Узел |
п/ст 1 |
п/ст 2 |
п/ст 3 |
п/ст 4 |
п/ст 5 |
с, руб. |
522,36 |
518,51 |
513,86 |
552,33 |
551,49 |
Р, МВт-ч |
45 |
45 |
45 |
20 |
45 |
Р3, МВт-ч |
45 |
45 |
45 |
20,0 |
45 |
Q, Мвар |
18 |
18 |
18 |
12 |
18 |
U, кВ |
111,2 |
111,7 |
113.9 |
105.6 |
105.7 |
Узел |
п/ст 6 |
п/ст 7 |
КЭС-1 |
ТЭЦ |
НОРЭМ |
с, руб |
564,6 |
520,92 |
510 |
540 |
513,03 |
Р, МВт-ч |
45 |
30 |
98.46 |
80,66 |
100 |
Р3, МВт-ч |
45 |
30 |
45, 58.46, 0 |
30, 30, 20,66 |
50, 50 |
Q, Мвар |
18 |
12 |
48 |
46 |
38,48 |
U, кВ |
103 |
111,7 |
115 |
108,1 |
229,8 |
Таблица 3
Результаты расчета аукциона по методу R-схемы
Но самый важный вывод заключается в том, что упрощенная расчетная модель дает практически аналогичные качественно и количественно результаты. Для дополнительной оценки возможностей метода R-схемы проведем расчет часовых узловых цен для одной из ценовых зон ЕЭС.
Расчет узловых цен по зонам свободного перетока ЕЭС России
Рассмотрим первую ценовую зону ЕЭС России в виде совокупности зон свободных перетока (ЗСП), в пределах которых электроэнергия и мощность могут быть замещены энергией и мощностью другого генерирующего оборудования в той же ЗСП, а замена электроэнергией и мощностью генерирующего оборудования, расположенного в иной ЗСП, может быть осуществлена только в пределах ограничений перетока между такими зонами. Изобразим на рис. 5 расчетную модель первой ценовой зоны ЕЭС России в виде графа, состоящего из 22 узлов, и обозначим предельные объемы поставки мощности между зонами свободного перетока (ЗСП).
В ходе решения рыночной задачи определяются равновесные цены ЗСП, хотя на самом деле в узлах, отнесенных к зоне свободного перетока, установится различный уровень узловых цен, обусловленный только потерями мощности в ЗСП, и эти значения будут лежать в некоторой окрестности равновесной цены ЗСП.
Апробируем предложенную методику в условиях ограниченной информации, выполнив прогноз конкурентного отбора ценовых заявок на рынке «на сутки вперед» (РСВ) для первой ценовой зоны ЕЭС России на 5 июня 2009 года при следующих допущениях:
-
• вся нагрузка является ценопринимающей;
-
• не учитывается импорт/экспорт электроэнергии;
-
• связь между ЗСП осуществляется линиями напряжением 500 и 220 кВ;
-
• доля ценопринимающего предложения в ценовых заявках тепловых электростанций составляет 33% и учитывает теплофикационный минимум и приоритетные объемы регулируемых договоров;
-
• для учета потерь внутри зоны свободного перетока увеличивается уровень потребления в ней на 5 %;
-
• готовность генерирующего оборудования к выработке и расходы на собственные нужды электростанций выражается коэффициентом использования Ки.
Для формирования ценовых заявок потребителей используем информацию о прогнозируемом потреблении, публикуемую на официальном сайте системного оператора [7, 8] (рис. 6).

Рис. 5. Расчетная модель с указанием допустимых перетоков мощности

Рис. 6. Суточный график нагрузки на 5.06.2009 г.
Результаты расчета конкурентного отбора ценовых заявок на 12 ч представлены на рис. 7.
Проанализируем сложившуюся ситуацию на рынке «на сутки вперед». ЗСП № 15 и № 28 имеют самые низкие значения равновесных цен вследствие большого объема ценопринимающего предложения от АЭС (15 - Балаковская АЭС, 28 - Кольская АЭС) и ГЭС и низкого уровня потребления внутри данных зон свободного перетока. Наибольший уровень равновесной цены приходится на ЗСП № 20 - 22, так как эти части энергосистемы России являются дефицитными и происходит запирание контролируемых сечений с соседними зонами, обеспечивающими приток более дешевой электроэнергии.
Сравним полученные в ходе пробного решения результаты с отчетом о конкурентном отборе ценовых заявок на РСВ, опубликованном на официальном сайте АТС на 12 ч 5.06.09 г., представив их в виде гистограммы (рис. 8). Разница полученных значений узловых цен по упрощенной методике и методике АТС обусловлена недостаточностью (неточностью) необходимой для решения информации (уровень потребления ЗСП, параметры контролируемых сечений, ценовые заявки) и упрощением модели.
В ходе пробного решения наиболее точные значения равновесных цен были получены для ЗСП, относящихся к объединенной энергосистеме Урала, для которых имелась более точная информация. Но несмотря на отклонения результатов прогноза по упрощенной методике от реальных значений, качественно динамика изменения узловых цен не нарушается, следовательно, данную

Рис. 7. Распределение узловых цен по ЗСП
■ Пробное решение а АТ С

7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
Рис. 8. Результаты прогноза узловых цен по ЗСП
методику можно использовать как практический инструмент прогнозирования узловых цен в отдельных региональных системах.
Результатом решения рыночной задачи являются индексы равновесных цен в узлах системы, по которым мощно судить о складывающейся ситуации на оптовом рынке электроэнергии (ОРЭ) и вырабатывать собственную стратегию поведения на рынке.
Список литературы Упрощенная модель определения узловых цен на рынке электроэнергии
- Регламенты рынка (НОРЭМ). -http://www.np-ats.ru/index.jsp?pid=524.
- Бартоломей, П.И. Оптимизация режимов энергосистем: учебное пособие/П.И. Бартоломей, Т.Ю. Паниковская. -Екатеринбург: УГТУ-УПИ, 2008. -164 с.
- Кетков, Ю.Л. MATLAB 7: программирование, численные методы: учебник/Ю.Л. Кетков, А.Ю. Кетков, М.М. Шулъц. -СПб.: БХВ-Петербург, 2005. -752 с.
- Ануфриев, И.Е. MATLAB 7: учебное пособие/И.Е. Ануфриев, А.Б. Смирнов, Е.Н. Смирнова. -СПб.: БХВ-Петербург, 2005. -1104 с.
- Арзамасцев, Д. А. АСУ и оптимизация режимов энергосистем: учебное пособие/Д.А. Арзамасцев, П.И. Бартоломей, A.M. Холян. -М.: Высшая школа, 1983. -272 с.
- Веников, В.А. Оптимизация режимов электростанций и энергосистем: учебник/В.А. Веников, В.Г. Журавлёв, Т.А. Филиппова. -М.: Энергоатомиздат, 1990. -487 с.
- План генерации и потребления объединенных энергосистем России. -http://www.soрorts.aspx?name=gen_cons_plan
- Результаты торговой сессии: Рынок на сутки вперед: Индексы равновесных цен по ЗСП. -http://www.atsenergo.ru/page/654.