[发明专利]基于集群的故障临界切除时间并行计算方法

摘要

基于集群的故障临界切除时间并行计算方法，根据故障的搜索区间长度分配计算核，确保搜索区间大的故障分到更多的核。如果故障的档位拆分步长小于计算精度，则按计算精度的要求来确定该故障所需要的核数。在采用绝对精度时，若档位拆分步长满足收敛精度，则采用较长的观察时间通过一次迭代直接求出暂态安全CCT和动态安全CCT；否则，先采用较短观察时间求取暂态安全CCT，再将暂态安全CCT作为上限，使用较长观察时间求取动态安全CCT。在采用相对精度时，先求取暂态安全CCT，再求取动态安全CCT。在考虑多种电力系统安全稳定裕度指标的情况下，求取功角稳定CCT、电压安全CCT、频率安全CCT和动态安全CCT中的最小者。

主权项

1.基于集群的故障临界切除时间并行计算方法，其特征在于包括下列步骤：1)确定暂态安全CCT和动态安全CCT的计算方式；2)根据故障i的故障切除时间变化量的上下限，获取其CCT搜索区间及搜索区间长度L_i；3)根据计算参数和搜索区间确定计算精度ξ_i；若计算精度采取相对值，则先计算出近似CCT，再根据收敛精度相对值来估计计算精度，筛选出搜索区间的1/2长度大于计算精度的故障，将各故障的方案数已确定标志置为0，标志为0的故障需要重新分配核，标志为1的故障方案数已确定；尚未分配的空闲核数为计算机集群当中的可用核数；4)将尚未分配的空闲核按故障搜索区间长度分配给方案数已确定标志为0的故障；设尚未分配的空闲核数为N_{idle_cpu}，故障数为M，故障的总搜索区间长度为则故障切除时间档位的理论拆分步长为故障i对应的方案数为如果Σi=1MNi,sch>Nidle_cpu,]]>则按搜索区间长度由小到大依次将N_i，sch减1至Σi=1MNi,sch=Nidle_cpu,]]>但N_i，sch不能为0以保证所有的故障的都参与计算，故障i的档位实际拆分步长为L_i/N_i，sch；5)如果存在L_i/(N_i，sch-1)＜2ξ_i的故障，则按2ξ_i的拆分步长确定该故障实际应分配的方案数N_i，sch，并将对应的方案数已确定标志置为1，尚未分配的空闲核数为空闲核数减去方案数已确定故障的方案数之和，执行步骤4)；6)根据各故障对应的N_i，sch将搜索区间等分，确定各计算方案对应的故障切除时间变化量，将不同的计算方案分发到不同的核进行并行计算，等待并收集稳定性评估结果；7)将裕度指标为正的方案号称为安全方案号，反之为不安全方案号；相邻的安全方案号和不安全方案号的对应的故障切除时间变化量为t_s和t_u，则该裕度指标对应的下一次迭代的搜索区间为[t_s，t_u]；如果不同裕度指标对应不同的[t_s，t_u]，则所关心的安全稳定CCT必然存在于[t_s，t_u]_min内，下一次CCT搜索的区间为[t_s，t_u]_min，同时给出其它裕度指标的近似CCT；如果所有故障搜索区间的1/2长度小于计算精度，则输出CCT计算结果；否则，执行步骤3)。