Автоматическое управление выработкой активной мощности ветроэлектростанций в условиях ограничений пропускной способности сети

background image

background image

28

Цыденов

 

Е

.

А

., 

ассистент

 

Отделения

 

электроэнергетики

 

и

 

электротехники

 

ФГАОУ

 

ВО

 «

НИ

 

ТПУ

»

Прохоров

 

А

.

В

., 

к

.

т

.

н

., 

доцент

 

Отделения

 

электроэнергетики

 

и

 

электротехники

 

ФГАОУ

 

ВО

 «

НИ

 

ТПУ

»

Ключевые

 

слова

автоматическое

 

регулирование

 

перетоков

 

активной

 

мощности

понижение

 

размерности

машинное

 

обучение

нейронные

 

сети

коэффициенты

 

долевого

 

участия

 

электростанций

расчет

 

установившегося

 

режима

возобновляемые

 

источники

 

энергии

ветровые

 

электростанции

в

о

з

о

б

н

о

в

л

я

е

м

а

я

 э

н

е

р

г

е

т

и

к

а

возобновляемая энергетика

Автоматическое управление 
выработкой активной мощности 
ветроэлектростанций в условиях 
ограничений пропускной 
способности сети

УДК

 621.316:621.311.24:004.855

В

 

работе

 

анализируются

 

ограничения

 

существующих

 

систем

 

автоматического

 

регулирования

 

перетоков

 

активной

 

мощно

-

сти

которые

 

в

 

большей

 

степени

 

проявляются

 

при

 

подключении

 

к

 

энергосистеме

 

источников

 

энергии

 

с

 

переменным

 

характером

 

работы

таких

 

как

 

солнечные

 

и

 

ветровые

 

электростанции

Для

 

преодоления

 

указанных

 

ограничений

 

авторами

 

предложена

 

методика

позволяющая

 

повысить

 

гибкость

 

и

 

селективность

 

управления

 

перетоками

 

активной

 

мощности

 

посредством

 

он

-

лайн

-

расчета

 

коэффициентов

 

долевого

 

участия

 

регулирующих

 

электростанций

Согласно

 

методике

расчет

 

коэффициентов

 

до

-

левого

 

участия

 

рассматривается

 

как

 

задача

 

оптимизации

реше

-

ние

 

которой

 

в

 

режиме

 

онлайн

 

позволяет

 

распределять

 

задание

 

внеплановой

 

мощности

 

между

 

регулирующими

 

электростанци

-

ями

 

с

 

учетом

 

их

 

характеристик

а

 

также

 

фактической

 

схемно

-

режимной

 

ситуации

 

в

 

энергосистеме

На

 

основе

 

предложенной

 

методики

 

разработан

 

программный

 

модуль

который

 

взаимодей

-

ствует

 

с

 

базой

 

данных

 

реального

 

времени

 

и

 

с

 

компонентом

 

для

 

расчета

 

установившегося

 

режима

что

 

позволяет

 

использовать

 

при

 

расчете

 

параметры

актуализируемые

 

в

 

темпе

 

обновления

 

телеметрии

В

 

работе

 

исследованы

 

различные

 

варианты

 

на

-

стройки

 

разработанного

 

модуля

позволяющие

 

реализовать

 

как

 

стандартные

так

 

и

 

более

 

сложные

 

стратегии

 

управления

необ

-

ходимые

 

для

 

сетей

 

с

 

ограниченной

 

пропускной

 

способностью

Кроме

 

того

авторами

 

предложено

 

использовать

 

модель

 

ма

-

шинного

 

обучения

а

 

именно

 

нейронную

 

сеть

 

с

 

полносвязными

 

слоями

для

 

расчета

 

перетоков

 

активной

 

мощности

Результаты

 

исследования

 

подтверждают

что

 

данное

 

решение

 

позволяет

 

обеспечить

 

повышение

 

скорости

 

расчета

 

коэффициентов

 

доле

-

вого

 

участия

необходимое

 

при

 

привлечении

 

источников

 

энергии

 

с

 

быстроизменяющимся

 

характером

 

работы

 

к

 

участию

 

во

 

вто

-

ричном

 

регулировании

Также

 

в

 

работе

 

рассматривается

 

приме

-

нение

 

алгоритма

 Lasso-

регрессии

 

для

 

сокращения

 

размерности

 

расчетной

 

модели

 

энергосистемы

 

в

 

соответствии

 

с

 

задачей

 

управления

 

и

 

повышения

 

быстродействия

 

расчета

 

перетоков

 

ал

-

горитмами

 

машинного

 

обучения

Предложенные

 

подходы

 

были

 

апробированы

 

в

 

рамках

 

расчетных

 

экспериментов

результаты

 

которых

 

подтвердили

 

преимущества

 

их

 

внедрения

 

в

 

существую

-

щие

 

системы

 

управления

 

перетоками

 

активной

 

мощности

.


background image

29

У

величение

 

доли

 

возобновляемых

 

источни

-

ков

 

энергии

  (

ВИЭ

в

 

составе

 

электроэнер

-

гетических

 

систем

 

приводит

 

к

 

появлению

 

новых

 

технологических

 

задач

связанных

 

со

 

стремлением

 

энергокомпаний

 

обеспечить

 

наиболее

 

эффективное

 

использование

 

их

 

ресурсов

 

в

 

услови

-

ях

 

ограничений

 

пропускной

 

способности

 

существую

-

щей

 

сетевой

 

инфраструктуры

Актуальность

 

проблемы

 

ограничения

 

мощности

 

ВИЭ

 

возрастает

 

с

 

увеличением

 

доли

 

ее

 

установлен

-

ной

 

мощности

 

в

 

энергосистеме

Согласно

 

опублико

-

ванным

 

отчетам

даже

 

при

 

относительно

 

небольшой

 

доле

 

ВИЭ

 

в

 

балансе

 

мощности

 

энергосистем

 

России

уже

 

проявляются

 

проблемы

приводящие

 

к

 

ограни

-

чению

 

выдачи

 

мощности

 

как

 

ветровых

так

 

и

 

солнеч

-

ных

 

электростанций

 [1]. 

К

 

основным

 

причинам

 

ввода

 

ограничений

 

относят

 [2]: 

– 

избыток

 

генерации

 

ВИЭ

 

в

 

часы

 

минимума

 

при

 

на

-

личии

 

тепловой

 

генерации

работающей

 

на

 

техно

-

логическом

 

минимуме

особенно

 

в

 

ночные

 

часы

;

– 

сетевые

 

ограничения

 (

значения

 

напряжения

 

ниже

 

допустимого

 

уровня

токовые

 

перегрузки

и

 

огра

-

ничения

 

по

 

динамической

 

устойчивости

 — 

низкая

 

инерция

).

Автоматизация

 

ограничения

 

выдачи

 

мощности

 

ВИЭ

 

позволяет

 

повысить

 

их

 

коэффициент

 

использо

-

вания

 

установленной

 

мощности

 (

КИУМ

), 

так

 

как

 

ввод

 

ограничений

 

по

 

команде

 

диспетчера

 

в

 

условиях

 

бы

-

стрых

 

и

 

слабо

 

прогнозируемых

 

изменений

 

мощности

 

ВИЭ

 

приводит

 

к

 

излишнему

 

увеличению

 

времени

 

действия

 

ограничений

 

и

 

их

 

избыточности

Поэтому

 

в

 

зарубежных

 

энергосистемах

где

 

доля

 

ВИЭ

 

значи

-

тельна

необходимость

 

принудительного

 

снижения

 

выработки

 

ВИЭ

 

является

 

актуальной

 

задачей

при

 

решении

 

которой

 

стараются

 

минимизировать

 

объ

-

ем

 

невыпускаемой

 

мощности

 

за

 

счет

 

автоматизации

 

оценки

 

схемно

-

режимных

 

условий

 

и

 

реализации

 

не

-

обходимых

 

управляющих

 

воздействий

 [2]. 

Для

 

этого

 

развиваются

 

программно

-

технические

 

комплексы

 

для

 

мониторинга

 

пропускной

 

способности

 

сети

по

-

добные

 

применяемой

 

в

 

России

 

системе

 

мониторинга

 

запасов

 

устойчивости

 [3], 

а

 

также

 

разрабатываются

 

стратегии

 

привлечения

 

ВИЭ

 

к

 

автоматическому

 

регу

-

лированию

 

частоты

 

и

 

перетоков

 

активной

 

мощности

 

(

АРЧМ

) [4, 5].

В

 

результате

помимо

 

роста

 

требований

 

к

 

элек

-

тростанциям

 

на

 

базе

 

ВИЭ

участвующим

 

во

 

вторич

-

ном

 

регулировании

данные

 

тенденции

 

определяют

 

дополнительные

 

требования

 

и

 

к

 

самим

 

системам

 

АРЧМ

такие

 

как

 

необходимость

 

более

 

гибкого

 

и

 

се

-

лективного

 

управления

 

перетоками

 

активной

 

мощ

-

ности

что

 

делает

 

актуальными

 

задачи

 

усовер

-

шенствования

 

существующих

 

систем

 

АРЧМ

 

для

 

их

 

соответствия

 

текущим

 

потребностям

 

отрасли

.

Для

 

повышения

 

качества

 

и

 

эффективности

 

АРЧМ

 

в

 

литературе

 

предлагаются

 

адаптивные

 

регуляторы

 

или

 

алгоритмы

 

адаптивной

 

настройки

 

традиционных

 

ПИД

-

регуляторов

которые

 

позволяют

 

улучшить

 

ди

-

намику

 

и

 

гибкость

 

регулирования

Однако

 

данные

 

подходы

 

либо

 

требуют

 

значительных

 

модификаций

 

существующих

 

систем

 

управления

либо

 

не

 

позволя

-

ют

 

повысить

 

их

 

селективность

 — 

способность

 

пере

-

распределять

 

перетоки

 

между

 

отдельными

 

связями

 

для

 

предотвращения

 

недопустимых

 

перегрузок

В

 

то

 

же

 

время

 

значительно

 

меньше

 

внимания

 

уделяет

-

ся

 

вопросам

 

распределения

 

управляющего

 

воздей

-

ствия

 

от

 

ПИД

-

регулятора

 

между

 

регулирующими

 

электростанциями

хотя

 

автоматизация

 

расчета

 

ко

-

эффициентов

 

долевого

 

участия

  (

КДУ

может

 

в

 

зна

-

чительной

 

степени

 

повысить

 

эффективность

 

управ

-

ления

.

Расчет

 

коэффициентов

 

долевого

 

участия

 — 

это

 

задача

 

оптимизации

где

 

целевая

 

функция

 

опреде

-

ляется

 

в

 

соответствии

 

с

 

целью

 

управления

Если

 

управление

 

осуществляется

 

для

 

регулирования

 

частоты

 

или

 

сальдо

 

перетока

значения

 

КДУ

 

могут

 

определяться

 

онлайн

 

с

 

помощью

 

программного

 

обе

-

спечения

которое

 

анализирует

 

вторичные

 

резервы

 

и

 

маневренность

 

регулирующих

 

электростанций

стоимость

 

регулирования

а

 

также

 

другую

 

телеме

-

трическую

 

и

 

рыночную

 

информацию

 [6–9]. 

В

 

этом

 

случае

 

АРЧМ

 

не

 

учитывает

 

влияние

 

отдельных

 

ре

-

гулирующих

 

электростанций

 

на

 

внутренние

 

связи

 

и

следовательно

задачи

 

селективного

 

управления

 

не

 

рассматриваются

.

С

 

другой

 

стороны

если

 

целью

 

АРЧМ

 

является

 

ре

-

гулирование

 

или

 

ограничение

 

перетоков

 

по

 

внутрен

-

ним

 

связям

  (

далее

 — 

регулирование

 

перетоков

), 

значения

 

КДУ

 

определяются

 

на

 

основе

 

анализа

 

эф

-

фективности

 

электростанций

который

 

заключается

 

в

 

их

 

ранжировании

 

по

 

способности

 

выполнять

 

регу

-

лирование

 

перетока

 

в

 

контролируемом

 

сечении

 

при

 

наименьшем

 

использовании

 

вторичных

 

резервов

 

и

 

обеспечивать

 

минимальное

 

влияние

 

на

 

смежные

 

сечения

 

для

 

предотвращения

 

их

 

перегрузки

 [10–12]. 

Данный

 

анализ

 

выполняется

 

несколько

 

раз

 

в

 

месяц

 

или

 

в

 

год

 

специалистами

 

диспетчерских

 

центров

а

 

его

 

методика

 

имеет

 

следующие

 

ограничения

:

1. 

Значения

 

КДУ

 

обновляются

 

редко

поэтому

 

не

 

от

-

ражают

 

эффективность

 

регулирования

 

электро

-

станций

 

в

 

соответствии

 

с

 

фактической

 

схемно

-

ре

-

жимной

 

ситуацией

 

в

 

энергосистеме

.

2. 

Существующий

 

анализ

 

эффективности

 

оценива

-

ет

 

влияние

 

электростанций

 

на

 

переток

 

при

 

пооче

-

редном

 

изменении

 

их

 

мощности

что

 

не

 

в

 

полной

 

мере

 

соответствует

 

реальному

 

процессу

 

управле

-

ния

где

 

изменение

 

загрузки

 

станций

 

производит

-

ся

 

одновременно

.

3. 

Сложность

 

и

 

глубина

 

анализа

 

при

 

расчете

 

КДУ

 

ограничена

 

возможностями

 

и

 

опытом

 

человека

поэтому

 

более

 

трудоемкие

 

задачи

такие

 

как

 

пре

-

дотвращение

 

перегрузки

 

смежных

 

связей

могут

 

быть

 

решены

 

не

 

полностью

 

или

 

не

 

решены

 

со

-

всем

.

Таким

 

образом

существующая

 

методика

 

расче

-

та

 

КДУ

 

для

 

задач

 

регулирования

 

частоты

 

и

 

сальдо

 

перетока

 

удовлетворяет

 

текущим

 

требованиям

 

в

 

от

-

расли

Однако

 

реализация

 

управления

 

для

 

решения

 

данных

 

задач

 

может

 

вызвать

 

перегрузки

 

внутренних

 

связей

 

энергосистемы

 

и

 

привести

 

к

 

срабатыванию

 

ограничителей

 — 

функций

 

предотвращения

 

пере

-

грузки

 

по

 

току

 

или

 

мощности

управляющие

 

воздей

-

ствия

 

которых

 

имеют

 

более

 

высокий

 

приоритет

чем

 

у

 

регуляторов

и

 

могут

 

быть

 

направлены

 

встречно

 

их

 

 4 (73) 2022


background image

30

ВОЗОБНОВЛЯЕМАЯ

ЭНЕРГЕТИКА

действиям

Чтобы

 

это

 

предотвратить

требуется

 

ана

-

лиз

 

параметров

 

режима

 

энергосистемы

 

при

 

опреде

-

лении

 

КДУ

Помимо

 

этого

недостатки

 

офлайн

-

рас

-

чета

 

КДУ

выполняемого

 

для

 

задач

 

регулирования

 

перетоков

снижают

 

эффективность

 

использования

 

вторичных

 

резервов

 — 

суммарное

 

изменение

 

мощ

-

ности

 

регулирующих

 

электростанций

 

значительно

 

превышает

 

ошибку

 

управления

 

и

 

более

 

сложные

 

стратегии

 

управления

 

остаются

 

нереализованными

Кроме

 

того

низкая

 

периодичность

 

обновления

 

зна

-

чений

 

КДУ

 

препятствует

 

эффективному

 

привлече

-

нию

 

ВИЭ

 

к

 

задачам

 

вторичного

 

регулирования

в

 

то

 

время

 

как

 

некоторые

 

отраслевые

 

стандарты

 

в

 

за

-

рубежных

 

странах

 

уже

 

регламентируют

 

требования

 

по

 

обеспечению

 

их

 

работы

 

под

 

управлением

 

систем

 

АРЧМ

 [5].

В

 

данной

 

статье

 

авторами

 

предложена

 

методи

-

ка

которая

 

объединяет

 

оптимизационный

 

алгоритм

 

с

 

компонентом

 

для

 

расчета

 

установившегося

 

режи

-

ма

 (

УР

и

 

позволяет

 

выполнять

 

оценку

 

КДУ

 

с

 

учетом

 

экономических

технологических

 

и

 

режимных

 

пара

-

метров

Предложенный

 

подход

 

легко

 

встраивается

 

в

 

существующую

 

архитектуру

 

систем

 

АРЧМ

позво

-

ляя

 

повысить

 

гибкость

 

и

 

селективность

 

регулиро

-

вания

Помимо

 

этого

авторы

 

разработали

 

и

 

проте

-

стировали

 

соответствующие

 

целевые

 

функции

 

для

 

реализации

 

как

 

стандартных

так

 

и

 

более

 

сложных

 

стратегий

 

управления

позволяющих

 

селективно

 

ре

-

гулировать

 

перетоки

 

в

 

сетях

 

с

 

ограниченной

 

пропуск

-

ной

 

способностью

.

Наличие

 

значительной

 

доли

 

источников

 

энергии

 

с

 

быстроизменяющимся

 

характером

 

работы

 

накла

-

дывает

 

повышенные

 

требования

 

к

 

скорости

 

расчета

 

коэффициентов

 

долевого

 

участия

В

 

большей

 

сте

-

пени

 

это

 

проявляется

 

при

 

привлечении

 

таких

 

источ

-

ников

 

ко

 

вторичному

 

регулированию

Поэтому

 

для

 

повышения

 

быстродействия

 

расчета

 

КДУ

 

в

 

работе

 

исследуется

 

возможность

 

применения

 

алгоритмов

 

машинного

 

обучения

 (

МО

), 

а

 

именно

 

искусственных

 

нейронных

 

сетей

 (

НС

). 

Данный

 

подход

 

рассматрива

-

ется

 

в

 

качестве

 

альтернативы

 

традиционным

 

алго

-

ритмам

 

расчета

 

УР

способной

 

в

 

значительной

 

сте

-

пени

 

ускорить

 

расчет

 

перетоков

 

активной

 

мощности

В

 

литературе

 

можно

 

найти

 

примеры

 

работ

анали

-

зирующих

 

применение

 

полносвязных

 

НС

 [13], 

свер

-

точных

 

НС

 [14] 

и

 

других

 

алгоритмов

 

МО

 

для

 

расчета

 

перетоков

 

мощности

Однако

 

несмотря

 

на

 

это

при

-

меров

 

использования

 

данных

 

алгоритмов

 

для

 

улуч

-

шения

 

существующих

 

систем

 

управления

 

реального

 

времени

 

не

 

найдено

.

РАЗРАБОТКА

 

МЕТОДИКИ

 

ОНЛАЙН

-

РАСЧЕТА

 

КОЭФФИЦИЕНТОВ

 

ДОЛЕВОГО

 

УЧАСТИЯ

 

ЭЛЕКТРОСТАНЦИЙ

Эффективность

 

систем

 

АРЧМ

 

может

 

быть

 

повы

-

шена

если

 

расчет

 

КДУ

 

будет

 

выполняться

 

онлайн

 

с

 

учетом

 

влияющих

 

факторов

значимых

 

для

 

рас

-

сматриваемой

 

задачи

 

управления

Для

 

этого

 

пред

-

лагается

 

использовать

 

алгоритм

структура

 

которого

 

изображена

 

на

 

рисунке

 1.

Согласно

 

рисунку

процесс

 

онлайн

-

расчета

 

КДУ

 

состоит

 

из

 

нескольких

 

этапов

Первоначально

  (

фи

-

олетовые

 

линии

из

 

базы

 

данных

 

реального

 

време

-

ни

 SCADA-

системы

 

диспетчерского

 

центра

 

посту

-

пает

 

актуальная

 

телеметрическая

 

информация

 

TM

 

в

 

модуль

 

оценки

 

состояния

Далее

 

оцененные

 

зна

-

чения

 

TM

'

 

передаются

 

в

 

модуль

 

для

 

расчета

 

уста

-

новившегося

 

режима

который

 

определяет

 

значения

 

режимных

 

параметров

  (

например

перетоков

 

через

 

контролируемое

 

сечение

 

и

 

передает

 

их

 

в

 

опти

-

мизационную

 

модель

куда

 

также

 

поступает

 

инфор

-

мация

 

об

 

объемах

 

вторичных

 

резервов

 

R

gen

текущих

 

значениях

 

мощности

 

регулирующих

 

электростанций

 

стоимости

 

выработки

 

C

gen

 

и

 

другая

 

оперативная

 

информация

.

Второй

 

этап

 — 

оптимизационный

 

процесс

  (

ро

-

зовые

 

линии

), 

в

 

котором

 

подбираются

 

мощности

 

регулирующих

 

электростанций

 

P

gen

 

таким

 

образом

чтобы

 

минимизировать

 

заранее

 

заданную

 

целевую

 

функцию

Для

 

рассматриваемой

 

оптимизационной

 

задачи

 

авторы

 

предлагают

 

использовать

 

следующий

 

общий

 

вид

 

целевой

 

функции

:

 

f

(

P

gen

) = 

i

w

 

f

i

(

P

gen

). (1)

Предложенная

 

целевая

 

функция

 

представляет

 

собой

 

взвешенную

 

сумму

 

оптимизационных

 

крите

-

риев

 

f

i

(

P

gen

), 

которые

 

определяются

 

в

 

соответствии

 

с

 

задачей

 

управления

 

и

 

могут

 

содержать

 

экономи

-

ческие

 (

например

C

gen

), 

технологические

 (

например

R

gen

и

 

режимные

 

параметры

 (

например

P

f

). 

При

 

этом

 

совместное

 

использование

 

нескольких

 

критериев

 

требует

 

нормализации

 

значений

чтобы

 

диапазоны

 

их

 

изменения

 

были

 

сопоставимы

 

между

 

собой

 [15]. 

Так

 

как

 

режимные

 

параметры

 

зависят

 

от

 

мощ

-

ности

 

регулирующих

 

электростанций

 

P

gen

их

 

необ

-

ходимо

 

пересчитывать

 

с

 

помощью

 

компонента

 

для

 

расчета

 

УР

 

на

 

каждой

 

итерации

 

оптимизационного

 

процесса

 

с

 

учетом

 

ряда

 

ограничений

:

1. 

Напряжения

 

в

 

каждом

 

узле

 

энергосистемы

 

долж

-

ны

 

находиться

 

в

 

допустимых

 

пределах

опреде

-

ляемых

 

устойчивостью

 

нагрузки

 

и

 

разрешенными

 

для

 

оборудования

 

уровнями

 

перенапряжений

.

2. 

Загрузка

 

генераторов

 

должна

 

определяться

 

с

 

уче

-

том

 

ограничений

 

на

 

передачу

 

мощности

 

по

 

усло

-

виям

 

статической

 

апериодической

 

устойчивости

.

3. 

Диапазон

 

изменения

 

мощности

 

участвующих

 

в

 

АРЧМ

 

станций

 

не

 

должен

 

превышать

 

фактиче

-

Рис

. 1. 

Структура

 

алгоритма

 

для

 

онлайн

-

расчета

 

коэффициентов

 

долевого

 

участия

 

электростанций

 

(

БДРВ

 — 

база

 

данных

 

реального

 

времени

ОМ

 — 

оп

-

тимизационная

 

модель

КДУ

 — 

блок

 

для

 

расчета

 

КДУ

ОС

 — 

модуль

 

оценки

 

состояния

УР

 — 

компонент

 

для

 

расчета

 

установившегося

 

режима

БД

 

АРЧМ

 — 

база

 

данных

 

системы

 

АРЧМ

)

TM

TM

'

PF

БД

РВ

ОС

УР

ОМ

КДУ

БД

АРЧМ

P

f

init

P

P

f

P

P

R

C

    ,

P


background image

31

ских

 

резервов

 

на

 

загрузку

 

и

 

на

 

разгрузку

.

По

 

окончании

 

оптимизацион

-

ного

 

процесса

  (

серые

 

линии

 

на

 

рисунке

 1) 

полученные

 

значения

 

мощностей

 

регулирующих

 

элек

-

тростанций

 

преобразуются

 

в

 

ко

-

эффициенты

 

долевого

 

участия

 

PF

Обновленные

 

значения

 

КДУ

 

сохраняются

 

в

 

базе

 

данных

 

АРЧМ

после

 

чего

 

цикл

 

расчета

 

начинает

-

ся

 

заново

 

с

 

использованием

 

акту

-

альных

 

значений

 

параметров

учи

-

тываемых

 

при

 

оптимизации

.

Таким

 

образом

предложенный

 

подход

 

может

 

быть

 

легко

 

интегрирован

 

в

 

существу

-

ющую

 

архитектуру

 SCADA/EMS 

приложений

В

 

то

 

же

 

время

 

оптимизационные

 

критерии

 

и

 

их

 

весовые

 

коэффициенты

 

w

i

 

в

 

функции

 (1) 

могут

 

быть

 

подобра

-

ны

 

с

 

учетом

 

особенностей

 

энергосистемы

 

и

 

опыта

 

ее

 

эксплуатации

.

РАСЧЕТ

 

КОЭФФИЦИЕНТОВ

 

ДОЛЕВОГО

 

УЧАСТИЯ

 

ДЛЯ

 

РАЗЛИЧНЫХ

 

ЗАДАЧ

 

УПРАВЛЕНИЯ

 

И

 

СОСТАВА

 

ГЕНЕРИРУЮЩЕГО

 

ОБОРУДОВАНИЯ

 1. 

Описание

 

исходных

 

данных

На

 

основе

 

предложенной

 

методики

 

с

 

помощью

 

языка

 

программирования

 Python 

был

 

разработан

 

программный

 

модуль

который

 

в

 

качестве

 

алгоритма

 

оптимизации

 

использует

 

метод

 

дифференциальной

 

эволюции

реализованный

 

в

 

библиотеке

 SciPy. 

Раз

-

работка

 

и

 

апробация

 

модуля

 

выполнены

 

на

 

персо

-

нальном

 

компьютере

 

со

 

следующими

 

параметра

-

ми

: Intel(R) Core(TM) i3-7100 3.90 

ГГц

, 8 

Гб

 DDR4 

2133 

МГц

 

ОЗУ

, 125 

Гб

 SSD.

Для

 

демонстрации

 

возможностей

 

разработанно

-

го

 

модуля

 

в

 

данном

 

разделе

 

приведены

 

результаты

 

расчета

 

КДУ

 

для

 

нескольких

 

задач

 

управления

При

 

проведении

 

экспериментов

 

была

 

подготовлена

 

мо

-

дель

 

части

 

энергосистемы

 

России

 (

рисунок

 2) 

с

 

теп

-

ло

и

 

гидроэлектростанциями

 

общей

 

мощностью

 

60 

ГВт

 (464 

узла

 

и

 742 

ветви

).

Чтобы

 

протестировать

 

онлайн

 

работу

 

модуля

авторы

 

использовали

 

архивную

 

телеметрическую

 

информацию

 

о

 

режиме

 

энергосистемы

 

за

 24 

часа

 

с

 

дискретностью

 5 

минут

При

 

этом

 

топология

 

сети

 

была

 

принята

 

постоянной

 

для

 

упрощения

 

расчетов

.

2. 

Результаты

 

расчетов

 

для

 

энергосистемы

 

с

 

тепло

и

 

гидрогенерацией

В

 

рассматриваемой

 

энергосистеме

 

АРЧМ

 

управ

-

ляет

 

двумя

 

гидроэлектростанциями

  (

ГЭС

1

 

и

 

ГЭС

2

для

 

регулирования

 

перетока

 

в

 

сечении

 F

1

Поэтому

 

разработанный

 

модуль

 

был

 

первоначально

 

про

-

тестирован

 

именно

 

для

 

этой

 

стандартной

 

задачи

 

управления

Для

 

расчета

 

установившегося

 

режима

 

был

 

выбран

 

метод

 

Ньютона

-

Рафсона

 (

Н

-

Р

), 

реализо

-

ванный

 

в

 

программном

 

комплексе

 RastrWin. 

Так

 

как

 

вторичные

 

резервы

 

гидроэлектростанций

 

в

 

большей

 

степени

 

определяются

 

правилами

 

водопользования

рельефом

 

и

 

сезонными

 

условиям

 — 

они

 

изменяют

-

ся

 

незначительно

 

в

 

течение

 

суток

 

и

 

могут

 

считаться

 

ВЭС

3

ВЭС

2

F

1

F

2

ГЭС

2

ГЭС

1

220/500

70 

ЭСт

База

ВЭС

1

постоянными

В

 

рассматриваемом

 

случае

 

авторы

 

приняли

 

объемы

 

резервов

 

на

 

загрузку

 

и

 

на

 

разгрузку

 

равными

 50 

МВт

 

для

 

каждой

 

ГЭС

.

Чтобы

 

настроить

 

модуль

алгоритм

 

работы

 

кото

-

рого

 

изображен

 

на

 

рисунке

 1, 

под

 

рассматриваемую

 

задачу

 

управления

 

перетоком

 

в

 

контролируемом

 

сечении

  F

1

предлагается

 

из

 

общего

 

вида

 

целевой

 

функции

 (1) 

сформулировать

 

следующую

:

 

f

(

P

gen

) = 

w

i

 

 f

i

(

P

gen

), (2)

 

f

1

(

P

gen

) = 

P

F1

 – (

 – 

)

, (3)

 

f

2

(

P

gen

) = 



k

R

 

˚

 

(

P

gen

 – 

).

 

n



1

. (4)

В

 

ходе

 

оптимизационного

 

процесса

 

алгоритм

 

стремится

 

так

 

подобрать

 

значения

 

мощностей

 

ре

-

гулирующих

 

электростанций

 

P

gen

чтобы

 

отклонение

 

перетока

 

в

 

контролируемом

 

сечении

 

P

F1

 

от

 

целево

-

го

 

значения

 

было

 

минимальным

Несмотря

 

на

 

то

что

 

КДУ

 

обновляются

 

онлайн

в

 

процессе

 

их

 

рас

-

чета

 

не

 

учитывается

 

фактическая

 

ошибка

 

управле

-

ния

Поэтому

 

для

 

выполнения

 

расчетов

 

необходи

-

мо

 

сформировать

 

эталонную

 

ошибку

 

управления

 

в

 

соответствии

 

с

 

рассматриваемой

 

задачей

 

управ

-

ления

Для

 

всех

 

экспериментов

представленных

 

в

 

работе

эталонная

 

ошибка

 

 — 

это

 

отклонение

 

перетока

 

в

 

контролируемом

 

объекте

 (

в

 

данном

 

слу

-

чае

в

 

сечении

  F

1

от

 

его

 

исходного

 

значения

 

Значение

 

 

принято

 

равным

 10 

МВт

 

и

 

при

 

экс

-

плуатации

 

может

 

быть

 

подобрано

 

в

 

соответствии

 

с

 

опытом

 

управления

 

энергосистемой

Второй

 

компонент

 (4) 

учитывает

 

резервы

 

электростанций

 

с

 

помощью

 

пропорциональных

 

коэффициентов

 

k

R

 

и

 

заставляет

 

алгоритм

 

оптимизации

 

миними

-

зировать

 

отклонение

 

мощностей

 

регулирующих

 

электростанций

 

от

 

их

 

начальных

 

значений

 

 

и

соответственно

использовать

 

при

 

регулирова

-

нии

 

только

 

наиболее

 

эффективные

 

электростан

-

ции

Величина

 

n

 

позволяет

 

распределять

 

команды

 

АРЧМ

 

так

чтобы

 

более

 

равномерно

 

расходовать

 

вторичные

 

резервы

 

мощности

Весовые

 

коэффи

-

циенты

 

w

1

 

и

 

w

2

 

определяют

 

важность

 

отдельных

 

критериев

 

и

 

могут

 

быть

 

настроены

 

при

 

необходи

-

мости

.

На

 

рисунке

 3 

отображены

 

результаты

 

расчета

 

КДУ

 

для

 

ГЭС

1

 

и

 

ГЭС

2

а

 

также

 

значения

 

КДУ

 

из

 

су

-

ществующей

 

системы

 

АРЧМ

 

в

 

рассматриваемый

 

промежуток

 

времени

Необходимо

 

отметить

что

 

на

 

Рис

. 2. 

Упрощенная

 

схема

 

тестовой

 

энергосистемы

 (

ЭСт

 — 

электростан

-

ция

ВЭС

 — 

ветроэлектростанция

ГЭС

 — 

гидроэлектростанция

)

 4 (73) 2022


background image

32

ВОЗОБНОВЛЯЕМАЯ

ЭНЕРГЕТИКА

рисунке

 3 

и

 

далее

 

значения

 

КДУ

 

приведены

 

с

 

учетом

 

влияния

 

регулирующих

 

электростанций

 

на

 

переток

 

в

 

контролируемом

 

сечении

 

или

 

связи

  (

знак

 

минуса

 

соответствует

 

действию

 

станции

 

на

 

разгрузку

а

 

знак

 

плюса

 — 

на

 

загрузку

).

Наибольшее

 

влияние

 

на

 

эффективность

 

регули

-

рующих

 

электростанций

 

оказывает

 

их

 

электрическая

 

удаленность

 

от

 

контролируемого

 

сечения

которая

 

определяется

 

топологией

 

сети

Так

 

как

 

в

 

данной

 

ра

-

боте

 

используется

 

постоянная

 

топология

значения

 

КДУ

 

на

 

рисунке

 3 

почти

 

не

 

изменяются

Также

 

по

 

рисунку

 

видно

что

 

рассчитанные

 

КДУ

 

схожи

 

с

 

уста

-

новленными

 

значениями

 

в

 

существующей

 

системе

 

АРЧМ

Небольшое

 

отличие

 

можно

 

объяснить

 

разни

-

цей

 

в

 

методике

 

расчета

 

и

 

анализируемых

 

схемно

-

ре

-

жимных

 

ситуациях

.

В

 

данном

 

эксперименте

 

каждый

 

оптимизацион

-

ный

 

цикл

выполняющий

 

расчет

 

КДУ

 

для

 5-

минутно

-

го

 

интервала

 

времени

в

 

течение

 

которого

 

значения

 

параметров

 

режима

 

принимаются

 

установившимися

выполняет

 

в

 

среднем

 1651 

расчет

 

УР

 

и

 

длится

 

около

 

23,1 

секунды

 

при

 

запуске

 

алгоритма

 

в

 

однопоточном

 

режиме

Также

 

стоит

 

отметить

что

 

время

 

расчета

 

од

-

ного

 

установившегося

 

режима

 

для

 

рассматриваемой

 

энергосистемы

 

составляет

 12 

мс

 

при

 

использовании

 

программного

 

комплекса

 

и

 

персонального

 

компьюте

-

ра

описанных

 

выше

.

3. 

Результаты

 

расчета

 

для

 

энергосистемы

 

с

 

тепло

-, 

гидро

и

 

ветрогенерацией

В

 

данном

 

разделе

 

представлены

 

результаты

 

апробации

 

модуля

 

при

 

наличии

 

ветровых

 

электро

-

станций

на

 

которых

 

размещены

 

вторичные

 

резер

-

вы

 

мощности

Для

 

реализации

 

описанной

 

ситуации

 

часть

 

существующей

 

генерации

 

была

 

замещена

 

вет

-

ропарками

 (

ВЭС

1

ВЭС

2

 

и

 

ВЭС

3

мощностью

 600, 600 

и

 300 

МВт

 

соответственно

 

и

 

подключена

 

к

 

системе

 

АРЧМ

.

В

 

левой

 

части

 

рисунка

 4 

отображена

 

выработка

 

активной

 

мощности

 

ВЭС

 

в

 

течение

 

суток

представ

-

ленная

 

в

 

процентах

 

от

 

установленной

 

мощности

 

и

 

принятая

 

для

 

иллюстрации

 

работы

 

модуля

 

соглас

-

но

 

данным

 

Национальной

 

лаборатории

 

по

 

изуче

-

нию

 

возобновляемой

 

энергии

 [16]. 

В

 

свою

 

оче

-

редь

величина

 

вторичного

 

резерва

 

ВЭС

 

задана

 

в

 

зависимости

 

от

 

располагаемой

 

мощности

Это

 

соответствует

 

методикам

описанным

 

в

 

работах

[2, 17], 

согласно

 

которым

 

резерв

 

на

 

загрузку

 

форми

-

руется

 

путем

 

принудительного

 

снижения

 

фактиче

-

ски

 

выдаваемой

 

мощности

 

на

 

некоторую

 

величину

 

относительно

 

располагаемой

в

 

то

 

время

 

как

 

резерв

 

на

 

разгрузку

 

может

 

задаваться

 

от

 

фактической

 

вы

-

даваемой

 

мощности

 

до

 

технологического

 

миниму

-

ма

 

станции

.

В

 

данной

 

работе

 

резервы

 

ВЭС

 

для

 

регулирова

-

ния

 

на

 

разгрузку

 

приняты

 

равными

 10% 

от

 

их

 

распо

-

лагаемой

 

мощности

 

в

 

каждый

 

момент

 

времени

Это

 

допущение

 

принято

так

 

как

 

практика

 

подключения

 

ВЭС

 

к

 

АРЧМ

 

в

 

данный

 

момент

 

еще

 

не

 

сформирова

-

лась

и

 

устоявшихся

 

нормативов

 

для

 

определения

 

вторичных

 

резервов

 

ВЭС

 

нет

.

В

 

данном

 

разделе

как

 

и

 

в

 

предыдущем

 

разде

-

ле

рассматривается

 

задача

 

определения

 

КДУ

 

для

 

управления

 

перетоком

 

в

 

сечении

  F

1

Результаты

 

расчета

представленные

 

в

 

правой

 

части

 

рисунка

 4, 

показывают

 

изменение

 

значений

 

КДУ

 

всех

 

регули

-

рующих

 

электростанций

 

вследствие

 

непостоянства

 

выработки

 

ВЭС

 

и

 

их

 

резервов

Это

 

доказывает

что

 

в

 

случае

 

использования

 

ВЭС

 

во

 

вторичном

 

регули

-

ровании

 

требуется

 

онлайн

-

расчет

 

КДУ

В

 

данных

 

экспериментах

 

оптимизационный

 

цикл

 

выполняет

 

в

 

среднем

 17 541 

расчет

 

УР

 

и

 

длится

 251,04 

се

-

кунды

.

4. 

Регулирование

 

перетока

 

с

 

учетом

 

загрузки

 

смежных

 

связей

Реализация

 

более

 

сложных

 

стратегий

 

управ

-

ления

 

требует

 

расширения

 

целевой

 

функции

 (2) 

дополнительными

 

критериям

 

оптимизации

На

-

пример

для

 

повышения

 

селективности

 

управле

-

ния

 

необходимо

 

учитывать

 

перетоки

 

мощности

 

Рис

. 3. 

Сравнение

 

рассчитанных

 

КДУ

 

со

 

значениями

 

в

 

существующей

 

системе

 

АРЧМ

Рис

. 4. 

Изменение

 

выработки

 

активной

 

мощности

 

ВЭС

 

в

 

течение

 

рассматриваемых

 

суток

 (

слева

), 

рассчитанные

 

значения

 

КДУ

 

регулирующих

 

ГЭС

 

и

 

ВЭС

 (

справа

)

–40

–45

–50

–55

–60

t

ч

КДУ

, %

0 3 6 9 12 15 

18 21 

24

ГЭС

1

ГЭС

2

ГЭС

АРЧМ

ГЭС

АРЧМ

0 3 6 9 12 15 18 21 24

t

ч

0 3 6 9 12 15 18 21 24

t

ч

100

75

50

25

0

0

–15

–30

–45

–60

P

, %

КДУ

, %

ГЭС

1

ВЭС

1

ВЭС

2

ВЭС

3

ГЭС

2


background image

33

в

 

смежных

 

сечениях

 

с

 

ограниченной

 

пропускной

 

способностью

чтобы

 

не

 

допустить

 

их

 

перегрузки

 

в

 

результате

 

регулирования

 

в

 

заданном

 

сечении

 

или

 

области

 

регулирования

Для

 

реализации

 

такого

 

способа

 

управления

 

предлагается

 

добавить

 

третий

 

компонент

 

f

3

(

P

gen

с

 

весовым

 

коэффициентом

 

w

3

 

в

 

це

-

левую

 

функцию

 (2):

 

f

(

P

gen

) = 

w

i

 

 

f

i

(

P

gen

), (5)

 

f

1

(

P

gen

) = 

P

F1

 – (

 – 

)

, (6)

 

f

2

(

P

gen

) = 



k

R

 

˚

 (

P

gen

 – 

).

 

n



1

, (7)

 

f

3

(

P

gen

) = 

P

F2

 – 

. (8)

Добавленный

 

компонент

 (8) 

определяет

 

штраф

 

для

 

алгоритма

 

оптимизации

 

за

 

отклонение

 

перето

-

ка

 

в

 

смежном

 

сечении

 

P

F2

 

от

 

его

 

исходного

 

значения

 

Таким

 

образом

подбирая

 

коэффициент

 

w

3

мож

-

но

 

добиться

 

меньшего

 

воздействия

 

на

 

смежное

 

се

-

чение

 F

2

вызванное

 

регулированием

 

перетока

 

в

 

кон

-

тролируемом

 

сечении

 F

1

.

В

 

данном

 

разделе

 

приведены

 

результаты

 

расчета

 

КДУ

 

с

 

тремя

 

различными

 

значениями

 

w

3

 — 

w

3

А

w

3B

w

3

С

для

 

которых

 

выполняется

 

условие

:

 0 

 

w

3

А

 

 

w

3B

 

 

w

3

С

. (9)

Результаты

 

расчетов

представленные

 

на

 

рисун

-

ке

 5, 

показывают

 

изменение

 

значений

 

КДУ

когда

 

ал

-

горитм

 

стремится

 

управлять

 

перетоком

 

в

 

контроли

-

руемом

 

сечении

 F

1

 

и

 

при

 

этом

 

минимизировать

 

вли

-

яние

 

на

 

смежное

 

сечение

 F

2

Электростанция

 

ВЭС

1

 

обладает

 

большим

 

влиянием

 

на

 

сечение

 F

2

поэтому

 

с

 

увеличением

 

w

3

 

значение

 

ее

 

КДУ

 

значительно

 

воз

-

растает

а

 

когда

 

объем

 

резервов

 

ВЭС

1

 

снижен

  (

ин

-

тервал

 

с

 5 

до

 12 

часов

 

на

 

рисунках

 4 

и

 5), 

переток

 

в

 

смежном

 

сечении

 

увеличивается

С

 

увеличением

 

коэффициента

 

w

3

 

влияние

 

на

 

переток

 

в

 

смежном

 

сечении

  F

2

 

снижается

что

 

по

-

зволяет

 

не

 

допустить

 

его

 

перегрузки

 

в

 

результате

 

регулирования

 

в

 

контролируемом

 

сечении

 F

1

Одна

-

ко

 

данная

 

стратегия

 

регулирования

 

приводит

 

к

 

до

-

полнительному

 

расходу

 

вторичных

 

резервов

 

на

 

раз

-

грузку

 — 

до

 3% 

в

 

экспериментах

представленных

 

на

 

рисунке

 5.

Таким

 

образом

выбор

 

оптимального

 

значения

 

коэффициента

 

w

3

с

 

одной

 

стороны

зависит

 

от

 

тре

-

бований

 

к

 

эффективному

 

использованию

 

вторичных

 

резервов

что

 

соответствует

 

стремлению

 

миними

-

зировать

 

ограничение

 

мощности

 

ВЭС

с

 

другой

 

сто

-

роны

 — 

от

 

требований

 

к

 

минимизации

 

влияния

 

на

 

переток

 

в

 

смежном

 

сечении

 

при

 

регулировании

 

в

 

ос

-

новном

 

сечении

.

Процесс

 

определения

 

КДУ

 

для

 

описанной

 

выше

 

задачи

 

в

 

среднем

 

занимает

 253,33 

секунды

 

и

 

требу

-

ет

 18 484 

расчета

 

установившегося

 

режима

.

Рис

. 5. 

Значения

 

КДУ

 

регулирующих

 

электростанций

отклонение

 

перетока

 

в

 

контролируемом

 (

P

F1

и

 

смежном

 

сечениях

 (

P

F2

), 

суммарное

 

изменение

 

мощности

 

регулирующих

 

электростанций

 (

P

gen

при

 

разных

 

значениях

 

w

3

 

(A — 0, B — w

3

B

, C — 

w

3

C

)

0

–25

–50

–75

–100

10

5

0

–5

–10

КДУ

, %

P

МВт

А

-1

А

-2

0

–25

–50

–75

–100

10

5

0

–5

–10

В

-1

В

-2

0 3 6 9 12 15 18 21 24

t

ч

0 3 6 9 12 15 18 

21 24

t

ч

0

–25

–50

–75

–100

10

5

0

–5

–10

С

-1

С

-2

ГЭС

1

ВЭС

1

ВЭС

2

ВЭС

3

ГЭС

2

P

F2

P

F1

P

gen

 4 (73) 2022


background image

34

ВОЗОБНОВЛЯЕМАЯ

ЭНЕРГЕТИКА

ПРИМЕНЕНИЕ

 

МАШИННОГО

 

ОБУЧЕНИЯ

 

ДЛЯ

 

ПОВЫШЕНИЯ

 

БЫСТРОДЕЙСТВИЯ

 

РАСЧЕТОВ

1. 

Описание

 

подхода

Представленные

 

ранее

 

результаты

 

экспери

-

ментов

 

доказывают

 

применимость

 

предложенной

 

методики

 

расчета

 

КДУ

 

для

 

различных

 

стратегий

 

регулирования

Однако

 

процесс

 

настройки

 

программ

-

ного

 

модуля

 

и

 

его

 

быстродействие

 

зависят

 

от

 

задачи

 

управления

 

и

 

структуры

 

энергосистемы

Первый

 

экс

-

перимент

 

с

 

двумя

 

регулирующими

 

электростанция

-

ми

 

показал

 

приемлемую

 

скорость

 

расчета

 

КДУ

но

 

с

 

увеличением

 

количества

 

электростанций

управля

-

емых

 

АРЧМ

процесс

 

расчета

 

стал

 

значительно

 

мед

-

леннее

Низкое

 

быстродействие

 

модуля

 

может

 

быть

 

существенным

 

недостатком

особенно

 

в

 

энергосис

-

темах

где

 

вторичные

 

резервы

 

предоставляют

 

ВИЭ

.

Требования

 

к

 

быстродействию

 

алгоритма

 

расчета

 

КДУ

 

определяются

 

по

 

большей

 

части

 

скоростью

 

из

-

менения

 

факторов

влияющих

 

на

 

результат

 

расчета

В

 

описанных

 

выше

 

экспериментах

 

используется

 

ин

-

формация

 

о

 

вторичных

 

резервах

 

и

 

эффективности

 

регулирующих

 

электростанций

В

 

этом

 

случае

 

для

 

задачи

 

регулирования

 

перетоков

 

в

 

энергосистемах

 

с

 

традиционным

 

составом

 

генерации

 

целесообразно

 

обновлять

 

значения

 

КДУ

 

по

 

факту

 

изменения

 

тополо

-

гии

так

 

как

 

именно

 

данный

 

фактор

 

в

 

основном

 

опре

-

деляет

 

эффективность

 

электростанций

  (

в

 

то

 

вре

-

мя

 

как

 

вторичные

 

резервы

 

изменяются

 

медленно

). 

В

 

случае

 

участия

 

ВЭС

 

во

 

вторичном

 

регулировании

 

важным

 

становится

 

обновление

 

КДУ

 

в

 

соответствии

 

с

 

изменением

 

доступной

 

мощности

 

ВЭС

которая

 

может

 

быть

 

измерена

 

или

 

приниматься

 

по

 

данным

 

прогноза

Если

 

учесть

что

 

точность

 

информации

 

о

 

прогнозируемой

 

мощности

 

ВЭС

 

и

соответствен

-

но

доступном

 

резерве

 

на

 

ее

 

разгрузку

 

увеличивает

-

ся

 

с

 

сокращением

 

периода

 

прогноза

то

 

в

 

качестве

 

предельного

 

времени

 

расчета

 

КДУ

 

можно

 

принять

 

часто

 

встречающийся

 

в

 

зарубежных

 

энергосистемах

 

интервал

 

прогноза

 — 5 

минут

который

 

обеспечивает

 

приемлемую

 

точность

 

прогноза

 

для

 

задач

 

регулиро

-

вания

 [18]. 

Однако

 

при

 

использовании

 

фактических

 

значений

 

мощности

 

ВЭС

 

расчет

 

КДУ

 

должен

 

выпол

-

няться

 

максимально

 

быстро

 

для

 

поддержания

 

наи

-

более

 

высокого

 

значения

 

КИУМ

 

путем

 

сокращения

 

длительности

 

ограничения

 

выработки

 

ВЭС

.

Одним

 

из

 

способов

 

повышения

 

быстродействия

 

разработанного

 

модуля

 

является

 

снижение

 

длитель

-

ности

 

каждого

 

шага

 

оптимизации

 

путем

 

ускорения

 

расчета

 

перетоков

 

активной

 

мощности

С

 

этой

 

целью

 

в

 

данной

 

работе

 

предлагается

 

вместо

 

традиционно

-

го

 

алгоритма

 

для

 

расчета

 

УР

 (

метода

 

Ньютона

-

Раф

-

сона

использовать

 

методы

 

машинного

 

обучения

В

 

отличие

 

от

 

метода

 

Н

-

Р

алгоритмы

 

МО

 

не

 

требуют

 

расчета

 

всех

 

режимных

 

параметров

 

и

 

могут

 

исполь

-

зоваться

 

для

 

оценки

 

только

 

тех

 

величин

которые

 

являются

 

значимыми

 

для

 

конкретной

 

задачи

напри

-

мер

перетоки

 

активной

 

мощности

 

в

 

контролируемом

 

и

 

смежном

 

сечениях

 

для

 

расчета

 

КДУ

.

Наиболее

 

распространенный

 

способ

 

описания

 

со

-

стояния

 

энергосистемы

 

с

 

постоянной

 

топологией

 — 

использование

 

значений

 

активных

 

и

 

реактивных

 

инъ

-

екций

 

мощности

 

в

 

каждом

 

узле

 

в

 

качестве

 

вектора

 

признаков

Однако

 

в

 

энергосистемах

 

могут

 

быть

 

тысячи

 

узлов

и

 

очевидно

что

 

инъекции

 

только

 

в

 

части

 

из

 

этих

 

узлов

 

могут

 

оказывать

 

влияние

 

на

 

контролируемые

 

па

-

раметры

Наличие

 

излишней

 

информации

 

в

 

векторе

 

признаков

 

может

 

привести

 

к

 

повышению

 

длительности

 

расчетов

 

и

 

снижению

 

точности

 

модели

 

МО

.

Для

 

снижения

 

размерности

 

признакового

 

про

-

странства

 

в

 

данной

 

работе

 

предлагается

 

использо

-

вать

 

модель

 Lasso-

регрессии

которая

 

представляет

 

собой

 

модель

 

линейной

 

регрессии

обученную

 

с

  L

1

-

регуляризацией

С

 

помощью

 

данного

 

подхода

 

процесс

 

понижения

 

размерности

 

выглядит

 

следующим

 

образом

:

1. 

Задается

 

коэффициент

 

регуляризации

.

2. 

Формируется

 

модель

 Lasso-

регрессии

 

с

 

исполь

-

зованием

 

всей

 

обучающей

 

выборки

.

3. 

Признаки

чей

 

вес

 

в

 

модели

 Lasso-

регрессии

 

больше

 

нуля

принимаются

 

как

 

значимые

 

для

 

рассматриваемой

 

задачи

 

и

 

используются

 

в

 

каче

-

стве

 

входных

 

параметров

 

модели

 

МО

 

для

 

расчета

 

перетоков

.

Применение

 Lasso-

регрессии

 

схоже

 

с

 

выпол

-

нением

 

анализа

 

эффективности

 

регулирующих

 

электростанций

в

 

рамках

 

которого

 

специалисты

 

из

-

меняют

 

значения

 

инъекций

 

мощности

 

в

 

генерирую

-

щих

 

узлах

 

и

 

отслеживают

 

реакцию

 

контролируемых

 

параметров

 

на

 

данные

 

изменения

Таким

 

образом

предложенный

 

подход

 

для

 

понижения

 

размерности

 

позволяет

 

выбрать

 

только

 

те

 

признаки

которые

 

ока

-

зывают

 

большее

 

влияние

 

на

 

контролируемые

 

пере

-

токи

 

и

 

другие

 

параметры

рассчитываемые

 

с

 

помо

-

щью

 

модели

 

МО

.

Предложенный

 

подход

 

для

 

расчета

 

перетоков

 

мощности

 

с

 

помощью

 

модели

 

машинного

 

обучения

 

схематично

 

изображен

 

на

 

рисунке

 6. 

Так

 

как

 

для

 

Рис

. 6. 

Предложенный

 

подход

 

для

 

расчета

 

перетоков

 

ак

-

тивной

 

мощности

 

с

 

помощью

 

модели

 

машинного

 

обуче

-

ния

 (

ОП

 — 

компонент

 

для

 

отбора

 

признаков

ФП

 — 

ком

-

понент

 

для

 

фильтрации

 

признаков

ФМ

 — 

компонент

 

для

 

формирования

 

моделей

 

машинного

 

обучения

МП

 — 

модель

 

машинного

 

обучения

 

для

 

расчета

 

перетоков

ОС

 

и

 

УР

 — 

модуль

 

для

 

оценки

 

состояния

 

и

 

компонент

 

для

 

расчета

 

установившегося

 

режима

как

 

на

 

рисунке

 1)

TM

'

ОП

ФП

ФП

ОС

ФМ

МП

P

fi

 lt

inj

P

fi

 lt

inj

Q

fi

 lt

inj

Q

fi

 lt

inj

P

inj, 

Q

inj

L

1

FL

МП

УР

Офлайн

Онлайн

P

inj, 

Q

inj

P

f

P

f

P

f

P

f

Об

учающая

 

выборк

а


background image

35

формирования

 

модели

 

используется

 

подход

 

обуче

-

ния

 

с

 

учителем

предварительно

 

должна

 

быть

 

подго

-

товлена

 

обучающая

 

выборка

 

с

 

известными

 

значени

-

ями

 

признаков

 (

P

inj

Q

inj

и

 

целевых

 

переменных

 (

P

f

). 

Первоначально

 (

зеленые

 

линии

обучающая

 

выбор

-

ка

 

передается

 

в

 

компонент

 

для

 

отбора

 

признаков

который

 

формирует

 

список

 

наиболее

 

значимых

 

для

 

рассматриваемой

 

задачи

 

параметров

  (

FL

). 

Далее

 

(

коричневые

 

линии

из

 

обучающей

 

выборки

 

выде

-

ляются

 

только

 

значимые

 

признаки

 (

), 

которые

 

впоследствии

 

передаются

 

в

 

компонент

 

формирова

-

ния

 

модели

 

для

 

расчета

 

перетоков

 

активной

 

мощно

-

сти

Данные

 

шаги

 

могут

 

быть

 

выполнены

 

несколько

 

раз

  (

голубые

 

линии

для

 

подбора

 

гиперпараметров

 

(

включая

 

коэффициент

 

регуляризации

 

на

 

этапе

 

от

-

бора

 

признаков

), 

позволяющих

 

достичь

 

необходи

-

мой

 

точности

 

модели

После

 

обучения

 

и

 

тестирова

-

ния

 

наилучшая

 

модель

 

и

 

соответствующий

 

перечень

 

значимых

 

признаков

 

используются

 

вместе

 

в

 

онлайн

-

режиме

 

как

 

компонент

 

для

 

расчета

 

установившегося

 

режима

 

из

 

рисунка

 1, 

позволяющий

 

определять

 

не

-

обходимые

 

для

 

оптимизации

 

параметры

 

режима

 

P

f

.

2. 

Применение

 

машинного

 

обучения

 

для

 

расчета

 

перетоков

 

активной

 

мощности

 

в

 

задаче

 

опре

-

деления

 

коэффициентов

 

долевого

 

участия

В

 

данном

 

разделе

 

рассматривается

 

та

 

же

 

задача

 

управления

что

 

и

 

в

 

разделе

 «

Результаты

 

расчета

 

для

 

энергосистемы

 

с

 

тепло

-, 

гидро

и

 

ветрогенерацией

», 

однако

 

алгоритм

 

на

 

рисунке

 1 

был

 

модифицирован

 

путем

 

замены

 

метода

 

Н

-

Р

 

на

 

модель

 

машинного

 

обу

-

чения

 

для

 

выполнения

 

расчета

 

перетоков

 

активной

 

мощности

Так

 

как

 

поиск

 

наилучшей

 

модели

 

МО

 

для

 

расчета

 

перетоков

 — 

это

 

задача

требующая

 

отдель

-

ного

 

исследования

в

 

данной

 

работе

 

для

 

упрощения

 

предлагается

 

использовать

 

модель

 

искусственной

 

нейронной

 

сети

 

с

 

полносвязными

 

слоями

 

как

 

алго

-

ритм

способный

 

аппроксимировать

 

любую

 

функ

-

цию

 [19].

Для

 

рассматриваемой

 

задачи

 

первоначально

 

был

 

составлен

 

вектор

 

признаков

который

 

содержит

 

инъ

-

екции

 

всех

 

узлов

кроме

 

балансирующего

 — 

всего

 

578 

значений

при

 

этом

 

рассчитываемым

 

значением

 

является

 

переток

 

в

 

контролируемом

 

сечении

 F

1

Обу

-

чающая

 

выборка

 

была

 

сформирована

 

путем

 

расчета

 

режима

 

методом

 

Ньютона

-

Рафсона

 

при

 

различных

 

комбинациях

 

значений

 

вектора

 

признаков

Диапазон

 

изменения

 

параметров

 

был

 

определен

 

с

 

помощью

 

архивной

 

телеметрической

 

информации

Получен

-

ная

 

выборка

 

составляет

 500 

тыс

режимов

 (60% — 

для

 

обучения

, 20% — 

для

 

валидации

, 20% — 

для

 

тестирования

).

Предложенная

 

модель

 

НС

 

была

 

реализована

 

с

 

помощью

 

библиотеки

 TensorFlow, 

а

 

ее

 

структура

 

была

 

выбрана

 

на

 

основе

 

подбора

 

таких

 

гиперпара

-

метров

как

 

функция

 

активации

количество

 

слоев

функция

 

потерь

 

и

 

т

.

д

Для

 

выбора

 

наилучшей

 

ком

-

бинации

 

гиперпараметров

включая

 

коэффициент

 

регуляризации

 Lasso-

регрессии

 

на

 

стадии

 

отбора

 

признаков

использовалась

 

методика

 

валидации

 

на

 

отложенной

 

выборке

Наилучшие

 

результаты

 

были

 

получены

 

моделью

 

с

 

двумя

 

скрытыми

 

слоя

-

ми

 

и

 

функцией

 

активации

 Leaky ReLU 

с

 

утечкой

 0,2, 

Рис

. 7. 

Значения

 

коэффициентов

 

долевого

 

участия

 

ГЭС

 

и

 

ВЭС

подключенных

 

к

 

АРЧМ

 (

сплошные

 

линии

 — 

ис

-

пользование

 

НС

 

для

 

расчета

 

перетоков

размытые

 

линии

 — 

результаты

 

из

 

раздела

 «

Результаты

 

расчета

 

для

 

энергосистемы

 

с

 

тепло

-, 

гидро

и

 

ветрогенера

-

цией

» 

с

 

использованием

 

метода

 

Н

-

Р

)

обученной

 

с

 

помощью

 

функции

 

потерь

 Huber 

и

 

оп

-

тимизатора

 Adam. 

В

 

результате

 

отбора

 

призна

-

ков

 

только

 89 

параметров

 

из

 578 

использовались

 

в

 

итоговой

 

модели

Средняя

 

абсолютная

 

ошибка

 

обученной

 

модели

 

на

 

тестовой

 

выборке

 

составила

 

менее

 0,03 

МВт

.

Результаты

 

расчета

 

КДУ

 

с

 

использованием

 

мо

-

дели

 

НС

 

для

 

расчета

 

перетоков

 (

рисунок

 7) 

схожи

 

с

 

полученными

 

выше

 

посредством

 

метода

 

Ньюто

-

на

-

Рафсона

Небольшие

 

отличия

 

можно

 

объяснить

 

точностью

 

модели

 

НС

.

Время

 

расчета

 

перетоков

 

активной

 

мощности

 

с

 

помощью

 

предложенной

 

модели

 

НС

 

составляет

 

1,1 

мс

что

 

в

 11 

раз

 

быстрее

 

метода

 

Ньютона

-

Раф

-

сона

Таким

 

образом

время

 

расчета

 

КДУ

 

с

 

исполь

-

зованием

 

НС

 

составило

 30,18 

с

Это

 

доказывает

 

применимость

 

моделей

 

МО

 

для

 

ускорения

 

расче

-

та

 

перетоков

которые

помимо

 

расчета

 

КДУ

могут

 

быть

 

использованы

 

в

 

других

 

онлайн

-

задачах

.

ЗАКЛЮЧЕНИЕ

 

И

 

ВЫВОДЫ

В

 

работе

 

показано

что

 

переменный

 

характер

 

ра

-

боты

 

ветровых

 

электростанций

особенно

 

в

 

случае

 

их

 

использования

 

во

 

вторичном

 

регулировании

требует

 

онлайн

-

расчета

 

коэффициентов

 

долевого

 

участия

 

для

 

повышения

 

эффективности

 

управле

-

ния

 

перетоками

 

активной

 

мощности

В

 

связи

 

с

 

этим

 

авторы

 

разработали

 

программный

 

модуль

который

 

легко

 

встраивается

 

в

 

существующую

 

архитектуру

 

систем

 

АРЧМ

и

 

протестировали

 

его

 

как

 

для

 

стан

-

дартных

так

 

и

 

для

 

более

 

сложных

 

стратегий

 

управ

-

ления

направленных

 

на

 

повышение

 

селективности

 

АРЧМ

 

в

 

сетях

 

с

 

ограниченной

 

пропускной

 

способ

-

ностью

.

В

 

ходе

 

исследований

 

было

 

определено

что

 

тра

-

диционные

 

методы

 

расчета

 

УР

 

позволяют

 

добить

-

ся

 

достаточного

 

быстродействия

 

расчета

 

КДУ

 

с

 

по

-

мощью

 

разработанной

 

методики

 

в

 

энергосистемах

 

с

 

тепло

и

 

гидрогенерацией

Однако

 

при

 

наличии

 

большого

 

числа

 

регулирующих

 

электростанций

 

или

 

при

 

привлечении

 

ВИЭ

 

ко

 

вторичному

 

регулирова

-

нию

 

скорости

 

обновления

 

КДУ

 

может

 

быть

 

недоста

-

точно

Для

 

повышения

 

быстродействия

 

алгоритма

 

0 3 6 9 12 15 18 21 24

t

ч

КДУ

, %

0

–15

–30

–45

–60

ГЭС

1

ВЭС

1

ВЭС

2

ВЭС

3

ГЭС

2

 4 (73) 2022


background image

36

ВОЗОБНОВЛЯЕМАЯ

ЭНЕРГЕТИКА

расчета

 

КДУ

 

авторы

 

предложили

 

использовать

 

и

 

адаптировали

 

к

 

рассматриваемой

 

задаче

 

модель

 

искусственной

 

нейронной

 

сети

 

с

 

полносвязными

 

слоями

 

как

 

альтернативу

 

стандартным

 

методам

 

расчета

 

УР

Описание

 

состояния

 

энергосистемы

 

произведено

 

с

 

помощью

 

вектора

 

инъекций

 

актив

-

ной

 

и

 

реактивной

 

мощности

 

в

 

узлах

размерность

 

которого

 

была

 

снижена

 

относительно

 

решаемой

 

за

-

дачи

 

управления

 

методом

 Lasso-

регрессии

.

Частота

 

обновления

 

значений

 

КДУ

 

определяет

-

ся

 

как

 

вычислительными

 

возможностями

так

 

и

 

эф

-

фективностью

 

управления

В

 

работе

 

авторы

 

показа

-

ли

как

 

изменяются

 

требования

 

к

 

быстродействию

 

расчетов

 

в

 

зависимости

 

от

 

типа

 

регулирующих

 

электростанций

Однако

 

вопрос

 

формулирования

 

более

 

точных

 

требований

 

к

 

скорости

 

расчета

 

КДУ

 

будет

 

рассмотрен

 

в

 

дальнейшем

Кроме

 

того

буду

-

щие

 

направления

 

работы

 

связаны

 

с

 

анализом

 

раз

-

личных

 

методов

 

машинного

 

обучения

 

и

 

с

 

исследо

-

ванием

 

других

 

подходов

 

к

 

понижению

 

размерности

 

вектора

 

признаков

 

для

 

рассматриваемых

 

задач

 

управления

Также

 

авторы

 

планируют

 

протестиро

-

вать

 

различные

 

оптимизационные

 

алгоритмы

 

для

 

повышения

 

точности

 

и

 

скорости

 

расчета

 

КДУ

ЛИТЕРАТУРА

 / REFERENCES

1. 

Информация

 

о

 

фактическом

 

ре

-

жиме

 

работы

 

объектов

 

ДПМ

 

ВИЭ

 

в

 2022 

году

АО

  «

СО

 

ЕЭС

» / 

Information on the actual operating 
mode of RES-operating facilities 
under the power delivery contract 
in 2022. JSC SO UPS. URL: https://
www.so-ups.ru/functioning/markets/
surveys/renewable/2022/. 

2.  Holttinen H., et al. Variable Renew-

able Energy Integration: Status 
Around the World. IEEE Power and 
Energy Magazine, 2021, vol. 19, 
no. 6, pp. 86-96.

3. 

Неуймин

 

В

.

Г

., 

Останин

 

А

.

Ю

., 

Тома

-

лев

 

А

.

А

Внедрение

 

системы

 

мони

-

торинга

 

запасов

 

устойчивости

 

при

 

планировании

 

и

 

управлении

 

элек

-

троэнергетическим

 

режимом

 

ОЭС

 

Сибири

 // 

Энергия

 

единой

 

сети

2019, 

 6(49). 

С

. 32–37 / Neuy-

min V.G., Ostanin A.Yu., Toma-
lev A.A. Introduction of the stability 
margin monitoring system at the 
stage of power mode planning and 
control of the United Power System 
of Siberia //

 Energiya yedinoy seti 

[United Grid Energy], 2019, no. 6(49), 
pp. 32-37. (In Russian)

4. 

Ключевые

 

вызовы

 

ЕЭС

 

России

 

в

 

процессе

 

глобального

 

энерго

-

перехода

 

АО

  «

СО

 

ЕЭС

» / Main 

challenges for UPS of Russia 
in the process of global energy 
transmission of JSC SO UPS.
URL: https://www.so-ups.ru/news/
press-release/press-release-view/
news/17180/.

5. Chernyakhovskiy I., et al. Grid-

Friendly Renewable Energy: Solar 
and Wind Participation in Automatic 
Generation Control Systems, 2019. 
URL: https://www.nrel.gov/docs/
fy19osti/73866.pdf.

6.  Bevrani H., Hiyama T. Intelli-

gent Automatic Generation Con-

trol. Intelligent Automatic Gen-
eration Control, CRC Press, 2017, 
290 p.

7.  Shankar R., Chatterjee K., Chatter-

jee T. K. A control strategy for load 
frequency control coordinating eco-
nomic load dispatch & load forecast-
ing via Kalman 

fi

 lter.  International 

Journal on Electrical Engineering 
and Informatics, 2012, vol. 4, no. 3, 
pp. 495-507.

8. Narvaez A. Implementation and 

tuning of the automatic genera-
tion control in the Energy Na-
tional Control Center of Ecua-
dor – CENACE. 2005 International 
Conference on Industrial Electron-
ics and Control Applications, 2005, 
pp. 1-4.

9. 

Жуков

 

А

.

В

., 

Сацук

 

Е

.

И

., 

Сафро

-

нов

 

А

.

Н

Развитие

 

системы

 

ав

-

томатического

 

вторичного

 

регу

-

лирования

 

частоты

 

и

 

перетоков

 

активной

 

мощности

 

в

 

ЕЭС

 

Рос

-

сии

 // 

Энергия

 

единой

 

сети

, 2014. 

 6(17). 

С

. 24–30 / Zhukov A.V., 

Satsuk E.I., Safronov A.N. Develop-
ment of the automated secondary 
regulation of frequency and active 
power 

fl

 ows in UPS of Russia //

 

Energiya yedinoy seti 

[United Grid 

Energy], 2014, no. 6(17), pp. 24-30. 
(In Russian)

10. Pengxiang W., Wenying L., Zhen-

gyi L. Application of Machine Learn-
ing Methods in Active Power Securi-
ty Correction of Power System. 2010 
International Conference on Electri-
cal and Control Engineering, 2010, 
pp. 476-479.

11. Zhang B., et al. Study on Fast Con-

trol Algorithm Based on Integra-
tive Sensitivity Analysis. 2009 15th 
International Conference on Intel-
ligent  System  Applications  to  Pow-
er Systems, 2009, pp. 1-4. URL: 

https://ieeexplore.ieee.org/docu-
ment/5352960.

12. Xiaobin L., et al. Active power 

fl

 ow 

adjustment based on sensitiv-
ity  analysis  of  DC  load 

fl

 ow  model. 

2012 Power Engineering and Auto-
mation Conference, 2012. pp. 1-4. 
URL: https://ieeexplore.ieee.org/
document/6612479.

13. Donnot  B.  Deep  learning  methods 

for predicting 

fl

 ows in power grids: 

novel architectures and algorithms. 
Université Paris Saclay (COmUE), 
2019. 

 2019SACLS060.

14. Du Y., et al. Achieving 100x Accel-

eration for N-1 Contingency Screen-
ing With Uncertain Scenarios Using 
Deep Convolutional Neural Net-
work. IEEE Transactions on Pow-
er Systems, 2019, vol. 34, no. 4, 
pp. 3303-3305.

15. Yang X.-S. Nature-Inspired Optimi-

zation Algorithms.Oxford: Elsevier, 
2014, pp. 197-211.

16. Draxl C., et al. The Wind Integra-

tion National Dataset (WIND) Tool-
kit. Applied Energy, 2015, vol. 151, 
pp. 355-366.

17. Zhang Y., et al. Investigating Power 

System Primary and Secondary Re-
serve Interaction under High Wind 
Power Penetration. United States, 
2016. URL: https://www.nrel.gov/
docs/fy17osti/64637.pdf.

18. Würth I., et al. Minute-Scale Fore-

casting of Wind Power – Results 
from the Collaborative Workshop 
of IEA Wind Task 32 and 36. Ener-
gies (Basel), 2019, vol. 12, no. 4, 
30 p.

19. Hornik K., Stinchcombe M., White H. 

Multilayer feedforward networks 
are universal approximators. Neu-
ral Networks, 1989, vol. 2, no. 5, 
pp. 359-366.


Оригинал статьи: Автоматическое управление выработкой активной мощности ветроэлектростанций в условиях ограничений пропускной способности сети

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

Читать онлайн

В работе анализируются ограничения существующих систем автоматического регулирования перетоков активной мощности, которые в большей степени проявляются при подключении к энергосистеме источников энергии с переменным характером работы, таких как солнечные и ветровые электростанции. Для преодоления указанных ограничений авторами предложена методика, позволяющая повысить гибкость и селективность управления перетоками активной мощности посредством онлайн-расчета коэффициентов долевого участия регулирующих электростанций. Согласно методике, расчет коэффициентов долевого участия рассматривается как задача оптимизации, решение которой в режиме онлайн позволяет распределять задание внеплановой мощности между регулирующими электростанциями с учетом их характеристик, а также фактической схемнорежимной ситуации в энергосистеме. На основе предложенной методики разработан программный модуль, который взаимодействует с базой данных реального времени и с компонентом для расчета установившегося режима, что позволяет использовать при расчете параметры, актуализируемые в темпе обновления телеметрии. В работе исследованы различные варианты настройки разработанного модуля, позволяющие реализовать как стандартные, так и более сложные стратегии управления, необходимые для сетей с ограниченной пропускной способностью. Кроме того, авторами предложено использовать модель машинного обучения, а именно нейронную сеть с полносвязными слоями, для расчета перетоков активной мощности. Результаты исследования подтверждают, что данное решение позволяет обеспечить повышение скорости расчета коэффициентов долевого участия, необходимое при привлечении источников энергии с быстроизменяющимся характером работы к участию во вторичном регулировании. Также в работе рассматривается применение алгоритма Lasso-регрессии для сокращения размерности расчетной модели энергосистемы в соответствии с задачей управления и повышения быстродействия расчета перетоков алгоритмами машинного обучения. Предложенные подходы были апробированы в рамках расчетных экспериментов, результаты которых подтвердили преимущества их внедрения в существующие системы управления перетоками активной мощности.

Поделиться:

«ЭЛЕКТРОЭНЕРГИЯ. Передача и распределение» № 5(86), сентябрь-октябрь 2024

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

Возобновляемая энергетика / Накопители
Щипицин П.И. Липилин М.С. Дворкин Д.В. Волков М.С. Чемборисова Н.Ш.
«ЭЛЕКТРОЭНЕРГИЯ. Передача и распределение» № 1(82), январь-февраль 2024

Система диагностики АКБ «Репей»

Энергоснабжение / Энергоэффективность Цифровая трансформация / Цифровые сети / Цифровая подстанция Возобновляемая энергетика / Накопители Диагностика и мониторинг
ООО НПП «Микропроцессорные технологии»
«ЭЛЕКТРОЭНЕРГИЯ. Передача и распределение»