2025-11-13T09:52:10.670865

Performance Comparison of Gate-Based and Adiabatic Quantum Computing for Power Flow Analysis

Kaseb, Moller, Palensky et al.
In this paper, we present the first direct comparison between gate-based quantum computing (GQC) and adiabatic quantum computing (AQC) for solving the AC power flow (PF) equations. Building on the Adiabatic Quantum Power Flow (AQPF) algorithm originally designed for annealing platforms, we adapt it to the Quantum Approximate Optimization Algorithm (QAOA). The PF equations are reformulated as a combinatorial optimization problem. Numerical experiments on a 4-bus test system assess solution accuracy and computational time. Results from QAOA are benchmarked against those obtained using D-Wave's Advantage system and Fujitsu's latest generation Digital Annealer, i.e., Quantum-Inspired Integrated Optimization software (QIIO). The findings provide quantitative insights into the performance trade-offs, scalability, and practical viability of GQC versus AQC paradigms for PF analysis, highlighting the potential of quantum algorithms to address the computational challenges associated with modern electricity networks in the Noisy Intermediate-Scale Quantum (NISQ).
academic

Сравнение производительности квантовых вычислений на основе вентилей и адиабатических квантовых вычислений для анализа потокораспределения

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

  • ID статьи: 2510.13378
  • Название: Performance Comparison of Gate-Based and Adiabatic Quantum Computing for Power Flow Analysis
  • Авторы: Zeynab Kaseb, Matthias Möller, Peter Palensky, Pedro P. Vergara
  • Классификация: quant-ph cs.NA cs.SY eess.SY math.NA
  • Дата публикации/конференция: 24-я конференция по вычислениям в энергосистемах (PSCC 2026)
  • Ссылка на статью: https://arxiv.org/abs/2510.13378

Аннотация

В данной статье впервые проводится прямое сравнение производительности квантовых вычислений на основе вентилей (GQC) и адиабатических квантовых вычислений (AQC) при решении уравнений потокораспределения переменного тока (PF). На основе алгоритма адиабатического квантового потокораспределения (AQPF), первоначально разработанного для платформ отжига, алгоритм адаптирован к квантовому приближённому алгоритму оптимизации (QAOA). Уравнения PF переформулированы как задача комбинаторной оптимизации и протестированы на 4-узловой тестовой системе для оценки точности решения и времени вычислений. Результаты QAOA сравниваются с результатами системы D-Wave Advantage и новейшего поколения цифрового отжига Fujitsu (QIIO). Исследование предоставляет количественные сведения о компромиссах производительности, масштабируемости и практической целесообразности парадигм GQC и AQC при анализе потокораспределения, подчёркивая потенциал квантовых алгоритмов в решении вычислительных задач современных электросетей в эпоху шумных квантовых вычислений среднего масштаба (NISQ).

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

Решаемая проблема

Анализ потокораспределения является фундаментальной задачей энергосистем, используемой для расчёта комплексных напряжений на всех шинах при заданных нагрузках, генерации и топологии сети. В сетях переменного тока анализ потокораспределения регулируется законами Кирхгофа, что приводит к системе нелинейных невыпуклых уравнений.

Значимость проблемы

  1. Ограничения традиционных методов: Классические итерационные численные методы (такие как Гаусс-Зейдель, Ньютон-Рафсон) могут отказать при крупномасштабных или плохо обусловленных случаях
  2. Вызовы современных электросетей: С увеличением распределённых энергоресурсов отказы сходимости влияют на надёжность и приводят к ошибочным решениям
  3. Вычислительная сложность: Современные электросети требуют алгоритмов потокораспределения, которые одновременно вычислительно эффективны и численно устойчивы

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

  • Гаусс-Зейдель: Сильно зависит от начальных значений, часто расходится при определённых режимах работы
  • Ньютон-Рафсон: Может не сходиться при сингулярности матрицы Якобиана, вычислительно затратен, показывает плохую производительность при высоких нагрузках или высокой доле возобновляемых источников энергии

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

Переформулирование задачи потокораспределения как задачи комбинаторной оптимизации путём дискретизации комплексных напряжений на шинах с использованием спиновых/двоичных переменных решения и преобразования в представление Изинга или QUBO открывает новые пути для использования преимуществ квантовых вычислений.

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

  1. Первая реализация: Предоставляет первую реализацию анализа комбинаторного потокораспределения с использованием QAOA
  2. Комплексное сравнение: Всестороннее сравнение GQC и AQC для анализа потокораспределения в эпоху NISQ
  3. Адаптация алгоритма: Успешная адаптация алгоритма AQPF, первоначально разработанного для платформ отжига, к QAOA
  4. Оценка производительности: Предоставляет количественные сведения о точности решения, времени вычислений и масштабируемости

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

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

Входные данные: Параметры энергосистемы (мощность генерации PG, мощность нагрузки PD, QD, матрица проводимости Y) Выходные данные: Комплексные напряжения всех шин V = μ + jω Ограничения: Удовлетворение установившихся уравнений баланса мощности

Модель комбинаторного анализа потокораспределения

1. Уравнения баланса мощности

Pi = PGi - PDi, ∀i ∈ {1, ..., N}     (1a)
Qi = QGi - QDi, ∀i ∈ {1, ..., N}     (1b)

2. Представление в прямоугольных координатах

Pi = Σk Gik(μiμk + ωiωk) + Bik(ωiμk - μiωk)     (2a)
Qi = Σk Gik(ωiμk - μiωk) - Bik(μiμk + ωiωk)     (2b)

3. Схема дискретизации

Применяется схема с одной спиновой переменной, где каждому μi и ωi назначается одна спиновая переменная решения:

μi := μ0i + sμi Δμi     (4a)
ωi := ω0i + sωi Δωi     (4b)

где sμi, sωi ∈ {±1} — спиновые переменные решения.

4. Построение модели Изинга

Преобразование задачи в минимизацию суммы квадратов остатков:

min s∈{±1}2N Σi (Pi - PGi + PDi)² + (Qi - QGi + QDi)²     (5)

5. Итеративная стратегия оптимизации

  • Приращения Δμi и Δωi постепенно уменьшаются на каждой итерации, обеспечивая переход от грубого поиска к точной оптимизации
  • Размер приращений контролируется функцией экспоненциального затухания

Реализация QAOA

Конструкция квантовой схемы

|ψ(γ,β)⟩ = ∏k=1^p e^(-iβkHM) e^(-iγkHC) |+⟩^⊗2N     (10)
  • HC: гамильтониан стоимости, специфичный для задачи
  • HM: смешивающий гамильтониан (операторы Паули-X)
  • γ, β: вариационные параметры

Процесс оптимизации

  1. Случайная инициализация параметров γ, β ∈ 0, 2π
  2. Обновление параметров с использованием оптимизатора Adam
  3. Оценка ожидаемой энергии ⟨HC⟩ через 1000 измерений
  4. Итеративная оптимизация до сходимости

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

Тестовая система

  • 4-узловая тестовая система: 1 балансирующий узел, 3 узла нагрузки
  • Размер системы среднего масштаба, подходящий для текущих ограничений оборудования NISQ

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

  1. QAOA: Симулятор вектора состояния lightning.qubit в PennyLane
  2. QA: Система D-Wave Advantage™ (примерно 5000 сверхпроводящих кубитов)
  3. QIIO: Программное обеспечение Fujitsu Quantum-inspired Integrated Optimization (поддерживает 100 000 полносвязных двоичных переменных)

Ключевые параметры

ПараметрQAOAQAQIIO
Количество переменных8 кубитов26 спиновых переменных20 переменных решения
Глубина схемыp=2--
Шаги оптимизации100--
Количество считываний10001000-
Порог сходимости1×10⁻³1×10⁻³1×10⁻³

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

Основные показатели производительности

РешательКол-во переменныхВремя компиляции сКоличество итерацийВремя на итерацию сОстаток
QA260.0032220.0155.18×10⁻⁴
QIIO200.025630.063.31×10⁻⁴
QAOA80.0330015.62.49×10⁻³

Сравнение точности решения

Результаты комплексных напряжений для 4-узловой системы в сравнении с эталонным решением Ньютона-Рафсона:

Методμ₁μ₂μ₃ω₁ω₂ω₃
NR0.9020.9160.890-0.092-0.080-0.104
QA0.9010.9150.889-0.093-0.080-0.105
QIIO0.9010.9150.889-0.092-0.080-0.105
QAOA0.9020.9160.890-0.089-0.078-0.099

Характеристики сходимости

  1. QIIO: Сходится быстрее всего (63 итерации), наивысшая точность
  2. QA: Требует 222 итерации, но в конечном итоге достигает высокой точности
  3. QAOA: Не достигает установленного порога в пределах 300 итераций, но качество решения приемлемо

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

  1. Стабильность оборудования: QA подвержена разрывам соединения на крупных системах, возникают ошибки "Remote end closed connection"
  2. Выбор системы: Advantage2_system1.5 на 20% быстрее, чем Advantage_system6.4, с лучшими результатами
  3. Масштабируемость: Методы AQC уже могут обрабатывать системы с 1354 узлами, тогда как GQC остаётся ограничена малыми масштабами

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

Применение квантовых вычислений в комбинаторной оптимизации

  • QAOA: Теоретически может достичь лучшего коэффициента приближения, чем классические эвристики, на задачах Max-Cut, разбиения графов и др.
  • Квантовый отжиг: Показывает перспективы в тестах оптимизации, имитационный отжиг на доступных машинах превосходит протоколы GQC

Квантовые вычисления в энергосистемах

  • Статья основана на предыдущих работах авторов, впервые применяя модель Изинга к анализу потокораспределения
  • Заполняет пробел в реализации GQC для анализа потокораспределения

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

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

  1. Проверка осуществимости: Все три квантовых метода способны восстановить решения потокораспределения, согласующиеся с классическим методом Ньютона-Рафсона
  2. Различия производительности: Методы AQC (QA, QIIO) превосходят методы GQC (QAOA) на текущем оборудовании
  3. Практическая применимость: QIIO показывает лучшую производительность по количеству итераций и точности

Ограничения

  1. Ограничения масштаба: Эксперименты проводились только на 4-узловой системе, что отражает текущие ограничения вычислительных затрат GQC
  2. Ограничения оборудования: QAOA ограничена количеством кубитов и верностью вентилей
  3. Ограничения NISQ: Шум и время когерентности текущих устройств ограничивают практическое применение

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

  1. Улучшение алгоритмов: Оптимизация параметров QAOA и конструкции схемы
  2. Развитие оборудования: Переоценка производительности с улучшением квантового оборудования
  3. Расширение масштаба: Тестирование масштабируемости алгоритмов на более крупных системах

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

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

  1. Новаторское исследование: Первое систематическое сравнение применения GQC и AQC в анализе потокораспределения
  2. Практическая ценность: Предоставляет важный справочник для применения квантовых вычислений в энергосистемах
  3. Полнота методологии: Полная цепь исследований от теоретического моделирования до экспериментальной проверки
  4. Объективная оценка: Честное сообщение о преимуществах, недостатках и ограничениях каждого метода

Недостатки

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

Влияние

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

Применимые области

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

Список литературы

Статья ссылается на 23 важные работы, охватывающие классические и передовые исследования в области анализа энергосистем, квантовых вычислений и комбинаторной оптимизации, обеспечивая прочную теоретическую базу для исследования.


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