[发明专利]基于低阶IGD-IRK的时滞电力系统小干扰稳定性分析方法

说明书

本发明公开了基于低阶IGD‑IRK的时滞电力系统小干扰稳定性分析方法，包括：建立时滞电力系统数学模型；根据时滞电力系统状态变量是否与时滞有关，时滞电力系统的微分方程改写为与时滞有关的部分及与时滞无关的部分；获得与重组后的时滞电力系统状态方程相对应的无穷小生成元；针对无穷小生成元基于隐式龙格‑库塔法低阶离散化，得到无穷小生成元低阶离散化近似矩阵；针对无穷小生成元低阶离散化近似矩阵进行位移逆变化，得到无穷小生成元低阶离散化近似逆矩阵；针对上述逆矩阵进行稀疏特征值计算，之后针对稀疏特征值再进行反变换和牛顿校验，得到时滞电力系统的精确特征值。

技术领域

本发明涉及电力系统技术领域，特别是涉及基于低阶IGD-IRK的时滞电力系统小干扰稳定性分析方法，IGD-IRK为英文“Infinitesimal Generator Discretization WithImplicit Runge-Kutta”的缩写，中文含义为：无穷小生成元隐式龙格-库塔离散化。

背景技术

当能源与电力的需求持续增加，建立跨区、跨国的大型互联电网是当代电力发展趋势。然而在互联初期，复杂的电气结构和薄弱的输电环节使得电力系统更容易出现区域间低频振荡现象。区间振荡的范围较广，振荡的机群之间的联系较为复杂，是制约系统小干扰稳定性的关键。以本地信号为输入的电力系统稳定器(Power System Stabilizer，PSS)能很好地抑制局部振荡模式，但是难以抑制区间振荡，从而难以保证大规模互联电网的稳定性。

广域测量系统(Wide-Area Measurement System，WAMS)的出现和发展，能够实时远距离采集动态参数，使得采用全局的、远端的信息实现电力系统稳定控制成为可能。广域测量系统主要由三部分组成，相量测量单元(Phasor Measurement Unit，PMU)，通信网络，以及监控系统。PMU同步测量系统各枢纽点的状态量，通过通信系统通道传送到监控系统，实现系统的保护控制与检测。经过分析处理的量测数据，不仅能检测系统低频振荡信息，还可以提升系统的阻尼水平和远距离大容量输电能力。

但是，广域量测数据在广域量测系统由不同的通信介质(如光纤、数字微波、电力线等)组成的通信系统中传输和处理时，存在明显的时延现象。时滞在控制系统中引入了一个滞后的相位，可能会使得系统的控制律失效，恶化系统的性能，甚至导致系统从小扰动稳定变为不稳定。因此分析时滞对于电力系统小干扰稳定性的影响，对于设计能够有效抑制区域间低频振荡的广域阻尼控制器，提升电力系统的安全、稳定运行水平，具有重要的理论意义和应用价值。

在现代电力系统中，时滞电力系统小干扰稳定性分析的关键就是关注系统的机电振荡问题。常用的时滞电力系统稳定性分析方法有函数变换法，时域法，预测补偿法等等。中国发明专利基于Padé近似的时滞电力系统特征值计算与稳定性判别方法，申请号为201210271783.8:[P].利用Padé近似多项式逼近时滞环节，进而计算系统最右侧的关键特征值，并判断系统的时滞稳定性。而特征分析法是电力系统小干扰稳定性分析的基本而有效的方法，已经形成了比较成熟和完善的理论。计算得到时滞电力系统的特征值，就可以运用经典的特征分析的思路和理论框架来分析系统小干扰稳定性，进行广域阻尼控制器的优化设计。

近年来，基于谱离散化的时滞系统部分特征值计算方法，开始用于时滞电力系统的稳定性分析。中国发明专利基于EIGD的大规模时滞电力系统特征值计算方法.201510055743.3[P].提出了一种基于显示IGD(Explicit Infinitesimal GeneratorDiscretization,EIGD)的大规模时滞电力系统特征值计算。利用计算得到时滞电力系统的特征值，进而判定系统的小干扰稳定性。

然而，上述时滞电力系统谱离散化计算方法生成的离散化矩阵维数很大，尽管采用稀疏算法进行计算，使其适用于大规模系统的特征值计算，但是矩阵固有的维数问题没有得到解决。

发明内容