Novel Approach Towards Global Optimality of Optimal Power Flow Using Quadratic Convex Optimization

Optimal Power Flow (OPF) can be modeled as a nonconvex Quadratically Constrained Quadratic Program (QCQP). Our purpose is to solve OPF to global optimality. To this end, we specialize the Mixed-Integer Quadratic Convex Reformulation method (MIQCR) to (OPF). This is a method in two steps. First, a Semi-Definite Programming (SDP) relaxation of (OPF) is solved. Then the optimal dual variables of this relaxation are used to reformulate OPF into an equivalent new quadratic program, where all the non-convexity is moved to one additional constraint. In the second step, this reformulation is solved within a branch-and-bound algorithm, where at each node a quadratic and convex relaxation of the reformulated problem, obtained by relaxing the non-convex added constraint, is solved. The key point of our approach is that the lower bound at the root node of the branch-and-bound tree is equal to the SDP relaxation value. We test this method on several OPF cases, from two-bus networks to more-than-a-thousand-buses networks from the MAT-POWER repository. Our first results are very encouraging.

[1]  Nikolaos V. Sahinidis,et al.  A polyhedral branch-and-cut approach to global optimization , 2005, Math. Program..

[2]  Jesse T. Holzer,et al.  Implementation of a Large-Scale Optimal Power Flow Solver Based on Semidefinite Programming , 2013, IEEE Transactions on Power Systems.

[3]  Santanu S. Dey,et al.  Matrix minor reformulation and SOCP-based spatial branch-and-cut method for the AC optimal power flow problem , 2017, Math. Program. Comput..

[4]  Paul A. Trodden,et al.  Local Solutions of the Optimal Power Flow Problem , 2013, IEEE Transactions on Power Systems.

[5]  Garth P. McCormick,et al.  Computability of global solutions to factorable nonconvex programs: Part I — Convex underestimating problems , 1976, Math. Program..

[6]  C. Josz Application of polynomial optimization to electricity transmission networks , 2016, 1608.03871.

[7]  Lorenz T. Biegler,et al.  Global optimization of Optimal Power Flow using a branch & bound algorithm , 2012, 2012 50th Annual Allerton Conference on Communication, Control, and Computing (Allerton).

[8]  Alain Billionnet,et al.  Extending the QCR method to general mixed-integer programs , 2010, Mathematical Programming.

[9]  R D Zimmerman,et al.  MATPOWER: Steady-State Operations, Planning, and Analysis Tools for Power Systems Research and Education , 2011, IEEE Transactions on Power Systems.

[10]  Jean Charles Gilbert,et al.  Application of the Moment-SOS Approach to Global Optimization of the OPF Problem , 2013, IEEE Transactions on Power Systems.

[11]  Daniel K. Molzahn,et al.  Moment/Sum-of-Squares Hierarchy for Complex Polynomial Optimization , 2015, 1508.02068.

[12]  Leo Liberti,et al.  Branching and bounds tighteningtechniques for non-convex MINLP , 2009, Optim. Methods Softw..

[13]  J. Lavaei,et al.  Network Topologies Guaranteeing Zero Duality Gap for Optimal Power Flow Problem , 2013 .

[14]  Jorge Nocedal,et al.  Knitro: An Integrated Package for Nonlinear Optimization , 2006 .

[15]  S. Low,et al.  Zero Duality Gap in Optimal Power Flow Problem , 2012, IEEE Transactions on Power Systems.

[16]  Jean B. Lasserre,et al.  Convexity in SemiAlgebraic Geometry and Polynomial Optimization , 2008, SIAM J. Optim..

[17]  Sourour Elloumi,et al.  Global solution of non-convex quadratically constrained quadratic programs , 2019, Optim. Methods Softw..

[18]  S. Low Convex relaxation of optimal power flow: A tutorial , 2013, 2013 IREP Symposium Bulk Power System Dynamics and Control - IX Optimization, Security and Control of the Emerging Power Grid.

[19]  Alain Billionnet,et al.  Exact quadratic convex reformulations of mixed-integer quadratically constrained problems , 2015, Mathematical Programming.

[20]  Jean Maeght,et al.  AC Power Flow Data in MATPOWER and QCQP Format: iTesla, RTE Snapshots, and PEGASE , 2016, 1603.01533.