二阶锥形式的配电网最优潮流问题（SOCP-OPF）推导过程

1.配电网潮流常规模型

一个简单的支路潮流示意图如下图所示，${\dot{V}}_i$、${\dot{V}}_i$表示节点$i$和节点$j$的复电压，$z_{ij}$、${\dot{I}}_{ij}$、${\tilde{S}}_{ij}$分别表示支路$i-j$上的复阻抗、复电流、复功率，${\tilde{S}}_j$表示节点$j$注入的复功率。

则潮流方程如下所示：
$${\dot{V}}_i-{\dot{V}}_j=z_{ij}{\dot{I}}_{ij}\text{\hspace{1em}(1)}$$
$${\tilde{S}}_{ij}={\dot{V}}_i{\dot{I}}_{ij}^{\ast }\text{\hspace{1em}(2)}$$
$${\tilde{S}}_j=\sum _{k\in \delta (j)}{\tilde{S}}_{jk}-\sum _{i\in \pi (j)}({\tilde{S}}_{ij}-z_{ij}|{\dot{I}}_{ij}{|}^2)\text{\hspace{1em}(3)}$$
其中，$k\in \delta (j)$表示功率从节点$j$流出的所有节点，$i\in \pi (j)$表示功率流向节点$j$的所有节点。

2.去相角松弛

式(1)、(2)、(3)中变量都为复数，首先将其松弛为实数。上式(3)可直接改写为如下形式：
$$p_j=\sum _{k\in \delta (j)}{p}_{jk}-\sum _{i\in \pi (j)}(p_{ij}-r_{ij}{I}_{ij}^2)\text{\hspace{1em}(4)}$$
$$q_j=\sum _{k\in \delta (j)}{q}_{jk}-\sum _{i\in \pi (j)}(q_{ij}-x_{ij}{I}_{ij}^2)\text{\hspace{1em}(5)}$$
将常规模型中的式(2)代入式(1)，得到下式
$${\dot{V}}_i-{\dot{V}}_j=z_{ij}\frac{{\tilde{S}}_{ij}^{\ast }}{{\dot{V}}_i^{\ast }}\Rightarrow {\dot{V}}_j={\dot{V}}_i-z_{ij}\frac{{\tilde{S}}_{ij}^{\ast }}{{\dot{V}}_i^{\ast }}$$
左右同乘以共轭，可得：
$${\dot{V}}_j{\dot{V}}_j^{\ast }=({\dot{V}}_i-z_{ij}\frac{{\tilde{S}}_{ij}^{\ast }}{{\dot{V}}_i^{\ast }})({\dot{V}}_i^{\ast }-z_{ij}^{\ast }\frac{{\tilde{S}}_{ij}}{{\dot{V}}_i})={\dot{V}}_i{\dot{V}}_i^{\ast }-z_{ij}{\tilde{S}}_{ij}^{\ast }-z_{ij}^{\ast }{\tilde{S}}_{ij}+z_{ij}{z}_{ij}^{\ast }\frac{{\tilde{S}}_{ij}{\tilde{S}}_{ij}^{\ast }}{{\dot{V}}_i{\dot{V}}_i^{\ast }}$$
将上式的实部虚部写开可得（其中$\text{j}$为虚数单位，${\text{j}}^2=-1$）
$$\begin{gathered} \Rightarrow V_j^2=V_i^2-[(r_{ij}+\text{j}{x}_{ij})(p_{ij}-\text{j}{q}_{ij})+(r_{ij}-\text{j}{x}_{ij})(p_{ij}+\text{j}{q}_{ij})]+|z_{ij}{|}^2{I}_{ij}^2 \\ =V_i^2-[(r_{ij}{p}_{ij}+x_{ij}{q}_{ij}-\text{j}{r}_{ij}{q}_{ij}+\text{j}{x}_{ij}{p}_{ij})+(r_{ij}{p}_{ij}+x_{ij}{q}_{ij}+\text{j}{r}_{ij}{q}_{ij}-\text{j}{x}_{ij}{p}_{ij})]+(r_{ij}^2+x_{ij}^2)I_{ij}^2 \end{gathered}$$
整理后可得：
$$V_j^2=V_i^2-2(r_{ij}{p}_{ij}+x_{ij}{q}_{ij})+(r_{ij}^2+x_{ij}^2)I_{ij}^2\text{\hspace{1em}(6)}$$
再对式(2)左右同乘以共轭
$${\tilde{S}}_{ij}{\tilde{S}}_{ij}^{\ast }={\dot{V}}_i{\dot{V}}_i^{\ast }{\dot{I}}_{ij}{\dot{I}}_{ij}^{\ast }$$
整理可得：
$$V_i^2{I}_{ij}^2=p_{ij}^2+q_{ij}^2\text{\hspace{1em}(7)}$$

3.二阶锥松弛

式(7)的最高次数为四次，首先降次。设$l_{ij}=I_{ij}^2$，$v_i=V_i^2$，则式(7)可改写为：
$$l_{ij}{v}_i=p_{ij}^2+q_{ij}^2$$
此时上式为二次等式约束，不便处理，进行二阶锥转化。
将等于号松弛为大于等于号，原最优潮流问题等价（证明详见参考文献[1]），可得下式：
$$l_{ij}{v}_i\ge p_{ij}^2+q_{ij}^2$$
左右两边乘4，加$(l_{ij}-v_i{)}^2$，可得
$$\begin{gathered} 4l_{ij}{v}_i+l_{ij}^2+v_i^2-2l_{ij}{v}_i\ge 4p_{ij}^2+4q_{ij}^2+(l_{ij}-v_i{)}^2 \\ \Rightarrow (l_{ij}+v_i{)}^2\ge (2p_{ij}{)}^2+(2q_{ij}{)}^2+(l_{ij}-v_i{)}^2 \end{gathered}$$
写成范数形式，得到下式：
$$l_{ij}+v_i\ge {\Vert \begin{array}{c}2p_{ij}\\ 2p_{ij}\\ l_{ij}-v_i\end{array}\Vert }_2\text{\hspace{1em}(8)}$$

4.二阶锥形式的最优潮流问题（SOCP-OPF）

综上，式(4)、(5)、(6)、(8)构成了二阶锥形式的潮流约束条件。则SOCP-OPF问题的形式如下所示：
$$\begin{gathered} \min f=P_{loss}=\sum _{(i,j)\in {\Omega }_{Network}}{l}_{ij}{r}_{ij} \\ \text{s.t.}{p}_j=\sum _{k\in \delta (j)}{p}_{jk}-\sum _{i\in \pi (j)}(p_{ij}-r_{ij}{l}_{ij}) \\ {q}_j=\sum _{k\in \delta (j)}{q}_{jk}-\sum _{i\in \pi (j)}(q_{ij}-x_{ij}{l}_{ij}) \\ {v}_j=v_i-2(r_{ij}{p}_{ij}+x_{ij}{q}_{ij})+(r_{ij}^2+x_{ij}^2)l_{ij} \\ {\Vert \begin{array}{c}2p_{ij}\\ 2p_{ij}\\ l_{ij}-v_i\end{array}\Vert }_2\le l_{ij}+v_i \\ {\underline{I}}_{ij}^2\le l_{ij}\le {\overline{I}}_{ij}^2 \\ {\underline{V}}_i^2\le v_i\le {\overline{V}}_i^2 \end{gathered}$$

参考文献

[1] M. Farivar和S. H. Low, 《Branch Flow Model: Relaxations and Convexification—Part I》, IEEE Trans. Power Syst., 卷 28, 期 3, 页 2554–2564, 8月 2013, doi: 10.1109/TPWRS.2013.2255317.
[2] 主动配电网SOCP_OPF学习笔记（2）推导SOCP_OPF模型