[发明专利]一种基于半监督谱聚类的最优解列断面搜索方法

权利要求书

1.一种基于半监督谱聚类的最优解列断面搜索方法，其特征在于，包括以下步骤：

步骤1：构建电力系统无向边权图G，并对其中所有的边进行赋权，构建图G的邻接权矩阵W，W的元素取为wij=|Pij|Dij,]]>其中，|Pij|=|pij|+|pji|2,]]>考虑网损影响，p_ij表示支路ij两端从节点i流向节点j的有功潮流，D_ij为支路ij两端节点之间的电气距离；

步骤2：计算邻接权矩阵W和度矩阵A，其元素分别为：

Wij=|Pij|Dij,(i,j)∈E00,(i,j)∉E0]]>    (式壹)

Aij=Σj=1,j≠inWij,i=j0,i≠j]]>    (式贰)

其中，E₀为图G的边集；

步骤3：执行第一次判断，判断系统是否发生大扰动？

若是，则执行第二次判断，判断系统是否已失稳？

若是，则顺序执行下述步骤4，第二次判断执行完毕；

若否，则回转执行所述的步骤3，第二次判断执行完毕；

若否，则执行第三次判断，判断是否到潮流采集时间？

若是，则回转执行所述的步骤1，第三次判断执行完毕；

若否，则回转执行所述的步骤3，第三次判断执行完毕；

步骤4：基于发电机分群信息，依据发电机同调约束或分离约束，利用式叁对邻接权矩阵W进行调整并更新度矩阵A：

W(nGi,nGj)=∞,(nGi,nGj)∈Must-link0,(nGi,nGj)∈Cannot-link]]>    (式叁)

其中，同调约束表示任意属于同调机群的两台发电机之间至少存在一条连通路径，分离约束表示任意两台属于非同调机群的发电机之间不存在任何连通途径，n_Gi和n_Gj表示不同的发电机节点，若发电机节点对(n_Gi,n_Gj)∈Must-link约束，则W(n_Gi,n_Gj)＝∞，表示任意两台同调的发电机间至少存在一条连通的路径；相反，若发电机节点对(n_Gi,n_Gj)∈Cannot-link约束，则W(n_Gi,n_Gj)＝0，表示任意两台非同调发电机间不存在任何连通路径；

步骤5：利用邻接权矩阵W和度矩阵A计算规范化拉普拉斯矩阵L_N＝A^-1L，其中，L为未规范化拉普拉斯矩阵，L＝A-W；

步骤6：求解方程L_NX＝λX的前k个最小的特征值所对应的特征向量X₁，X₂，…，X_k；

步骤7：令X_i(i＝1,2,…,k)∈R^n×k构成矩阵Y，其行向量y_j∈R^k(j＝1,2,…,n)对应系统的各节点的编号；

步骤8：利用改进的PAM算法对行向量y_j∈R^k(j＝1,2,…,n)进行优化聚类以获得最优解列断面，并以此为解列具体方案将原系统解列为k个相互独立的孤岛：G₁，G₂，…，G_k，实现式肆所示的最小复合有功潮流冲击，从而确定最优主动解列断面；

P~sum=minV1,V2,...,Vk⋐V(12Σi∈Vs,j∈Vt|Pij|Dij)(s,t=1,2,...,k,s≠t)]]>    (式肆)

式中：为最小复合有功潮流冲击；V₁，V₂，…，V_k表示电网遭受大干扰后形成k个同调机群。

2.根据权利要求1所述的基于半监督谱聚类的最优解列断面搜索方法，其特征在于：步骤8中所述的改进的PAM算法，包括以下子步骤：

步骤8.1：计算对象y_j∈R^k(j＝1,2,…,n)相互间的欧式距离矩阵S并建立存储索引表，设置迭代次数N，对于每一个对象y_j∈R^k(j＝1,2,…,n)，计算其对应的v_j，表达式为：

vj=Σp=1nSpjΣl=1nSpl(j=1,2,...,n)]]>   (式伍)

以k个最小的v_j对应的节点作为初始中心点；

步骤8.2：任选一个非中心点作为新的中心点，通过查询距离矩阵S，将所有对象分配到距离最小的中心点归为一类，并为所有对象对当前和新的中心点计算平方误差和；

步骤8.3：判断前后两次平方误差和大小，若后者较小，则用新的中心点替换当前中心点；

步骤8.4：判断是否达到迭代次数N？

若是，则形成最优的K个中心点，平方误差和保持不变，完成聚类；

若否，则返回执行所述的步骤8.2。