2025-11-14T17:49:11.423932

Bilevel Optimization for Improved Flexibility Aggregation Models of Electric Vehicle Fleets

Härtel, von Bonin
Electric vehicle (EV) fleets are expected to become an increasingly important source of flexibility for power system operations. However, accurately capturing the flexibility potential of numerous and heterogeneous EVs remains a significant challenge. We propose a bilevel optimization formulation to enhance flexibility aggregations of electric vehicle fleets. The outer level minimizes scheduling deviations between the aggregated and reference EV units, while the inner level maximizes the aggregated unit's profits. Our approach introduces hourly to daily scaling factor mappings to parameterize the aggregated EV units. Compared to simple aggregation methods, the proposed framework reduces the root-mean-square error of charging power by 78~per cent, providing more accurate flexibility representations. The proposed framework also provides a foundation for several potential extensions in future work.
academic

Двухуровневая оптимизация для улучшенных моделей агрегирования гибкости автопарков электромобилей

Основная информация

  • ID статьи: 2506.04843
  • Название: Bilevel Optimization for Improved Flexibility Aggregation Models of Electric Vehicle Fleets
  • Авторы: Philipp Härtel (Fraunhofer IEE), Michael von Bonin (Fraunhofer IEE)
  • Классификация: math.OC cs.SY eess.SY
  • Дата публикации: 5 июня 2025 г. (препринт arXiv)
  • Ссылка на статью: https://arxiv.org/abs/2506.04843

Аннотация

Автопарки электромобилей (ЭМ) ожидаются в качестве все более важного ресурса гибкости при эксплуатации электроэнергетических систем. Однако точное отражение потенциала гибкости многочисленных неоднородных электромобилей остается значительной проблемой. В данной работе предложена формулировка двухуровневой оптимизации для улучшения агрегирования гибкости автопарков ЭМ. Верхний уровень минимизирует отклонение расписания между агрегированным блоком ЭМ и эталонным блоком ЭМ, нижний уровень максимизирует прибыль агрегированного блока. Метод вводит отображение масштабирующих коэффициентов от часовой к суточной шкале для параметризации агрегированного блока ЭМ. По сравнению с простыми методами агрегирования предложенная схема снижает среднеквадратическую ошибку мощности зарядки на 78%, обеспечивая более точное представление гибкости.

Предпосылки исследования и мотивация

Определение проблемы

С массовым распространением электромобилей автопарки ЭМ станут важным ресурсом гибкости электроэнергетической системы. Однако точное моделирование гибкости заряда-разряда крупномасштабных неоднородных автопарков ЭМ сталкивается со значительными проблемами:

  1. Вычислительная сложность: Отдельное моделирование каждого ЭМ вычислительно неосуществимо при крупномасштабном планировании электроэнергетических систем
  2. Проблемы неоднородности: Различные ЭМ имеют различные режимы вождения, емкость батареи, поведение зарядки и другие характеристики
  3. Переоценка гибкости: Существующие простые методы агрегирования часто переоценивают фактический потенциал гибкости автопарка

Значимость исследования

Прогноз немецкого рынка электроэнергии на 2045 год показывает, что медианное потребление электроэнергии ЭМ остается ниже 10 ГВт в течение года, но верхний квартиль часто превышает 40 ГВ, демонстрируя огромную изменчивость. Введение возможностей взаимодействия транспорта и сети (V2G) еще больше усиливает потребность в точном представлении гибкости.

Ограничения существующих методов

  1. Наивное агрегирование: Оценка гибкости автопарка путем простого суммирования доступности зарядки и потребностей в вождении отдельных ЭМ, игнорирующая взаимодействие и зависимости между транспортными средствами
  2. Агрегирование виртуальной батареи: Моделирование автопарка ЭМ как виртуальной системы накопления энергии, но по-прежнему зависит от вручную настраиваемых универсальных масштабирующих коэффициентов
  3. Методы кластеризации: Группировка ЭМ с аналогичными режимами зарядки, но отсутствие оптимизированной стратегии агрегирования

Основные вклады

  1. Предложена схема двухуровневой оптимизации: Впервые применена двухуровневая оптимизация к агрегированию гибкости автопарков ЭМ, верхний уровень оптимизирует масштабирующие коэффициенты для минимизации отклонения расписания, нижний уровень оптимизирует прибыль рыночного участия агрегированного блока
  2. Разработан механизм отображения масштабирующих коэффициентов: Введено отображение масштабирующих коэффициентов временной шкалы от часовой к суточной, способное отражать характеристики гибкости различных временных гранулярностей
  3. Значительное повышение точности агрегирования: По сравнению с простыми методами агрегирования среднеквадратическая ошибка мощности зарядки снижена на 78%, обеспечивая более точное представление гибкости
  4. Предоставлена масштабируемая схема: Обеспечена инфраструктура для будущих расширений (таких как V2G, различные типы пользователей, быстрая зарядка и т.д.)

Подробное описание метода

Определение задачи

Входные данные: Доступность заряда-разряда отдельных блоков ЭМ, траектории состояния заряда (SOC), цены на электроэнергию, потребности в вождении Выходные данные: Оптимизированные масштабирующие коэффициенты агрегированного блока ЭМ (AEV) и стратегия расписания Цель: Минимизация отклонения расписания между AEV и эталонным блоком отдельного ЭМ

Архитектура двухуровневой оптимизации

Верхний уровень (задача верхнего уровня)

Минимизация отклонения расписания между агрегированным блоком и эталонным блоком ЭМ:

minΞOLtT[γC(xu,tCX^Vu,tC)2+γD(xu,tDX^Vu,tD)2+γS(xu,tSX^Vu,tS)2]\min_{\Xi^{OL}} \sum_{t \in T} \left[ \gamma^C (x^C_{u,t} - \hat{X}^C_{V_u,t})^2 + \gamma^D (x^D_{u,t} - \hat{X}^D_{V_u,t})^2 + \gamma^S (x^S_{u,t} - \hat{X}^S_{V_u,t})^2 \right]

Ограничения включают:

  • Отображение масштабирующих коэффициентов: xu,tC=κu,τCvVuXv,tCx^C_{u,t} = \kappa^C_{u,\tau} \sum_{v \in V_u} X^C_{v,t}, где τ=fnC(t)\tau = f^C_n(t)
  • Неотрицательность масштабирующих коэффициентов: κu,τC,κu,τCR+\kappa^C_{u,\tau}, \overline{\kappa}^C_{u,\tau} \in \mathbb{R}_+

Нижний уровень (задача нижнего уровня)

Максимизация прибыли рыночного участия агрегированного блока:

minΞILtTΠt(xu,tCxu,tD)\min_{\Xi^{IL}} \sum_{t \in T} \Pi_t (x^C_{u,t} - x^D_{u,t})

Ограничения:

  • Ограничения мощности заряда-разряда: xu,tCxu,tCxu,tC\underline{x}^C_{u,t} \leq x^C_{u,t} \leq \overline{x}^C_{u,t}
  • Непрерывность SOC: xu,t+1S=ρuxu,tSΦu,tDRΦu,tTH+ηuCxu,tC1ηuDxu,tDx^S_{u,t+1} = \rho_u x^S_{u,t} - \Phi^{DR}_{u,t} - \Phi^{TH}_{u,t} + \eta^C_u x^C_{u,t} - \frac{1}{\eta^D_u} x^D_{u,t}

Механизм отображения масштабирующих коэффициентов

Разработана периодическая функция отображения: fnC(t)=fnD(t)=fnS(t)=tmod(7×24)nf^C_n(t) = f^D_n(t) = f^S_n(t) = \left\lfloor \frac{t \bmod (7 \times 24)}{n} \right\rfloor

Это отображение обеспечивает циклирование в течение недели (168 часов), при этом каждые n часов в цикле группируются, что позволяет обрабатывать каждый рабочий день.

Переформулировка условий KKT

Благодаря выпуклости задачи нижнего уровня двухуровневая задача переформулируется как однослойное смешанное целочисленное линейное программирование:

  • Введены двойственные переменные λu,t,μu,t()\lambda_{u,t}, \mu^{(\cdot)}_{u,t}
  • Добавлены условия оптимальности KKT
  • Использован метод большого M для переформулировки условий дополняющей нежесткости

Экспериментальная установка

Набор данных

  • Временной диапазон: 3-недельный период с часовым разрешением (данные Германии за январь 2012 г.)
  • Конфигурация ЭМ: 900 профилей маршрутов коммутирующих пассажиров на основе исторических данных обследования поездок
  • Цены на электроэнергию: Прогнозируемые цены на электроэнергию в Германии на 2035 год
  • Спецификации транспортных средств: Включая емкость батареи, эффективность вождения, мощность домашней зарядки-разряда, ограничения SOC

Показатели оценки

  • Основной показатель: Среднеквадратическая ошибка (RMSE) мощности зарядки
  • Показатели оптимизации: Значение целевой функции, оптимальный разрыв
  • Измерения сравнения: Траектории SOC, отклонение расписания мощности зарядки

Методы сравнения

  • Simple Aggregation (SA): Эвристический метод агрегирования на основе виртуального накопления энергии
  • Варианты AEV: Различные отображения масштабирующих коэффициентов (1ч, 2ч, 4ч, 6ч, 24ч)

Детали реализации

  • Решатель: Gurobi Optimizer
  • Упрощенная установка: V2G не рассматривается, внимание сосредоточено на проверке основного метода
  • Типы пользователей: Рассматриваются только коммутирующие пассажиры для снижения сложности неоднородности

Результаты экспериментов

Основные результаты

МетодОптимальное значениеОптимальная нижняя границаОтносительный разрывRMSE (МВт)
AEV 1ч2891.691957.9532.29%2.926
AEV 2ч3734.222921.8921.75%3.125
AEV 4ч4130.893709.8910.19%3.164
AEV 6ч4300.004162.433.20%3.256
AEV 24ч5666.595609.961.00%3.827
SA---18.084

Ключевые выводы:

  • Метод AEV достиг снижения RMSE на 78% по сравнению с SA
  • Более тонкозернистые отображения (1ч, 2ч) производят меньшие отклонения расписания
  • Отображение 24ч достигло оптимального разрыва 1%

Анализ масштабирующих коэффициентов

  • Агрегированная максимальная доступность мощности зарядки снижена на 0%-18%
  • Максимальная траектория SOC снижена до 60%
  • Диапазон минимальной траектории SOC составляет 0%-15% от агрегированного максимума SOC
  • Наблюдаются закономерности изменений между различными рабочими днями

Анализ поведения расписания

Из рисунка 4 можно наблюдать:

  • Расписание блока AEV высоко согласуется с расписанием эталонного отдельного ЭМ
  • Метод SA производит более выраженные пики зарядки в периоды низких цен
  • Метод AEV лучше отражает фактические ограничения гибкости

Связанные работы

Методы агрегирования ЭМ

  1. Агрегирование виртуальной батареи: Моделирование автопарка ЭМ как виртуальной системы накопления энергии с учетом ограничений заряда-разряда и емкости батареи
  2. Методы кластеризации: Группировка ЭМ с аналогичными режимами зарядки для избежания смешивания неоднородных транспортных средств
  3. Эквивалентные модели: Заимствование концепций составных или эквивалентных моделей из систем гидроэнергетики

Приложения двухуровневой оптимизации

  • Использование в системах гидроэнергетики для определения точных упрощенных моделей
  • Расчет рынков электроэнергии и планирование емкости
  • Оптимизация реагирования на спрос и накопления энергии

Заключение и обсуждение

Основные выводы

  1. Эффективность метода: Схема двухуровневой оптимизации значительно повышает точность агрегирования гибкости автопарков ЭМ
  2. Вычислительная осуществимость: Эффективное решение достигнуто путем переформулировки KKT
  3. Практическая ценность: Обеспечено точное представление гибкости ЭМ для крупномасштабного планирования электроэнергетических систем

Ограничения

  1. Упрощающие предположения: Текущая версия не рассматривает V2G, зарядку на рабочем месте, быструю зарядку и другие сценарии
  2. Неоднородность пользователей: Рассматриваются только коммутирующие пассажиры, не охватываются различные типы пользователей
  3. Вычислительная сложность: Некоторые экземпляры не достигают глобального оптимума (относительный разрыв > 1%)

Направления будущих исследований

  1. Расширение функциональности: Интеграция возможностей V2G, множественных типов пользователей, инфраструктуры быстрой зарядки
  2. Анализ чувствительности: Различные ценовые сценарии, масштабирующие коэффициенты, специфичные для рабочих дней
  3. Роль агрегатора: Использование двухуровневой структуры для анализа механизмов стимулирования коммерческих агрегаторов гибкости
  4. Приложения в реальном времени: Расширение на расписание в реальном времени и участие на рынке

Глубокая оценка

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

  1. Инновационность метода: Впервые систематически применена двухуровневая оптимизация к агрегированию автопарков ЭМ с прочной теоретической базой
  2. Полнота экспериментов: Сравнение множественных отображений временной шкалы с убедительными результатами
  3. Высокая практическая ценность: Повышение точности на 78% имеет важное значение для планирования электроэнергетических систем
  4. Расширяемость схемы: Обеспечена хорошая инфраструктура для будущих исследований

Недостатки

  1. Ограничение сценариев: Текущие экспериментальные сценарии относительно упрощены, отсутствует проверка важных функций, таких как V2G
  2. Вычислительная эффективность: Некоторые экземпляры имеют значительный разрыв решения, что может повлиять на практическое применение
  3. Чувствительность параметров: Отсутствует анализ чувствительности к ключевым параметрам (таким как веса отклонения γ)
  4. Проверка масштабируемости: Недостаточная проверка производительности метода на автопарках большего размера

Влияние

  1. Академический вклад: Предоставлена новая парадигма оптимизации для области агрегирования ЭМ
  2. Инженерная ценность: Может быть непосредственно применена в инструментах планирования электроэнергетических систем
  3. Поддержка политики: Обеспечена техническая основа для разработки механизмов рынка гибкости ЭМ

Применимые сценарии

  1. Планирование электроэнергетических систем: Оценка гибкости ЭМ при принятии долгосрочных инвестиционных решений
  2. Проектирование рынков: Проектирование механизмов участия коммерческих агрегаторов ЭМ на рынках электроэнергии
  3. Оптимизация эксплуатации: Координированное управление автопарками ЭМ в распределительных сетях

Библиография

Статья цитирует 16 связанных работ, охватывающих в основном:

  • Методы агрегирования ЭМ: агрегирование виртуальной батареи, методы кластеризации
  • Теория двухуровневой оптимизации: условия KKT, смешанное целочисленное переформулирование
  • Эквивалентное моделирование систем гидроэнергетики: составные модели, двухуровневые приложения
  • Немецкие данные о поездках: обследование MiD, статистика транспорта

Общая оценка: Это высококачественная статья с важным вкладом в область агрегирования гибкости автопарков ЭМ. Метод инновационен, эксперименты полны, результаты значительны, обеспечивая эффективное решение ключевых технических проблем интеграции крупномасштабных ЭМ. Несмотря на некоторые ограничения, его инновационная ценность и практический потенциал делают его важным прогрессом в этой области.