一种基于多支路切除和背离路径的潮流转移识别算法的制作方法

本发明涉及电网公司电力系统保护与安全控制领域，具体地说，是一种基于多支路切除和背离路径的潮流转移识别算法。

背景技术：

电网中由故障支路切除引起的潮流转移是近年来世界各国发生输电线路连锁切除事故的重要原因，严重危害着系统的安全稳定运行。在此背景下，目前电网调度中心严格按照《电力系统安全稳定导则》要求安排电网运行方式，运行中严格执行“n-1”准则，如设备检修，机组调控等操作，在系统运行状态发生变化时及时消除安全隐患，这样可以使由单一元件开断引发支路连锁跳闸的几率越来越小；但是，国内大部分电网还没有达到能够抵御同一断面多条支路切除严重故障的强度，这主要是因为一方面校验计算复杂，另一方面若要满足多支路切除安全校验需提高线路的热稳定极限功率。由于电网缺乏对多支路切除情况的应对措施，一旦发生此类故障时任由后备保护切除线路，大大增加了输电线路连锁跳闸的几率。所以，对多条支路切除后的输电断面的行为分析，即快速找出潮流转移较大的支路并判断其上潮流是否越限，对预防线路连锁跳闸具有重要意义。

技术实现要素：

本发明的目的是针对现有方法的不足，提出一种基于多支路切除和背离路径的潮流转移识别算法，以适应于发生多支路切除情况下潮流转移的识别，在减少分析计算量的同时有效避免了输电断面漏选问题。

为了达到上述目的，本发明采用的技术方案是：

一种基于多支路切除和背离路径的潮流转移识别算法，该方法根据故障支路切除后的潮流转移特征，首先深入分析了转移到非故障线路上的潮流组成，其次在一定的节点深度内搜索切除支路两端节点间的前k条最短路径集，最终识别出受潮流转移影响较大的线路构成输电断面。本发明在切除支路的潮流转移估算、k值选取和断面判据方面作了重要改进，以适应于发生多支路切除情况下潮流转移的识别；重新建立了一种k的取值方案，并采用背离路径算法确定最短路径集，在减少分析计算量的同时有效避免了输电断面漏选问题。ieee30节点系统的仿真验证了本发明算法的有效性。。

一种基于多支路切除和背离路径的潮流转移识别算法，所述方法包括以下步骤：

步骤1：对转移到非故障支路上的潮流组成进行分析，发生多支路切除时，支路k上的转移潮流可看成多条以等效潮流为初始潮流的支路单独开断后转移到支路k上潮流的叠加，将多支路切除简化为多条单一支切除情况，避免了直流法中对节点阻抗矩阵的大量修正计算，并实现了对潮流转移的快速估算；

步骤2：依据图搜索中背离路径算法，在一定的节点深度内搜索切除支路两端节点间的前k条最短路径，确定潮流转移较大的区域范围。文中结合了多支路切除的潮流转移特征，重新建立一种k的取值方案，并采取边搜索边修正k值的策略，避免搜索范围设置不当而漏选支路的问题；

步骤3：考虑电网的实时运行状态，给出了输电断面的选取原则，最终识别出实际潮流增加较大的支路集，且能成功找出潮流反向增大的线路。

作为进一步描述，步骤1中提到若电网中发生支路切除，由切除支路l转移到支路k上的潮流量(以初始潮流方向为正)可通过l对k的支路开断分布因子dk-l表示：

δpk＝dk-lpl

其中：pl为切除支路l的初始有功潮流；δpk为转移到支路k上的潮流。

由上式可知，若已知分布因子即可计算切除支路对任意支路的潮流转移量，而dk-l可由系统的节点阻抗矩阵计算得到，其计算式如下：

其中：xk-l、xl-l、xk、xl只与电网的拓扑参数有关。对于一个确定的网络，支路l对任意支路k的分布因子dk-l随之确定。

作为进一步描述，步骤2提到由电路基本方程可知：当发生支路切除时，可将切除支路用附加电流源(与支路切除前的电流等量、反向)替代后的等值无源网络中的潮流分布来近似切除支路的潮流转移。若对附加电源搜索其前k条最短闭合路径回路，则这些回路包含该等值网络中绝大多数电流较大的输电线路，也即包含了受潮流转移影响较大的支路，理论上说明了前k最短路径法搜索输电断面的可行性。

多支路切除时转移潮流的分布规律与单一支路切除情况差异较大，因此需重新设计k的取值方案。

进一步地，步骤3中提到目前已积累许多特点各异的前k最短路径算法，其中背离路径算法通过求取第k-1条路径的背离路径来确定其第k条最短路径，与本发明中k值的修正过程相符，且每派生一条背离路径只需执行一次dijstra算法，算法的时间复杂度较小，更易满足搜索的实时性要求。综合上述优点，本发明中采用背离路径算法对发生多支路切除时受潮流转移影响严重的线路进行识别。

与现有方法相比，本发明的有益效果是：本发明将路径搜索引入多支路切除的潮流转移识别算法中，是一种新的尝试。本发明中通过分析计算，在获取潮流估算式的同时得到经修正后的路径数，进而结合路径搜索与实时运行状态识别出潮流增加较大的支路。所提算法具有如下特点：潮流转移分析模块采用了极限的思想，保持了初始网络的拓扑结构，避免对节点阻抗矩阵进行修改，大大减少矩阵运算量，提高算法效率；考虑了输电线路的初始潮流，从而能够识别得到实时的潮流增加较大支路集，所得结果更加符合实际电网运行情况；采用边搜索边修正k值的方法，避免搜索范围设置不当而漏选支路的问题；给出了反映潮流耦合强度的比例系数，通过计算该系数的取值可粗略判断多支路切除事故的严重程度。

附图说明

图1为：本发明一实施例的多支路切除的潮流转移模型。

图2为：本发明一实施例的潮流转移分析流程图。

图3为：本发明一实施例的路径搜索流程图。

图4为：本发明一实施例的潮流转移识别流程图。

图5为：本发明一实施例的ieee30节点系统简化拓扑图。

具体实施方式

下面通过附图和实施例对本发明的技术方案作进一步详细说明。

本发明所述的一种基于多支路切除和背离路径的潮流转移识别算法，包括如下步骤：

步骤1：对转移到非故障支路上的潮流组成进行分析，发生多支路切除时，支路k上的转移潮流可看成多条以等效潮流为初始潮流的支路单独开断后转移到支路k上潮流的叠加，将多支路切除简化为多条单一支切除情况，避免了直流法中对节点阻抗矩阵的大量修正计算，并实现了对潮流转移的快速估算；

假设初始潮流分别pa与pb的支路a、b发生故障切除，支路k为网络中除支路a、b外的任意一条支路。多支路切除潮流转移可以分解为切除支路间无穷次潮流的来回折返过程。首先，支路a、b的初始潮流分别依据分布因子dk-a、dk-b转移到支路k上，此外pa按照分布因子db-a转移到支路b，pb按照分布因子da-b转移到支路a，完成一次支路间潮流交换过程；其次，上述支路间交换的潮流替代了切除支路的初始潮流，继续对网络进行转移，依此类推，直至潮流全部转移到网络中其他支路上，如图1所示，图中：表示支路b第n次转移到支路a上的潮流；表示支路a第n次转移到支路b上的潮流。

由图1可知：pa和依据分布因子dk-a向支路k进行转移，pb和依据分布因子dk-b向支路k进行转移。将他们转移到支路k上的分量进行累加，即可得到转移到支路k上的总潮流为：

可分解为以下两部分：

1)pb＝0时

此时有如下关系式：

由上两式可推得：

2)pa＝0时

此时有如下关系式：

由上两式推得：

将两部分中推导出的最终两式叠加可求得潮流为：

同理可得为：

将上两式代入第一个式子中可得：

δpk＝dk-ap′a+dk-bp′b

其中：p′a、p′b分别为支路a、b的等效潮流，其意义：将等效潮流替代初始潮流后，可以消除开断支路间的潮流耦合关系，从而直接对支路k进行潮流转移。由上述推导可知，当断线支路确定后，p′a与p′b可通过初始条件计算得到。

由最后一个式子可得：发生多支路切除时，支路k上的转移潮流可看成多条以等效潮流为初始潮流的支路单独开断后转移到支路k上潮流的叠加，将多支路切除简化为多条单一支切除情况，避免了直流法中对节点阻抗矩阵的大量修正计算，并实现了对潮流转移的快速估算。

步骤2：依据图搜索中背离路径算法，在一定的节点深度内搜索切除支路两端节点间的前k条最短路径，确定潮流转移较大的区域范围。文中结合了多支路切除的潮流转移特征，重新建立一种k的取值方案，并采取边搜索边修正k值的策略，避免搜索范围设置不当而漏选支路的问题；

多支路切除时转移潮流的分布规律与单一支路切除情况差异较大，因此需重新设计k的取值方案。转移到支路k的潮流为独立的两部分：一部分为支路a的等效潮流p′a转移到支路k上的分量δpk1；另一部分为支路b的等效潮流p′b转移到支路k上的分量δpk2。

δpk＝δpk1+δpk2＝dk-ap′a+dk-bp′b

本文采用的取值方案：首先由δpk1确定支路a的基本路径数ka，其次把δpk2作为辅助判据，修正ka取值，将经修正后的ka作为支路a的所求路径数。其步骤如下：

步骤1：确定基本路径数ka

δpk1是支路a直接对支路k进行转移的潮流，与支路a发生单支路切除的转移潮流分布规律是一致的，可类比于单支路切除时的k值求取方法：

其中：β为潮流转移系数的阈值，考虑算法的效率，一般取0.3。

此时，支路a已用等效潮流p′a代替其初始潮流pa。则转移到其余支路上的潮流是支路a单独开断时的p′a/pa倍，令a＝p′a/pa，则ka应满足：

结合上两式可得：

ka＝|[-ak]|

其中[·]算子表示不超过括号内的最小整数。

步骤2：修正ka

由最短路径可知，除去前ka条最短路径外，转移到支路a两端节点间的第ka+1条最短路径上的潮流是其余最短路径中的最大值。因此，只需判断第ka+1条最短路径中加入分量δpk2后，其上转移潮流是否越限即可，修正过程如下所述：

①搜索第ka+1条最短路径。

②找出该路径中与前ka条路径不关联的线路构成新增线路集，计算新增线路集中任意线路t对支路a的潮流转移系数δpt/pa，若满足下式，则令ka＝ka+1。

其中δpt为转移到支路t上的潮流量，ε为潮流转移系数的上限系数。当δpt/pa大于ε时，说明δpk2占δpt的比例较大，即支路t有更大的几率包含在支路b两端节点间的前若干条最短路径中，为了避免重复搜索，故不修正ka。经过大量仿真验证，当ε＝0.5时，能够有效抑制ka的增大且搜索结果包含所有受潮流转移影响较大的支路。

③重复执行①②直至第ka+1条最短路径上没有满足上式的新增线路。

同理可得支路b的路径数kb的取值。

综合上述改进方案：考虑了多支路切除潮流转移特征，并合理的在单支路切除时k值选取原则的基础上进行改进，易于实现；所有参数均能通过初始条件求得，说明本方案可行；一般a值不会太大，满足快速性要求。

步骤3：考虑电网的实时运行状态，给出了输电断面的选取原则，最终识别出实际潮流增加较大的支路集，且能成功找出潮流反向增大的线路。

识别算法流程按照下述模块化进行设计：

1)潮流转移分析模块

该模块的功能是提供网络拓扑结构与切除支路的等效潮流，为路径搜索提供初始数据，其流程如图2所示。

2)路径搜索模块

该模块的功能是搜索切除支路两端节点间的若干条最短路径，对受多支路切除潮流转移影响较大的区域进行界定。由于电网中受潮流转移影响较大的元件不超过全网的20％，所以没必要进行全网分析，只需在一定的节点深度t内进行搜索即可，进而提高算法效率，搜索流程如图3所示。

3)潮流转移识别模块

该模块功能是结合电网的实时运行状态，识别出路径集中潮流实际增加较大的线路；将潮流增量大于某一特定值的线路加入输电断面。

δpk.0＞γpi.min

其中：δpk.0为支路k的潮流增量；pi.min为最小的切除支路潮流；γ为有功增加系数阈值，本文取0.3。

然而，电网中潮流实际增加较大输电线路存在两种情况：1)其上转移潮流与初始潮流方向一致；2)其上转移潮流与初始潮流方向相反。对于不同情况，线路上的潮流增量是不同的。

对于情况1：线路上的潮流增量即为转移到其上的潮流量，即δpk.0＝δpk。代入上式可得输电断面的判据：

δpk＞γpi.min

对于情况2：转移到线路的潮流δpk＜0，仅仅反应潮流减小，为了反应潮流增量，给出了该情况下潮流增加的判据：

|δpk|-2pk＞0

其中pk为支路k的初始潮流。

此时线路潮流增加量δpk.0＝|δpk|-2pk，带入第一个式子可得：

|δpk|-2pk＞γpi.min

因此，该模块的流程设计如图4所示。

对ieee30节点系统进行仿真验证，将不能与线路构成闭合回路的发电机与负荷节点从网络中剔除，得到的系统简化拓扑如图5所示。k等于4是相对保守的估计，一般第k条最短路径中的新增支路的潮流转移系数达不到1/k。因此，基准值k取3不会影响潮流转移的识别结果。

算例1：开断支路a、b分别为l18-19与l16-17。经计算后支路a的等效潮流p′a＝4.59mw，a＝1.51，基本路径数ka＝|[-ak]|＝5；同理可得kb＝5。路径搜索数据如表1、表2所示。

表1中基本路径数为5条，第6条路径中新出现支路l6-10的潮流转移系数为0.48，满足约束，说明l16-17的等效潮流转移到l6-10的分量使得其潮流转移系数超过阈值，则对ka进行修正；继续搜索第7条路径，然而第7条路径中新增支路l4-2与l2-6的潮流转移系数分别为0.13与0.14，不满足约束，不需要修正ka。同理表2中kb的修正过程与表1类似，最终经修正后ka＝6、kb＝5。

将路径集中的线路依据断面选取原则进行识别，其结果如表3所示。

算例2：开断支路a、b分别为l28-27与l12-15，经计算后支路a的等效潮流p′a＝18.38mw，a＝1.10，基本路径数ka＝|[-ak]|＝4；同理可得kb＝4。路径搜索结果与识别结果如表4、表5和表6所示。

由结果可知，算例1、算例2中分别搜索到10条、15条潮流增量较大的线路，相较于精确潮流计算所得的有功增加量，说明本发明对转移到非故障支路上的潮流计算方法能有效地近似线路的实际有功增加值，为进一步判断支路切除后输电断面内线路是否过载提供参考依据。并且有效找出潮流反向增大的支路，如表6中线路l21-22与l25-27。

对比分析算例1和算例2：从等效潮流与初始潮流的比值a可以看出，算例1中潮流间的耦合强度相较于算例2而言更为严重，导致算例1经修正后的路径数大于算例2的数值。在极端情况下，如开断线路为l15-23与l27-29时，切除线路间的开断分布因子为0，比值a＝1，则开断潮流间没有耦合关系，此时转移潮流在网络中的分布和单支路切除时是一致的。因此，系数a可以反映该支路与其他断线支路间的潮流耦合强度，并能够进一步预判此切除事故的严重程度。a值越接近1，切除支路的开断潮流在网络中的分布越接近单支路切除情况，随着a值增大，其开断潮流的分布与单支路切除情况时的差异也越大，事故也越严重。

表1为本发明一实施例的对l18-19的路径搜索数据。

表2为本发明一实施例的对l16-17的路径搜索数据。

表3为本发明一实施例的l18-19和l16-17切除事故的识别结果。

表4为本发明一实施例的对l28-27的路径搜索数据。

表5为本发明一实施例的对l12-15的路径搜索数据。

表6为本发明一实施例的l28-27和l12-15切除事故的识别结果。

表1

表2

表3

表4

表5

表6

从上述结果可知，本发明将路径搜索引入多支路切除的潮流转移识别算法中，是一种新的尝试。本发明中通过分析计算，在获取潮流估算式的同时得到经修正后的路径数，进而结合路径搜索与实时运行状态识别出潮流增加较大的支路。所提算法具有如下特点：潮流转移分析模块采用了极限的思想，保持了初始网络的拓扑结构，避免对节点阻抗矩阵进行修改，大大减少矩阵运算量，提高算法效率；考虑了输电线路的初始潮流，从而能够识别得到实时的潮流增加较大支路集，所得结果更加符合实际电网运行情况；采用边搜索边修正k值的方法，避免搜索范围设置不当而漏选支路的问题；给出了反映潮流耦合强度的比例系数，通过计算该系数的取值可粗略判断多支路切除事故的严重程度。