CN102891485A - 基于序分量法的弱环配电网三相解耦潮流计算方法 - Google Patents

Abstract

本发明公开了一种基于序分量法的弱环配电网三相解耦潮流计算方法，包括以下步骤：（1）基于序分量法对弱环配电网进行三相潮流建模，并实现三序解耦；（2）在步骤（1）的基础上，建立弱环配电网三相解耦潮流计算方法，所述弱环配电网三相解耦潮流算法是基于道路的回路分析法。因此，该算法保留了面向支路的前推回推法的优点，而其收敛性不受回路多少的影响，具有很好的收敛性。同时，由于采用三序解耦，减少***潮流计算节点数目，降低了潮流计算量，提高了计算效率。最后，通过6母线和33母线测试算例验证了本发明方法的正确性和良好的收敛性，以及较强的回路处理能力，迭代次数随着闭合环路的增加而减少。

Description

基于序分量法的弱环配电网三相解耦潮流计算方法

技术领域

本发明属于电力***分析与计算领域，具体涉及一种基于序分量法的弱环配电网三相解耦潮流计算方法。

背景技术

配电网潮流计算是配电***分析的一项重要内容，它是对配电***规划设计和运行方式的合理性、可靠性及经济性进行定量分析的重要依据。国内外学者根据配电网的特点提出了各种配电网潮流算法，如前推回代法、隐式Zbus高斯法、回路阻抗法、改进牛顿法及快速解耦法等。其中前推回代法由于充分利用了配电网的结构特点，并且其具有物理概念明晰、编程简单、没有大矩阵计算、计算速度快、收敛性较好、非常适合于求解辐射状配电网潮流等优点而被广泛应用。但这些方法处理环网能力较弱。

相分量法对三相直接建模，物理概念清晰，易处理不对称元件和负荷。但是，对称元件三相不能解耦，导致潮流计算时节点数目较多，计算量较大，计算效率不高。同时，平衡节点和PV节点的处理也存在困难。

因此，需要一种新的弱环配电网三相潮流计算方法以解决上述问题。

发明内容

发明目的：本发明针对现有技术中配电网三相潮流计算方法存在的缺陷，提供一种基于序分量法的弱环配电网三相解耦潮流计算方法。

技术方案：为解决上述技术问题，本发明的基于序分量法的弱环配电网三相解耦潮流计算方法采用如下技术方案：

一种基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，包括以下步骤：

(1)基于序分量法对弱环配电网进行三相潮流建模，并实现三序解耦；

(2)在步骤(1)的基础上，建立弱环配电网三相解耦潮流计算方法，所述弱环配电网三相解耦潮流算法是基于道路的回路分析法。

发明原理：采用序分量法时，可以对三相对称元件实现三序解耦，从而减少***潮流计算节点数目，提高计算效率。因此基于序分量法的三相潮流建模方法在处理线路空间对称、负荷三相不平衡的配电网具有很好的效果和计算优势。可以将一组不对称的“a”、“b”、“c”三相分量分解为三组三相对称的序分量，那么，三相潮流计算就变成了计算三组三相对称序分量中的一相。因此，配电网三相不平衡潮流计算的计算量会减少2/3，在保持较好的收敛性的情况下，可以给配电网三相潮流计算带来更快的计算速度。本发明针对现有配电网三相潮流计算方法存在的一些不足，提出了一种基于序分量法的弱环配电网三相解耦潮流计算方法，该方法首先基于序分量法进行配电***三相潮流建模，实现三序解耦。

有益效果：本发明针对现有配电网三相潮流计算方法存在的一些不足，提出了一种基于序分量法的弱环配电网三相解耦潮流计算方法，该方法首先基于序分量法进行配电***三相潮流建模，实现三序解耦。接着，在配电网三序解耦的情况下，基于道路的回路分析法提出一种改进的弱环配电网三序解耦潮流算法。该算法保留了面向支路的前推回推法的优点，而其收敛性不受回路多少的影响，具有很好的收敛性。同时，由于采用三序解耦，减少***潮流计算节点数目，降低了潮流计算量，提高了计算效率。最后，通过6母线和33母线测试算例验证了本发明方法的正确性和良好的收敛性，以及较强的回路处理能力，迭代次数随着闭合环路的增加而减少。

附图说明

图1为本发明流程结构示意图；

图2为6母线三相不平衡弱环配电网***；

图3为33母线三相不平衡弱环配电网***；

图4为算例2在表4中的7种运行方式下各节点C相电压幅值分布图。

具体实施方式

下面结合附图、具体实施方式和算例，进一步阐明本发明，应理解这些算例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

请参阅图1所示，本发明的基于序分量法的弱环配电网三相解耦潮流计算方法，包括以下步骤：(1)基于序分量法对弱环配电网进行三相潮流建模，并实现三序解耦；(2)在步骤(1)的基础上，建立弱环配电网三相解耦潮流计算方法，弱环配电网三相解耦潮流算法基于道路的回路分析法。

其中，步骤(1)弱环配电网三相潮流建模包括以下步骤：

A、确定弱环配电网络的树，树具有节点、树支和连支，分别给节点、树支和连支进行编号；

B、确定弱环配电网拓扑结构参数，弱环配电网拓扑结构参数包括节点数N、支路数b和独立回路数l；根节点是电源且作为参考节点，独立节点个数为n＝N-1，支路数b＝n+l，其中，弱环配电网络的树包括树支和连支，树支条数与独立节点个数相等，连支条数与独立回路数相等，即树支n条，连支l条，支路数等于树支条数与连支条数的和。

C、获取网络参数，设定参考节点三相电压

各节点三相电压相量矩阵为

在弱环配电***三序网络中，电源节点(即参考节点)的三序电压矩阵为

各节点三序电压矩阵为

其中，a＝e^j2π/3， $A=\frac{1}{3}\left[\begin{array}{ccc}1& 1& 1\\ 1& a& a^2\\ 1& a^2& a\end{array}\right],$ $A^{-1}=\left[\begin{array}{ccc}1& 1& 1\\ 1& a^2& a\\ 1& a& a^2\end{array}\right].$

D、计算三序网络参数

三序网络参数

为基于支路i的序阻抗

形成的对角阵(n×n阶)，其中，上标s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；支路i的三相阻抗为则有

其中， $Z_{\mathrm{bi}}^{\mathrm{0,1,2}}=\left[\begin{array}{ccc}{Z}_{\mathrm{bi}}^0& 0& 0\\ 0& Z_{\mathrm{bi}}^1& 0\\ 0& 0& Z_{\mathrm{bi}}^2\end{array}\right],$ $Z_{\mathrm{bi}}^{a,b,c}=\left[\begin{array}{ccc}{Z}_{\mathrm{iaa}}& Z_{\mathrm{iab}}& Z_{\mathrm{iac}}\\ Z_{\mathrm{iba}}& Z_{\mathrm{ibb}}& Z_{\mathrm{ibc}}\\ Z_{\mathrm{ica}}& Z_{\mathrm{icb}}& Z_{\mathrm{icc}}\end{array}\right].$

步骤(2)弱环配电网三相解耦潮流计算方法包括以下步骤：

E、计算三序网络道路矩阵T_s(n×b阶)和回-支关联矩阵B_s(l×b阶)，另对于注入序电流为零的节点，在三序网络道路矩阵T_s中把该节点所对应行删去后形成新矩阵为T_sg(m×b阶)，其中，下标s＝0、1、2，分别表示三序网络模型中的零序、正序和负序网络模型，m为节点注入序电流不为零的节点个数。

F、计算回路序阻抗矩阵

以及其逆矩阵

其中，s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；

G、给弱环配电网各节点三相电压赋初始值

其中，E_n＝[E,E,…,E]^T，共n个E，E为3×3单位矩阵；

H、计算k次迭代时节点i注入的各相电流

其中，

是节点i各相注入功率，

是节点i各相并联导纳之和，p＝a,b,c，i＝1,2,…,m；m为节点注入序电流不为零的节点个数，k为迭代次数变量；

I、计算k次迭代时节点i注入的各序电流

i＝1,2,…,m；m为节点注入序电流不为零的节点个数，k为迭代次数变量；

J、计算k次迭代时的

其中，

为k次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵(m×1阶)，m为节点注入序电流不为零的节点个数；

为回路序阻抗矩阵的逆矩阵，即

且

s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；

K、计算k次迭代时的其中，1_n＝[1,1,…,1]^T，共n个1；s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；

L、基于逆变换计算k次迭代时节点i三相电压相量i＝1,2,…,n；

M、判断和

幅值之差是否满足收敛精度要求，若满足，则结束迭代；否则转步骤H。

请参阅图1所示，本发明的基于序分量法的弱环配电网三相解耦潮流计算方法，具体包括如下步骤：

1)确定弱环配电网络的树，给节点编号，规定树中的根节点的编号最小，设为“0”，其余节点按其离根节点的远近来编号，离根节点越远的节点编号越大。而树支的编号则规定为取两端节点编号中的大者，连支的编号依次从树支的最大编号后加1。

2)确定弱环配电网拓扑结构参数，包括节点数，支路数，独立回路数。不妨设三相弱环配电网有N个节点、b条支路和l个独立回路，假设首节点(即根节点)是电源且作为参考节点，则独立节点个数为n＝N-1，支路数b＝n+l。其中，弱环配电网络的树包括树支和连支，树支条数与独立节点个数相等，连支条数与独立回路数相等，即，树支n条，连支l条；支路数等于树支条数与连支条数的和。

3)设根节点是电源且作为参考节点，并设定电源节点三相电压

各节点三相电压相量矩阵为

在弱环配电***三序网络中，可以得出电源节点的三序电压为

各节点三序电压矩阵为其中，令a＝e^j2π/3， $A=\frac{1}{3}\left[\begin{array}{ccc}1& 1& 1\\ 1& a& a^2\\ 1& a^2& a\end{array}\right],$ $A^{-1}=\left[\begin{array}{ccc}1& 1& 1\\ 1& a^2& a\\ 1& a& a^2\end{array}\right].$

4)计算三序网络参数

为基于支路i的序阻抗

形成的对角阵(n×n阶)，其中，上标s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；支路i的三相阻抗为

则有

其中， $Z_{\mathrm{bi}}^{\mathrm{0,1,2}}=\left[\begin{array}{ccc}{Z}_{\mathrm{bi}}^0& 0& 0\\ 0& Z_{\mathrm{bi}}^1& 0\\ 0& 0& Z_{\mathrm{bi}}^2\end{array}\right],$ $Z_{\mathrm{bi}}^{a,b,c}=\left[\begin{array}{ccc}{Z}_{\mathrm{iaa}}& Z_{\mathrm{iab}}& Z_{\mathrm{iac}}\\ Z_{\mathrm{iba}}& Z_{\mathrm{ibb}}& Z_{\mathrm{ibc}}\\ Z_{\mathrm{ica}}& Z_{\mathrm{icb}}& Z_{\mathrm{icc}}\end{array}\right].$

5)计算三序网络道路矩阵T_s(n×b阶)和回-支关联矩阵B_s(l×b阶)；另对于注入序电流为零的节点，在三序网络道路矩阵T_s中把该节点所对应行删去后形成新矩阵为T_sg(m×b阶)。其中，下标s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型，m为节点注入序电流不为零的节点个数。

6)计算回路序阻抗矩阵以及其逆矩阵

其中，s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型。

7)给弱环配电网各节点三相电压赋初始值其中，E_n＝[E,E,…,E]^T，共n个E，E为3×3单位矩阵。

8)计算k次迭代时节点i注入的各相电流

其中，是节点i各相注入功率，

是节点i各相并联导纳之和，p＝a,b,c，i＝1,2,…,m；m为节点注入序电流不为零的节点个数，k为迭代次数变量。

9)计算k次迭代时节点i注入的各序电流

i＝1,2,…,m。

10)计算k次迭代时的

其中，

为k次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵(m×1阶)，m为节点注入序电流不为零的节点个数；

为回路序阻抗矩阵的逆矩阵，即且

s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型。

步骤10)的公式推导如下：

针对具有N个节点、b条支路和l个独立回路的三相弱环配电网，假设首节点(根节点)是电源且作为参考节点，则独立节点个数为n＝N-1，支路数b＝n+l。其中，树支n条，连支l条。

对于连通图中一颗选定的树，由于基本回路中仅包含一条连支，基本回路数等于连支数，基本回路-支路(下面简称“回-支”)关联矩阵B描述基本回路、树支、连支之间的联系。其中回-支关联矩阵B是一个l×b阶矩阵，假定连支支路的正方向都是从大号节点指向小号节点，基本回路的正方向与连支支路的正方向相同，如果支路k在回路j内，且二者方向相同，则B(j,k)＝1，如果支路k在回路j内，且二者方向相反，则B(j,k)＝-1，如果支路k不在回路j内，则B(j,k)＝0。

一个节点的道路是指节点沿树到根所经过的路径上的支路集合，节点的道路强调的是路径上的支路，对于一个给定的树，节点的道路是唯一的，节点的道路只由树支支路组成，用道路矩阵T描述道路。其中道路矩阵T是一个n×b阶矩阵，假定道路的正方向都是从电源点(即根节点)指向各节点，各树支支路正方向与道路正方向相同，如果支路k在道路i上，则T(i,k)＝1，反之T(i,k)＝0。道路矩阵T是一个稀疏矩阵，利用稀疏技术可以降低内存需求。

在弱环配电序网络中，定义节点注入序电流向量矩阵为

支路序电流向量矩阵为连支序电流(也即回路序电流)向量矩阵为在解耦的各序网模型电路中，可以获得各序网络的道路矩阵和回-支关联矩阵分别为T₀、T₁、T₂和B₀、B₁、B₂，并依据KCL电流定律，支路序电流

与节点注入序电流

回路序电流

满足如下等式：

${\stackrel{\·}{I}}_b^s=T_s^T{\stackrel{\·}{I}}_n^s+B_s^T{\stackrel{\·}{I}}_l^s---\left(1\right)$

其中，s＝0,1,2，分别表示序网络模型中的零序、正序和负序网络模型。上式表明支路序电流由两部分组成：一是节点注入序电流

对在道路上的支路的贡献，用

表示，二是回路序电流

对在回路上的支路的贡献，用表示，则有

${\stackrel{\·}{I}}_b^s={\stackrel{\·}{I}}_{b1}^s+{\stackrel{\·}{I}}_{b2}^s---\left(2\right)$

${\stackrel{\·}{I}}_{b1}^s=T_s^T{\stackrel{\·}{I}}_n^s---\left(3\right)$

${\stackrel{\·}{I}}_{b2}^s=B_s^T{\stackrel{\·}{I}}_l^s---\left(4\right)$

式(3)给出了

之间的关联，但是，在实际***中不是每个节点都有注入序电流，对于注入序电流为零的节点，在各序网的道路矩阵T_s中把该节点所对应行删去后形成新矩阵为T_sg，此时式(3)变为

${\stackrel{\·}{I}}_{b1}^s=T_{\mathrm{sg}}^T{\stackrel{\·}{I}}_g^s---\left(5\right)$

式(5)中

为去除注入序电流为零的节点后形成的新注入序电流矩阵(m×1阶)，m为节点注入序电流不为零的节点个数。

对任一弱环配电***序分量电路模型中，基于欧姆定律，支路特性约束为

${\stackrel{\·}{U}}_b^s=Z_b^s{\stackrel{\·}{I}}_b^s---\left(6\right)$

依据KVL电压定律，基本回路约束为

$B_s{\stackrel{\·}{U}}_b^s=0---\left(7\right)$

其中，

为配电网支路序电压矩阵(n×1阶)；为基于支路i的序阻抗

形成的对角阵(n×n阶)，s＝0,1,2，分别表示序网络模型中的零序、正序和负序网络模型。

由式(2)、(4)、(5)、(6)代入式(7)可得：

$B_s{Z}_b^s{B}_s^T{\stackrel{\·}{I}}_l^s+B_s{Z}_b^s{T}_{\mathrm{sg}}^T{\stackrel{\·}{I}}_g^s=0---\left(8\right)$

令

为回路序阻抗矩阵，其逆矩阵为

则有

${\stackrel{\·}{I}}_l^s=-Y_l^s{B}_s{Z}_b^s{T}_{\mathrm{sg}}^T{\stackrel{\·}{I}}_g^s---\left(9\right)$

${\stackrel{\·}{I}}_b^s=T_{\mathrm{sg}}^T{\stackrel{\·}{I}}_g^s-B_s^T{Y}_l^s{B}_s{Z}_b^s{T}_{\mathrm{sg}}^T{\stackrel{\·}{I}}_g^s---\left(10\right)$

${\stackrel{\·}{U}}_b^s=Z_b^s{\stackrel{\·}{I}}_b^s=(Z_b^s-Z_b^s{B}_s^T{Y}_l^s{B}_s{Z}_b^s)T_{\mathrm{sg}}^T{\stackrel{\·}{I}}_g^s---\left(11\right)$

设电源节点三相电压相量矩阵为

各节点三相电压相量矩阵为

在配电序网络中，可以得出电源节点的三序电压矩阵为

各节点三序电压矩阵为

那么，在各序网络模型中，可知任一节点与电源节点的序电压差等于从此节点开始沿着该节点的道路到达电源节点所经支路的支路序电压之和，即(设1_n＝[1,1,…,1]^T，共n个1；s＝0,1,2，分别表示序网络模型中的零序、正序和负序网络模型)：

$\mathrm{\Δ}{\stackrel{\·}{U}}_n^s=1_n{\stackrel{\·}{U}}_0^s-{\stackrel{\·}{U}}_n^s=T_s{\stackrel{\·}{U}}_b^s=T_s(Z_b^s-Z_b^s{B}_s^T{Y}_l^s{B}_s{Z}_b^s)T_{\mathrm{sg}}^T{\stackrel{\·}{I}}_g^s---\left(12\right)$

${\stackrel{\·}{U}}_n^s=1_n{\stackrel{\·}{U}}_0^s-\mathrm{\Δ}{\stackrel{\·}{U}}_n^s---\left(13\right)$

11)计算k次迭代时的

其中，1_n＝[1,1,…,1]^T，共n个1；s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型。

12)基于逆变换计算k次迭代时节点i三相电压相量

i＝1,2,…,n。

13)判断

和

幅值之差是否满足收敛精度要求。若满足，则结束迭代；否则转步骤8)。

算例分析

算例1

图2为6母线三相不平衡弱环配电网，变压器Y_n-y_n接线方式，对***进行了一些调整，线路空间对称，即线路相分量阻抗矩阵完全对称，而三相负荷不平衡，其单线图如图2所示。

下面就开环和闭环2种运行方式进行讨论，2种运行方式下，基于本文算法的潮流收敛情况如表1所示，其中收敛精度为10^-6。

表1算例1的2种运行方式及潮流收敛迭代次数

从表1可以看出，方式2的潮流计算比方式1的收敛更快，这是因为：在图1的***中，环路闭合后负荷离电源的电气距离缩短，负荷电流的变化引起负荷节点的电压波动变小，所以算法的迭代次数减少了，程序收敛更快了。

在方式1情况下，基于本发明计算方法的潮流计算结果如表2所示。

表2在方式1情况下潮流计算结果

在方式2情况下，基于本发明计算方法的潮流计算结果如表3所示。

表3在方式2情况下潮流计算结果

比较表2(方式1)和表3(方式2)中对应节点的三相电压也可以看出，测试***在闭环运行(即方式2)时，整个***的三相电压水平都得到了较大的提升，因为环路闭合后改善了整个***的功率流向，对电压有很好的提升作用，所以***的节点电压水平都得到了极大的改善，同时，也提高了算法的收敛性。这些在算例1的仿真中都已经得到了验证。

算例2

图3为33母线***，对负荷参数进行了调整，扩展为三相，但三相负荷不平衡，假定线路空间对称，即线路相分量阻抗矩阵完全对称，5条环路，其单线图如图3所示。

下面就表4所示的7种运行方式进行讨论，在7种运行方式下，各节点C相电压幅值分布如图4所示。

表4算例2的7种运行方式及潮流收敛迭代次数

由表4可以看出，Case1时，5条环路全部打开，算法的迭代次数为7次，迭代次数最多，收敛最慢。当Case7时，5条环路全部闭合，算法的迭代次数为4次，迭代次数最少，收敛最快。就算法收敛性的总体趋势来看，算法的迭代次数随着闭合环路的增加而减少。

图4为算例2在表4中的7种运行方式下各节点C相电压幅值分布图，从图4可以看出，Case1时，5条环路全部打开，节点电压分布最差，这也充分的反映了配电***负荷三相不平衡的特性，同时，也表明了表4中迭代次数最多的原因所在。随着环路的闭合，节点电压水平都得到了提高。特别是在Case7时，5条环路全部闭合，节点电压水平提高到最好。这说明：随着环路闭合的越多，由于***负荷三相不平衡给配电网带来的节点电压水平低问题都会得到更好的改善和提升，而且算法的收敛性也得到了提高。

Claims (7)

1.一种基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，包括以下步骤：

(1)基于序分量法对弱环配电网进行三相潮流建模，并实现三序解耦；

(2)在步骤(1)的基础上，建立弱环配电网三相解耦潮流计算方法，所述弱环配电网三相解耦潮流算法是基于道路的回路分析法。

2.如权利要求1所述的基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，步骤(1)所述的弱环配电网三相潮流建模包括以下步骤：

A、确定弱环配电网络的树，所述树具有节点、树支和连支，分别给所述节点、树支和连支进行编号；

B、确定弱环配电网拓扑结构参数，所述弱环配电网拓扑结构参数包括节点数N、支路数b和独立回路数l；

C、获取网络参数，设定电源节点三相电压

各节点三相电压相量矩阵为

D、计算三序网络参数

其中，上标s＝0、1、2，分别表示三序网络模型中的零序、正序和负序网络模型。

3.如权利要求2所述的基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，所述电源节点是根节点且作为参考节点，独立节点个数为n＝N-1，支路数b＝n+l，其中，弱环配电网络的树包括树支和连支，树支条数与独立节点个数相等，连支条数与独立回路数相等，即树支n条，连支l条，支路数等于树支条数与连支条数的和。

4.如权利要求2所述的基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，在弱环配电***三序网络中，电源节点的三序电压矩阵为

各节点三序电压矩阵为

其中，a＝e^j2π/3， $A=\frac{1}{3}\left[\begin{array}{ccc}1& 1& 1\\ 1& a& a^2\\ 1& a^2& a\end{array}\right],$ $A^{-1}=\left[\begin{array}{ccc}1& 1& 1\\ 1& a^2& a\\ 1& a& a^2\end{array}\right].$

5.如权利要求2所述的基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，步骤D中所述三序网络参数

为基于支路i的序阻抗

形成的对角阵(n×n阶)，其中，上标s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；支路i的三相阻抗为

则有

其中， $Z_{\mathrm{bi}}^{\mathrm{0,1,2}}=\left[\begin{array}{ccc}{Z}_{\mathrm{bi}}^0& 0& 0\\ 0& Z_{\mathrm{bi}}^1& 0\\ 0& 0& Z_{\mathrm{bi}}^2\end{array}\right],$ $Z_{\mathrm{bi}}^{a,b,c}=\left[\begin{array}{ccc}{Z}_{\mathrm{iaa}}& Z_{\mathrm{iab}}& Z_{\mathrm{iac}}\\ Z_{\mathrm{iba}}& Z_{\mathrm{ibb}}& Z_{\mathrm{ibc}}\\ Z_{\mathrm{ica}}& Z_{\mathrm{icb}}& Z_{\mathrm{icc}}\end{array}\right].$

6.如权利要求1所述的基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，步骤(2)所述的弱环配电网三相解耦潮流计算方法包括以下步骤：

E、计算三序网络道路矩阵T_s(n×b阶)和回-支关联矩阵B_s(l×b阶)，其中，下标s＝0、1、2，分别表示三序网络模型中的零序、正序和负序网络模型；

F、计算回路序阻抗矩阵以及其逆矩阵

其中，s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；

G、给弱环配电网各节点三相电压赋初始值

其中，E_n＝[E,E,…,E]^T，共n个E，E为3×3单位矩阵；

H、计算k次迭代时节点i注入的各相电流

其中，

是节点i各相注入功率，

是节点i各相并联导纳之和，p＝a,b,c，i＝1,2,…,m；m为节点注入序电流不为零的节点个数，k为迭代次数变量；

I、计算k次迭代时节点i注入的各序电流

i＝1,2,…,m；m为节点注入序电流不为零的节点个数，k为迭代次数变量；

J、计算k次迭代时的

其中，

为k次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵(m×1阶)，m为节点注入序电流不为零的节点个数；

为回路序阻抗矩阵的逆矩阵，即

且

s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；

K、计算k次迭代时的其中，1_n＝[1,1,…,1]^T，共n个1；s＝0,1,2，分别表示三序网络模型中的零序、正序和负序网络模型；

L、基于逆变换计算k次迭代时节点i三相电压相量

i＝1,2,…,n；

M、判断

和

幅值之差是否满足收敛精度要求，若满足，则结束迭代；否则转步骤H。

7.如权利要求1所述的基于序分量法的弱环配电网三相解耦潮流计算方法，其特征在于，步骤J中，对于注入序电流为零的节点，在三序网络道路矩阵T_s中把该节点所对应行删去后形成新矩阵为T_sg(m×b阶)，其中，下标s＝0、1、2，分别表示三序网络模型中的零序、正序和负序网络模型，m为节点注入序电流不为零的节点个数。

Images