[发明专利]一种最小切负荷计算方法

权利要求书

1.一种最小切负荷计算方法，其特征在于，所述方法包括以下步骤：

(1)利用蒙特卡洛法在搜索空间中选择随机初始点，在所述随机初始点中任取一个样本x₁作为第一个类别Γ₁的聚类中心c₁，x₁对应的增长方向b₁作为所述第一个类别Γ₁的增长方向矢量t₁；

(2)对第i个样本x_i，计算b_i与各已存在类别增长方向矢量t_i之间的夹角θ_ij(j＝1，2，...M)，M为已存在类别数目；

(3)确定与所述bi夹角最小的类别Γ_p及最小夹角θ_p；

(4)判断所述最小夹角θ_p是否满足最小相似度，如果是，执行步骤(5)；如果否，执行步骤(6)；

(5)将所述样本x_i归入所述Γ_p类，同时按照第一公式更新聚类中心c_p和增长方向矢量t_p后，执行步骤(2)，当所有样本计算完毕，结束聚类过程，执行步骤(7)；

(6)M＝M+1，将所述样本x_i作为一个新的类别Γ_M，t_M＝b_i，执行步骤(2)，当所有样本计算完毕，结束聚类过程，执行步骤(7)；

(7)以M个聚类中心组成的样本空间Φ₂＝{c₁，c₂，...，c_M}调用连续潮流来计算电压稳定临界点，获取所述电压稳定临界点；

(8)通过所述电压稳定临界点处的灵敏度信息以及分岔原理求取局部切平面法向量；

(9)将所述电压稳定临界点依据所述局部切平面法向量的相似性聚类成K类，Ψ＝{Ψ₁，Ψ₂，...，Ψ_K}；对每一类聚类后的所述电压稳定临界点分别使用最小二乘法进行超平面的拟合，获取稳定域超平面；

(10)定义一个与所述稳定域超平面平行，且沿负荷增长方向b与所述稳定域超平面距离为λ_cr的安全超平面，求解所述安全超平面的系数；

(11)根据所述安全超平面的系数计算最小切负荷量(Δy)；

其中，所述第一公式具体为：cp=1NpΣxj∈Γpxj,]]>tp=1NpΣxj∈Γpbj,]]>

N_p为纳入x_i后Γ_p类别总样本数。