牛顿拉夫逊潮流计算中检测雅可比矩阵奇异性和网络孤岛的新方法.docVIP

牛顿-拉夫逊潮流计算中检测雅可比矩阵奇异性和网络孤岛的新方法 高新技术—ZUJ.|—I~U.UZChinaNewTechnologiesandProducts国团圜圈囫盈●■■囡■■■-■_■ 牛顿一拉夫逊潮流计算中检测雅可比矩阵 奇异性和网络孤岛的新方法 张勇明王涛周敏 (重庆市电力公司綦南供电局,重庆401420) 摘要:本文提出一种在牛顿一拉夫逊潮流计算中检测雅可比矩阵奇异性和网络孤岛的新方法.即:利用分解后的H矩阵,进4~5,1 (行)缌陛相关性的检测,判断奇异性,回代计算列相关系数,根据对相关系数值的检测,判断那些母线在同一个孤岛中. 关键词:网络孤岛;矩阵奇异;计算 中图分类号:TP30文献标识码:A 1.问题的提出 电力系统静态安全分析,可靠性评估,都要 进行潮流计算,如果采用牛顿一拉夫逊方法,经 常会遇到雅可比矩阵奇异的现象.其实,在常规 牛顿一拉夫:’赴{十算潮流计算时,有时,也会遇到 这种情况,由于这种现象对所研究的问题结果 没有产生特别的影响,因此,人们往往不太注 意,一笔代过.其实,在某些情况下,这种现象对 所研究的问题有意想不到的影响,对此,很有必 要进行分析.我们知道,引起雅可比矩阵出现奇 异的原因很多,网络孤岛(无平衡节点的网络) 是其中的—个原因.那么,怎样检测雅可比矩阵 奇异性,怎样识别网络孤岛,是本文所要解决问 题. 2目前研究的现状 雅可比矩阵奇异现象,在牛顿一拉夫逊计 算潮流中经常会遇到,这是由多种原因造成的. 网络孤岛是雅可比矩阵出现奇异性的—个重要 原因.许多文献由此提出了许多识别网络孤岛 拓朴论基础上的方法.例如:基于Boolean乘法 关联矩阵的判别方法,一种简单的树类型寻找 的方法,深度优先搜索方法等. 拓朴论基础上搜索网络孤岛至少有以下三 个缺陷.第一,网络孤岛并不是雅可比矩阵出现 奇异的唯一原因;第二,即便是一个分裂的网 络,只要孤岛中存在平衡母线,潮流问题也是可 以解决的;第三,拓朴论上的方法,不论网络是 否连通,在潮流计算前都被执行,比较费时. 参考文献【2】中方法,克服了上述缺陷只需 对分解后的稚可比矩阵进行判断,如果奇异,仅 需对分解后的雅可比矩阵进行回代运算,计算 列相关系数,判断网络孤岛.但在该文中,用NM 方法比用FD方法速度慢得多. 基本原理 A是nxn阶矩阵,无论A是否奇异,都可 经过如下规格化,消去过程,分解成LDU的形 式: aik,=aik/aii(Vkgt;(1a) ahk=ahk--a~aik/a(Vhgt;(1b) 由式(1a),(1b问’明显看出: 当aii=0时,规格化,消去过程不能进行下 去.对此,作如下处理: (i)如果aik=0(Vkgt;,就绕过此行继续 进行分解过程. (ii)如果Vkgt;i)不全是零,就推迟此行 的规格化,消去过程. A矩阵经过(1a),(1喊的运算后,可表示为: A=LDU=LU’=L 其中:L是nxn阶下三角矩阵,对角元素全 为1; U是n×n阶上三角矩阵,对角元素全为1; D是nxn阶对角矩阵. 由线性代数知识,可引出如下定理: 定理一矩阵奇异当且仅当矩阵分解后, (2a)或(2嘁中一个成立: aft=0且ai0(Vkgt;(2a) ai.=0且=0(Vkgt;(2b) 即:u中有一零行或中有一零列,那么 矩阵奇异. 证明:A=LDU=LU=LUf31 必要性:lAI=lLDUI=lL1lDlIU l=lLIIul=ILl1UI 由于U中有一零行或L中有一零列,则l ul:O或lLI=0 所以,IAI=0,矩阵奇异. 充分陛:矩阵A奇异,那么IAl=0 由(3)式可得:I【=0由于D是对角矩 阵,因此,至少有一对角元素为0. 因为U=DU,LVLD 所以(2a)或(2b)工弋有一个成立,即U中有一 零行或中有一零列. u中行为零,这是行相关隋况;L丰中列为 零,这是列相关隋况. A矩阵奇异,那么A矩阵行向量,列向量线 性相关,即: ∑ici=0 (4a) (4b) 其中,ri和c分别是A矩阵的行向量和列向 量,6i和.是不全为零的系数. 另一个有用的结论在下面给出: 搬台一当行相关或列相关发生时,即(4a) 或(4b1式成立时,对应最后个不为零的8i和 . 是u的零行或L的零列. 证明:设rk和C是最后一个8k≠0(≠ 0)对应的行或列.显然,由前k一1行(列)向量确 定的向量空间的维数不因增加第k行(列)而增 加.因为行(列)空间的维数等于u1或中非零 行(列)的个数,这可推出uf中第k行(列) 必是—个零行(列o 网络孤岛的识别方法 4.1原理 系统中若存在孤立节点,在形成导纳矩阵 就可检测出,因为孤立节点在导纳矩阵中所在 的行和列全是零.因此,在本文中,不考虑孤立 节点的情况.由于平衡节点不参与雅可