[发明专利]动态拓扑电网中直流潮流的分析方法

说明书

技术领域

动态拓扑电网中直流潮流的分析方法属于电网分析和电网安全技术领域。

在大规模电网中，当某些传输线因发生故障而断开时，电网的拓扑就会发生变化，从而引起整个网络中的潮流发生重新分配。如果潮流重新分配以后部分传输线上潮流过载，那些线路也会断开，从而进一步引起电网拓扑的变化。如此下去就有可能造成电网的连锁故障，可能导致大面积停电事故。对于可能发生传输线故障的电网来说，其拓扑结构是动态变化的，因此动态拓扑下的潮流分析对于防止电网连锁故障具有重要意义。

背景技术

直流潮流分析仅考虑潮流中的有功功率部分，与交流潮流分析相比具有分析速度快的优势，经常用于对精度要求不很高的潮流分析中。传统的直流潮流分析方法是针对连通电网设计的，直流潮流分析时需要在电网中选取一个参考节点并规定其相角取值为零，然后才能求解其它节点的相角和传输线上的直流潮流。如果电网并非连通的而是包含若干个孤岛，传统的直流潮流分析方法必须首先通过拓扑分析找出同属于一个孤岛的节点，然后将每个孤岛当作一个连通系统分别进行求解。当系统具有动态拓扑时，系统连通与否以及系统中包含的孤岛数目都会造成直流潮流表达式形式及分析过程的改变。因此，传统的直流潮流分析方法不适合应用于动态拓扑系统的潮流分析。

本发明是针对传统直流潮流分析方法应用于动态拓扑系统时的困难提出的一种新的直流潮流分析方法。该方法在计算前不需要对系统的拓扑结构进行分析，也不需要在系统中选取参考节点，采用广义逆求出直流潮流结果。可以证明，该方法对于连通与不连通的系统都适用，且不论系统是否连通及系统中包含多少个孤岛，该方法的直流潮流表达形式和分析过程都是相同的。如果将系统中部分传输线的连通与否作为二值变量，该方法求得的带有变量的潮流结果将是系统拓扑结构的函数，可以反映不同的系统拓扑结构对直流潮流的影响。因此，相对于传统的直流潮流分析方法，本发明在应用于动态拓扑系统的潮流分析问题时具有明显的优势。

发明内容

本发明提供了一种动态拓扑电网中直流潮流的分析方法。该方法在计算时不需要对电网拓扑的连通性进行分析，不需要在电网中选取参考节点，采用广义逆求解潮流方程。

对于一个包含n个网络节点和m条传输线的电网，所述方法是在计算机中依次按以下步骤实现的：

步骤(1)，向计算机输入一个关联矩阵D，该关联矩阵的维数为m×n，n为母线节点数，m为所述传输线的条数，每一列对应于所述电网中的每一个已经编好顺序的节点，每一行对应于所述电网中已经编好顺序的两个节点之间的一条传输线，元素值为1表示该节点为该传输线的起点，元素值为0表示该节点为该传输线的终点，从而得到传输线的参考方向，其它元素的值均为0，因而，该关联矩阵D中的列数表示所述电网中节点的总数，行数表示所述电网中传输线的总数；

步骤(2)，向计算机输入一个传输线导纳矩阵C，该传输线导纳矩阵C的维数为m×m，是一个对角矩阵，其m个对角线元素分别是m条传输线的电抗的倒数，其余元素都为0，所述传输线导纳矩阵C中传输线的排列顺序是与所述关联矩阵D中传输线的排列顺序相一致的；

步骤(3)，该计算机测定所述电网中每个节点上的净注入功率，得到电网节点的功率向量P，P＝[p₁，…，p_n]^T，所述每个节点上的净注入功率等于该节点上发电机的发电量与该节点上负载耗电量之差，所述向量P无论对于连通的电网或者不连通电网中的每一个孤岛，其净注入功率之和都为0；

步骤(4)，该计算机按下式计算所述电网的加权拉普拉斯矩阵(Laplacian)L：

L＝D^TCD；

步骤(5)，该计算机用奇异值分解法求出步骤(4)得到的矩阵L的广义逆矩阵L⁺；

步骤(6)，该计算机按下式计算所述电网中的直流潮流向量f：

f＝CDL⁺P，

向量f的分量为每一条传输线上的有功功率流，f的分量为正表示潮流方向与传输线的参考方向一致，f的分量为负表示潮流方向与传输线的参考方向相反，所述传输线的排列顺序与所述关联矩阵D和传输线导纳矩阵C的顺序一致。

相对于传统的潮流分析方法，所述方法具有如下优点：

第一，该方法的上述步骤对任意拓扑系统均适用。不论系统是否连通，步骤(4)中求得的系统加权拉普拉斯矩阵(Laplacian)的广义逆都是存在且唯一的，系统的直流潮流都可以使用步骤(5)中的计算式得到。