配电网分区形成多微网恢复负荷策略

本文梳理了两篇论文的内容，主题都是在配电网中进行微网划分。第一部分的第一篇文章比较简略，部分过程不够详细，而第二部分第二篇文章比较详细，且增加了新的内容，并对其进行了进一步的阐述。本文使用的变量和参数名称在第一部分和第二部分并不一致，主要是为了和原文保持一致。

灾后配电网分区形成多微网恢复负荷(A New Model for Resilient Distribution Systems by Microgrids Formation)

用于解决配电网分区形成多微网恢复负荷的问题。参考文章[1]Ding T, Lin Y, Li G, et al. A new model for resilient distribution systems by microgrids formation[J]. IEEE Transactions on Power Systems, 2017, 32(5): 4145-4147.

模型约束

运行约束

$$\{\begin{array}{l}{P}_{DG,j}-P_{L,j}={\sum }_{s\in \delta (j)}{H}_{js}-{\sum }_{i\in \pi (j)}{H}_{ij}\\ Q_{DG,j}-Q_{L,j}={\sum }_{s\in \delta (j)}{G}_{js}-{\sum }_{i\in \pi (j)}{G}_{ij}\end{array},\forall j\in \mathcal{V}(1)$$
公式(1)显示了功率的平衡，其中，$P_{DG,j}$和$Q_{DG,j}$分别表示节点$j$的DG所发出的有功功率和无功功率。$P_{L,j}$和$Q_{L,j}$分别表示节点$j$的负荷所消耗的有功功率和无功功率。$H_{js}$和$G_{js}$分别表示节点$j$与节点$s$之间的有功功率和无功功率。$\delta (j)$和$\pi (j)$分别表示节点$j$的子节点和父节点。$\mathcal{V}$表示所有节点集合。

$$\{\begin{array}{l}-M\left(1-c_{ij}\right)+\left(r_{ij}{H}_{ij}+x_{ij}{G}_{ij}\right)/U_0\le U_j-U_i\\ \le M\left(1-c_{ij}\right)+\left(r_{ij}{H}_{ij}+x_{ij}{G}_{ij}\right)/U_0\\ -S_{ij}^{\max }{c}_{ij}\le H_{ij}\le S_{ij}^{\max }{c}_{ij},-S_{ij}^{\max }{c}_{ij}\le G_{ij}\le S_{ij}^{\max }{c}_{ij}\end{array},\forall \left(i,j\right)\in \mathcal{E}(2)$$
公式(2)显示了DistFlow方程。具体来说，如果该支路处于闭合状态，则该支路的电压差受潮流约束，支路潮流应被限制；否则，电压差是任意的，支路流必须为零。其中，M是一个足够大的正数，$c_{ij}$表示支路$ij$之间是否闭合，$c_{ij}=1$为闭合，$c_{ij}=0$为断开，$r_{ij}$和$x_{ij}$分别表示支路$ij$之间的电阻和电抗，$U_0$表示基准电压，$U_j$表示节点$j$的电压幅值，$S_{ij}^{\max }$表示支路$ij$之间的最大潮流。$\mathcal{E}$表示所有支路集合。

这个公式说明了，当该支路处于闭合状态时，即$c_{ij}=1$，其潮流应被限制，第一个不等式退化为一个等式，即$\left(r_{ij}{H}_{ij}+x_{ij}{G}_{ij}\right)/U_0=U_j-U_i$；如果没有闭合，即$c_{ij}=0$，其潮流可以是任意的，第一个不等式变为$-M\le U_j-U_i\le M$，即$U_j-U_i$可以是任意的。第二第三个不等式同理。

U j m i n ≤ U j ≤ U j m a x   , j ∈ V \ { D G } ( 3 ) U j = U j 0   , j ∈ { D G } ( 4 ) U_{j}^{\mathrm{min}}\leq U_{j}\leq U_{j}^{\mathrm{max}}\:,\quad j\in\mathcal{V}\backslash\{DG\}\quad(3)\\ U_{j}=U_{j}^{0}\:,\quad j\in\{DG\}\quad(4) Ujmin​≤Uj​≤Ujmax​,j∈V\{DG}(3)Uj​=Uj0​,j∈{DG}(4)

$\mathcal{V}$表示所有节点集合， { D G } \{DG\} {DG}表示DG节点集合。公式(3)表示除了DG节点之外的节点$j$的电压幅值应在一个上下限之间，其中，$U_j^{\min }$和$U_j^{\max }$分别表示节点$j$的电压幅值的下限和上限。公式(4)表示DG节点的电压幅值应等于指定的电压。

P D G , j min ⁡ ≤ P D G , j ≤ P D G , j max ⁡   ,   Q D G , j min ⁡ ≤ Q D G , j ≤ Q D G , j max ⁡   , j ∈ { D G }   ( 5 ) 0 ≤ P L , j ≤ P L , j 0 , j ∈ V ( 6 ) P_{DG,j}^{\min}\leq P_{DG,j}\leq P_{DG,j}^{\max}\:,\:Q_{DG,j}^{\min}\leq Q_{DG,j}\leq Q_{DG,j}^{\max}\:,\quad j\in\left\{DG\right\}\:(5)\\ 0\leq P_{L,j}\leq P_{L,j}^{0},\quad j\in\mathcal{V}\quad(6) PDG,jmin​≤PDG,j​≤PDG,jmax​,QDG,jmin​≤QDG,j​≤QDG,jmax​,j∈{DG}(5)0≤PL,j​≤PL,j0​,j∈V(6)

公式(5)表示DG节点的有功和无功功率应在一个上下限之间，其中，$P_{DG,j}^{\min }$和$P_{DG,j}^{\max }$分别表示DG节点的有功功率的下限和上限，无功功率同理。公式(6)表示节点$j$的负荷所消耗的有功功率应小于正常情况下的负荷需求。

拓扑约束

需要满足放射状拓扑约束，经过精简之后为：

• 每个子图为连通图；
• 分支数等于节点数减去给定子图数。

第一个约束可以考虑为：
设计一个具有相同拓扑结构的虚拟网络，在这个网络中，每个子图只有一个电源，而所有其他节点都有单位的负荷需求，由于虚拟网络与原始电力网络具有相同的拓扑结构，因此它们具有相同的连通性。因此，在虚拟网络中每个节点的功率平衡意味着在源节点和所有其他节点之间至少存在一条路径，从而使该子图必须是连通的。源节点是可以任意选取的，在这里我们选择每个微电网中的分布式发电（DG）节点作为源节点。因此，我们有：

∑ s ∈ δ ( j ) F j s − ∑ i ∈ π ( j ) F i j = − 1 ,   j ∈ V \ { D G } ( 7 ) ∑ s ∈ δ ( j ) F j s − ∑ i ∈ π ( j ) F i j = W j ,   j ∈ { D G } ( 8 ) \sum_{s\in\delta(j)}F_{js}-\sum_{i\in\pi(j)}F_{ij}=-1,\:j\in\mathcal{V}\backslash\{DG\}\quad(7)\\ \sum_{s\in\delta(j)}F_{js}-\sum_{i\in\pi(j)}F_{ij}=W_{j},\:j\in\{DG\}\quad(8) s∈δ(j)∑​Fjs​−i∈π(j)∑​Fij​=−1,j∈V\{DG}(7)s∈δ(j)∑​Fjs​−i∈π(j)∑​Fij​=Wj​,j∈{DG}(8)

其中，$F_{ij}$表示虚拟网络中节点$i$和$j$之间的流量，$W_j$表示虚拟网络中源节点$j$所供给的总功率。$\mathcal{V}$表示所有节点集合， { D G } \{DG\} {DG}表示DG节点集合。公式(7)表示除了DG节点之外的节点$j$都会消耗单位的功率，公式(8)表示DG节点的总流量为$W_j$。

− M c i j ≤ F i j ≤ M c i j ∀ i j ∈ E ( 9 ) − M ( 2 − c i j ) ≤ F i j ≤ M ( 2 − c i j ) ∀ i j ∈ E ( 10 ) W j ≥ 1 j ∈ { D G } ( 11 ) -Mc_{ij}\leq F_{ij}\leq Mc_{ij}\quad\forall ij\in\mathcal{E}\quad(9)\\ -M\left(2-c_{ij}\right)\leq F_{ij}\leq M\left(2-c_{ij}\right)\quad\forall ij\in\mathcal{E}\quad(10)\\ W_j\geq1\quad j\in\{DG\}\quad(11) −Mcij​≤Fij​≤Mcij​∀ij∈E(9)−M(2−cij​)≤Fij​≤M(2−cij​)∀ij∈E(10)Wj​≥1j∈{DG}(11)

公式(9)和(10)表示$F_{ij}$的上下限。公式(11)表示DG节点的总流量至少为1。

第二个约束可以考虑为：
∑ i j ∈ E c i j = ∣ V ∣ − ∣ { D G } ∣ ( 12 ) \sum_{ij\in\mathcal{E}}c_{ij}=\begin{vmatrix}\mathcal{V}\end{vmatrix}-\begin{vmatrix}\{DG\}\end{vmatrix}\quad(12) ij∈E∑​cij​= ​V​ ​− ​{DG}​ ​(12)
该约束与文字描述的约束相同。由于每个微电网只包含一个DG，因此微电网的数量与DG的数量相同。 ∣ { D G } ∣ \begin{vmatrix}\{DG\}\end{vmatrix} ​{DG}​ ​是一个预先设定好的参数。

目标函数

令$w_j$是节点$j$所连接的负荷的权重，$w_j$可以根据节点$j$的重要性或其他因素进行设置。目标函数可以定义为：

$$\begin{gathered} \max \sum _{j\in \mathcal{V}}{w}_i{P}_{L,i}(13-a) \\ \text{s.t. Operation Constraints: (1)-(6)}(13-b) \\ \text{Topology Constraints: (7)-(12)}(13-c) \end{gathered}$$

也就是最大化所有节点的负荷需求。

变量总结

本模型包含的变量数有如下表所示：

变量名	变量类型	所属集合
$c_{ij}$	0-1变量	$\mathcal{E}$
$F_{ij}$	连续变量	$\mathcal{E}$
$U_j$	连续变量	$\mathcal{V}$
$P_{DG,j}$	连续变量	∣ { D G } ∣ \begin{vmatrix}\{DG\}\end{vmatrix} ​{DG}​ ​
$Q_{DG,j}$	连续变量	∣ { D G } ∣ \begin{vmatrix}\{DG\}\end{vmatrix} ​{DG}​ ​
$P_{L,j}$	连续变量	$\mathcal{V}$
$Q_{L,j}$	连续变量	$\mathcal{V}$
$H_{js}$	连续变量	$\mathcal{E}$
$G_{js}$	连续变量	$\mathcal{E}$
$W_j$	连续变量	∣ { D G } ∣ \begin{vmatrix}\{DG\}\end{vmatrix} ​{DG}​ ​

进一步考虑每个微网应存在至少一个可控分布式电源

本文考虑到了包含连结开关的拓扑重构，以及在微网中对于分布式电源使用主从控制。即：

1. 考虑主从分布式发电机的弹性微网形成模型，其中每个孤岛只有一个主发电机，以保证系统的自给自足。
2. 整个系统的拓扑结构可以通过分段开关和联络开关进行重构，从而使一个馈线上的负荷可以转移到微网形成模型中的另一个馈线上，以恢复更多负荷。
3. 另外，采用混合整数二阶锥规划 (MISOCP) 松弛方法来求解所提出的模型。

参考文献[2]Ding T, Lin Y, Bie Z, et al. A resilient microgrid formation strategy for load restoration considering master-slave distributed generators and topology reconfiguration[J]. Applied energy, 2017, 199: 205-216.

约束条件

配电系统条件约束

令$\mathcal{G}=(\mathcal{V},\mathcal{E})$为配电系统的节点和支路集合，其中$\mathcal{V}$表示节点集合，$\mathcal{E}$表示支路集合。灾害之后，会有一些线路被破坏而无法使用，即处于打开状态，则令该状态的线路集合为${\mathcal{E}}_0$，可见存在关系${\mathcal{E}}_0\subset \mathcal{E}$。设0-1决定变量$y_{ij}$，代表线路$ij$是否处于打开状态，即$y_{ij}=1$表示线路$ij$处于闭合状态，$y_{ij}=0$表示线路$ij$处于打开状态。则该变量存在以下约束：

$$y_{ij}=0,\forall ij\in {\mathcal{E}}_0(1)$$

主从控制模式表明，每一个微网中只会有一个DG充当主控制单元，其他的DG都充当从属的控制单元。所以微网数量等于主DG控制单元的数量，令$\Theta$表示DG主控制节点集合，那么$|\Theta |$就是DG主控制节点的数量。

网络拓扑重构下辐射状的约束

辐射状约束可以总结为：

1. 闭合分支的数量等于母线数量减去子图的数量；
2. 保证每个子图的连通性；

则第一个可以用下面约束表示：
$$\sum _{ij\in \mathcal{E}}{y}_{ij}=|\mathcal{V}|-|\Theta |-N-R(2)$$

其中，$|\mathcal{V}|$是节点的数量，$|\Theta |$是DG主控制节点的数量，$N$是变电站的数量，$R$是负荷孤岛的数量，而负荷孤岛就是其所连接的节点既没有变电站也没有DG，仅仅只有用电负荷。

由下图可见，$N$实际上是一个提前设定好的参数，而$R$和故障线的位置、联结开关和DG的位置有关：

想要确定负荷孤岛的数量$R$，相当于找到只包含负载母线的连通分量的数量。无向图的连通分量是该图的连通子图，如果无向图是连通图，则只有一个连通分量。如果有多个连通分量，则可以使用深度优先搜索或广度优先搜索等遍历算法来找到无向图的连通分量。从顶点开始遍历后，将访问从该顶点可达的所有顶点，如果还有其他连通分量，则在这一次遍历完成后仍然会有未访问的顶点。从这些未访问顶点中的一个开始，可以找到另一个连通分量。如果继续此过程，直到访问所有顶点，就可以识别出所有连通分量。这个过程的算法如下所示：

from collections import deque

def find_connected_components(graph):
    visited = [False] * len(graph)
    compNum = 0
    
    for i in range(len(graph)):
        if not visited[i]:
            compNum += 1
            q = deque()
            q.append(i)
            visited[i] = True
            
            while q:
                w = q.popleft()
                
                for k in graph[w]:
                    if not visited[k]:
                        visited[k] = True
                        q.append(k)
                        
    return compNum

# 示例图（邻接表形式）
graph = [
    [1, 2],
    [0, 3],
    [0, 4],
    [1],
    [2],
    [6],
    [5]
]

print(find_connected_components(graph))  

可见，在输入邻接表形式的示例图之后，得到的结果为2，即该图有两个连通分量。

第二个情况可以用下面约束表示，原因在上一章已经阐述过：

$$\begin{gathered} \sum _{s\in \delta (j)}{F}_{js}-\sum _{i\in \pi (j)}{F}_{ij}=-1,j\in \mathcal{V}\setminus \Pi (3) \\ \sum _{s\in \delta (j)}{F}_{js}-\sum _{i\in \pi (j)}{F}_{ij}=W_j,j\in \Pi (4) \\ -M{y}_{ij}⩽F_{ij}⩽M{y}_{ij},\forall ij\in \mathcal{E}(5) \\ -M(2-y_{ij})⩽F_{ij}⩽M(2-y_{ij}),\forall ij\in \mathcal{E}(6) \\ {W}_j⩾1,j\in \Pi (7) \end{gathered}$$

其中，$F_{ij}$是在虚构的配电线ij上的流量，$W_j$是虚构的网络中，电源所共计的总功率。$\Pi$是全体DG的集合。$\delta (j)$是节点j的子节点集合，$\pi (j)$是节点j的父节点集合。$M$是一个足够大的正数。

微网形成的约束

微网形成问题可以视为一个图划分问题。将$\mathcal{G}$分割成$N_P$个子图${\mathcal{G}}_i=({\mathcal{V}}_i,{\mathcal{E}}_i)i=1,2\dots N_P$，则有$\mathcal{G}={\mathcal{G}}_1\cup {\mathcal{G}}_2\dots \cup {\mathcal{G}}_{N_p}$，且${\mathcal{G}}_i\cap {\mathcal{G}}_j=\varnothing ,\forall i\ne j$。同样可以表示为，$\mathcal{V}={\mathcal{V}}_1\cup {\mathcal{V}}_2\dots \cup {\mathcal{V}}_{N_p},{\mathcal{V}}_i\ne \varnothing$，并且${\mathcal{V}}_i\cap {\mathcal{V}}_j=\varnothing ,\forall i\ne j$。

那么通过切割之后，所有的微网的节点形成了一个割集$({\mathcal{V}}_1,{\mathcal{V}}_2,\dots ,{\mathcal{V}}_{N_p})$，定义支路集合${\mathcal{E}}_b\subseteq \mathcal{E}$，集合中的所有支路一边连接的节点是${\mathcal{V}}_i$，另一边是${\mathcal{V}}_j$，且$i\ne j$。

对于每个微网，都需要满足功率流约束。这组方程式被称为分支潮流模型，具体如下：

$$\{\begin{array}{l}{P}_{DG,j}-P_{L,j}={\sum }_{k\in \delta (j)}{H}_{jk}-{\sum }_{i\in \pi (j)}(H_{ij}-r_{ij}{l}_{ij}),\forall j\in \mathcal{V}\\ Q_{DG,j}-Q_{L,j}={\sum }_{k\in \delta (j)}{G}_{jk}-{\sum }_{i\in \pi (j)}(G_{ij}-x_{ij}{l}_{ij})+b_{s_j}{u}_j,\forall j\in \mathcal{V}\\ u_j=u_i-2(r_{ij}{H}_{ij}+x_{ij}{G}_{ij})+(r_{ij}^2+x_{ij}^2)l_{ij},\forall (i,j)\in \mathcal{E}\\ H_{ij}^2+G_{ij}^2=l_{ij}{u}_i,\forall (i,j)\in \mathcal{E}\end{array}(8)$$

其中，$P_{DG,j}$和$Q_{DG,j}$分别表示节点$j$的DG所发出的有功功率和无功功率，$P_{L,j}$和$Q_{L,j}$分别表示节点$j$的负荷所消耗的有功功率和无功功率，$H_{ij}$和$G_{ij}$分别表示节点$i$和节点$j$之间的有功功率和无功功率，$u_i$表示节点$i$的电压幅值平方，$l_{ij}$表示支路$ij$上的电流平方，$b_{s_j}$连接到母线$j$的充电电容大小，$r_{ij}$和$x_{ij}$分别表示支路$ij$上的电阻和电抗。

直观地说，如果支路$ij$断开，那么分支电流和分支潮流都应该为零，即$H_{ij}$,$G_{ij},l_{ij}$都为零，且这两个节点之间的电压幅值不受到公式(8)中的第三个公式的约束。可见，属于割集的边应该被断开，因此，每个微网功率流方程的分割约束可以表示为：

$$\begin{gathered} \{\begin{array}{l}{P}_{DG,j}-P_{L,j}=\sum _{k\in \delta (j)}{H}_{jk}-\sum _{i\in \pi (j)}(H_{ij}-r_{ij}{I}_{ij}),\forall j\in \mathcal{V}\\ Q_{DG,j}-Q_{L,j}=\sum _{k\in \delta (j)}{G}_{jk}-\sum _{i\in \pi (j)}(G_{ij}-x_{ij}{I}_{ij})+b_s{u}_j,\forall j\in \mathcal{V}\end{array}(9) \\ {u}_j-u_i=\{\begin{array}{ll}-2(r_{ij}{H}_{ij}+x_{ij}{G}_{ij})+(r_{ij}^2+x_{ij}^2)I_{ij}& \text{if}(i,j)\in \mathcal{E}\cap y_{ij}=1\\ \text{arbitrary}& \text{otherwise}\end{array}(10) \\ {G}_{ij}\text{or}{H}_{ij}:\{\begin{array}{ll}\in [-S_{ij}^{\text{max}},S_{ij}^{\text{max}}]& \text{if}(i,j)\in \mathcal{E}\cap y_{ij}=1\\ =0& \text{otherwise}\end{array}(11) \\ {l}_{ij}:\{\begin{array}{ll}\in [0,l_{ij}^{\text{max}}]& if(i,j)\in \mathcal{E}\cap y_{ij}=1\\ =0& \text{otherwise}\end{array}(12) \end{gathered}$$

其中，$S_{ij}^{\text{max}}$是支路$ij$上的最大允许潮流，$l_{ij}^{\text{max}}$是支路$ij$上的最大允许电流平方。

然而，约束(10)-(12)是“if-else”类型的约束，意味着一次只能有一个约束成立。这些约束可以使用大M方法等价地转换为一系列仿射约束，如下所示：
$$\begin{gathered} \{\begin{array}{l}-M(1-y_{ij})\le u_j-u_i+2(r_{ij}{H}_{ij}+x_{ij}{G}_{ij})-(r_{ij}^2+x_{ij}^2)l_{ij}\le M(1-y_{ij})\\ -S_{ij}^{\max }{y}_{ij}\le H_{ij}\le S_{ij}^{\max }{y}_{ij},-S_{ij}^{\max }{y}_{ij}\le G_{ij}\le S_{ij}^{\max }{y}_{ij}\\ 0\le l_{ij}\le (l_{ij}^{\max }{)}^2{y}_{ij}\end{array},\forall (i,j)\in \mathcal{E}(13) \\ {H}_{ij}^2+G_{ij}^2=l_{ij}{u}_i,\forall (i,j)\in \mathcal{E}(14) \end{gathered}$$

如果$y_{ij}=1$，则约束(10)-(12)成立的第一个约束成立，否则约束(10)-(12)的第二个约束成立。

配电系统物理约束

$$\begin{gathered} \left\{(P_{DG,j},Q_{DG,j})|\begin{array}{rl}& |Q_{DG,j}|\le P_{DG,j}\tan ({\cos }^{-1}{\theta }_j)\\ & P_{DG,j}^2+Q_{DG,j}^2\le (S_{DG,j}^{\max }{)}^2\\ & 0\le P_{DG,j}\le P_{DG,j}^0\end{array}\right\},\forall j\in \Pi (15) \\ (U_j^{\min }{)}^2\le u_j\le (U_j^{\max }{)}^2,\forall j\in \mathcal{V}\setminus \Theta (16) \\ {u}_j=(U_j^0{)}^2,\forall j\in \Theta (17) \\ \{\begin{array}{l}0\le P_{L,j}\le P_{L,j}^0\\ Q_{L,j}=P_{L,j}\tan ({\cos }^{-1}{\alpha }_j)\end{array},\forall j\in \mathcal{V}(18) \end{gathered}$$

其中，${\theta }_j$表示第$j$个DG的最大功率因数角，$S_{DG,j}^{\max }$表示第$j$个DG的最大容量，$P_{DG,j}^0$表示第$j$个DG的正常情况输出有功功率，$U_j^0$表示第$j$个母线上的主控DG的给定电压幅值，$U_j^{\min },U_j^{\max }$分别表示第$j$个母线上的电压幅值的上下限，$P_{L,j}^0$表示第$j$个负载需求的最大有功功率，${\alpha }_j$表示第$j$个负载需求的功率因数角。

其中，约束(15)描述了每个DG的可行域；约束(16)和(17)规定了电压幅值的限制；约束(18)暗示了有功和无功功率应同时削减，并且应保持每个负荷的功率因数。

目标函数

令$w_j$是负荷的优先级，则目标函数为：
$$\begin{gathered} \max \sum _{j\in \mathcal{V}}{w}_j{P}_{L,j}(19.\text{a}) \\ \text{s.t.}\text{Constraints}(2)-(7),(9),(13)-(18)(19.\text{b}) \end{gathered}$$

凸松弛求解方法

上述问题是一个混合整数非线性非凸的规划(MINNP)问题，关键在于约束(14)是非线性的。可以使用凸松弛方法来找到最优潮流中一个较好的近似解。这可以通过将问题重新表述为连续变量的凸二次锥规划来实现，整个问题可以表述为一个混合整数二阶锥规划 (MISOCP)。

对于二次等式，可以将其松弛为不等式来得到锥松弛，如下所示：
$$H_{ij}^2+G_{ij}^2⩽l_{ij}{u}_i,\forall (i,j)\in \mathcal{E}(20)$$

虽然约束被松弛了，但原有的等式约束仍然存在在其中，所以不会损失这一部分约束。

进一步，上述公式可以被重新化为一个二阶锥公式，如下所示：
$${\Vert \left[\begin{array}{c}2H_{ij}\\ 2G_{ij}\\ l_{ij}-u_i\end{array}\right]\Vert }_2\le l_{ij}+u_i,\forall (i,j)\in \mathcal{E}(21)$$

也就是将约束(21)中左侧的三维向量，求取其二阶范数，可得：
$$\sqrt{(2H_{ij}{)}^2+(2G_{ij}{)}^2+(l_{ij}-u_i{)}^2}\le l_{ij}+u_i$$

两边都是正数，将两边同时进行平方，可得：
$$\begin{gathered} (2H_{ij}{)}^2+(2G_{ij}{)}^2+(l_{ij}-u_i{)}^2\le (l_{ij}+u_i{)}^2 \\ \Rightarrow 4H_{ij}^2+4G_{ij}^2+l_{ij}^2-2l_{ij}{u}_i+u_i^2\le l_{ij}^2+2l_{ij}{u}_i+u_i^2 \end{gathered}$$

经过整理之后就可见约束(21)与约束(20)是等价的。

那么，最终的模型就可以得到如下所示：
$$\begin{gathered} \max \sum _{j\in \mathcal{V}}{w}_j{P}_{L,j}(22.\text{a}) \\ \text{s.t.}\text{Constraints}(2)-(7),(9),(13),(15)-(18)\text{and}(21)(22.\text{b}) \end{gathered}$$

二阶锥松弛评估

为了评估二阶锥松弛的准确性，一个松弛间隙误差定义为：
$$\text{Gap}=\underset{\forall (i,j)\in \mathcal{E}}{\max }|H_{ij}^2+G_{ij}^2-l_{ij}{u}_i|$$

实际上就是将约束(21)中的两侧项移到一边，然后计算计算其差的最大值，如果其间隙差十分接近于0，那么就可以认为该松弛是很精确的。

经过论文[2]的测试，得到其最大间隙为$10^{-4}$左右，十分接近于0，这说明二阶锥松弛和原始的非凸模型几乎完全匹配，从而保证了最优性。

变量与参数总结

变量表如下

变量名	变量说明
$y_{ij}$	二元变量，节点$i$到节点$j$的连接状态
$u_i$	节点$i$的电压幅值平方
$H_{ij}$	配电线路$ij$上的有功潮流
$G_{ij}$	配电线路$ij$上的无功潮流
$l_{ij}$	支路$ij$上的电流平方
$P_{DG,j}$	DG节点$j$的输出有功功率
$Q_{DG,j}$	DG节点$j$的输出无功功率
$P_{L,j}$	节点$j$的负荷所消耗的有功功率
$Q_{L,j}$	节点$j$的负荷所消耗的无功功率
$F_{ij}$	配电线路$ij$上的虚构潮流
$W_j$	虚构网络中源节点所供给的功率

参数表如下

参数名	参数说明	参数名	参数说明
$i,j,s,k$	节点的编号	$ij$	节点$i$和节点$j$之间的支路编号
$\mathcal{E}$	网络中所有支路的集合	$\Pi$	网络中所有DG节点的集合
${\mathcal{E}}_o$	网络中所有断开支路的集合	$\mathcal{V}$	网络中所有节点的集合
$\Theta$	主控DG节点的集合	$\delta (j)$	节点j的子节点集合
$\pi (j)$	节点j的父节点集合	$|\mathcal{E}|$	网络中全部支路的数量
$|\mathcal{V}|$	网络中全部节点的数量	$|\Pi |$	网络中全部DG节点的数量
$|\Theta |$	网络中全部主控DG节点的数量	$r_{ij}$	支路$ij$上的电阻
$x_{ij}$	支路$ij$上的电抗	$b_{s_j}$	节点$j$的充电电容大小
$P_{L,j}^0$	节点$j$的负荷所消耗的正常有功功率	$Q_{L,j}^0$	节点$j$的负荷所消耗的正常无功功率
$P_{DG,j}^0$	第$j$个DG的正常情况输出有功功率	$U_j^0$	第$j$个母线上的主控DG的给定电压幅值
${\theta }_j$	第$j$个DG的最大功率因数角	$S_{DG,j}^{\max }$	第$j$个DG的最大容量
$S_{ij}^{\max }$	支路$ij$上的最大容量	$I_{ij}^{\max }$	支路$ij$上的最大电流
$N$	变电站的数量	$R$	负荷孤岛的数量
$w_j$	负荷的优先级	$M$	一个足够大的正数
${\alpha }_j$	第$j$个负载需求的功率因数角