直接计算静态电压稳定临界点的新方法研究.pdfVIP

直接计算静态电压稳定临界点的新方法呋 2 Pei 江伟1，王成山1，余贻鑫1，ZHANG (1．天津大学电气与自动化工程学院，天津市南开区300072： 2．EPRI U．S．，PaloAlto，CA，USA，94303) 摘要：文中提出一种降阶求解静态电压稳定l临界点的新技术．其特点是原直接法的(2n+1)维牛顿迭代方程 组可降阶为(肘1)维线性方程组求解．与解高维方程的传统直接法相比该方法计算量小、易于采用 稀疏技术实现，适合于大规模电力系统电压稳定I临界点的在线求解．另外，还提出一种利用负荷参 数的二阶导数进行临界点预测的新方法，解决了直接法各状态变量及右特征向量的初值难于确定的 问题。在1000母线系统上的计算表明该文方法具有良好的收敛性。 基金项目：国家自然科学基金重大项目；美国电力科学研究院合作科研基金项目(EP-P11543／C5729)． NationalNaturalScienceFoundationof ProjectSupportedby China；EPRI 关键词t静态电压稳定；鞍结分岔：牛顿法；二阶导数；矩阵降阶；右特征向量 1 引言 近年来大系统电压崩溃事故频频发生，电压稳定【l之】问题又一次成为专家学者们研究的热 点。在电压稳定研究中，静态电压稳定临界点(也称为潮流方程鞍型分岔点，Fold分岔点，二 次转折分岔点，以下简称为SNB点)的计算一直是一项基础性课题。精确的SNB点不仅提供 诸如电压稳定裕度AmⅨ的重要信息，还能为预防控制、紧急控制提供必要的灵敏度信息[3叫。研 究、设计一个快速而精确的SNB点计算方法，对电压稳定在线安全评价具有重要意义。 在SNB点处，潮流方程雅克比矩阵奇异，造成潮流方程难于求解。为了克服这一困难，早 期计算SNB点的方法有电纳法‘51，目前较常用的有连续法[硒1和直接法【9-121等。连续法通过追踪 PV曲线来获取近似的SNB点。虽然连续法被认为是比较可靠的一种方法，但是计算速度较慢 的缺点使其较适合离线计算而不适于在线应用。直接法通过求解非线性特征方程组直接求取 SNB点。目前直接法的应用不如连续法广泛，其原因主要是直接法需求解的非线性方程组比连 续潮流方程组的维数增加一倍，所需内存空间及计算量增加许多，也不易采用稀疏技术求解。 另外，SNB点初值的选取较为困难，常导致直接法牛顿迭代难于收敛。 本文在文献【13】的基础上提出一种降阶求解SNB点的新方法。其特点是通过变量替换，把 直接法的(2n+1)维牛顿迭代方程组转变为四个同系数矩阵的(斛1)维线性方程组求解。由 于所要求解的线性方程组维数降低一半，而且每次迭代只需进行一次LU分解，因此改进后的 直接法的计算效率大大提高。另外本文还提出一种利用负荷参数的二阶导数进行SNB点初值预 估的新方法。最后在Ontario1000母线系统算例上验证了本文算法的高效性。 Hydro 2 SNB点特征方程的求解 2．1 SNB点的特征方程 用f(x，A)=0表示以负荷为参数的潮流方程，则采用右特征向量形式的SNB点处的特 征方程为【131： f(x，旯)=0 (1) L(x，∞·v=O (2)’ ％一l=0 (3) 其中，式(1)是潮流方程，，：X×震呻y，X，l，∈R”，A∈R代表有功负荷参数；式(2)是潮流 雅可比矩阵奇异性条件方程，，∈R”；式(3)是规范化方程，％表示y的第p个分量。 在SNB点处正具有与简单零特征值对应的左特征向量60(行向量)与右特征向量y。并 有： 60·厶≠0 ∞·(厶·，)·p≠0 解非线性方程组(