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

直接计算静态电压稳定临界点的新方法 2 Pei 江伟’，王成山1，余贻鑫1，ZHAN6 (1．天津大学电气与自动化工程学院，天津市南开区300072：2．EPRIU．S．，Palo Alto，CA，USA，94303) 肇要t文中提出一种降阶求解静态电压稳定临界点的新技术．其特点是原直接法的(2n+1)维牛顿迭 代方程组可降阶为(斛．1)维线性方程组求解．与解高维方程的传统直接法相比该方法计算量小、易于采用稀 疏技术实现。适合于大规模电力系统电压稳定临界点的在线求解。另外，还提出一种剩用负荷参数的二阶导 数进行临界点预测的新方法，解决了直接法各状态变量及右特征向量的初值难于确定的问题。在1000母线系 统上的计算表明该文方法具有忘好的收敛性． 关键词：静态电压稳定，鞍结分岔；牛顿法；=阶导数；矩阵降阶；右特征向量 1目l言 近年来大系统电压崩溃事故频频发生，电压稳定‘”1问题又一次成为专家学者们研究的热 点。在电压稳定研究中，静态电压稳定临界点(也称为潮流方程鞍型分岔点，Fold分岔点， 二次转折分岔点，以下简称为SNB点)的计算一直是一项基础性课题。精确的SNB点不仅提 供诸如电压稳定裕度厶。的重要信息，还能为预防控制、紧急控制提供必要的灵敏度信息【卅。 研究、设计一个快速而精确的SNB点计算方法，对电压稳定在线安全评价具有重要意义． 在SNB点处，潮流方程雅克比矩阵奇异，造成潮流方程难于求解。为了克服这一困难， 早期计算SNB点的方法有电纳法嘲，目前较常用的有连续法哪】和直接法[9-”1等。连续法通过 追踪Py曲线来获取近似的SNB点。虽然连续法被认为是比较可靠的一种方法，但是计算速 度较慢的缺点使其较适合离线计算而不适于在线应用。直接法通过求解非线性特征方程组直 接求取SNB点。目前直接法的应用不如连续法广泛，其原因主要是直接法需求解的非线性方 程组比连续潮流方程组的维数增加一倍，所需内存空间及计算量增加许多，也不易采用稀疏 技术求解。另外，SNB点初值的选取较为困难，常导致直接法牛顿迭代难于收敛。 本文在文献B3]的基础上提出一种降阶求解SNB点的新方法。箕特点是通过变量替换， 把直接法的(2n+1)维牛顿迭代方程组转变为四个同系数矩阵的(肘1)维线性方程组求解。 由于所要求解的线性方程组维数降低一半，而且每次迭代只需进行一次LU分解，因此改进后 的直接法的计算效率大大提高。另外本文还提出一种利用负荷参数的二阶导数进行SNB点初 值预估的新方法。最后在Ontario1000母线系统算例上验证了本文算法的高效性。 Hydro 2 SNB点特征方程的求解 2．1 SNB点的特征方程 用f(x，A)=O表示以负荷为参数的潮流方程。则采用右特征向量形式的SNB点处的特 征方程为㈣： ，o，A)=0 (1) ￡(x，A)·V；0 (2) ％一1=0 (3) 雅可比矩阵奇异性条件方程，v∈置。；式(3)是规范化方程，匕表示v的第P个分量。 在SNB点处￡具有与简单零特征值对应的左特征向量∞(行向量)与右特征向量v。并 有： ∞·(丘·V)·y≠0 解非线性方程组(1邮)的牛顿迭代方程为 I厶·，正 I△’I—I正。，I(4) 01 陪纠0‘ 0阱嘲0∽L从J 【 j l J 式中勺6F，表示除第p个元素为1外其他均为0的单位向量(如无特殊注明，本文所指 向量均为列向量)；厶‘pE矗“。·其第‘行、第J列元素定义为∑k=l仁UX生I：O}Xk。唯)口 可以证明．在SNB点处方程(4)左端的系数矩阵非奇异‘…。r’ 直接求解牛顿迭代方程(4)具有一定的困难。方程(4)