一种基于实时动态等值的电力系统快速暂态稳定仿真方法
【技术领域】
[0001] 本发明属于电力系统自动化，特别涉及了一种适用于在线计算的电力系统暂态稳 定仿真方法。
【背景技术】
[0002] 电网的不断扩大和电力市场的出现使得电力系统的运行环境更加复杂，对电网安 全稳定运行的要求也越来越高。尽管我国在稳定控制系统领域已经有相当多的研究成果， 但是还存在诸多不足。稳定控制中"离线计算、实时匹配"的方式计算量大，对运行方式和 网络结构变化适应能力差，很容易出现失配情况；"在线预决策、实时匹配"测量误差和传送 丢失都有可能造成所确定的系统运行方式与实际运行方式失配，从而造成预决策失误；"实 时决策、实时控制"则是最理想的稳定控制手段，它要求根据检测到的故障信息，按当时接 线方式和潮流方式，超实时计算并实施控制，完全避免运行工况和故障的失配问题。但相应 的技术难度最大，要求有能够对受扰系统做出快速准确预测和控制的良好算法，如何开发 出超实时的算法解决这类系统的暂态稳定性问题值得研究。
[0003] 现代电力系统在线实时的安全稳定控制越来越重要，而常规的稳定分析方法花费 时间长，无法应用于在线实时应用场合。为解决电力系统计算中的简化问题提出了动态等 值的方法，主要使用在以下三种场合：大规模电力系统的离线暂态稳定分析、离线动态稳定 分析及在线动态安全分析。与这三种状态相适应，动态等值方法也可分为三大类：同调等值 法、模态等值法及估计等值法。前两种均用于离线分析，而在线分析中常需要对系统作在线 的动态等值，而由于系统的结构、运行工况多变，一般不可能离线先做等值然后在线调用， 因此只有估计等值法适用于在线分析。
[0004] 传统的估计等值法同样存在一些不足：在选取外部系统的等值模型时具有经验 性，每次迭代修正参数时步长要确保参数快速收敛，并且很多情况下需要考虑噪声影响并 作滤波及相关分析，数学处理较为复杂。
【发明内容】
[0005] 本发明目的在于针对现有技术的不足，提出一种基于实时动态等值的电力系统快 速暂态稳定仿真方法，该方法可克服电力系统暂态稳定分析计算中，现有的数值积分方法 计算量大，计算速度不能满足电力系统在线计算要求的缺点。本发明的方法是一种以在线 计算为目标的超实时暂态稳定仿真算法，采用一种等值系统节点导纳矩阵快速辨识方法， 求取等值系统网络参数，显著简化了电力系统网络，减少了大规模电力系统暂态稳定仿真 的计算量。
[0006] 本发明是通过以下技术方案实现的：一种基于实时动态等值的电力系统快速暂态 稳定仿真方法，包括以下步骤：
[0007] 步骤1 :对原始系统进行暂态稳定数值积分计算，只仿真到等值时刻tz;
[0008] 步骤2 :若未发现无界角度间隙则确认是稳定的情况，若发现了无界角度间隙则 执行以下的步骤；
[0009] 步骤3 :根据等值时刻前的发电机功角曲线对系统进行分群；
[0010] 步骤4 :根据步骤3的分群结果，对每一群的发电机进行模型聚合获得该群的等值 发电机，并应用等值系统节点导纳矩阵快速辨识方法获得等值系统的节点导纳矩阵；
[0011] 步骤5 :在等值系统上完成暂态稳定数值积分计算。
[0012] 上述技术方案中，步骤2中的通过判断是否出现无界角度间隙以排除极安全情 况，所述的无界角度间隙是：在数值积分的每个积分步结束处将各发电机功角从大到小排 列一次，每两个相邻功角之间构成一个角度间隙，如果存在某角度间隙对于任意给定的正 值β，总可以找到一个时间te使得t>t e时该角度间隙必大于β，则称该角度间隙为无界 角度间隙，在电力系统实际应用中，往往按经验取有限值为门限值β，而将大于β值的角 度间隙成为无界角度间隙。
[0013] 在原始系统中进行数值积分至等值时刻，等值时刻大于故障切除时刻。从开始时 刻至等值时刻观察原始系统各发电机的功角轨迹，如果发现了无界角度间隙则进一步进行 暂态稳定分析，反之，则确认为是极安全（极为稳定）的情况。为了尽可能减少在原始系统 上的数值积分时间，本发明公开的仿真算法可选取较小的等值时刻，所以在选取角度间隙 门限值时也取相对较小的值。
[0014] 步骤3中的所述的对系统分群的方法如下：
[0015] 若在分析开始至等值时刻的时段内发现了无界角度间隙，则在等值时刻对发电机 进行分群。将某两台发电机划为同群的标准为在等值时刻前它们的功角摇摆曲线具有相同 的运动趋势。本发明的方法中计算两台发电机等值时刻前所有积分步的功角差值形成一组 数据，求取该组数据的方差值作为这两台发电机功角摇摆曲线同趋势的评价指标。设定方 差判定阈值ε，满足如下条件的两台发电机判断为同群：
[0019] 式中νζ是至等值时刻总的数值积分步数；δ ki、δ k]分别为第k时步第i台机与第 j台机的功角值；^:为t e [0, tz]时间段内第i台机与第j台机功角差值平均值；#为 t e [0, τζ]时间段内第i台机与第j台机功角差值的方差；ε为给定的方差判定阈值。
[0020] 步骤4中的发电机模型聚合的方法如下：
[0021] 应用静态eeac的模型聚合方法进行模型聚合，经典模型下第i群的等值发电机全 部参数计算公式见如下公式：
[0031] 其中111是系统的分群个数，1 = 1，2*"111;1111是第1群的发电机个数；1';1、01、乂 /,1、 e' p ρηι、p& ixl、iyl、s j别表示模型聚合后第i群等值发电机的惯性时间常数、阻尼 系数、直轴暂态电抗、暂态电势、机械功率、电磁功率、机端电流实部、机端电流虚部及功角； γ4、difc、；4fc、私、表示第i群第k台发电机的惯性时间常数、阻尼系数、直轴 暂态电抗、暂态电势及机械功率；、li、·分别表示等值时刻第i群第k台 发电机的电磁功率、机端电流实部、机端电流虚部及功角。
[0032] 步骤4中所述的等值系统节点导纳矩阵快速辨识方法如下：
[0033] 若等值时刻将系统分成了 m群，则在原始系统上从等值时刻继续做m-1步数值积 分计算，积分步长为δ t，通过数值积分计算可以获得在τζ，τζ δ t，…，tz (m_l) δ t总共m 个时步内，每一台原始发电机在每一时步的机端电流实部ix、机端电流虚部、以及功角δ。
[0034] 发电机经典模型下，应用如下公式计算获得第i群等值发电机第s时步机端电流 实部ixls及机端电流虚部i yls
[0037] 其中，lp、/yi>·分别表示第i群第k台发电机第s时步的机端电流实部、虚部。
[0038] 通过如下公式计算第i群等值发电机第s时步的机端电压实部vxls及机端电压虚 部 vyls
[0040] 再通过如下公式计算第i群等值发电机第s时步的虚拟注入电流实部r xls及虚 拟注入电流虚部"yls
[0042] 上述两个公式中有：
[0045] 其中ral、x' dl、δ pe'汾别表示该计算时步第i群等值发电机绕组阻抗、交轴暂 态电抗、功角以及暂态电势。
[0046] 经过上述步骤获得了 m个时步的所有等值发电机机端电压与虚拟注入电流，通过 如下公式计算得到等值系统的节点导纳矩阵：
[0048] 其中//s、泛&分别为第i群等值发电机第s时步的机端虚拟注入电流及机端电压 的复向量，即：
[0050] 其中，j为虚数单位。
[0051] 本发明的有益效果：本发明针对暂态稳定分析在线实时应用的目标，在时域仿真 法的基础上利用动态等值的思想，不划分研究系统与外部系统，避免了外部系统划分的经 验性；在原始系统上进行短时间的数值积分计算，根据功角摇摆曲线将发电机分成若干群， 分群时保留功角曲线的多群特性；应用静态eeac的模型聚合方法对每一群发电机组进行 等值，节省了估计等值法中复杂的数学处理及噪声滤波等步骤；应用提出的等值系统节点 导纳矩阵快速辨识方法获得等值系统网络参数，最后在等值系统上完成数值积分计算。由 于通常分群数不会太大，辨识等值系统导纳矩阵的耗时不长，使本发明提出的方法不仅具 备常规暂态稳定数值积分法精度高、模型适应性强、收敛性好等优点，而且大幅度简化了计 算系统，减少了计算量。由于动态等值是针对具体系统的给定运行方式及给定扰动下完成 的，适用于在线的快速暂态稳定分析。
【具体实施方式】
[0052] 以下结合具体算例对本发明作进一步说明。
[0053] 本发明提出的一种基于实时动态等值的快速暂态稳定仿真方法，包括以下步骤
[0054] 步骤1 :对原始系统进行暂态稳定数值积分计算，只仿真到等值时刻tz;
[0055] 步骤2 :若未发现无界角度间隙则确认是稳定的情况，若发现了无界角度间隙则 执行以下的步骤；
[0056] 步骤3 :根据等值时刻前的发电机功角曲线对系统进行分群。
[0057] 步骤4 :根据步骤3的分群结果，对每一群的发电机进行模型聚合获得该群的等值 发电机，并应用等值系统节点