牛顿法潮流计算的高效综合稀疏技术 integrated sparse technique of newton power flow calculation.pdfVIP

牛顿法潮流!计算的高效综合稀疏技术 颜伟1，黄正波1，余娟·，张海兵z，项 (1．重庆大学输配电装备及系统安全与新技术国家重点实验室，重庆400030； 2．重庆市电力公司电网检修分公司，重庆400010) 摘要：提出了一种快速实现潮流计算的牛顿法综合潮流稀疏技术。雅可比矩阵的节点分块结构能有效提 高矩阵的形成、修正与线性方程组求解效率。基于此创建了一种由十字链表层和二叉链表层构成的二层链 表结构．十字链表层存储雅可比矩阵，二叉链表层存储节点导纳矩阵，两者之间的对应元素通过指针直接 关联．在雅可比矩阵形成与修正过程中，通过两层链表之间的关联结构可直接从二叉链表层中提取导纳信 息形成或修正十字链表层中的雅可比矩阵．避免消元操作引入的注入元对原始雅可比矩阵结构的破坏所带 来的影响。十字链表层可直接应用于分块雅可比线性方程组求解操作，同时，通过保留链表结构等措施进 一步提高线性方程组求解速度。通过IEEE57到波兰2746节点等5个网络的潮流计算表明：所提出的潮流 综合稀疏技术相对于流行的稀疏技术。效率优势明显。 关键词：牛顿法潮流计算；稀疏技术；雅可比矩阵修正；分块矩阵；线性方程组；动态链表 744 中图分类号：TM 文献标志码：A 文章编号：1004．9649(2010)07．0019-05 接定位待修正元素地址，降低了修正速度。 0引言 为此．本文分析认为节点分块雅可矩阵结构可 有效提高牛顿潮流计算效率。基于此．创建了一种二 稀疏技术(1_2】对提高潮流计算效率具有重要意层链表结构．可以快速定位待修正的雅可比矩阵单 义．并取得了大量研究成果[3--10)。相关的研究成果一元。完成矩阵的修正；同时进行快速的分块矩阵的高 般涉及：常用的节点编号技术有半动态编号法【引、近 斯消元。上述以二层链表为核心的技术统称为综合 似最小度法[”、压缩带宽法…等：常用的存储技术有 潮流稀疏技术。 十字链表法[，驯、单链表法【8J、三元素法…；常用的消 元技术有高斯消元法和LU分解[10】等．而潮流稀疏1 节点分块雅可比矩阵结构在牛顿潮流计 技术便是上述3类操作的组合策略。另外。通过预留 算中的优势 注入元空间、保留链表对角元地址、引入因子表以保 留矩阵高斯消元路径呻]等技术可进一步改进上述3 本文的分析均基于直角坐标形式的雅可比矩 类操作。 阵。但是。文献[10]指出，将极坐标形式的雅可比矩 雅可比矩阵可看作以节点为基础的由2x2阶 阵留出(n—m)行空行和(，l—rn)列空列后，再以节点 子阵[鸶玺]、[≈≈】构成的分块雅可比矩阵c节为单元进行分块．则极坐标形式的雅可比矩阵与直 角坐标形式的节点分块雅可比矩阵具有相同结构。 点分块雅可比矩阵)．但现有研究并没有充分考虑利 因此．本文对于直角坐标形式的雅可比矩阵的推论 用这些特点对算法进行改进。针对雅可比矩阵的形 与应用同样适用于极坐标形式的雅可比矩阵潮流 成和修正．现有研究一般将雅可比矩阵和节点导纳 计算。 矩阵分别存储在2个十字链表之中．雅可比矩阵的 1．1 在雅可比矩阵形成和修正中的效率优势 形成和修正需要在2个链表之间进行数据查找操 节点数为Ⅳ的网络．节点分块雅可比矩阵的2x 作。在第一次迭代后，由于雅可比链表中填入了高斯 消元操作引入的注入元单元．2类链表便具有不同 2阶单元[≥乞】，[笔鸶】的修正公式如下c以％ 的结构。导致后续的雅可比矩阵修正过程中．无法直 为例，％、五和毛与前者具有相似的结构)： 收稿日期：2009一10-11；修回日期：2010-05-04 作者简介：颜伟(1968一)，男，四川乐山人。教授，博士生导师，从事电力系统