[发明专利]一种直角坐标牛顿法潮流计算方法

摘要

本发明公开了一种直角坐标牛顿法潮流计算方法，该方法从直角坐标牛顿法潮流计算的基本原理出发，在分析其基本修正方程的特点基础上提出了一种直角坐标牛顿法潮流计算的雅可比矩阵构成方法来改善潮流计算收敛性。该方法包括以下步骤：原始数据输入和电压初始化；形成节点导纳矩阵；形成雅可比矩阵J；解方程及修正电压实部e、虚部f；节点及支路数据输出。由于该方法对雅可比矩阵J的部分对角元素计算公式进行了修正，解决了直角坐标牛顿法潮流计算在分析含有小阻抗支路系统时的收敛性问题，本发明迭代次数、内存占有量和运算速度都与常规直角坐标牛顿法潮流计算处理不含小阻抗系统时相当，计算量不但没有增加，反而略有减少。

主权项

1、一种直角坐标牛顿法潮流计算方法，包括以下步骤：A、原始数据输入和电压初始化；B、形成节点导纳矩阵；C、形成雅可比矩阵J；D、解方程及修正电压实部e、虚部f；E、节点及支路数据输出；其特征在于：所述的雅可比矩阵J的部分对角元素计算公式如下：$\frac{\∂\mathrm{\Δ}{P}_i}{\∂e_i}=-G_{\mathrm{ii}}{e}_i-B_{\mathrm{ii}}{f}_i$$\frac{\∂\mathrm{\Δ}{P}_i}{\∂f_i}=B_{\mathrm{ii}}{e}_i-G_{\mathrm{ii}}{f}_i$$\frac{\∂\mathrm{\Δ}{Q}_i}{\∂e_i}=B_{\mathrm{ii}}{e}_i-G_{\mathrm{ii}}{f}_i$$\frac{\∂\mathrm{\Δ}{Q}_i}{\∂f_i}=G_{\mathrm{ii}}{e}_i+B_{\mathrm{ii}}{f}_i.$