Метод увеличенного лагранжиана на GPU для безопасности-ограниченного оптимального потока мощности переменного тока

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

• ID статьи: 2510.13333
• Название: An Augmented Lagrangian Method on GPU for Security-Constrained AC Optimal Power Flow
• Авторы: François Pacaud, Armin Nurkanović, Anton Pozharskiy, Alexis Montoison, Sungho Shin
• Категория: math.OC (Оптимизация и управление)
• Дата публикации: 15 октября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2510.13333

Аннотация

В данной работе предложен новый алгоритм для решения крупномасштабной задачи безопасности-ограниченного оптимального потока мощности переменного тока (AC-SCOPF). Метод основан на методе увеличенного лагранжиана с нелинейными ограничениями (NCL) с использованием внутренних точечных методов для решения подзадач. NCL имеет два ключевых преимущества для крупномасштабных задач SC-OPF: во-первых, NCL может обрабатывать сложные задачи, такие как неразрешимые задачи или модели с дополнительными ограничениями; во-вторых, член увеличенного лагранжиана естественным образом регуляризует ньютоновскую линейную систему во внутренних точечных методах, что позволяет решать ньютоновскую систему без поворота, что может быть эффективно распараллелено на GPU. Численные результаты показывают, что MadNCL может полностью решить AC-SCOPF с 500 узлами и 256 отказами на GPU менее чем за 3 минуты, в то время как Knitro требует более 3 часов для поиска эквивалентного решения.

Исследовательский контекст и мотивация

Описание проблемы

В сетях передачи электроэнергии оптимальное планирование обычно вычисляется путем решения задачи безопасности-ограниченного оптимального потока мощности (SCOPF). Это планирование минимизирует заданный критерий (стоимость или потери в сети) при учете физических ограничений (потоки мощности, ограничения потока линий) и емкости генераторов. Кроме того, планирование должно оставаться допустимым при ряде аварийных ситуаций, соответствующих отказам линий или генераторов (критерий безопасности N-1).

Важность проблемы

SCOPF обычно формулируется как крупномасштабная задача линейного программирования, называемая DC-SCOPF, размер которой растет линейно с количеством отказов. Однако это достигается за счет линеаризации нелинейных физических ограничений, что влияет на точность решения. Тем не менее, решение AC-SCOPF с использованием исходной нелинейной формулировки остается открытой проблемой.

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

Нелинейная формулировка имеет две проблемы:

1. AC-SCOPF является огромной нелинейной задачей программирования, размер которой превышает возможности передовых нелинейных оптимизационных решателей, таких как Ipopt или Knitro
2. Дополнительные ограничения в модели AC-SCOPF численно сложны в обработке и требуют специализированных алгоритмов

Исследовательская мотивация

Характеристики крупномасштабной AC-SCOPF подталкивают алгоритм к границам, поскольку количество дополнительных ограничений растет линейно с количеством отказов. Для решения этой проблемы авторы предлагают использовать метод увеличенного лагранжиана на основе метода увеличенного лагранжиана с нелинейными ограничениями (NCL) для решения AC-SCOPF.

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

1. Первое применение метода увеличенного лагранжиана: Впервые применены алгоритмы на основе увеличенного лагранжиана для решения корректирующей SCOPF с дополнительными ограничениями
2. Реализация с ускорением на GPU: Разработан MadNCL, реализация NCL с поддержкой ускорения на GPU, использующая NVIDIA cuDSS для линейного решения
3. Обработка дополнительных ограничений: Продемонстрировано, что MadNCL хорошо обрабатывает дополнительные ограничения в AC-SCOPF и эффективно обнаруживает неразрешимые задачи
4. Значительное повышение производительности: Достигнуто значительное ускорение по сравнению с традиционными методами, версия GPU более чем в 6 раз быстрее версии CPU

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

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

Задача безопасности-ограниченного оптимального потока мощности переменного тока (AC-SCOPF) определяется как:

$\min_{x,u} f(x^0, u^0)$

при ограничениях:

• Базовый случай: $g^0(x^0, u^0) = 0$, $h^0(x^0, u^0) \leq 0$
• Для каждого отказа $k \in \{1, \cdots, K\}$:
  • $g^k(u^0, x^k, u^k) = 0$
  • $h^k(x^k, u^k) \leq 0$
  • $0 \leq G(x^k, u^k) \perp H(x^k, u^k) \geq 0$

где дополнительные ограничения возникают из:

1. Автоматического регулирования выработки (AGC): регулирование активной мощности для стабилизации частоты
2. Переключения PV/PQ: управление напряжением и ограничения реактивной мощности

Архитектура модели

Переформулировка MPCC

Переформулировка AC-SCOPF как математического программирования с дополнительными ограничениями (MPCC):

$\min_{w \in \mathbb{R}^n} \phi(w) \quad \text{s.t.} \quad \begin{cases} c(w) = 0, \quad w_0 \geq 0 \\ 0 \leq w_1 \perp w_2 \geq 0 \end{cases}$

Алгоритм NCL

NCL работает на двух уровнях:

• Внешняя итерация: обновление параметра штрафа $\rho^{(n)}$ и оценок множителей $(λ^{(n)}, ν_0^{(n)})$
• Внутренняя итерация: решение ограниченной нелинейной подзадачи:

$\min_{w,r,t} L_ρ(w, r, t, λ^{(n)}, ν_0^{(n)})$

при ограничениях: $c(w) - r = 0, \quad W_1W_2e \leq t, \quad (w_0, w_1, w_2) \geq 0$

Структура ньютоновской системы

Ньютоновская система подзадачи имеет следующую структуру:

$\begin{bmatrix} A & B^⊤ \\ B & -C \end{bmatrix} \begin{bmatrix} Δw \\ Δy \end{bmatrix} = \begin{bmatrix} r_1 \\ r_2 \end{bmatrix}$

где регуляризация, обеспечиваемая членом увеличенного лагранжиана, позволяет использовать разложение без поворота.

Технические инновации

1. Естественная регуляризация: член увеличенного лагранжиана естественным образом регуляризует ньютоновскую линейную систему, сохраняя невырожденность системы даже когда строгая дополнительность не выполняется
2. Разложение без поворота: регуляризация позволяет использовать методы без поворота, такие как символическое разложение Холецкого, которые могут быть эффективно распараллелены на GPU
3. Обнаружение неразрешимости: когда задача неразрешима, NCL автоматически переходит к задаче допустимости путем увеличения параметра штрафа $ρ^{(n)}$ до бесконечности

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

Наборы данных

Использованы примеры из библиотеки MATPOWER:

• 118ieee, ACTIVSg200, 300ieee, ACTIVSg500
• 1354pegase, ACTIVSg2000, 2869pegase
• Количество отказов варьируется от 2 до 256

Метрики оценки

• Время решения: общее время решения и время на одну итерацию
• Количество итераций: количество итераций внутренних точечных методов
• Значение целевой функции: значение целевой функции оптимального решения
• Допустимость: способность обнаруживать неразрешимые отказы

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

• Knitro: передовой оптимизационный решатель с поддержкой MPCC, использующий метод точного штрафа $\ell_1$
• MadNCL-CPU: версия CPU с использованием HSL MA57
• MadNCL-GPU: версия GPU с использованием NVIDIA cuDSS

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

• Язык программирования: Julia 1.11
• Допуск сходимости: 1e-6
• Минимальный параметр барьера: $μ_{min} = 10^{-7}$
• Оборудование: процессор AMD EPYC 7430, GPU NVIDIA A30 (24GB памяти)

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

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

Производительность фильтрации отказов

При фильтрации отказов MadNCL значительно превосходит Knitro:

MadNCL как минимум в 10 раз быстрее на примерах с более чем 300 узлами.

Решение крупномасштабной AC-SCOPF

Для примера ACTIVSg500 с увеличением количества отказов:

Абляционные эксперименты

Производительность GPU vs CPU

Повышение производительности MadNCL-GPU по сравнению с MadNCL-CPU:

• Для K≥64 версия GPU примерно в 6 раз быстрее версии CPU
• Для K≥64 версия GPU более чем в 20 раз быстрее Knitro

Анализ времени на одну итерацию

С увеличением количества отказов время MadNCL-GPU на одну итерацию растет наиболее медленно, демонстрируя хорошую масштабируемость.

Экспериментальные выводы

1. Масштабируемость: MadNCL демонстрирует отличную масштабируемость, способен обрабатывать задачи с почти миллионом переменных
2. Робастность: NCL может автоматически обнаруживать и обрабатывать неразрешимые задачи
3. Эффективность параллелизма: реализация GPU полностью использует преимущества параллельных вычислений
4. Численная стабильность: регуляризация увеличенного лагранжиана повышает численную стабильность

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

Основные направления исследований

1. Методы решения MPCC: включая прямые методы, методы регуляризации и методы штрафа
2. Оптимизация электроэнергетических систем: различные стратегии решения DC-SCOPF и AC-SCOPF
3. Ускорение на GPU: перенос алгоритмов оптимизации на платформу GPU

Вклад данной работы

По сравнению с существующими работами, в данной статье впервые применен метод увеличенного лагранжиана к AC-SCOPF с дополнительными ограничениями и реализовано эффективное ускорение на GPU.

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

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

1. MadNCL может эффективно решать крупномасштабные задачи AC-SCOPF, обрабатывая почти миллион переменных
2. Версия с ускорением на GPU достигает десятикратного повышения производительности по сравнению с традиционными CPU решателями
3. Метод увеличенного лагранжиана обеспечивает робастное решение для обработки дополнительных ограничений

Ограничения

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

Будущие направления

1. Решение проблемы плохой обусловленности ньютоновской системы во внутренних точечных методах
2. Расширение на более крупные примеры (десятки тысяч узлов, сотни отказов)
3. Улучшение методов предобусловливания для повышения численной стабильности

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

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

1. Инновационность метода: впервые применен NCL к AC-SCOPF, новый технический подход
2. Качество реализации: высокое качество реализации на GPU, полное использование преимуществ параллельных вычислений
3. Полнота экспериментов: комплексная экспериментальная оценка, включая тесты масштабируемости и робастности
4. Практическая ценность: значительное повышение производительности делает возможным крупномасштабные приложения в реальном времени

Недостатки

1. Теоретический анализ: отсутствует анализ сходимости NCL на задачах SCOPF
2. Численная стабильность: остаются проблемы численной стабильности на примерах максимального размера
3. Универсальность: применимость метода в основном ограничена оптимизацией электроэнергетических систем

Влияние

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

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

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

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

Статья цитирует 31 связанную работу, охватывающую моделирование SCOPF, методы решения MPCC, теорию увеличенного лагранжиана и оптимизацию на GPU, обеспечивая прочную теоретическую основу для исследования.