CN114172155A - 多重故障开断的潮流计算方法、系统、设备及存储介质 - Google Patents

多重故障开断的潮流计算方法、系统、设备及存储介质 Download PDF

Info

Publication number
CN114172155A
CN114172155A CN202111488863.4A CN202111488863A CN114172155A CN 114172155 A CN114172155 A CN 114172155A CN 202111488863 A CN202111488863 A CN 202111488863A CN 114172155 A CN114172155 A CN 114172155A
Authority
CN
China
Prior art keywords
branch
branches
disconnection
flow
tidal current
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111488863.4A
Other languages
English (en)
Inventor
王淼
於益军
吕闫
孙博
陈之栩
史沛然
黄宇鹏
史普鑫
于钊
郑晓雨
吴建斌
罗卫华
张国芳
邹鹏
石上丘
李劲松
罗雅迪
韩巍
王伟
王少芳
马晓忱
孙略
郎燕生
李理
李媛媛
陶蕾
宋旭日
黄运豪
刘升
门德月
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanxi Electric Power Co Ltd
North China Grid Co Ltd
State Grid Sichuan Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Shanxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanxi Electric Power Co Ltd
North China Grid Co Ltd
State Grid Sichuan Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Shanxi Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Shanxi Electric Power Co Ltd, North China Grid Co Ltd, State Grid Sichuan Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Shanxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202111488863.4A priority Critical patent/CN114172155A/zh
Publication of CN114172155A publication Critical patent/CN114172155A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供一种多重故障开断的潮流计算方法、系统、设备及存储介质，方法包括：根据电网基态潮流断面，计算N‑1支路开断分布因子；对多重故障所有支路依次进行N‑1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N‑1开断，直到满足收敛要求；将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。该方法能够有效的减少计算规模，提高了计算速度。

Description

多重故障开断的潮流计算方法、系统、设备及存储介质
技术领域
本发明属于电力调度领域，涉及一种多重故障开断的潮流计算方法、系统、设备及存储介质。
背景技术
电力市场环境下，区域安全校核基于电网未来运行断面，采用预想事故分析的方法来评估电网的安全风险，以便及早采取相应的措施，防范电网事故的发生，目前面临着电网设备更多、电网结构更紧密、电源类型更复杂、新能源消纳需求更迫切、电网运行风险更大的问题。
预想故障分析已实现了《电力系统安全稳定导则》中要求的N-1故障开断分析，N-1安全分析的实质是电力系统运行的稳态分析问题，即潮流问题。但随着市场环境下负荷密度的不断增大，电网结构越来越紧密,N-n故障的发生率也随之增加，对电网只进行“N-1”的检验已无法满足实际电网安全稳定运行的要求，在N-1基础上进行N-n多重故障分析成为电网安全可靠性评估的关键问题之一。但大规模电力系统元件众多，N-n多重故障组合数量很大，对多重故障都进行详细拓扑分析和潮流计算工作量很大，且计算速度慢。
发明内容
为解决现有技术中存在的问题，本发明提供了一种多重故障开断的潮流计算方法、系统、设备及存储介质，该方法能够有效的减少计算规模，提高了计算速度。
为达到上述目的，本发明采用了以下技术方案：
一种多重故障开断的潮流计算方法，包括以下步骤：
根据电网基态潮流断面，计算N-1支路开断分布因子；
对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
作为本发明的进一步改进，所述根据电网基态潮流断面，计算N-1支路开断分布因子，包括以下步骤：
当电网支路开断引起其他支路的潮流变化，变化量计算如下：
Figure BDA0003397656720000021
其中，Pl为开断支路l的功率，
Figure BDA0003397656720000022
为支路l开断后所引起的支路k的潮流变化量，Dk-l为支路开断分布因子，计算公式如下所示：
Figure BDA0003397656720000023
其中，xk为支路k的电抗值，xl为支路l的电抗值,i,j为支路l的首末节点，m,n为支路k的首末节点，Xmi为节点阻抗矩阵的第m行、第i列元素。
作为本发明的进一步改进，所述对多重故障支路组依次进行N-1开断，计算开断支路后所引起的其他支路潮流变化，具体为：
所述对于n条故障支路组，从第1条到第n条依次开断，计算第k条支路的潮流变化，具体包括：
开断多重故障的第1条支路，则第1条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000024
其中，
Figure BDA0003397656720000025
为第a条支路在第1步的潮流值，
Figure BDA0003397656720000026
为第a条支路的潮流初值，Da-1为第1条开断支路对第a条支路的转移因子,
Figure BDA0003397656720000027
为第k条支路在第1步的潮流变化值，n为多重故障的支路数，a＝2,…，n；
开断多重故障的第2条支路，则第2条支路潮流为0，其潮流转移到其他支路上，也有一部分潮流转移到了第1条支路上：
Figure BDA0003397656720000031
第n步，开断多重故障的第n条支路，则第n条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000032
此时，a＝1,2,…，n-1
第1轮n条支路开断计算完成后，计算第k条支路在所有开断步数下的潮流变化量ΔPk1：
Figure BDA0003397656720000033
其中ΔPks为第s轮开断后，第k条支路的潮流变化量,s＝1,2,……,M,M为总的开断轮次数。
作为本发明的进一步改进，所述判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求，具体为：
在每一轮支路组开断结束后，判断
Figure BDA0003397656720000034
是否小于收敛门槛，如果大于收敛门槛，返回重新从第1条到第n条依次开断，进行新的一轮开断计算；
第M轮开断计算完成后，各开断支路潮流小于收敛门槛，则完成计算。
作为本发明的进一步改进，所述对多重故障所有支路依次进行N-1开断步骤，采用多线程并行计算。
作为本发明的进一步改进，所述将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流，具体包括：
将所有开断步数下的第k条支路潮流变化量相加，得到第k条支路总的潮流变化量：
Figure BDA0003397656720000041
一种多重故障开断的潮流计算系统，包括：
分布因子计算模块，用于根据电网基态潮流断面，计算N-1支路开断分布因子；
开断潮流变化计算模块，用于对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
叠加潮流计算模块，用于将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
作为本发明的进一步改进，所述分布因子计算模块，具体用于：
当电网支路开断引起其他支路的潮流变化，变化量计算如下：
Figure BDA0003397656720000042
其中，Pl为开断支路l的功率，
Figure BDA0003397656720000043
为支路l开断后所引起的支路k的潮流变化量，Dk-l为支路开断分布因子，计算公式如下所示：
Figure BDA0003397656720000044
其中，xk为支路k的电抗值，xl为支路l的电抗值,i,j为支路l的首末节点，m,n为支路k的首末节点，Xmi为节点阻抗矩阵的第m行、第i列元素。
作为本发明的进一步改进，所述开断潮流变化计算模块，具体用于：
所述对于n条故障支路组，从第1条到第n条依次开断，计算第k条支路的潮流变化，具体包括：
开断多重故障的第1条支路，则第1条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000051
其中，
Figure BDA0003397656720000052
为第a条支路在第1步的潮流值，
Figure BDA0003397656720000053
为第a条支路的潮流初值，Da-1为第1条开断支路对第a条支路的转移因子,
Figure BDA0003397656720000054
为第k条支路在第1步的潮流变化值，n为多重故障的支路数，a＝2,…，n；
开断多重故障的第2条支路，则第2条支路潮流为0，其潮流转移到其他支路上，也有一部分潮流转移到了第1条支路上：
Figure BDA0003397656720000055
第n步，开断多重故障的第n条支路，则第n条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000056
此时，a＝1,2,…，n-1
第1轮n条支路开断计算完成后，计算第k条支路在所有开断步数下的潮流变化量ΔPk1：
Figure BDA0003397656720000057
其中ΔPks为第s轮开断后，第k条支路的潮流变化量,s＝1,2,……,M,M为总的开断轮次数。
作为本发明的进一步改进，所述叠加潮流计算模块，具体用于：
将所有开断步数下的第k条支路潮流变化量相加，得到第k条支路总的潮流变化量：
Figure BDA0003397656720000061
一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述多重故障开断的潮流计算方法的步骤。
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述多重故障开断的潮流计算方法的步骤。
与现有技术相比，本发明的有益效果体现在：
本发明提出了提出一种适用于N-n多重故障开断的快速潮流计算方法，利用基态潮流断面的支路开断分布因子进行N-x多重开断快速扫描，将N-x多重故障分解为多个N-1故障，将N-n多重故障开断分解为多个N-1开断，通过开断后支路潮流变化的叠加，得到各支路的潮流，通过依次开断多重故障的每一条支路，根据支路开断分布因子计算N-1开断后的支路潮流转移，直到多重故障的每一条支路潮流小于收敛门槛，分析其他支路的潮流情况。能够有效的减少计算规模，提高了计算速度，本发明也适用于N-2的开断分析，适用范围更广。
附图说明
图1为本发明多重故障开断的潮流计算方法流程简图；
图2为本发明多重故障开断的潮流计算方法流程详细图；
图3为本发明多重故障开断的潮流计算系统示意图；
图4为本发明一种电子设备结构示意图。
具体实施方式
下面将参考附图并结合实施例来详细说明本发明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
以下详细说明均是示例性的说明，旨在对本发明提供进一步的详细说明。除非另有指明，本发明所采用的所有技术术语与本申请所属领域的一般技术人员的通常理解的含义相同。本发明所使用的术语仅是为了描述具体实施方式，而并非意图限制根据本发明的示例性实施方式。
相关术语解释
N-1故障开断分析：对电网支路逐个开断，本支路的潮流变为零，潮流转移到其他支路，计算其他支路的潮流变化情况；
多重故障：多个支路组进行故障开断分析；
开断分布因子：电网支路开断后，其他支路的潮流变化量与本支路的潮流值的比值，小于1。
本发提出一种适用于N-n多重故障开断的快速潮流计算方法，基于N-1的支路开断分布因子，通过反复支路开断和计算支路的转移潮流，直到支路开断的潮流都转移到其他支路。本发明将N-x多重故障开断分解为多个N-1开断，通过支路开断后其他支路潮流变化的叠加，得到其他各支路的潮流，能够有效的解决多重开断下的潮流快速计算的问题，提升故障扫描计算效率。
如图1所示，本发明提出了提出一种适用于N-n多重故障开断的快速潮流计算方法，包括以下步骤：
首先，根据电网基态潮流断面，计算N-1支路开断分布因子；
其次，对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
最后，将开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
本发明将N-n多重故障开断分解为多个N-1开断，通过支路组开断后其他支路潮流变化的叠加，得到其他各支路的潮流，能够有效的减少计算规模，提高了计算速度。尤其是将支路开断分布因子用于N-n多重故障分析中，基于基态的潮流和拓扑结构，根据N-1条件下的开断分布因子，快速计算支路组开断后其他支路的潮流变化，分析设备越限情况，筛选出严重故障，然后针对严重的故障，采用详细潮流计算，得到更为准确的电网安全风险评估结果。
如图2所示，本发明的方法具体包括以下步骤：
步骤(1)：根据电网基态潮流断面，计算N-1支路开断分布因子；所述步骤(1)，计算N-1支路开断分布因子主要包括以下步骤：
当电网支路开断引起其他支路的潮流变化，其变化量可用以下公式表示：
Figure BDA0003397656720000081
其中，Pl为支路开断l的功率，
Figure BDA0003397656720000082
为支路l开断后所引起的支路k的潮流变化量，Dk-l为支路开断分布因子，其计算公式如下所示：
Figure BDA0003397656720000083
其中，xk为支路k的电抗值，xl为支路l的电抗值,i,j为支路l的首末节点，m,n为支路k的首末节点，Xmi为节点阻抗矩阵的第m行、第i列元素。
步骤(2)：对多重故障支路组依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化；对于n条故障支路组，从第1条到第n条依次开断，计算第k条支路的潮流变化，主要包括以下步骤：
2-1.第1步，开断多重故障的第1条支路，则第1条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000084
其中，
Figure BDA0003397656720000091
为第a条支路在第1步的潮流值，
Figure BDA0003397656720000092
为第a条支路的潮流初值，Da-1为第1条支路开断对第a条支路的转移因子,
Figure BDA0003397656720000093
为第k条支路在第1步的潮流变化值，n为多重故障的支路数，a＝2,…，n。
2-2.第2步，开断多重故障的第2条支路，则第2条支路潮流为0，其潮流转移到其他支路上，也有一部分潮流转移到了第1条支路上：
Figure BDA0003397656720000094
2-3.第n步，开断多重故障的第n条支路，则第n条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000095
此时，a＝1,2,…，n-1
2-4.第1轮n条支路开断计算完成后，计算第k条支路在所有开断步数下的潮流变化量ΔPk1：
Figure BDA0003397656720000096
其中ΔPks为第s轮开断后，第k条支路的潮流变化量,s＝1,2,……,M,M为总的开断轮次数。
在每一轮支路组开断结束后，判断
Figure BDA0003397656720000097
是否小于收敛门槛，如果大于收敛门槛，则回到2-1，进行新的一轮开断计算。
步骤(3)：将开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流；所述步骤(3)，具体步骤如下：
第M轮开断计算完成后，各支路开断潮流小于收敛门槛，则完成计算，将所有开断步数下的第k条支路潮流变化量相加，得到第k条支路总的潮流变化量：
Figure BDA0003397656720000101
由于N-n多重开断的规模巨大，但各个开断故障之间计算独立，可采用多线程并行计算加快计算速度。
如图3所示，本发明还提供一种多重故障开断的潮流计算系统，包括：
分布因子计算模块，用于根据电网基态潮流断面，计算N-1支路开断分布因子；
开断潮流变化计算模块，用于对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
叠加潮流计算模块，用于将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
具体的，所述分布因子计算模块，具体用于：
当电网支路开断引起其他支路的潮流变化，变化量计算如下：
Figure BDA0003397656720000102
其中，Pl为开断支路l的功率，
Figure BDA0003397656720000103
为支路l开断后所引起的支路k的潮流变化量，Dk-l为支路开断分布因子，计算公式如下所示：
Figure BDA0003397656720000104
其中，xk为支路k的电抗值，xl为支路l的电抗值,i,j为支路l的首末节点，m,n为支路k的首末节点，Xmi为节点阻抗矩阵的第m行、第i列元素。
具体的，所述开断潮流变化计算模块，具体用于：
所述对于n条故障支路组，从第1条到第n条依次开断，计算第k条支路的潮流变化，具体包括：
开断多重故障的第1条支路，则第1条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000111
其中，
Figure BDA0003397656720000112
为第a条支路在第1步的潮流值，
Figure BDA0003397656720000113
为第a条支路的潮流初值，Da-1为第1条开断支路对第a条支路的转移因子,
Figure BDA0003397656720000114
为第k条支路在第1步的潮流变化值，n为多重故障的支路数，a＝2,…，n；
开断多重故障的第2条支路，则第2条支路潮流为0，其潮流转移到其他支路上，也有一部分潮流转移到了第1条支路上：
Figure BDA0003397656720000115
第n步，开断多重故障的第n条支路，则第n条支路潮流为0，其潮流转移到其他支路上：
Figure BDA0003397656720000116
此时，a＝1,2,…，n-1
第1轮n条支路开断计算完成后，计算第k条支路在所有开断步数下的潮流变化量ΔPk1：
Figure BDA0003397656720000117
其中ΔPks为第s轮开断后，第k条支路的潮流变化量,s＝1,2,……,M,M为总的开断轮次数。
具体的，所述叠加潮流计算模块，具体用于：
将所有开断步数下的第k条支路潮流变化量相加，得到第k条支路总的潮流变化量：
Figure BDA0003397656720000121
如图4所示，本发明第三个目的是提供一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述多重故障开断的潮流计算方法的步骤。
所述多重故障开断的潮流计算方法包括以下步骤：
根据电网基态潮流断面，计算N-1支路开断分布因子；
对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
本发明第四个目的是提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述多重故障开断的潮流计算方法的步骤。
所述多重故障开断的潮流计算方法包括以下步骤：
根据电网基态潮流断面，计算N-1支路开断分布因子；
对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (12)

1.一种多重故障开断的潮流计算方法，其特征在于，包括以下步骤：
根据电网基态潮流断面，计算N-1支路开断分布因子；
对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
2.根据权利要求1所述的方法，其特征在于，
所述根据电网基态潮流断面，计算N-1支路开断分布因子，包括以下步骤：
当电网支路开断引起其他支路的潮流变化，变化量计算如下：
Figure FDA0003397656710000011
其中，Pl为开断支路l的功率，
Figure FDA0003397656710000012
为支路l开断后所引起的支路k的潮流变化量，Dk-l为支路开断分布因子，计算公式如下所示：
Figure FDA0003397656710000013
其中，xk为支路k的电抗值，xl为支路l的电抗值,i,j为支路l的首末节点，m,n为支路k的首末节点，Xmi为节点阻抗矩阵的第m行、第i列元素。
3.根据权利要求1所述的方法，其特征在于，
所述对多重故障支路组依次进行N-1开断，计算开断支路后所引起的其他支路潮流变化，具体为：
所述对于n条故障支路组，从第1条到第n条依次开断，计算第k条支路的潮流变化，具体包括：
开断多重故障的第1条支路，则第1条支路潮流为0，其潮流转移到其他支路上：
Figure FDA0003397656710000021
其中，
Figure FDA0003397656710000022
为第a条支路在第1步的潮流值，
Figure FDA0003397656710000023
为第a条支路的潮流初值，Da-1为第1条开断支路对第a条支路的转移因子,
Figure FDA0003397656710000024
为第k条支路在第1步的潮流变化值，n为多重故障的支路数，a＝2,…，n；
开断多重故障的第2条支路，则第2条支路潮流为0，其潮流转移到其他支路上，也有一部分潮流转移到了第1条支路上：
Figure FDA0003397656710000025
第n步，开断多重故障的第n条支路，则第n条支路潮流为0，其潮流转移到其他支路上：
Figure FDA0003397656710000026
此时，a＝1,2,…，n-1
第1轮n条支路开断计算完成后，计算第k条支路在所有开断步数下的潮流变化量ΔPk1：
Figure FDA0003397656710000027
其中ΔPks为第s轮开断后，第k条支路的潮流变化量,s＝1,2,……,M,M为总的开断轮次数。
4.根据权利要求3所述的方法，其特征在于，
所述判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求，具体为：
在每一轮支路组开断结束后，判断
Figure FDA0003397656710000031
是否小于收敛门槛，如果大于收敛门槛，返回重新从第1条到第n条依次开断，进行新的一轮开断计算；
第M轮开断计算完成后，各开断支路潮流小于收敛门槛，则完成计算。
5.根据权利要求1所述的方法，其特征在于，
所述对多重故障所有支路依次进行N-1开断步骤，采用多线程并行计算。
6.根据权利要求1所述的方法，其特征在于，
所述将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流，具体包括：
将所有开断步数下的第k条支路潮流变化量相加，得到第k条支路总的潮流变化量：
Figure FDA0003397656710000032
7.一种多重故障开断的潮流计算系统，其特征在于，包括：
分布因子计算模块，用于根据电网基态潮流断面，计算N-1支路开断分布因子；
开断潮流变化计算模块，用于对多重故障所有支路依次进行N-1开断，计算支路开断后所引起的其他支路潮流变化，判断多重故障支路所有支路潮流是否满足收敛要求，如果不满足，继续进行N-1开断，直到满足收敛要求；
叠加潮流计算模块，用于将支路开断引起的潮流变化叠加到其他支路的基础潮流上，得到所有支路的潮流。
8.根据权利要求7所述的系统，其特征在于，
所述分布因子计算模块，具体用于：
当电网支路开断引起其他支路的潮流变化，变化量计算如下：
Figure FDA0003397656710000033
其中，Pl为开断支路l的功率，
Figure FDA0003397656710000041
为支路l开断后所引起的支路k的潮流变化量，Dk-l为支路开断分布因子，计算公式如下所示：
Figure FDA0003397656710000042
其中，xk为支路k的电抗值，xl为支路l的电抗值,i,j为支路l的首末节点，m,n为支路k的首末节点，Xmi为节点阻抗矩阵的第m行、第i列元素。
9.根据权利要求7所述的系统，其特征在于，
所述开断潮流变化计算模块，具体用于：
所述对于n条故障支路组，从第1条到第n条依次开断，计算第k条支路的潮流变化，具体包括：
开断多重故障的第1条支路，则第1条支路潮流为0，其潮流转移到其他支路上：
Figure FDA0003397656710000043
其中，
Figure FDA0003397656710000044
为第a条支路在第1步的潮流值，
Figure FDA0003397656710000045
为第a条支路的潮流初值，Da-1为第1条开断支路对第a条支路的转移因子,
Figure FDA0003397656710000046
为第k条支路在第1步的潮流变化值，n为多重故障的支路数，a＝2,…，n；
开断多重故障的第2条支路，则第2条支路潮流为0，其潮流转移到其他支路上，也有一部分潮流转移到了第1条支路上：
Figure FDA0003397656710000047
第n步，开断多重故障的第n条支路，则第n条支路潮流为0，其潮流转移到其他支路上：
Figure FDA0003397656710000051
此时，a＝1,2,…，n-1
第1轮n条支路开断计算完成后，计算第k条支路在所有开断步数下的潮流变化量ΔPk1：
Figure FDA0003397656710000052
其中ΔPks为第s轮开断后，第k条支路的潮流变化量,s＝1,2,……,M,M为总的开断轮次数。
10.根据权利要求7所述的系统，其特征在于，
所述叠加潮流计算模块，具体用于：
将所有开断步数下的第k条支路潮流变化量相加，得到第k条支路总的潮流变化量：
Figure FDA0003397656710000053
11.一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1-6任一项所述多重故障开断的潮流计算方法的步骤。
12.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1-6任一项所述多重故障开断的潮流计算方法的步骤。
CN202111488863.4A 2021-12-07 2021-12-07 多重故障开断的潮流计算方法、系统、设备及存储介质 Pending CN114172155A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111488863.4A CN114172155A (zh) 2021-12-07 2021-12-07 多重故障开断的潮流计算方法、系统、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111488863.4A CN114172155A (zh) 2021-12-07 2021-12-07 多重故障开断的潮流计算方法、系统、设备及存储介质

Publications (1)

Publication Number Publication Date
CN114172155A true CN114172155A (zh) 2022-03-11

Family

ID=80484158

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111488863.4A Pending CN114172155A (zh) 2021-12-07 2021-12-07 多重故障开断的潮流计算方法、系统、设备及存储介质

Country Status (1)

Country Link
CN (1) CN114172155A (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104092210A (zh) * 2014-07-14 2014-10-08 国家电网公司 一种具备n-2开断快速扫描功能的在线静态安全分析方法
CN107959284A (zh) * 2016-10-14 2018-04-24 华北电力大学（保定） 一种基于多支路切除和背离路径的潮流转移识别算法
CN110661265A (zh) * 2019-09-19 2020-01-07 吉林大学 一种基于支路开断分布因子的安全约束最优潮流计算方法
WO2020078109A1 (zh) * 2018-10-17 2020-04-23 中国电力科学研究院有限公司 识别电网薄弱断面的方法、装置及存储介质

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104092210A (zh) * 2014-07-14 2014-10-08 国家电网公司 一种具备n-2开断快速扫描功能的在线静态安全分析方法
CN107959284A (zh) * 2016-10-14 2018-04-24 华北电力大学（保定） 一种基于多支路切除和背离路径的潮流转移识别算法
WO2020078109A1 (zh) * 2018-10-17 2020-04-23 中国电力科学研究院有限公司 识别电网薄弱断面的方法、装置及存储介质
CN110661265A (zh) * 2019-09-19 2020-01-07 吉林大学 一种基于支路开断分布因子的安全约束最优潮流计算方法

Similar Documents

Publication Publication Date Title
CN103995215B (zh) 一种基于多层次反馈调整的智能电网故障诊断方法
CN110635479B (zh) 一种限制短路电流运行方式的智能辅助决策方法及系统
CN101562340B (zh) 一种求解电力系统故障后临界近似潮流的方法
CN112347716B (zh) 基于q学习的电网脆弱点检测方法、系统、设备及介质
CN106875127B (zh) 统一潮流控制器可靠性建模及其接入电网可靠性评估方法
CN111900724B (zh) 一种电力系统暂态稳定紧急控制在线决策方法、装置及存储介质
CN115453356B (zh) 一种动力设备运行状态监测分析方法、系统、终端及介质
CN112462198A (zh) 一种基于自编码器的电网故障线路判定方法及系统
CN113468745B (zh) 一种基于历史故障的配电网可靠性快速评估方法及系统
CN111769545B (zh) 一种在相继的多重开断场景下潮流评估方法及装置
KR101568873B1 (ko) 대규모 계통 적용을 위한 고속화 방법
Huang et al. Vulnerability analysis of electrical cyber physical systems using a simulation platform
CN111864744B (zh) 高比例水电系统调速器控制模式在线切换方法及系统
CN114172155A (zh) 多重故障开断的潮流计算方法、系统、设备及存储介质
Gusev et al. Recursion based contingency analysis of an electrical power system
CN110021936B (zh) 含环配电网运行状态的确定方法及装置
CN112769124A (zh) 一种基于潮流转移和追踪的电力系统快速运行风险评估方法
CN116541748A (zh) 基于人工智能技术的电网故障诊断方法和系统
CN115954956A (zh) 一种配电网分布式电源接入容量评估方法及系统
CN113300367B (zh) 电网关键断面在线快速识别方法、输电限额制定方法及装置
CN115587459A (zh) 一种基于图最大密度的配电网辐射状约束建模方法及系统
CN113985207A (zh) 一种电网运行设备故障监测方法、系统、装置及存储介质
CN113283703A (zh) 电网预防控制决策生成方法、系统、存储介质及计算设备
KR101736256B1 (ko) 고장전류 기여율 산정 장치 및 방법
CN113945860A (zh) 考虑高速开关动作时序的短路电流在线校核方法、装置及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20220311