Unbalanced optimal power flow refers to a class of optimization problems subject to the steady state physics of three-phase power grids with nonnegligible phase unbalance. Significant progress on this problem has been made on the mathematical modelling side of unbalanced OPF, however there is a lack of information on implementation aspects as well as data sets for benchmarking. One of the key problems is the lack of definitions of current and voltage bounds across different classes of representations of the power flow equations. Therefore, this tutorial-style paper summarizes the structural features of the unbalanced (optimal) power problem for three-phase systems. The resulting nonlinear complex-value matrix formulations are presented for both the bus injection and branch flow formulation frameworks, which typically cannot be implemented as-is in optimization toolboxes. Therefore, this paper also derives the equivalent real-value formulations, and discusses challenges related to the implementation in optimization modeling toolboxes. The derived formulations can be re-used easily for continuous and discrete optimization problems in distribution networks for a variety of operational and planning problems. Finally, bounds are derived for all variables involved, to further the development of benchmarks for unbalanced optimal power flow, where consensus on bound semantics is a pressing need. We believe benchmarks remain a cornerstone for the development and validation of scalable and reproducible optimization models and tools. The soundness of the derivations is confirmed through numerical experiments, validated w.r.t. OpenDSS for IEEE test feeders with 3x3 impedance matrices.
NomenclatureThis article depends on the definition of a variety of scalar, vector and matrix parameters and variables related to grid buses and branches (see Tables 1-5). The core variables are current, voltage and power, whereas parameters are mainly impedance and variable bounds. Fig. 1 summarizes the variables and parameters defined in the fundamental 3×3 branch model for which the well known The Π-equivalent model used. With respect to balance networks, both series and shunt elements are represented by full complex-valued matrices including the mutual impedance coupling between the conductors. Using the Π-equivalent model, the branch current can be split into a series component s and a shunt component sh , respectively. All circuit element voltages are defined w.r.t. (local) ground voltage , = 0 . Table 1 illustrates typography and mathematical notation used throughout; Table 2 defines sets and indices; Table 3 defines parameters; Table 4 defines typical engineering variables; Table 5 defines lifted variables. Black and red colors indicate real-valued variables and parameters, respectively. Blue and brown colors are used for complex-valued variables and parameters, respectively.