CN102593823B - 一种基于叠加原理的配电网在线潮流计算方法 - Google Patents

Abstract

本发明公开了一种基于叠加原理的配电网在线潮流计算方法。该方法利用离线潮流计算结果，事先计算得出各负荷节点单位负荷电流在各支路引起的电压降/电流向量；进行在线潮流计算时，根据负荷节点负荷电流以及单位负荷电流引起的支路电压降/电流向量，利用线性系统的叠加原理计算支路电压降/电流向量，进而计算出节点电压、支路电流等。本发明可以很好的解决配电网在线潮流计算问题，无需对节点进行重新编号，直接利用离线潮流计算的节点编号即可。本发明对开环、闭环潮流计算均可适用，无需对环网潮流进行特殊处理；既可用于单相潮流计算，也可用于三相潮流计算；具有计算速度快、占用内存少、稳定性好的特点。

Description

一种基于叠加原理的配电网在线潮流计算方法

技术领域：

本发明涉及电力系统运行与安全分析技术领域，具体涉及一种基于叠加原理的配电网在线潮流计算方法。

背景技术：

潮流计算是电力系统中应用最为广泛、最基本和最重要的电气计算，它的任务是根据给定的网络结构及运行条件，计算整个网络的潮流分布。

潮流计算的结果，无论是对于现有系统运行方式的分析研究，还是对规划阶段设计方案的分析比较，都是必不可少的。电力系统静态、暂态稳定计算也需要利用潮流计算的结果，这些都属于离线计算的范畴。

随着现代化的调度控制中心的建立，为了对电力系统进行实时安全监控，需要根据实时数据库提供的信息，判断系统当前的运行状态并对预想事故进行安全分析，这就需要进行在线潮流计算。

配电网具有以下特点：网络拓扑一般呈辐射状、但也有短时环网运行状态，R/X比值较大，支路数和节点数十分庞大，三相不平衡等。这些特点导致网络雅克比矩阵出现不同程度的病态特征，传统的牛顿拉夫逊法、快速解耦法等算法对配电网不再适用。专家学者们研究提出了适用于配电网特点的潮流计算方法，如：隐式高斯法、前推回代法、回路阻抗法等。这些算法可直接利用节点及支路参数，无需利用雅克比矩阵，编程简单，算法收敛性好，但都要用到复杂的节点编号方法，节点编号对算法实现有很大影响。

配电网可以看作是线性网络，线性系统的叠加原理对配电网同样适用。部分专家学者在处理弱环配电网潮流计算中，提出将环网运行的配电网转化为辐射型配电网，利用叠加原理，将环网潮流视为环网前的辐射网潮流与环网两端电压向量差引起的潮流的叠加。

配电网各个支路的电压降和支路通过的电流也适用叠加原理，可以看作是各个负荷电流共同作用的结果。当网络结构和运行方式维持恒定时，单位负荷电流在各个支路引起的电压降以及电流向量是恒定的，可以利用离线潮流计算结果计算得出。

在线潮流计算主要用于运行中电力系统的监视和实时控制，除了可靠的收敛特性外，计算速度快、占用内存少也是最基本的要求。常规的潮流计算方法，不能很好的满足配电网在线潮流计算对计算速度和收敛性的要求。

发明内容：

本发明提出一种基于叠加原理的配电网在线潮流计算方法。该方法利用离线潮流计算结果，事先计算得出各负荷节点单位负荷电流在各支路引起的电压降/电流向量；进行在线潮流计算时，根据负荷节点负荷电流以及单位负荷电流引起的支路电压降/电流向量，利用线性系统的叠加原理计算支路电压降/电流向量，进而计算出节点电压、支路电流等。

该方法的技术原理如下：

对于一个实际配电网，当网络结构和运行方式一定时，各支路电压降向量适用叠加原理，即：一个支路的电压降向量等于各个负荷节点在该支路引起的电压降向量之和，某负荷节点负荷电流在一个支路引起的电压降向量与该负荷电流成正比。利用离线潮流计算结果，计算出各负荷节点单位负荷电流在各支路引起的电压降向量；进行在线潮流计算时，根据各负荷节点负荷电流及上述计算出的各负荷节点单位负荷电流在各支路引起的电压降向量，利用线性系统的叠加原理计算出各支路电压降向量，进而计算出各节点电压、各支路损耗、各节点/支路功率等。

上述求解方法对支路电流的求解同样适用。对于一个实际配电网，当网络结构和运行方式一定时，各支路电流向量适用叠加原理，即：一个支路的电流向量等于各个负荷节点在该支路引起的电流向量之和，某负荷节点负荷电流在一个支路引起的电流向量与该负荷电流成正比。利用离线潮流计算结果，计算出各负荷节点单位负荷电流在各支路引起的电流向量；进行在线潮流计算时，根据各负荷节点负荷电流及上述计算出的各负荷节点单位负荷电流在各支路引起的电流向量，利用线性系统的叠加原理计算出各支路通过的电流向量，进而计算出各支路电压降、各节点电压、各支路损耗、各节点/支路功率等。

本发明提供的一种基于叠加原理的配电网在线潮流计算方法，其改进之处在于，利用离线潮流计算的结果和线性系统的叠加原理进行在线潮流计算；

配电网的网络结构和运行方式一定时，利用离线潮流计算结果，事先计算出各负荷节点单位负荷电流在各支路引起的电压降向量，进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电压降向量，利用线性系统的叠加原理计算出各支路电压降向量，进而计算出各节点电压、各支路损耗和各节点/支路功率；或者利用离线潮流计算结果，事先计算出各负荷节点单位负荷电流在各支路引起的电流向量，进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电流向量，利用线性系统的叠加原理计算出各支路通过的电流向量，进而计算出各支路电压降、各节点电压、各支路损耗和各节点/支路功率。

其中，所述利用离线潮流计算结果，事先计算出各负荷节点单位负荷电流在各支路引起的电压降向量；进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电压降向量，利用线性系统的叠加原理计算出各支路电压降向量，进而计算出各节点电压、各支路损耗和各节点/支路功率，包括以下步骤：

(1)依据离线潮流计算结果，列写电压降落矩阵：

$\mathrm{\ΔU}=\left[\begin{array}{cccccc}{\mathrm{\ΔU}}_{11}& {\mathrm{\ΔU}}_{12}& ...& {\mathrm{\ΔU}}_{1j}& ...& {\mathrm{\ΔU}}_{1n}\\ {\mathrm{\ΔU}}_{21}& {\mathrm{\ΔU}}_{22}& ...& {\mathrm{\ΔU}}_{2j}& ...& {\mathrm{\ΔU}}_{2n}\\ .& .& & .& & .\\ .& .& ...& .& ...& .\\ .& .& & .& & .\\ {\mathrm{\ΔU}}_{i1}& {\mathrm{\ΔU}}_{i2}& ...& {\mathrm{\ΔU}}_{\mathrm{ij}}& ...& {\mathrm{\ΔU}}_{\mathrm{in}}\\ .& .& & .& & .\\ .& .& ...& .& ...& .\\ .& .& & .& & .\\ {\mathrm{\ΔU}}_{m1}& {\mathrm{\ΔU}}_{m2}& ...& {\mathrm{\ΔU}}_{\mathrm{mj}}& ...& {\mathrm{\ΔU}}_{\mathrm{mn}}\end{array}\right]$

式中：ΔU为m*n阶矩阵；m为负荷节点数；n为支路数，ΔU_ij为在节点i注入单位电流，其余节点注入电流为0时，支路j的电压降向量。

(2)根据各负荷节点负荷电流，计算各支路电压降向量：

${\mathrm{\ΔU}}_j=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{I}_i*{\mathrm{\ΔU}}_{\mathrm{ij}}$

式中：ΔU_j为支路j的电压降向量；I_i为负荷节点i的负荷电流。

(3)从馈线首端至末端，依次计算各个节点电压：

U_i＝U_k-ΔU_j

式中：k，i分别为支路j的首端节点和末端节点，U_k为首端节点k的电压，U_i为末端节点i的电压。

(4)计算各支路电流：

${\mathrm{\ΔI}}_j=\frac{{\mathrm{\ΔU}}_j}{Z_j}$

式中：ΔI_j表示通过支路j的电流向量，Z_j表示支路j的阻抗参数。

(5)计算功率分布及功率损耗，输出结果。

其中，利用离线潮流计算结果，事先计算出各负荷节点单位负荷电流在各支路引起的电流向量；进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电流向量，利用线性系统的叠加原理计算出各支路通过的电流向量，进而计算出各支路电压降、各节点电压、各支路损耗和各节点/支路功率，包括如下步骤：

(1)依据离线潮流计算结果，列写电流矩阵：

$\mathrm{\ΔI}=\left[\begin{array}{cccccc}{\mathrm{\ΔI}}_{11}& {\mathrm{\ΔI}}_{12}& ...& {\mathrm{\ΔI}}_{1j}& ...& {\mathrm{\ΔI}}_{1n}\\ {\mathrm{\ΔI}}_{21}& {\mathrm{\ΔI}}_{22}& ...& {\mathrm{\ΔI}}_{2j}& ...& {\mathrm{\ΔI}}_{2n}\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ {\mathrm{\ΔI}}_{i1}& {\mathrm{\ΔI}}_{i2}& ...& {\mathrm{\ΔI}}_{\mathrm{ij}}& ...& {\mathrm{\ΔI}}_{\mathrm{in}}\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ {\mathrm{\ΔI}}_{m1}& {\mathrm{\ΔI}}_{m2}& ...& {\mathrm{\ΔI}}_{\mathrm{mj}}& ...& {\mathrm{\ΔI}}_{\mathrm{mn}}\end{array}\right]$

式中：ΔI为m*n阶矩阵；m为负荷节点数；n为支路数；ΔI_ij表示在负荷节点i注入单位电流，其余负荷节点注入电流为0时，支路j通过的电流向量。

(2)根据各负荷节点负荷电流，计算各支路电流向量：

${\mathrm{\ΔI}}_j=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{I}_i*{\mathrm{\ΔI}}_{\mathrm{ij}}$

式中：ΔI_j为支路j的电流向量；I_i为负荷节点i的负荷电流。

(3)计算各支路电压降向量：

ΔU_j＝ΔI_j*Z_j

式中：ΔU_j为支路j的电压降向量；Z_j为支路j的阻抗参数。

(4)从馈线首端至末端，依次计算各个节点电压：

U_i＝U_k-ΔU_j

式中：U_k表示支路j的首端节点k的电压，U_i表示支路j的末端节点i的电压。

(5)计算功率分布及功率损耗，输出结果。

其中，当网络结构和运行方式发生变化时，采用如下方式对ΔU进行修正：

(a)对于闭环设计、开环运行的网络，增加一条支路时ΔU的修正：

从原有负荷节点i引出一条支路，同时增加一个负荷节点，编号n+1，由于负荷节点数增加，支路数增加，电压降落矩阵ΔU增加一行一列：(矩阵中(m+1)行和(n+1)列为增加的一行一列)

式中：ΔU′为修正后的电压降落矩阵，为(m+1)*(n+1)维。

对于j＝1、2…n时：

ΔU′_j＝ΔU_j+I_(m+1)*ΔU_(m+1)

式中：ΔU′_j为修正后的支路j的电压降向量；

对于j＝n+1时：

${{\mathrm{\ΔU}}_{n+1}}^{\′}=\underset{i=1}{\overset{m+1}{\mathrm{\Σ}}}{I}_i*{\mathrm{\ΔU}}_{i(n+1)}$

(b)对于闭环设计、开环运行的网络，切除一条支路时ΔU的修正：

将支路k断开，同时减少一个负荷节点k，由于负荷节点数减少，支路数减少，电压降落矩阵ΔU减少一行一列：(矩阵中虚线代表减少的一行一列)

式中：ΔU′为修正后的电压降落矩阵，为(m-1)*(n-1)维。

对于j＝1、2…n，(j≠k)时：

ΔU′_j＝ΔU_j-I_k*ΔU_kj

(c)对于多电源供电网络，切除其中一个电源时ΔU的修正：

可用一个与该电源出口电流大小相等、方向相反的电流源来等效，即可看作增加一条支路、增加一个负荷节点：

对于j＝1、2…n时：

ΔU′_j＝ΔU_j+I_(m+1)*ΔU_(m+1)j

式中：I_(m+1)与被切除电源原有电流方向相反。

其中，所述配电网在线潮流计算方法用于单相潮流计算或三相潮流计算。

当利用离线潮流计算结果，事先计算出各负荷节点单位负荷电流在各支路引起的电流向量；进行在线潮流计算时，当网络结构和运行方式发生变化，参考对ΔU进行修正的方法对ΔI进行修正。

本发明中提到的负荷电流是指有功电流，对无功电流同样适用，如果为复电流向量，则功率因数应维持恒定值。

本发明可以很好的解决配电网在线潮流计算问题，与现有技术比，本发明的有益效果为：

1、无需对节点进行重新编号，直接利用离线潮流计算的节点编号即可；

2、对开环、闭环潮流计算均可适用，无需对环网潮流进行特殊处理；

3、既可用于单相潮流计算，也可用于三相潮流计算；

4、计算速度快、占用内存少、稳定性好。

附图说明：

图1为本发明的总体流程图。

图2为本发明提供的根据单位负荷电流引起的支路电压降进行在线潮流计算的流程图。

图3为本发明提供的根据单位负荷电流引起的支路电流进行在线潮流计算的流程图。

具体实施方式：

下面结合附图对本发明的具体实施方式作进一步的详细说明。

配电网可以看作是线性网络，线性系统的叠加原理对配电网同样适用。对于一个实际配电网，当网络结构和运行方式一定时，各支路电压降向量适用叠加原理，即：一个支路的电压降向量等于各个负荷节点负荷在该支路引起的电压降向量之和，某负荷节点负荷在一个支路引起的电压降向量与该负荷电流成正比。

同理，对于一个实际配电网，当网络结构和运行方式一定时，各支路电流向量适用叠加原理，即：一个支路的电流向量等于各个负荷节点负荷在该支路引起的电流向量之和，某负荷节点负荷在一个支路引起的电流向量与该负荷电流成正比。

当网络结构和运行方式维持恒定时，单位负荷电流在各个支路引起的电压降以及电流向量是恒定的，可以利用离线潮流计算结果事先计算得出。进行在线潮流计算时，可直接利用上述计算结果和负荷电流，应用叠加原理进行求解。本发明的总体思路如图1所示。

当节点负荷电流已知时，根据单位负荷电流在各支路引起的电压进行在线潮流计算流程如图2所示，具体包括如下步骤：

(1)依据离线潮流计算结果，列写电压降落矩阵：

$\mathrm{\ΔU}=\left[\begin{array}{cccccc}{\mathrm{\ΔU}}_{11}& {\mathrm{\ΔU}}_{12}& ...& {\mathrm{\ΔU}}_{1j}& ...& {\mathrm{\ΔU}}_{1n}\\ {\mathrm{\ΔU}}_{21}& {\mathrm{\ΔU}}_{22}& ...& {\mathrm{\ΔU}}_{2j}& ...& {\mathrm{\ΔU}}_{2n}\\ .& .& & .& & .\\ .& .& ...& .& ...& .\\ .& .& & .& & .\\ {\mathrm{\ΔU}}_{i1}& {\mathrm{\ΔU}}_{i2}& ...& {\mathrm{\ΔU}}_{\mathrm{ij}}& ...& {\mathrm{\ΔU}}_{\mathrm{in}}\\ .& .& & .& & .\\ .& .& ...& .& ...& .\\ .& .& & .& & .\\ {\mathrm{\ΔU}}_{m1}& {\mathrm{\ΔU}}_{m2}& ...& {\mathrm{\ΔU}}_{\mathrm{mj}}& ...& {\mathrm{\ΔU}}_{\mathrm{mn}}\end{array}\right]$

式中：ΔU为m*n阶矩阵；m为负荷节点数；n为支路数，ΔU_ij为在节点i注入单位电流，其余节点注入电流为0时，支路j的电压降向量。

(2)根据各负荷节点负荷电流，计算各支路电压降向量：

${\mathrm{\ΔU}}_j=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{I}_i*{\mathrm{\ΔU}}_{\mathrm{ij}}$

式中：ΔU_j为支路j的电压降向量；I_i为负荷节点i的负荷电流。

(3)从馈线首端至末端，依次计算各个节点电压：

U_i＝U_k-ΔU_j

式中：k，i分别为支路j的首端节点和末端节点，U_k为首端节点k的电压，U_i为末端节点i的电压。

(4)计算各支路电流：

${\mathrm{\ΔI}}_j=\frac{{\mathrm{\ΔU}}_j}{Z_j}$

式中：ΔI_j表示通过支路j的电流向量，Z_j表示支路j的阻抗参数。

(5)计算功率分布及功率损耗，输出结果。

对于支路j，功率损耗可用下述公式计算：

ΔP_j＝ΔU_jΔI_j

支路j的首端节点k和末端节点i的功率分别为：

P_k＝U_kΔI_j

P_i＝U_iΔI_j

当节点负荷电流已知时，根据单位负荷电流在各支路引起的电流进行在线潮流计算，计算流程如图3所示，具体步骤如下：

(I)依据离线潮流计算结果，列写电流矩阵：

$\mathrm{\ΔI}=\left[\begin{array}{cccccc}{\mathrm{\ΔI}}_{11}& {\mathrm{\ΔI}}_{12}& ...& {\mathrm{\ΔI}}_{1j}& ...& {\mathrm{\ΔI}}_{1n}\\ {\mathrm{\ΔI}}_{21}& {\mathrm{\ΔI}}_{22}& ...& {\mathrm{\ΔI}}_{2j}& ...& {\mathrm{\ΔI}}_{2n}\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ {\mathrm{\ΔI}}_{i1}& {\mathrm{\ΔI}}_{i2}& ...& {\mathrm{\ΔI}}_{\mathrm{ij}}& ...& {\mathrm{\ΔI}}_{\mathrm{in}}\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ {\mathrm{\ΔI}}_{m1}& {\mathrm{\ΔI}}_{m2}& ...& {\mathrm{\ΔI}}_{\mathrm{mj}}& ...& {\mathrm{\ΔI}}_{\mathrm{mn}}\end{array}\right]$

式中：ΔI为m*n阶矩阵；m为负荷节点数；n为支路数；ΔI_ij表示在负荷节点i注入单位电流，其余负荷节点注入电流为0时，支路j通过的电流向量。

(II)根据各负荷节点负荷电流，计算各支路电流向量：

${\mathrm{\ΔI}}_j=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{I}_i*{\mathrm{\ΔI}}_{\mathrm{ij}}$

式中：ΔI_j为支路j的电流向量；I_i为负荷节点i的负荷电流。

(III)计算各支路电压降向量：

ΔU_j＝ΔI_j*Z_j

式中：ΔU_j为支路j的电压降向量；Z_j为支路j的阻抗参数。

(IV)从馈线首端至末端，依次计算各个节点电压：

U_i＝U_k-ΔU_j

式中：U_k表示支路j的首端节点k的电压，U_i表示支路j的末端节点i的电压。

(V)计算功率分布及功率损耗，输出结果。

该方法既可用于单相潮流计算，也可用于三相潮流计算；应用本发明对三相分别进行计算，即可得到配电网三相潮流计算结果。

当网络结构和运行方式发生变化时，ΔU和ΔI也发生变化，可以事先计算得出各种不同运行方式下的ΔU和ΔI，也可对现有ΔU和ΔI进行修正。当根据单位负荷电流引起的支路电压降进行在线潮流计算时，本实施例选用下述方法对ΔU进行修正：

(a)对于闭环设计、开环运行的网络，增加一条支路时ΔU的修正：

从原有负荷节点i引出一条支路，同时增加一个负荷节点，编号n+1，由于负荷节点数增加，支路数增加，电压降落矩阵ΔU增加一行一列：

式中：ΔU′为修正后的电压降落矩阵，为(m+1)*(n+1)维。

对于j＝1、2…n时：

ΔU′_j＝ΔU_j+I_(m+1)*ΔU_(m+1)

式中：ΔU′_j为修正后的支路j的电压降向量；

对于j＝n+1时：

${{\mathrm{\ΔU}}_{n+1}}^{\′}=\underset{i=1}{\overset{m+1}{\mathrm{\Σ}}}{I}_i*{\mathrm{\ΔU}}_{i(n+1)}$

(b)对于闭环设计、开环运行的网络，切除一条支路时ΔU的修正：

将支路k断开，同时减少一个负荷节点k，由于负荷节点数减少，支路数减少，电压降落矩阵ΔU减少一行一列：

式中：ΔU′为修正后的电压降落矩阵，为(m-1)*(n-1)维。

对于j＝1、2…n，(j≠k)时：

ΔU′_j＝ΔU_j-I_k*ΔU_kj

(c)对于多电源供电网络，切除其中一个电源时ΔU的修正：

可用一个与该电源出口电流大小相等、方向相反的电流源来等效，即可看作增加一条支路、增加一个负荷节点：

对于j＝1、2…n时：

ΔU′_j＝ΔU_j+I_(m+1)*ΔU_(m+1)j

式中：I_(m+1)与被切除电源原有电流方向相反。

本实施例中提到的负荷电流是指有功电流，对无功电流同样适用，如果为复电流向量，则功率因数应维持恒定值。

上述说明均是基于负荷电流已知的情况下进行的，对于负荷节点已知量为有功功率、无功功率的情形，也可用该方法进行计算，计算步骤类似，负荷电流可利用功率计算得出：

$I_i={\left(\frac{S_i}{U_i}\right)}^{*}$

式中：S_i为节点i的功率，U_i为节点i的电压(可取根节点电压)，I_i为节点i的负荷电流。

最后应该说明的是：结合上述实施例仅说明本发明的技术方案而非对其限制。所属领域的普通技术人员应当理解到：本领域技术人员可以对本发明的具体实施方式进行修改或者等同替换，但这些修改或变更均在申请待批的权利要求保护范围之中。

Claims (5)

1.一种基于叠加原理的配电网在线潮流计算方法，其特征在于，利用离线潮流计算的结果和线性系统的叠加原理进行在线潮流计算；

配电网的网络结构和运行方式一定时,利用离线潮流计算结果，先计算出各负荷节点单位负荷电流在各支路引起的电压降向量，进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电压降向量，利用线性系统的叠加原理计算出各支路电压降向量，进而计算出各节点电压、各支路损耗和各节点/支路功率；或者利用离线潮流计算结果，先计算出各负荷节点单位负荷电流在各支路引起的电流向量，进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电流向量，利用线性系统的叠加原理计算出各支路通过的电流向量，进而计算出各支路电压降、各节点电压、各支路损耗和各节点/支路功率。

2.如权利要求1所述的方法，其特征在于，所述利用离线潮流计算结果，先计算出各负荷节点单位负荷电流在各支路引起的电压降向量；进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电压降向量，利用线性系统的叠加原理计算出各支路电压降向量，进而计算出各节点电压、各支路损耗和各节点/支路功率，包括以下步骤：

（1）依据离线潮流计算结果，列写电压降落矩阵：

$\mathrm{\ΔU}=\left[\begin{array}{cccccc}\mathrm{\Δ}{U}_{11}& \mathrm{\Δ}{U}_{12}& ...& \mathrm{\Δ}{U}_{1j}& ...& \mathrm{\Δ}{U}_{1n}\\ \mathrm{\Δ}{U}_{21}& \mathrm{\Δ}{U}_{22}& ...& \mathrm{\Δ}{U}_{2j}& ...& \mathrm{\Δ}{U}_{2n}\\ .& .& & .& & .\\ .& .& ...& .& ...& .\\ .& .& & .& & .\\ \mathrm{\Δ}{U}_{i1}& \mathrm{\Δ}{U}_{i2}& ...& \mathrm{\Δ}{U}_{\mathrm{ij}}& ...& \mathrm{\Δ}{U}_{\mathrm{in}}\\ .& .& & .& & .\\ .& .& ...& .& ...& .\\ .& .& & .& & .\\ \mathrm{\Δ}{U}_{m1}& \mathrm{\Δ}{U}_{m2}& ...& \mathrm{\Δ}{U}_{\mathrm{mj}}& ...& \mathrm{\Δ}{U}_{\mathrm{mn}}\end{array}\right]$

式中：ΔU为m*n阶矩阵；m为负荷节点数；n为支路数，ΔU_ij为在节点i注入单位电流，其余节点注入电流为0时，支路j的电压降向量；

（2）根据各负荷节点负荷电流，计算各支路电压降向量：

$\mathrm{\Δ}{U}_j=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{I}_i*\mathrm{\Δ}{U}_{\mathrm{ij}}$

式中：ΔU_j为支路j的电压降向量；I_i为负荷节点i的负荷电流；

（3）从馈线首端至末端，依次计算各个节点电压：

U_i＝U_k-ΔU_j

式中：k，i分别为支路j的首端节点和末端节点，U_k为首端节点k的电压，U_i为末端节点i的电压；

（4）计算各支路电流：

$\mathrm{\Δ}{I}_j=\frac{\mathrm{\Δ}{U}_j}{Z_j}$

式中：ΔI_j表示通过支路j的电流向量，Z_j表示支路j的阻抗参数；

（5）计算功率分布及功率损耗，输出结果。

3.如权利要求1所述的方法，其特征在于，利用离线潮流计算结果，先计算出各负荷节点单位负荷电流在各支路引起的电流向量；进行在线潮流计算时，根据各负荷节点负荷电流和所述计算出的各负荷节点单位负荷电流在各支路引起的电流向量，利用线性系统的叠加原理计算出各支路通过的电流向量，进而计算出各支路电压降、各节点电压、各支路损耗和各节点/支路功率，包括如下步骤：

（1）依据离线潮流计算结果，列写电流矩阵：

$\mathrm{\ΔI}=\left[\begin{array}{cccccc}\mathrm{\Δ}{I}_{11}& \mathrm{\Δ}{I}_{12}& ...& \mathrm{\Δ}{I}_{1j}& ...& \mathrm{\Δ}{I}_{1n}\\ \mathrm{\Δ}{I}_{21}& \mathrm{\Δ}{I}_{22}& ...& \mathrm{\Δ}{I}_{2j}& ...& \mathrm{\Δ}{I}_{2n}\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ \mathrm{\Δ}{I}_{i1}& \mathrm{\Δ}{I}_{i2}& ...& \mathrm{\Δ}{I}_{\mathrm{ij}}& ...& \mathrm{\Δ}{I}_{\mathrm{in}}\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ .& .& .& .& .& .\\ \mathrm{\Δ}{I}_{m1}& \mathrm{\Δ}{I}_{m2}& ...& \mathrm{\Δ}{I}_{\mathrm{mj}}& ...& \mathrm{\Δ}{I}_{\mathrm{mn}}\end{array}\right]$

式中：ΔI为m*n阶矩阵；m为负荷节点数；n为支路数；ΔI_ij表示在负荷节点i注入单位电流，其余负荷节点注入电流为0时，支路j通过的电流向量；

（2）根据各负荷节点负荷电流，计算各支路电流向量：

$\mathrm{\Δ}{I}_j=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{I}_i*\mathrm{\Δ}{I}_{\mathrm{ij}}$

式中：ΔI_j为支路j的电流向量；I_i为负荷节点i的负荷电流；

（3）计算各支路电压降向量：

ΔU_j＝ΔI_j*Z_j

式中：ΔU_j为支路j的电压降向量；Z_j为支路j的阻抗参数；

（4）从馈线首端至末端，依次计算各个节点电压：

U_i＝U_k-ΔU_j

式中：U_k表示支路j的首端节点k的电压，U_i表示支路j的末端节点i的电压；

（5）计算功率分布及功率损耗，输出结果。

4.如权利要求1所述的方法，其特征在于，当网络结构和运行方式发生变化时，采用下述方式对ΔU进行修正：

（a）对于闭环设计、开环运行的网络，增加一条支路时ΔU的修正：

从原有负荷节点i引出一条支路，同时增加一个负荷节点，编号n+1，电压降落矩阵ΔU增加一行一列：

式中：ΔU′为修正后的电压降落矩阵，为（m+1）*(n+1)维；

对于j=1、2…n时：

ΔU'_j＝ΔU_j+I_(m+1)*ΔU_(m+1)

式中：ΔU'_j为修正后的支路j的电压降向量；

对于j=n+1时：

$\mathrm{\Δ}{U_{n+1}}^{\′}=\underset{i=1}{\overset{m+1}{\mathrm{\Σ}}}{I}_i*\mathrm{\Δ}{U}_{i(n+1)}$

（b）对于闭环设计、开环运行的网络，切除一条支路时ΔU的修正：

将支路k断开，同时减少一个负荷节点k，电压降落矩阵ΔU减少一行一列：

式中：ΔU′为修正后的电压降落矩阵，为（m-1）*（n-1）维；

对于j=1、2…n，（j≠k）时：

ΔU'_j＝ΔU_j-I_k*ΔU_kj

（c）对于多电源供电网络，切除其中一个电源时ΔU的修正：

用一个与该电源出口电流大小相等、方向相反的电流源等效，看作增加一条支路、增加一个负荷节点：

对于j=1、2…n时：

ΔU'_j＝ΔU_j+I_(m+1)*ΔU_(m+1)j

式中：I_(m+1)与被切除电源原有电流方向相反。

5.如权利要求1所述的方法，其特征在于，所述配电网在线潮流计算方法用于单相潮流计算或三相潮流计算。

Images