[论文解读] Quadratically constrained quadratic programs on acyclic graphs with application to power flow

Subhonmesh Bose, Dennice F. Gayme+2 more|arXiv (Cornell University)|Mar 26, 2012

Vehicle Routing Optimization Methods被引用 33

一句话总结

该论文证明了当约束矩阵满足一个技术性条件时，具有无环图结构的非凸二次约束二次规划（QCQPs）可在多项式时间内求解。该结果被应用于径向配电网中的最优潮流（OPF）问题，表明可通过半定规划松弛（SDR）与秩恢复获得全局最优解；当该条件不成立时，提出了基于扰动的启发式方法。

ABSTRACT

This paper proves that non-convex quadratically constrained quadratic programs can be solved in polynomial time when their underlying graph is acyclic, provided the constraints satisfy a certain technical condition. When this condition is not satisfied, we propose a heuristic to obtain a feasible point. We demonstrate this approach on optimal power flow problems over radial networks.

研究动机与目标

• 识别一类可在非凸性条件下于多项式时间内求解的无环图上的非凸QCQPs。
• 通过识别秩一恢复的充分条件，将半定规划松弛（SDR）的适用范围扩展至非凸QCQPs。
• 将理论结果应用于径向配电网中的最优潮流（OPF）问题，其中底层图结构为无环图。
• 为不满足充分条件的QCQPs开发基于扰动的启发式方法，以确保获得可行的近似最优解。
• 通过仿真表明，该启发式方法能可靠地生成径向网络中OPF问题的近似最优可行解。

提出的方法

• 定义一个使用复变量和厄米特矩阵的QCQP，其中目标函数和约束均为二次型。
• 根据约束矩阵和目标矩阵的稀疏模式构建无向图，变量作为顶点，矩阵元素非零的位置作为边。
• 证明：若所得图是一棵树（即无环且连通），且约束矩阵满足与最小半正定秩相关的技术条件，则可通过SDR在多项式时间内求解该QCQP。
• 引入一种扰动技术：在目标矩阵中添加一个较小的正定项，使SDR变为严格凸，从而确保在极限情况下可恢复秩一解。
• 使用一系列扰动问题，其中扰动参数δ逐渐减小，对每个问题求解SDR，得到一系列收敛于原始QCQP全局最优解的解。
• 对于秩一条件不成立的情况，采用固定小δ的扰动方法，以在多项式时间内获得目标值与最优值相差不超过ε的可行解，确保多项式时间近似。

实验结果

研究问题

• RQ1 在何种条件下，无环图上的非凸QCQPs可被多项式时间求解？
• RQ2 当约束矩阵非半正定时，无环图上QCQPs的半定规划松弛（SDR）能否获得全局最优解？
• RQ3 最小半正定秩条件在确保从SDR中实现秩一恢复方面起什么作用？
• RQ4 当充分条件不满足时，如何利用基于扰动的方法获得近似最优的可行解？
• RQ5 该方法能否有效应用于径向配电网中的最优潮流（OPF）问题？

主要发现

• 若约束矩阵满足与最小半正定秩相关的技术条件，则无环图上的非凸QCQPs可在多项式时间内求解。
• 当无环图上QCQPs的SDR产生秩一解时，可精确恢复原始QCQP的最优解。
• 对于SDR解非秩一的QCQPs，扰动技术可确保在多项式时间内计算出与全局最优解相差不超过ε的可行解。
• 该扰动方法可保证收敛至一个可行点，其目标值可任意接近真实最优值，且误差受用户定义容差的约束。
• 该理论框架成功应用于径向网络中的最优潮流（OPF）问题，其中无环结构使得通过SDR与秩恢复实现高效求解成为可能。
• 仿真结果表明，即使充分条件不满足，该启发式方法始终能为径向网络中的OPF问题生成近似最优的可行解。