The convex approximation of Flexibility Regions at TSO-DSO Interfaces with SOCP Relaxation: Introduce losses to the objective function
Problem Statement
\[ \mathcal{R} = \{(P_0, Q_0) \in \mathbb{R}^2 : \text{power flow equation, and operational constraints}\} \]
The boundary scanning method from [1] allows us to approximate the flexibility region from the outside by solving a series of OPF problems.
However, the flexibility region is not convex in general. Without guarantees of convexity, it is difficult to use the flexibility region and derive cost functions.
A Case Study of the Convexity
IEEE 33-bus System
System configuration:
- 8 curtailable RES generators
- 2 battery storage systems that could work in 4 quadrants
The left and bottom edges of the flexibility region are convex.
Convexity of the left and bottom edges of the flexibility region with the injection limits
Let:
- \(x\) the active power injection at the slack bus
- \(y\) the reactive power injection at the slack bus
- \(\lambda\) the vector of all other injection variables, \(\Lambda\) a convex feasible set of \(\lambda\)
- \(F(x,y,\lambda) = 0\) the AC power flow equations (non-convex)
The bottom boundary is given by minimizing \(y\) for a given \(x\): \[ f(x) = \inf_{y, \lambda} \{ y \mid (x, y, \lambda) \in \mathcal{S} \} \] where \(\mathcal{S} = \{ (x, y, \lambda) \mid F(x,y,\lambda) = 0, \lambda \in \Lambda \}\) and the flexibility region is given by \(\mathcal{R} = \{ (x,y) \mid \exists \lambda: (x, y, \lambda) \in \mathcal{S} \}\).
The left boundary is given by minimizing \(x\) for a given \(y\): \[ g(y) = \inf_{x, \lambda} \{ x \mid (x, y, \lambda) \in \mathcal{S} \} \]
Because the objective function is linear and therefore convex, \(f(x)\) and \(g(y)\) will be convex if the set \(\mathcal{S}\) is convex.
With SOCP relaxation: \[ f_\text{relaxed}(x) = \min_{y, \lambda} \{ y \mid F_\text{relaxed}(x,y,\lambda) \le 0, \lambda \in \Lambda \} \] Then \(F_\text{relaxed}\) represents a convex set (the intersection of affine hyperplanes and second-order cones), \(f_\text{relaxed}(x)\) and \(g_\text{relaxed}(y)\) are convex functions.
If the SOCP relaxation is exact for every boundary point, then \(F_\text{relaxed}(x,y,\lambda) \le 0\) is equivalent to \(F(x,y,\lambda) = 0\), and the convexity of \(f(x)\) and \(g(y)\) is ensured.
\(g(x)\) and \(f(y)\), i.e., the left and bottom envelopes of the flexibility region with the injection limits, are convex functions if the SOCP relaxation of the AC power flow equations is exact along the boundaries.
AC Power Flow Problems in the Lifted Space
From [2], we learn that the AC power flow equations can be reformulated in a lifted space \(\Psi\).
Starting from the radial AC network, define the lifted variables collected into \(\Psi\):
\[ w_i = |V_i|^2, \]
\[ c_{ik} = V_iV_k\cos(\theta_i-\theta_k), \]
\[ s_{ik} = V_iV_k\sin(\theta_i-\theta_k). \]
The flexibility outputs are linear functions of the lifted variables:
\[ y= \begin{bmatrix} P_s\\ Q_s \end{bmatrix} = C_s\Psi. \]
Let \(u\) represent injections at non-slack buses, all operational quantities become linear:
\[ u=C_{ns}\Psi. \] with injection limits: \[ Au\le b. \]
Operational constraints (voltage limits and line flow limits) are represented as
\[ M\Psi\le m. \]
The only remaining non-convexity is just the rank-1 condition for the \(\Psi\) space:
\[ c_{ik}^2+s_{ik}^2=w_iw_k. \]
Thus the AC flexibility problem becomes
\[ \begin{aligned} &y = C_s\Psi,\\ &u = C_{ns}\Psi,\\ &Au \le b,\\ &M\Psi \le m,\\ &c_{ik}^2+s_{ik}^2 = w_iw_k. \end{aligned} \]
The goal is to characterize a larger region of \(y \in \mathbb{R}^2\).
SOCP Approximation of the Flexibility Region in the Lifted Space
SOCP Relaxation
Replace
\[ c_{ik}^2+s_{ik}^2=w_iw_k \]
with
\[ c_{ik}^2+s_{ik}^2\le w_iw_k. \]
For radial networks this corresponds to relaxing the rank-1 condition of the AC power flow equations.
The feasible set becomes convex:
\[ \mathcal X \subseteq \mathcal X_{\rm soc}. \]
Its projection onto
\[ (P_s,Q_s) \]
is therefore also convex.
Explicit Losses-aware Flexibility Region
The key observation is that the problematic portions of the flexibility boundary are associated with inefficient operating points that trying to maximize losses (the right and top boundaries we observe in Figure 2).
Since total losses are also linear in the lifted variables,
\[ P_{\rm loss} = l_P^T\Psi \sum_{(i,k)} = g_{ik} (w_i+w_k-2c_{ik}), \]
an explicit losses budget can be imposed:
\[ P_{\rm loss}\le \eta_p \quad Q_{\rm loss}\le \eta_q \]
\(\eta_p\) and \(\eta_q\) can be calculated using NLP-based methods (solving 2 OPF problems to estimate).
The resulting feasible set is
\[ \mathcal X_{\rm soc,\eta} \subseteq \mathcal X_{\rm soc}. \]
This step intentionally sacrifices a small amount of flexibility in exchange for:
- avoidance of inefficient operating points,
- improve relaxation tightness,
- preservation of convexity.
If the exactness of the SOCP relaxation is guaranteed for all boundary points, then \(\mathcal X_{\rm soc,\eta}\) is an inner convex approximation of the original flexibility region \(\mathcal X\).
- We cannot guarantee that the SOCP relaxation is exact for all points on the boundary.
- The maximization problem will lead to a maximization of losses, which gives a huge relaxation gap.
- We assume that, with a small losses budget, the relaxation gap will be small / acceptable and only lead to a small neighborhood of the feasible states. (Need to be investigated further)
- We could check the exactness numerically as a post-processing step by solving the original AC power flow equations for the boundary points.
Because \(\mathcal X_{\rm soc,\eta}\) is convex, its projection \(\mathcal F_{\rm soc,\eta} = \{y=C_s\Psi\}\) is also convex.
Unlike the traditional boundary tracing, the flexibility region can therefore be represented by:
- the SOCP model itself as a convex constraint,
- explicit polyhedral approximations (e.g. boundary points → convex hull, supporting hyperplanes),
- inner ellipsoids,
- …
Introduce Losses to the Objective Function
To express the FR explicitly, we solve a series of support-function problems: \[ \min d_\gamma^Ty = \min \cos(\gamma) P_s + \sin(\gamma) Q_s \]
However, as shown here, the explicit losses upper bounds usually lead to numerical instability when trying to maximize the consumption at the slack bus, and thus the resulting flexibility region may lose 50% of its original size and be conservative.
Instead of using hard losses bounds, we can also introduce the losses into the objective function: \[ \min d_\gamma^Ty + l^T\Psi = \min \cos(\gamma) P_s + \sin(\gamma) Q_s + P_{\rm loss} + Q_{\rm loss} \]
For radial networks (shown in the numerical examples below), this method is effective, since the augmented objective is strictly increasing with respect to every branch-current variable involved in the SOCP relaxation:
\[ F_\gamma = d_\gamma^Ty + l^T\Psi = d_\gamma^Ty + \sum_{(i,k)} (r_{ik} + x_{ik}) l_{ik} \] \[ \frac{\partial F_\gamma}{\partial l_{ik}} = r_{ik}(1+\cos\gamma) + x_{ik}(1+\sin\gamma) > 0 \]
Let \(y_\text{aug}^*\) denote the optimal solution of the augmented objective, and \(y^*\) denote the optimal solution of the original AC objective. If the SOCP relaxation is exact at the optimum of the augmented problem and therefore AC feasible, we have: \[ d_\gamma^Ty_\text{aug}^* \ge d_\gamma^Ty^* \] therefore, the augmented boundary does not extend outside the true region along the direction of \(\gamma\).
Numerical Examples
| Case | Time NLP [s] | Memory NLP [MB] | Time SOCP [s] | Memory SOCP [MB] | Area Ratio |
|---|---|---|---|---|---|
| 4-node | 0.07 | 6.22 | 0.01 | 1.02 | 0.99 |
| 12-node | 0.09 | 15.94 | 0.02 | 2.76 | 0.99 |
| 33-node | 0.18 | 50.91 | 0.06 | 7.96 | 0.99 |
| 69-node | 0.39 | 123.10 | 0.14 | 16.97 | 0.99 |
| 74-node | 0.50 | 152.23 | 0.12 | 19.24 | 0.99 |
| 85-node | 0.44 | 181.85 | 0.18 | 21.39 | 0.98 |
| 94-node | 0.52 | 220.12 | 0.18 | 24.88 | 0.98 |
| 118-node | 0.58 | 311.04 | 0.34 | 30.85 | 0.99 |
| 1197-node | 17.85 | 20077 | 5.07 | 607.73 | 0.98 |
| Swiss-PDGs | 0.46 | 188.48 | 0.41 | 22.02 | 0.99 |
Current Challenges
SOCP Exactness
The current approach relies on post-processing verification of the relaxation gap, which is not theoretically strict.
Meshed networks are even more challenging …
Next Step: try AROPF [3] to guarantee exactness of the SOCP relaxation.
Flexibility cost function
The value function provides a natural definition of the optimal flexibility cost: \[ f(P_s,Q_s) = \min_{\Psi} J(\Psi) \]
Currently, the cost function is defined by point-wise evaluation of the value function within the estimated flexibility region.