发明内容
本发明的目的是提供一种基于SOM-FAT的电力系统连锁故障预测方法,能避免发生系统范围的大停电事件。
本发明所采用的技术方案是,一种基于SOM-FAT的电力系统连锁故障预测方法,包括:建立初始故障支路集,检测初始故障支路集中的支路是否为故障支路,若为故障支路,记录该支路路径,并计算发生概率,若为正常支路,计算该支路与其他支路间的关联性大小,并建立下级故障支路集合;依次检测下级故障支路集合中的支路是否为故障支路,若为故障支路,记录该支路路径,并计算发生概率,若为正常支路,计算该支路与其他支路间的关联性大小;完成对初始故障支路集中的每条支路的检测,得到故障支路集。
本发明的特点还在于:
包括以下步骤:
步骤1、设计故障集评价指标,选择故障集评价指标较大的支路组成初始故障支路集;
步骤2、断开初始故障支路集中的支路,判断系统是否解列失稳,如果解列失稳,则记录该条支路路径,并计算发生概率;若系统稳定,则根据上下级支路关联性指标,计算该支路与其他支路间的关联性大小;
步骤3、对上下级支路关联性指标进行分类,将上下级支路关联性指标较高的支路作为下级故障支路集合,并计算每条下级支路发生故障的概率;
步骤4、依次断开下级故障支路集合中的支路,判断系统是否解列失稳,如果解列失稳,则记录该条支路路径,并计算发生概率;若系统稳定,则根据上下级支路关联性指标,计算该支路与其他支路间的关联性大小,直至搜索完下级故障支路集合中的每条支路;
步骤5、重复步骤2-4,完成初始故障支路集中的每条支路检测,得到故障支路集。
步骤1中故障集评价指标El为:
El=η1Cl×η2Tl×η3λl (1);
上式中,Cl为节点扰动对支路的潮流冲击性,Tl为支路在电力系统中的重要性,λl为支路自身故障概率,η1、η2、η3均为阈值。
节点扰动对支路的潮流冲击性Cl计算过程如下:
系统处于平衡状态时支路的潮流为Pl0,当某节点a负荷突增,支路1的潮流改变为Pla,则节点a对支路l的潮流冲击Ela为:
ΔEla=Pla-Pl0 (2);
则节点a对系统的潮流冲击ΔEa为:
上式中:S为系统支路数量;
则支路l的潮流冲击率ηla为:
结合公式1-3,节点a在系统的潮流分布熵HD(a)为:
所以,节点a扰动对支路1的潮流冲击Cla为:
步骤2中上下级支路关联性指标Ml为:
Ml=Sl×Dl×Bl (9);
上式中,Bl为上级支路潮流转移熵对下级支路冲击性,Sl为上级支路故障切除后对下级支路的影响度,表示如下:
上式中,Fl j为支路l在故障发生后的传输功率,Fl j-1为支路1在故障发生前的传输功率,为支路1-1在故障发生前的传输功率;
Dl为支路在电网运行中的负载率指标,表示如下:
上式中,Flmax为支路l的最大传输功率。
上级支路潮流转移熵对下级支路冲击性Bl的计算过程如下:
支路l-1断开对支路l传输裕度影响的潮流转移熵Hl,l-1为:
Hl,l-1=-ηl,l-1lnηl,l-1 (12);
上式中,ηl,l-1为支路l-1开断对支路l的传输裕度影响值;
则支路l-1潮流转移熵对于支路l的冲击性Bl为:
步骤3具体包括:
步骤3.1、采用生成随机数的方式对输入层的上下级支路关联性指标Ml和竞争层节点的值Wj(j=1,2,…,N)进行初始化,随机数的取值范围为[0,1],如下式:
并对学习率η(0)以及迭代次数进行初始化,学习率的取值范围为0<η(0)<1;
步骤3.2、寻找获胜节点;
步骤3.3、确定优胜邻域,并通过下式调整获胜节点及优胜邻域内节点的权值:
上式中,dj表示优胜邻域内节点与获胜节点的欧氏距离,N(dj)为邻域函数,
步骤3.3、循环学习直至迭代结束,获得聚类结果;
步骤3.4、将聚类结果中上下级支路关联性指标最高的一类支路作为下级故障支路集合,并计算每条下级支路发生故障的概率。
优胜邻域为:
r(t)=Ce-t/T (19);
上式中,r(t)表示t次迭代时邻域范围的半径,C为竞争层维度的1/2,T为总迭代次数。
步骤3.2寻找获胜节点的方法为欧氏距离法或余弦法。
步骤3.2中采用自适应的方式来调整学习率,学习率与迭代次数成反比:
η(t+1)=η(0)e-t/T (22)。
本发明的有益效果是:
本发明的基于SOM-FAT的电力系统连锁故障预测方法,在传统故障树预测的基础上,以节点扰动对支路的潮流冲击性、支路在电力系统中的重要性、支路自身故障概率为初始故障集评价指标,以上级支路潮流转移熵对下级支路冲击性、上级支路故障切除后对下级支路的影响度及支路在电网运行中的负载率为上下级支路关联性指标,能更准确的确定初始故障集、下级故障支路,得到故障支路集,以避免发生系统范围的大停电事件;对SOM聚类算法进行改进,提出用自适应的方式来解决算法中学习率难以确定的问题,使得聚类结果更加准确。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种基于SOM-FAT的电力系统连锁故障预测方法,包括:
建立初始故障支路集,检测初始故障支路集中的支路是否为故障支路,若为故障支路,记录该支路路径,并计算发生概率,若为正常支路,计算该支路与其他支路间的关联性大小,并建立下级故障支路集合;依次检测下级故障支路集合中的支路是否为故障支路,若为故障支路,记录该支路路径,并计算发生概率,若为正常支路,计算该支路与其他支路间的关联性大小;完成对初始故障支路集中的每条支路的检测,得到故障支路集。
如图1所示,具体包括以下步骤:
步骤1、设计故障集评价指标,建立初始故障支路集;
故障集评价指标El为:
El=η1Cl×η2Tl×η3λl (1);
上式中,Cl为节点扰动对支路的潮流冲击性,Tl为支路在电力系统中的重要性,λl为支路自身故障概率,λl可查阅《电力系统安全运行导则》得到,η1、η2、η3分别为三个指标的阈值,可根据支路实际运行情况,设定适当的阈值。综合三个指标因素,选取故障集评价指标El较大的支路组成初始故障支路集。
节点扰动对支路的潮流冲击性Cl,该指标考虑电网节点由于负荷突然变化造成的扰动对于电网的冲击问题,计算过程如下:
系统处于平衡状态时支路的潮流为Pl0,当某节点a负荷突增,引起系统扰动,支路l的潮流改变为Pla,节点a受到扰动后,则节点a对支路l的潮流冲击Ela(即支路l的潮流增加量)为:
ΔEla=Pla-Pl0 (2);
则节点a对系统的潮流冲击ΔEa为:
上式中:S为系统支路数量;
则采用支路l的潮流冲击率ηla,表示支路l在节点a对系统的潮流冲击中所占的比例:
结合公式1-3,节点a在系统的潮流分布熵HD(a),HD(a)反映了节点在发生扰动后系统所受到的潮流冲击的分布特性,表示如下:
所以,支路l在系统节点受到扰动后,自身所承担的潮流冲击的大小Cla(即节点a扰动对支路l的潮流冲击)为:
支路在电力系统中的重要性Tl计算过程如下:
一般情况下,在发电机Gm向负荷Ln传送功率的过程中,支路Pmn承担的部分越多,表明该支路在电力系统中占有比较重要的位置,表明其在电力系统中的重要性比较高,定义单一发电负荷对(m,n)的支路潮流介数为:
式中:Pij(m,n)为总传输功率Pm,n在支路i,j上的分量,Pij(m,n)/Pmn为支路i,j对发电负荷对(m,n)的输电贡献比例因子,min(Sm,Sn)为该发电负荷对的权重,该权重取实际发电量和实际负荷消耗量中的较小值,表征了输电支路所传输的最大可用功率。
由于电网中每条支路都承担了不同发电负荷对之间的功率传送任务,从而一条支路可能被多个发电负荷对的功率传输路径经过,其被经过的次数也反映了该支路在电网当前运行方式下的重要程度,即电网中所有发电负荷对的累加效果,定义支路的重要性指标如下:
式中:G为支路提供功率的所有发电机节点集合,L为消耗支路功率的所有负荷节点集合。
步骤2、断开初始故障支路集中的支路,判断系统是否解列失稳,如果解列失稳,则记录该条支路路径,并计算发生概率;若系统稳定,则根据上下级支路关联性指标,计算该支路与其他支路间的关联性大小;
上下级支路关联性指标Ml为:
Ml=Sl×Dl×Bl (9);
上式中,Bl为上级支路潮流转移熵对下级支路冲击性,Sl为上级支路故障切除后对下级支路的影响度,即在上级故障支路l-1发生故障前后,下级支路l的传输功率变化与上级支路l-1在故障前的传输功率的比值表示如下:
上式中,Fl j为支路l在故障发生后的传输功率,Fl j-1为支路l在故障发生前的传输功率,为支路l-1在故障发生前的传输功率;
Dl为支路在电网运行中的负载率指标,定义为支路l在支路l-1发生故障且将要传播到下一级时,其传输功率与其所承担的最大传输功率的比值,表示如下:
上式中,Flmax为支路l的最大传输功率。
上级支路潮流转移熵对下级支路冲击性Bl的计算过程如下:
支路l-1断开对支路l传输裕度影响的潮流转移熵Hl,l-1为:
Hl,l-1=-ηl,l-1lnηl,l-1 (12);
上式中,ηl,l-1为支路l-1开断对支路l的传输裕度影响值,表示如下:
则支路l-1潮流转移熵对于支路l的冲击性Bl(即支路潮流转移熵冲击性)为:
步骤3、对上下级支路关联性指标进行分类,将上下级支路关联性指标较高的支路作为下级故障支路集合,并计算每条下级支路发生故障的概率;
步骤3.1、采用生成随机数的方式对输入层的上下级支路关联性指标Ml和竞争层节点的值Wj(j=1,2,…,N)进行初始化,随机数的取值范围为[0,1],如下式:
并对学习率η(0)以及迭代次数进行初始化,学习率的取值范围为0<η(0)<1;
步骤3.2、寻找获胜节点;
输入层的一个输入向量为一个样本数据。对于当前输入的样本数据,在竞争层中寻找与之最匹配的节点(即获胜节点,通常使用相近性来度量);采用欧氏距离法或余弦法寻找获胜点,欧氏距离法采用两个向量间的欧氏距离来度量,距离越小,相近性越高,欧氏距离法表示如下;
余弦法使用两个向量间的夹角来度量,夹角越小,余弦越大,相近性越高,表示如下:
竞争层的所有节点中除最匹配的节点输出为1,其余节点的输出均为0:
步骤3.3、确定优胜邻域(即优胜节点的一个邻域范围),该邻域范围的原则是邻域半径与迭代次数成反比,以在保证相邻节点差异性的同时保证一定的相似性,优胜邻域为:
r(t)=Ce-t/T (19);
上式中,r(t)表示t次迭代时邻域范围的半径,C为竞争层维度的1/2,T为总迭代次数;
并通过下式调整获胜节点及优胜邻域内节点的权值:
上式中,dj表示优胜邻域内节点与获胜节点的欧氏距离;N(dj)为邻域函数,
采用自适应的方式来调整学习率,学习率与迭代次数成反比:
η(t+1)=η(0)e-t/T (22);
步骤3.3、当前样本数据学习过程结束后,输入层传入下一个样本数据,迭代上述过程,直至t=T;或者定义最小的学习率min(η),当η(t)<min(η)时,结束迭代,聚类结束,获得聚类结果;
步骤3.4、将聚类结果中上下级支路关联性指标最高的一类支路作为下级故障支路集合,并计算每条下级支路发生故障的概率。
步骤4、依次断开下级故障支路集合中的支路,判断系统是否解列失稳,如果解列失稳,则记录该条支路路径,并计算发生概率;若系统稳定,则根据上下级支路关联性指标,计算该支路与其他支路间的关联性大小,直至搜索完下级故障支路集合中的每条支路;
步骤5、重复步骤2-4,完成初始故障支路集中的每条支路检测,得到故障支路集,预测结束。
通过以上方式,本发明的基于SOM-FAT的电力系统连锁故障预测方法,在传统故障树预测的基础上,以节点扰动对支路的潮流冲击性、支路在电力系统中的重要性、支路自身故障概率为初始故障集评价指标,以上级支路潮流转移熵对下级支路冲击性、上级支路故障切除后对下级支路的影响度及支路在电网运行中的负载率为上下级支路关联性指标,能更准确的确定初始故障集、下级故障支路;对SOM聚类算法进行改进,提出用自适应的方式来解决算法中学习率难以确定的问题,使得聚类结果更加准确。