一种适合研究目的使用的牛顿法潮流计算方法

摘要

一种适合研究目的使用的牛顿法潮流计算方法，为以牛顿法潮流计算为基础进行进一步研究的科研人员提供一个易于修改和维护的牛顿法潮流计算算法。采用极坐标牛顿法潮流计算的基本原理步骤，在分析潮流计算各个组成部分的基础上提出了一种通过简单逻辑判断来避免不必要运算以提高潮流计算计算速度的方法。本发明的技术方案是在形成雅可比矩阵模块和修正方程高斯消去模块中通过简单判断来避免不必要的运算。采用本发明对一个445节点实际大型电网进行了计算，计算时间为0.484s，计算速度完全能够满足科研需要。

主权项

一种适合研究目的使用的牛顿法潮流计算方法，采用极坐标牛顿法潮流计算的基本原理步骤，其特征在于：在形成雅可比矩阵模块和修正方程高斯消去模块中通过简单判断来避免不必要的运算；其中所述的形成雅可比矩阵模块包括以下步骤：步骤1：设置行号i＝1；步骤2：判断节点i是否为平衡节点，如果是平衡节点转至步骤15；步骤3：设置列号j＝1；步骤4：判断节点j是否为平衡节点，如果是平衡节点转至步骤13；步骤5：判断节点i和节点j之间的导纳实部Gij和虚部Bij是否都为0，如果都为0转至步骤13；步骤6：计算雅可比矩阵元素Ji+i‑1，j+j‑1，如果j≠i，根据式(3)计算，如果j＝i，根据式(11)计算；步骤7：判断节点j是否为PQ节点，如果不是PQ节点转至步骤9；步骤8：计算雅可比矩阵元素Ji+i‑1，j+j，如果j≠i，根据式(4)计算，如果j＝i，根据式(12)计算；步骤9：判断节点i是否为PQ节点，如果不是PQ节点转至步骤13；步骤10：计算雅可比矩阵元素Ji+i，j+j‑1，如果j≠i，根据式(5)计算，如果j＝i，根据式(13)计算；步骤11：判断节点j是否为PQ节点，如果不是PQ节点转至步骤13；步骤12：计算雅可比矩阵元素Ji+i，j+j，如果j≠i，根据式(6)计算，如果j＝i，根据式(14)计算；步骤13：令j＝j+1；步骤14：判断j是否大于n，如果j不大于n，则返回到步骤4；步骤15：令i＝i+1；步骤16：判断i是否大于n，如果i不大于n，则返回到步骤2；否则结束；所述的修正方程高斯消去模块的步骤是：步骤1：设置当前行号i＝1；步骤2：判断矩阵元素aii是否为0，如果aii为0，则转至步骤17；步骤3：设置k＝1；步骤4：判断k是否小于i，如果k不小于i，则转至步骤12；步骤5：判断矩阵元素aik是否为0，如果aik为0，则转至步骤11；步骤6：设置当前列号j＝k+1；步骤7：判断j是否大于n，如果j大于n，则转至步骤10；步骤8：根据式(18)对系数矩阵进行消去运算；步骤9：令j＝j+1，返回到步骤7；步骤10：根据式(20)对右端常数项进行消去运算；步骤11：令k＝k+1，返回到步骤4；步骤12：设置当前列号j＝i+1；步骤13：判断j是否大于n，如果j大于n，则转至步骤16；步骤14：根据式(19)对系数矩阵进行规格化运算；步骤15：令j＝j+1，返回到步骤13；步骤16：根据式(21)对右端常数项进行规格化运算；步骤17：令i＝i+1；步骤18：判断i是否大于n，如果i不大于n，则返回到步骤2；否则结束。