一种配电网灾前故障场景模拟预测方法及装置与流程

1.本发明涉及配电网韧性恢复技术领域，尤其涉及一种配电网灾前故障场景模拟预测方法及装置。


背景技术：

2.电力系统韧性是指电网在遭受重大灾害、人为攻击等极端情况下，改变自身状态减少故障损失，并尽快恢复正常供电水平的能力。配电网需要针对灾害做出提前部署与预防，合理考虑灾害不确定性进行灾前部署决策。
3.韧性研究关注极端自然灾害，这类小概率时间很难通过历史数据得到未来灾害造成的影响，不能简单采用概率模型衡量灾害结果。在这种情况下，场景抽样是一种常见的分析场景不确定性的方法，通过对场景聚类得出各场景的发生概率，从而为基于场景进行的优化决策提供基础。然而现有的场景聚类方法往往应用于连续变量的不确定性(如风电、光伏出力的不确定性)，鲜有针对配电网故障场景的聚类方法，因此有必要建立一种配电网故障场景聚类方法，为配电网灾前韧性决策提供依据。


技术实现要素：

4.本发明提供一种配电网灾前故障场景模拟预测方法及装置，用以解决现有技术中缺少配电网故障场景聚类研究的缺陷，实现基于配电网故障场景进行聚类，为配电网灾前韧性决策提供依据。
5.本发明提供一种配电网灾前故障场景模拟预测方法，包括：
6.建立配电网故障概率模型；
7.基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合；
8.在所述配电网故障场景集合中进行拓扑重构建立故障拓扑重构模型；
9.在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度；
10.根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景；
11.在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景；
12.根据所述第二场景聚类结果得到配电网灾前故障场景模拟预测结果。
13.根据本发明提供的一种配电网灾前故障场景模拟预测方法，所述配电网为拓扑结构，包括多条输电线路；所述基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合，具体包括：
14.将预设灾害信息输入所述配电网故障概率模型，得到所述配电网故障概率模型输
出的配电网各条输电线路故障的概率，其中，所述预设的灾害信息包括灾害持续时间；
15.根据各输电线路故障概率，对各输电线路进行灾害持续时间范围内的一次故障场景随机生成，得到一个故障场景；
16.重复进行随机生成的过程，得到预设生成场景总数的故障场景，形成所述配电网故障场景集合。
17.根据本发明提供的一种配电网灾前故障场景模拟预测方法，其特征在于，所述在所述配电网故障场景集合中进行拓扑重构建立故障拓扑重构模型，具体包括：所述故障拓扑重构模型的目标函数为：
[0018][0019][0020][0021]
其中，表示配电网中第i节点t时刻的电荷恢复量，ω
i
表示第i节点的负荷重要度，表示配电网中所有节点集合，t表示故障场景中灾害持续时间。
[0022]
根据本发明提供的一种配电网灾前故障场景模拟预测方法，所述在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度，具体包括：
[0023]
配电网在单个故障场景中单个节点的负荷损失量为：
[0024][0025]
其中，ω
i
表示第i节点的负荷重要度，表示配电网中第i节点t时刻的电荷恢复量，表示配电网中第i节点电荷恢复量上限，表示配电网中所有节点集合；
[0026]
单个故障场景下故障损失表示为节点故障损失向量p
loss
：
[0027][0028]
根据本发明提供的一种配电网灾前故障场景模拟预测方法，所述在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度，具体包括：
[0029]
应用树编辑距离(ted,tree edit distance)衡量配电网故障场景的拓扑相似度；
[0030]
两故障场景t1,t2之间的ted为：
[0031][0032]
其中，dist(t1,t2)表示t1,t2之间的ted，z＝{z1，z2,
…
,z
n
}表示将故障场景从t1转换为t2的一系列树编辑操作，|z|表示序列的操作总数，表示可完成故障场景t1,t2转换的操作序列集合，c(z
i
)表示操作z
i
的代价。
[0033]
根据本发明提供的一种配电网灾前故障场景模拟预测方法，所述根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景，具体包括：
[0034]
将通过配电网拓扑重构可消除的故障场景以及无故障场景归为一类第一场景聚类结果；
[0035]
在剩余的故障场景中选取k个聚类中心为初始聚类中心，k为大于零的正整数；
[0036]
以各初始聚类中心为中心分别进行至少一次故障场景聚集，将故障场景聚集后的各结果设为中间聚类结果；
[0037]
计算各中间聚类结果的中心：当各中间聚类结果的中心分别与对应的初始聚类中心一致时，停止聚类，各中间聚类结果分别为所述第一场景聚类结果；当各中间聚类结果的中心分别与对应的初始聚类中心不一致时，继续聚类，直至各中间聚类结果的中心分别与对应的初始聚类中心一致。
[0038]
根据本发明提供的一种配电网灾前故障场景模拟预测方法，所述在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景，具体包括：
[0039]
确定第一场景聚类结果中各故障场景的类间距离；
[0040]
确定第三数量的数值；
[0041]
按照第一场景聚类结果中各故障场景每次合并形成类簇时，选取故障场景合并后的类簇与其他类簇类间距离最大的方式进行聚类，形成第三数量的第二场景聚类结果。
[0042]
本发明还一种配电网灾前故障场景模拟预测装置，包括：
[0043]
故障概率建模单元，用于建立配电网故障概率模型；
[0044]
故障场景生成单元，用于基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合；
[0045]
故障拓扑重构单元，用于在所述配电网故障场景集合中对每个故障场景建立故障拓扑重构模型；
[0046]
故障损失及场景拓扑相似度处理单元，用于在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度；
[0047]
第一聚类单元，用于根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景；
[0048]
第二聚类单元，用于在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景；
[0049]
输出单元，用于根据所述第二场景聚类结果得到配电网灾前故障场景模拟预测结果。
[0050]
本发明还提供一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述配电网灾前故障场景模拟预测方法的步骤。
[0051]
本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述配电网灾前故障场景模拟预测方法的步骤。
[0052]
本发明提供的配电网灾前故障场景模拟预测方法及装置，针对配电网韧性恢复中灾害不确定性问题，本发明通过场景生成
‑
场景聚类方法，通过建立配电网元件在灾害下的故障概率模型，根据灾害预报信息进行场景生成，并利用场景的故障损失和拓扑相似度进行场景聚类，为基于场景进行的优化决策提供基础，以在灾前分析灾害对配电网的影响。
附图说明
[0053]
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0054]
图1是本发明提供的配电网灾前故障场景模拟预测方法流程图；
[0055]
图2是图1中步骤120的具体流程图；
[0056]
图3是图1中步骤150的具体流程图；
[0057]
图4是图1中步骤160的具体流程图；
[0058]
图5是应用本发明方法进行仿真实验的配电网节点分布结构示意图；
[0059]
图6是应用本发明方法进行仿真实验的拓扑聚类后仿真实验效果示意图；
[0060]
图7是发明提供的配电网灾前故障场景模拟预测装置结构示意图；
[0061]
图8是本发明提供的电子设备的结构示意图。
具体实施方式
[0062]
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0063]
在韧性相关的应用中往往涉及自然灾害这一不确定性因素，无法在事前获得灾害造成的配电网故障信息。为了减小灾害损失，配电网需要针对灾害做出提前部署与预防，合理考虑灾害不确定性进行灾前部署决策。本发明实施例中，通过建立配电网元件在灾害下的故障概率模型，根据灾害预报信息进行场景生成，并利用场景的故障损失和拓扑相似度
进行场景聚类，为后文基于场景进行的优化决策提供基础。
[0064]
韧性研究关注极端自然灾害，这类小概率时间很难通过历史数据得到未来灾害造成的影响，不能简单采用概率模型衡量灾害结果。本发明实施例建立一种配电网故障场景聚类方法，为灾前韧性决策提供依据。配电网故障场景聚类主要有以下难点：
[0065]
(1)故障场景涉及拓扑变化，需要考虑拓扑相似度对场景聚类的影响；
[0066]
(2)考虑故障损失和拓扑相似度的情况下，如何依据多维度的场景距离进行聚类。
[0067]
为解决以上问题，本发明实施例提供的场景聚类方法包含场景生成、场景聚类两个步骤。场景生成中，假设配电网元件的故障概率模型已知；在场景聚类中，引入树编辑距离衡量拓扑相似度，并在聚类过程中采用先依据故障损失聚类、再依据拓扑相似度聚类的两阶段方式进行聚类。
[0068]
如图1所示，本发明实施例提供一种配电网灾前故障场景模拟预测方法，包括：
[0069]
步骤110：建立配电网故障概率模型；
[0070]
由于配电网故障概率模型的建模在现有技术中已有研究，本发明实施例中，在现有故障概率模型的基础上进行一定简化，因此本发明实施例中，主要考虑台风这类灾害天气对配电网产生的影响。配电网元件中，假设故障由配电网线路中的杆塔造成，线路上的杆塔受风速和结冰影响而故障。
[0071]
基于以上假设和现有技术中模型，可以得到如下配电网线路故障概率模型：
[0072][0073][0074][0075]
(1)式中p
f,pole,t
表示t时刻某杆塔在风速w
t
(m/s)下的损坏概率，(2)式中表示t时刻某线路因风速导致的损坏概率，其由线路内的杆塔损坏概率和杆塔数量决定，(3)式中p
f,line,t
表示综合风力和结冰影响导致的损坏概率。
[0076]
本发明实施例中应用的故障概率模型中，为简化起见，配电网内的杆塔参数设定一致，线路的损坏概率差异只受线路上杆塔数量影响。在其他实施例中，可根据需要替换更精确的配电网故障概率模型。
[0077]
步骤120：基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合；
[0078]
本步骤中，根据配电网在灾害下的故障概率模型进行场景生成，获得一定数量的灾害故障场景。
[0079]
本步骤中，在配电网故障概率模型确定后，应用蒙特卡洛模拟方法生成故障场景。
[0080]
步骤130：在所述配电网故障场景集合中进行拓扑重构建立故障拓扑重构模型；
[0081]
衡量不同配电网故障场景的相似度是进行故障场景聚类的先决条件。针对配电网的灾害场景进行了聚类，其中故障场景的相似度以配电网各个节点的故障损失作为依据，相似的故障场景造成的各个节点的损失也相似。这一出发点是从故障影响角度出发，能一定程度上反映场景相似度。然而，在配电网的故障场景中，仅从故障损失角度进行聚类存在
不足。配电网恢复时，根据故障情况会产生不同的拓扑重构结果，而拓扑结构影响电源与负荷供电路径，故障严重时会划分为孤岛，因而拓扑结构对具体的恢复决策影响很大。不同的拓扑重构结果可能导致相似的故障损失，而这会产生本质上不同的配电网韧性恢复决策，因而本发明实施例中，在进行聚类前需要先对故障场景集合建立配电网拓扑重构结构。
[0082]
步骤140：在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度；
[0083]
本发明实施例从场景故障损失与场景拓扑相似度两方面衡量场景的相似度，以场景的故障损失及场景间的拓扑相似度作为基础进行场景聚类。
[0084]
步骤150：根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景；
[0085]
本发明实施例中，进行故障损失聚类时，先将可通过配电网重构消除故障影响或无故障的场景归为一类，减少后续聚类的计算量。
[0086]
步骤160：在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景；
[0087]
本发明实施例中，拓扑相似度聚类采用层次聚类方法，层次聚类法适合在给定样本间距离的情况下使用，其基本思想是，首先合并相似度最高或距离最近的两个样本为一类，并按照此原则不断对剩余样本进行合并，直至最终所有样本合为一类或达到终止条件。
[0088]
步骤170：根据所述第二场景聚类结果得到配电网灾前故障场景模拟预测结果。
[0089]
本发明实施例中，配电网为拓扑结构，包括多条输电线路，多个节点。
[0090]
具体地，如图2所示，步骤120具体包括：
[0091]
步骤121：将预设灾害信息输入所述配电网故障概率模型，得到所述配电网故障概率模型输出的配电网各条输电线路故障的概率，其中，所述预设的灾害信息包括灾害持续时间；
[0092]
灾害信息包括台风强度、线路结冰概率、持续时间等。
[0093]
步骤122：根据各输电线路故障概率，对各输电线路进行灾害持续时间范围内的一次故障场景随机生成，得到一个故障场景；
[0094]
步骤123：重复进行随机生成的过程，得到预设生成场景总数的故障场景，形成所述配电网故障场景集合。
[0095]
根据预设生成场景总数进行与预设生成场景总数数量一致的k次模拟生成，一次模拟生成的场景是灾害持续时间t的时长范围内的各线路的故障情况；
[0096]
将每次灾害持续时间t划分为t个时段，在每个时段内根据故障概率模型计算配电网每条线路的故障概率，t≥0，t≥0；在每个时段内对每条线路生成0
‑
1范围的随机数，若随机数小于该时段计算出的线路故障概率，且该时段之前该线路没有故障，则认为此时段该线路故障；
[0097]
重复上述每个时段的故障场景生成过程t次，完成每次灾害持续时间t的一次故障场景生成过程；
[0098]
重复上述一次故障场景生成过程k次，完成预设生成场景总数数量的故障场景生成。
[0099]
本发明实施例中，故障拓扑重构模型使用来自主网的电力进行自主恢复，因此步骤130中故障拓扑重构模型的目标函数为：
[0100][0101]
其中，表示配电网中第i节点t时刻的电荷恢复量，ω
i
表示第i节点的负荷重要度，表示配电网中所有节点集合，t表示故障场景中灾害持续时间。
[0102]
本发明实施例中，步骤140中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，具体包括：
[0103]
在步骤130建立的故障拓扑重构模型进行求解，得到配电网在单个故障场景中单个节点的负荷损失量为：
[0104][0105]
其中，ω
i
表示第i节点的负荷重要度，表示配电网中第i节点t时刻的电荷恢复量，表示配电网中第i节点电荷恢复量上限，表示配电网中所有节点集合；
[0106]
单个故障场景下故障损失表示为节点故障损失向量p
loss
：
[0107][0108]
本发明实施例中，步骤140中得到各故障场景之间的场景拓扑相似度，具体包括：
[0109]
为了衡量配电网故障重构后的拓扑相似度，本发明实施例应用树编辑距离(ted，tree edit distance)进行计算。这一概念指一个树状结构的图转化为另一树状结构的图的最小操作代价。在配电网领域，将这一概念应用在配电网规划领域，本发明实施例中，应用ted衡量配电网故障场景的拓扑相似度。
[0110]
两故障场景t1，t2之间的ted为：
[0111][0112]
其中，dist(t1，t2)表示t1，t2之间的ted，z＝{z1，z2，
…
，z
n
}表示将故障场景从t1转换为t2的一系列树编辑操作，|z|表示序列的操作总数，表示可完成故障场景t1，t2转换的操作序列集合，c(z
i
)表示操作z
i
的代价，优选地每次操作的代价均设置为1。
[0113]
对于任意树t1，t2，t3，ted应满足如下性质：
[0114]
(1)dist(t1，t2)≥0且dist(t1，t1)＝0；
[0115]
(2)dist(t1，t2)＝dist(t2，t1)；
[0116]
(3)dist(t1，t3)≤dist(t1，t2) dist(t2，t3)
[0117]
本步骤中在进行拓扑相似度计算时，配电网某一场景的拓扑结构通过步骤130建立的故障拓扑重构模型计算得到。本发明实施例中，在故障拓扑重构模型构建中引入了“虚拟节点”概念，即使配电网经过重构后形成多个孤岛，孤岛的根节点也需与“虚拟节点”相连，因此所有配电网重构结果均可生成一个以“虚拟节点”为根节点的树，方便进行拓扑相
似度计算。因此本发明实施例中，采用zhang
‑
shasha算法进行ted计算，该算法可以通过在python平台调用zss库使用。
[0118]
本发明实施例中，通过步骤140的计算结果进行场景聚类的过程，分别包括按照计算的故障损失进行第一阶段聚类及按照拓扑相似度进行第二阶段聚类。
[0119]
第一阶段聚类，如图3所示，步骤150是具体包括：
[0120]
步骤151：将通过配电网拓扑重构可消除的故障场景以及无故障场景归为一类第一场景聚类结果；
[0121]
本步骤中，故障损失聚类时，首先将可通过配电网重构消除故障影响或无故障的场景归为一类，减少后续计算量。
[0122]
对剩余的故障场景聚类时本发明实施例采用k均值(k
‑
means)聚类算法：
[0123]
步骤152：在剩余的故障场景中选取k个聚类中心为初始聚类中心，k为大于零的正整数；
[0124]
步骤153：以各初始聚类中心为中心分别进行至少一次故障场景聚集，将故障场景聚集后的各结果设为中间聚类结果；
[0125]
步骤154：计算各中间聚类结果的中心：当各中间聚类结果的中心分别与对应的初始聚类中心一致时，停止聚类，各中间聚类结果分别为所述第一场景聚类结果；当各中间聚类结果的中心分别与对应的初始聚类中心不一致时，继续聚类，直至各中间聚类结果的中心分别与对应的初始聚类中心一致。
[0126]
在步骤152
‑
154中进行聚类时，各场景中的距离指的是欧式距离(euclidean metric)，这一距离要求输入数据为连续数据，故障损失符合k
‑
means对输入数据为连续数据的要求。本发明实施例中，k
‑
means聚类算法中的k值需要提前设定，在步骤160中需要应用ted计算场景拓扑距离，因此希望第一场景聚类结果中包括故障场景每类场景的数量较小，以减小后续计算量，因此在选取k值时根据不同k值的聚类效果，在较好的聚类结果中选取一个较大的k值进行聚类。在应用k
‑
means聚类算法时，可以在python平台调用scikit
‑
learn库使用。
[0127]
第二阶段聚类，如图4所示，步骤160具体包括：
[0128]
步骤161：确定第一场景聚类结果中各故障场景的类间距离；
[0129]
通过ted计算k
‑
means聚类算法后第一场景聚类结果内场景两两之间的距离，可以得到一个类内的距离矩阵ed：
[0130][0131]
(8)式中n为第二数量，矩阵中的元素为对应位置的灾后场景间的ted，根据ted性质，该矩阵为一对称矩阵，实际只需计算n(n
‑
1)/2次即可。
[0132]
得到距离矩阵后，本发明实施例中的拓扑相似度聚类采用层次聚类方法。层次聚类法适合在给定样本间距离的情况下使用，其基本思想是，首先合并相似度最高或距离最近的两个样本为一类，并按照此原则不断对剩余样本进行合并，直至最终所有样本合为一
类或达到终止条件。进行层次聚类时需要选取出差异尽量大的若干代表场景，因此终止条件设置为：给定聚层次聚类数，使用ward法进行层次聚类。
[0133]
步骤162：确定第三数量的数值；
[0134]
确定第三数量即为在每个第一景聚类结果中要得到第三数量类的第二场景聚类结果。
[0135]
步骤163：按照第一场景聚类结果中各故障场景每次合并形成类簇时，选取故障场景合并后的类簇与其他类簇类间距离最大的方式进行聚类，形成第三数量的第二场景聚类结果。
[0136]
本步骤应用的ward法层次聚类的类间距离与常用的single，complete，average等方法不同，其引入下式定义类间距离，作为判断类簇合并的依据：
[0137][0138]
(9)式中u为类簇s和t合并后的类簇，v为其他未使用的类簇。|*|表示类簇内的样本个数，s＝|v| |s| |t|。每次合并类簇时，选取合并后与其他类簇距离最大的类簇。这本质上是在衡量合并时误差平方和(ess,error sum of squares)的变化情况，其定义如下式：
[0139][0140]
(10)式中x
i
为样本的值，所有样本整体的ess为样本各类簇的ess之和。根据ward法进行层次聚类时，每次合并的类簇是使得整体ess增加最少的两个类簇，这实质上是在要求合并的两个类簇合并后类内样本要尽量相似，类且与其他类簇尽量不相近。这一方法能有效得到互相之间差异最大的类簇，体现场景的代表性。
[0141]
层次聚类终止时的类簇数量对于聚类效果影响很大，在单次聚类时，可根据需求和经验设置类簇数量，而在本发明实施例中，需对故障损失聚类后每个类进行层次聚类，对每个类单独设置类簇数量并不现实。因此本发明实施例中，层次聚类的聚类数使用l法确定：
[0142]
通常认为较好的聚类效果较好时，聚类数量和聚类的类间距离在拐点附近效果较好。l法对“类簇数
‑
类间距离”图上靠近左侧和右侧的点分别进行线性回归进行拟合，综合拟合误差最小时得到最终的类簇数。综合拟合误差定义如下式：
[0143][0144]
c
*
＝argmin
c
rmse
c
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12)
[0145]
(11)式中c表示“类簇数
‑
类间距离”图进行左侧和右侧划分的类簇数；b表示最大类簇数；rmse(l
c
),rmse(r
c
)分别表示左侧直线和右侧直线拟合的根均方误差，其中直线l
c
包含类数为2至c的点，直线r
c
包含类数c 1至b的点。在求(12)式中的最优聚类数c
*
时，分别将c设置为3至b
‑
2进行计算，即可得到最终的c
*
。
[0146]
本发明实施例中，步骤170具体包括：在步骤160得到的各个第二场景聚类结果中，随机选取一个故障场景作为该类的故障场景，并将该类第二场景聚类结果中的故障场景数量与预设生成场景总数的比值作为该类第二场景聚类结果的生成概率。
[0147]
为验证聚类方法效果，设置以下算例进行测试。本算例使用的配电网拓扑如图5所示，节点1连接电网主网，图中虚线为配电网联络线，配电网中，节点15、节点20有电动公交充电站，节点9、节点18、节点29有分布式电源。
[0148]
在场景生成阶段，假定台风灾害持续72小时，模拟台风过境：风力从2级变化至15级，再从15级变化回4级。共生成场景1000个。配电网韧性恢复持续6个时段，每个时段1h。
[0149]
故障损失聚类：
[0150]
对生成的1000个场景依次计算故障损失，得到其中216个场景重构后无损失，可进行合并；另外784个场景需要进行故障损失聚类。取k
‑
means算法聚类数为11，得到聚类结果。为了检验聚类效果，使用t
‑
sne方法对聚类结果进行可视化。t
‑
sne是一种非线性降维算法，可以有效将高维数据降维至2
‑
3维空间，观察聚类效果。
[0151]
拓扑聚类：
[0152]
根据故障损失聚类结果，对故障损失每一类进行层次聚类。以故障损失第9类场景的层次聚类结果为例，结果如图6所示。
[0153]
如图6可看出，第9类故障场景共有30个场景，并根据l法最终聚为4类，从每类中随机选取一个场景，得到第9类的聚类结果。图中每类代表场景对应概率自左向右分别为1.1％，0.6％，0.3％，1.0％，由聚类场景数除以场景总数1000得到。其他类别的故障场景也进行同样的聚类处理。
[0154]
最终，1000个场景聚类结果如下：无故障损失场景216个，聚类为1个场景；有故障损失场景784个，聚类为100个场景，每一类的聚类结果汇总与表1中。
[0155]
表1有故障损失场景聚类结果汇总
[0156][0157]
由表1可看出1000个场景最终聚类为101个场景。
[0158]
下面对本发明实施例提供的配电网灾前故障场景模拟预测装置进行描述，下文描述的配电网灾前故障场景模拟预测装置与上文描述的配电网灾前故障场景模拟预测方法可相互对应参照，如图7所示，本发明实施例提供的一种配电网灾前故障场景模拟预测装置，包括：
[0159]
故障概率建模单元710，用于建立配电网故障概率模型；
[0160]
由于配电网故障概率模型的建模在现有技术中已有研究，本发明实施例中，在现有故障概率模型的基础上进行一定简化，因此本发明实施例中，主要考虑台风这类灾害天气对配电网产生的影响。配电网元件中，假设故障由配电网线路中的杆塔造成，线路上的杆塔受风速和结冰影响而故障。
[0161]
故障场景生成单元720，用于基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合；
[0162]
在配电网故障概率模型确定后，应用蒙特卡洛模拟方法生成故障场景。
[0163]
故障拓扑重构单元730，用于在所述配电网故障场景集合中对每个故障场景建立故障拓扑重构模型；
[0164]
故障损失及场景拓扑相似度处理单元740，用于在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度；
[0165]
衡量不同配电网故障场景的相似度是进行故障场景聚类的先决条件。针对配电网的灾害场景进行了聚类，其中故障场景的相似度以配电网各个节点的故障损失作为依据，相似的故障场景造成的各个节点的损失也相似。这一出发点是从故障影响角度出发，能一定程度上反映场景相似度。然而，在配电网的故障场景中，仅从故障损失角度进行聚类存在不足。配电网恢复时，根据故障情况会产生不同的拓扑重构结果，而拓扑结构影响电源与负荷供电路径，故障严重时会划分为孤岛，因而拓扑结构对具体的恢复决策影响很大。不同的拓扑重构结果可能导致相似的故障损失，而这会产生本质上不同的配电网韧性恢复决策，因而本发明实施例中，在进行聚类前需要先对故障场景集合建立配电网拓扑重构结构。
[0166]
第一聚类单元750，用于根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景；
[0167]
本发明实施例中，进行故障损失聚类时，先将可通过配电网重构消除故障影响或无故障的场景归为一类，减少后续聚类的计算量。
[0168]
第二聚类单元760，用于在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景；
[0169]
本发明实施例中，拓扑相似度聚类采用层次聚类方法，层次聚类法适合在给定样本间距离的情况下使用，其基本思想是，首先合并相似度最高或距离最近的两个样本为一类，并按照此原则不断对剩余样本进行合并，直至最终所有样本合为一类或达到终止条件。
[0170]
输出单元770，用于根据所述第二场景聚类结果得到配电网灾前故障场景模拟预测结果。
[0171]
本发明实施例中，故障场景生成单元720具体包括：
[0172]
预设信息确定子单元，用于将预设灾害信息输入所述配电网故障概率模型，得到所述配电网故障概率模型输出的配电网各条输电线路故障的概率，其中，所述预设的灾害信息包括灾害持续时间；
[0173]
单次场景生成子单元，用于根据各输电线路故障概率，对各输电线路进行灾害持续时间范围内的一次故障场景随机生成，得到一个故障场景；
[0174]
循环子单元，用于控制单次场生成子单元重复进行随机生成的过程，得到预设生成场景总数的故障场景，形成所述配电网故障场景集合。
[0175]
本发明实施例中，第一聚类单元750具体包括：
[0176]
无故障聚类子单元，用于将通过配电网拓扑重构可消除的故障场景以及无故障场
景归为一类第一场景聚类结果；
[0177]
聚类中心确定子单元，用于在剩余的故障场景中选取k个聚类中心为初始聚类中心，k为大于零的正整数；
[0178]
中间聚类子单元，用于以各初始聚类中心为中心分别进行至少一次故障场景聚集，将故障场景聚集后的各结果设为中间聚类结果；
[0179]
聚类终止子单元，用于计算各中间聚类结果的中心：当各中间聚类结果的中心分别与对应的初始聚类中心一致时，停止聚类，各中间聚类结果分别为所述第一场景聚类结果；当各中间聚类结果的中心分别与对应的初始聚类中心不一致时，继续聚类，直至各中间聚类结果的中心分别与对应的初始聚类中心一致。
[0180]
本发明实施例中，第二聚类单元760具体包括：
[0181]
类间距确定子单元，用于确定第一场景聚类结果中各故障场景的类间距离；
[0182]
类簇数确定子单元，用于确定第三数量的数值；
[0183]
类簇聚类子单元，用于按照第一场景聚类结果中各故障场景每次合并形成类簇时，选取故障场景合并后的类簇与其他类簇类间距离最大的方式进行聚类，形成第三数量的第二场景聚类结果。
[0184]
下面结合图8描述本发明实施例提供的一种电子设备的实体结构示意图，如图8所示，该电子设备可以包括：处理器(processor)810、通信接口(communications interface)820、存储器(memory)830和通信总线840，其中，处理810，通信接口820，存储830通过通信总线840完成相互间的通信。处理器810可以调用存储器830中的逻辑指令，以执行配电网灾前故障场景模拟预测方法，该方法包括：建立配电网故障概率模型；基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合；在所述配电网故障场景集合中进行拓扑重构建立故障拓扑重构模型；在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度；根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景；在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景；根据所述第二场景聚类结果得到配电网灾前故障场景模拟预测结果。
[0185]
此外，上述的存储器830中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0186]
另一方面，本发明实施例还提供一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法所提供的配电网灾前故障场景模
拟预测方法，该方法包括：建立配电网故障概率模型；基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合；在所述配电网故障场景集合中进行拓扑重构建立故障拓扑重构模型；在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度；根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景；在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景；根据所述第二场景聚类结果得到配电网灾前故障场景模拟预测结果。
[0187]
又一方面，本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各提供的配电网灾前故障场景模拟预测方法：建立配电网故障概率模型；基于所述配电网故障概率模型，根据预设灾害信息及预设生成场景总数，生成配电网故障场景集合；在所述配电网故障场景集合中进行拓扑重构建立故障拓扑重构模型；在所述配电网故障场景集合中根据所述故障拓扑重构模型处理得到每个故障场景的故障损失，以及各故障场景之间的场景拓扑相似度；根据各个故障场景的故障损失按照第一预设规则进行聚类，得到第一数量的第一场景聚类结果，其中每个所述第一场景聚类结果中包括第二数量个故障场景；在每个第一场景聚类结果中根据各故障场景之间的场景拓扑相似度按照第二预设规则进行聚类，得到第三数量的第二场景聚类结果，其中每个所述第二场景聚类结果包括第四数量个故障场景；根据所述第二场景聚类结果得到配电网灾前故障场景模拟预测结果。
[0188]
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
[0189]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
[0190]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。