基于PIGD-IRK的大规模时滞电力系统特征值高效计算方法.pdfVIP

第 47 卷 第 12 期 2023 年 6 月 25 日 Vol . 47 No . 12 Jun . 25 ，2023 DOI ： 10. 7500/AEPS20220624006 基于 PIGD -IRK 的大规模时滞电力系统特征值高效计算方法 1 1 1 2 贾小凡 ，叶 华 ，刘玉田 ，刘牧阳 （1. 电网智能化调度与控制教育部重点实验室（山东大学），山东省济南市250061 ； 2 . 新疆大学电气工程学院，新疆维吾尔自治区乌鲁木齐市830046 ） 摘要 ：电力系统控制回路中固有的时滞会影响控制器的性能，甚至危害系统的稳定运行。为此，提 出了一种基于无穷小生成元隐式龙格-库塔部分离散化（PIGD -IRK ）的特征值计算方法，用以分析 时滞对大规模电力系统小干扰稳定性的影响。首先,分析了制约无穷小生成元隐式龙格-库塔离散 化（IGD-IRK ）方法计算效率的瓶颈；然后,应用部分谱离散化思想大幅度降低其构造的无穷小生成 元离散化矩阵的维数；接着，利用离散化矩阵的结构特点将其从阶梯型分块矩阵相似变换为常数矩 阵的逆和分块上三角矩阵的乘积，从而可在特征值计算过程中充分利用离散化矩阵和系统增广状 态矩阵的稀疏特性。最后，通过四机两区域系统和2 个实际电网的计算结果验证了所提方法的准 确性和高效性。 关键词 ：时滞；小干扰稳定；谱离散化；无穷小生成元；隐式龙格-库塔法 0 引言 求解得到时滞电力系统在复平面上右侧的所有关键 特征值，且计算结果的准确度不受时滞大小的影响。 在电力系统的控制回路中，反馈信号和控制信 在IGD 类方法中，文献［11 ，14 ］所提出的无穷 号在采集、路由、传输、处理和执行过程中普遍存在 一定的时滞。例如，基于同步量测的广域阻尼控制 小生成元隐式龙格-库塔离散化（IGD with implicit ［1 ］ Runge-Kutta ，IGD-IRK ）方法计算得到的系统关键 系统的通信时滞在几十到数百毫秒范围内变化 ， 频率控制系统的通信时滞和采样诱导时滞高达数 特征值具有高精度，但计算效率较低。由于IGD- ［2 ］ IRK 生成的离散化矩阵为阶梯型分块矩阵且维数较 秒 ，基于模块化多电平换流器的高压直流 （modular multilevel converter based high -voltage 高，导致在特征值计算过程中无法充分利用离散化 direct current ，MMC-HVDC ）输电系统的链路延时 矩阵和系统增广状态矩阵的稀疏特性求解其逆矩阵 一般为400~ 600 μs［3 ］。在实际的工程系统中，时滞 与向量的乘积 （matrix inversion-vector product ， 导致理论上设计良好的控制器的实际性能不佳甚至 MIVP ）运算。IGD-IRK 方法应用于大规模时滞电 ［4-7 ］ 失效，成为电力系统稳定运行的潜在威胁 。 力系统的高效性有待进一步改善。 与传统电力系统不同，时滞电力系统的特征方 鉴于此，本文提出了一种基于无穷小生成元隐 程为含有指数项的超越方程，给系统的小干扰稳定 式龙格- 库塔部分离散化（partial IGD -IRK ，PIGD - 性分析增加很大的难度。为了避免直接求解超越方 IRK ）的大规模时滞电力系统特征值高效、准确计算 程 ，文献 ［8 ］和文献