一种搜索圆半径可控的电力系统小干扰特征值计算方法

摘要

本发明公开了属于电力系统运行与控制技术领域的一种搜索圆半径可控的电力系统小干扰特征值计算方法。设置搜索圆圆心、搜索圆半径等参数；生成电力系统小干扰分析的系数矩阵A，进行位移逆变换后得到矩阵T_u；对矩阵T_u进行k步的Arnoldi分解；将k步的Arnoldi分解扩展到m步的Arnoldi分解；筛选H_m的特征值，经隐式位移QR处理后，重新得到k步的Arnoldi分解；重启动迭代得到搜索圆内的所有特征值。本发明提出的搜索圆半径可控的电力系统小干扰特征值计算方法，通过动态增加特征值计算个数，将可控参数由特征值个数转换为搜索圆半径，能够自动计算出指定位置和大小的搜索圆内部所有特征值。

主权项

一种搜索圆半径可控的电力系统小干扰特征值计算方法，其特征在于，包括以下步骤：1)设置搜索圆圆心u、搜索圆半径r、特征值个数k的初始值、Krylov子空间维数m、特征值个数动态增加量Δk、收敛标准tol、初始向量v₀；2)电力系统经潮流计算获取稳态运行信息，初始化电力系统小干扰分析动态元件参数，利用电力系统小干扰分析动态元件参数生成电力系统小干扰分析的系数矩阵A，将电力系统小干扰分析的系数矩阵A进行位移逆变换，得到矩阵T_u；3)对矩阵T_u进行k步的Arnoldi分解，k步的Arnoldi分解为：T_uV_k＝V_kH_k+f_ke_k^T   (1)式(1)中，T_u为系数矩阵A位移逆变换后的矩阵；矩阵V_k∈C^n╳k，其列向量为Krylov子空间Span{v₀，Av₀，A²v₀，…，A^k‑1v₀}的一组正交基，C^n╳k表示n行k列的矩阵，v₀为初始向量，v₀∈C^k╳1＝(0,0,…,0,1)^T，C^k╳1表示k行1列的矩阵；H_k为Hessenberg矩阵，H_k＝V_k^TAV_k∈C^k╳k，C^k╳k表示k行k列的矩阵；f_k是维数为k的残差向量；e_k是第k个元素为1的k维标准正交基；T表示转置；4)将k步的Arnoldi分解扩展到m步的Arnoldi分解，T_uV_m＝V_mH_m+f_me_m^T   (2)式(2)中，T_u为系数矩阵A位移逆变换后的矩阵；矩阵V_m∈C^n╳m，其列向量为Krylov子空间Span{v₀，Av₀，A²v₀，…，A^m‑1v₀}的一组正交基，C^n╳m表示n行m列的矩阵，v₀为初始向量，v₀∈C^k╳1＝(0,0,…,0,1)^T，C^k╳1表示k行1列的矩阵；H_m为Hessenberg矩阵，H_m＝V_m^TAV_m∈C^m╳m，C^m╳m表示m行m列的矩阵；f_m是维数为m的残差向量；e_m是第m个元素为1的m维标准正交基；T表示转置；5)筛选Hessenberg矩阵H_m的m‑k个模值最小的特征值，经隐式位移QR处理后，重新得到k步的Arnoldi分解；6)判断残差||f⁺_k||是否小于tol，如果残差||f⁺_k||小于tol，执行步骤7)，否则，返回步骤4)重启动；7)将矩阵T_u的特征值变换成电力系统小干扰分析的系数矩阵A的特征值；8)计算实际搜索圆半径，判断实际搜索圆半径是否大于搜索圆半径r，若实际搜索圆半径大于搜索圆半径r，系数矩阵A的特征值即为搜索圆内的所有特征值，计算结束；否则，令k’＝k+Δk作为新的特征值个数，m’＝m+Δm作为新的Krylov子空间维数，返回步骤4)重启动。