一种提取极坐标潮流方程雅可比矩阵的方法与流程

本发明属于电力潮流计算方法技术领域，具体涉及一种提取极坐标潮流方程雅可比矩阵的方法。

背景技术：

电力系统潮流计算是分析电力系统稳态运行情况的一种基本计算。数学原理上，它是求解多元非线性方程组的问题，应用最为广泛的求解方法是牛顿法。使用牛顿法进行求解的计算量，很大程度上在于形成雅可比矩阵这一步骤，降低这一步骤的复杂度，成为提升算法的计算效率，减少算法对计算机资源占用的重要手段。

技术实现要素：

本发明的目的在于针对现有技术中的上述问题，提供一种提取极坐标潮流方程雅可比矩阵的方法，可降低雅可比矩阵形成过程的复杂度，提升潮流算法的计算速度。

为实现上述发明目的，本发明采用了如下技术方案：

一种提取极坐标潮流方程雅可比矩阵的方法，包括如下步骤：

s1.根据复数的节点功率方程按节点电压的极坐标形式分别对v和δ求偏导，得：

其中，表示节点功率；表示节点功率增量；diag表示取对应矢量的对角矩阵；表示电压向量；表示的共轭复数；表示节点导纳矩阵；表示的共轭复数；v表示电压幅值；δ表示电压相角；

s2.展开步骤s1中得到的式(1)，得到复功率不平衡量表达式：

式中，δv为电压幅值增量，δδ为电压相角增量；

s3.把步骤s2中得到的式(2)的δv项先乘以diag(v)后，再除去diag(v)，得：

s4.对步骤s3中得到的式(3)进行符号化替换，令

将式(3)变为：

s5.从步骤s4得到的式(4)中分别提取矩阵d、矩阵j的实部和虚部，得：

式中，re表示提取对应矩阵的实部；im表示提取对应矩阵的虚部；dr表示d的实部；jr表示j的实部；di表示d的虚部；ji表示j的虚部；

s6.将步骤s5中的式(5)代入步骤s4中的式(4)中，得：

s7.将步骤s6得到的式(6)整理为矩阵形式，其系数矩阵即为雅可比矩阵：

式中，δp表示有功功率增量；δq表示无功功率增量。

相比于现有技术，本发明的优势在于：

本发明所提供的一种提取极坐标潮流方程雅可比矩阵的方法，在现有使用牛顿法求解潮流问题形成雅可比矩阵的过程中，从节点复功率方程出发，直接提取雅可比矩阵，简化了雅可比矩阵的形成过程；同时，此过程不存在三角函数运算，大大降低了计算量及出错概率。本发明的方法极大地降低了雅可比矩阵形成过程的复杂度，提升大电网潮流计算速度，减少算法对计算机cpu、内存的占用。

附图说明

图1是本发明一种提取极坐标潮流方程雅可比矩阵的方法的流程图。

具体实施方式

以下结合实施例及其附图对本发明技术方案作进一步非限制性的详细说明。

如图1所示，一种提取极坐标潮流方程雅可比矩阵的方法，包括如下步骤：

s1.从节点复功率方程出发，根据复数的节点功率方程按节点电压的极坐标形式分别对v和δ求偏导，得：

其中，表示节点复功率；表示节点复功率增量；diag表示取对应矢量的对角矩阵；表示电压向量；表示的共轭复数；表示节点导纳矩阵；表示的共轭复数；v表示电压幅值；δ表示电压相角；为实数向量；

s2.展开步骤s1中得到的式(1)，得到复功率不平衡量表达式：

式中，δv为电压幅值增量，δδ为电压相角增量；

s3.把步骤s2中得到的式(2)的δv项先乘以diag(v)后，再除去diag(v)，得：

s4.对步骤s3中得到的式(3)进行符号化替换，令

将式(3)变为：

s5.从步骤s4得到的式(4)中分别提取矩阵d、矩阵j的实部和虚部，得：

式中，re表示提取对应矩阵的实部；im表示提取对应矩阵的虚部；dr表示d的实部；jr表示j的实部；di表示d的虚部；ji表示j的虚部；j为虚部符号；

s6.将步骤s5中的式(5)代入步骤s4中的式(4)中，得：

s7.将步骤s6得到的式(6)整理为矩阵形式，其系数矩阵即为雅可比矩阵：

式中，δp表示有功功率增量；δq表示无功功率增量。

实施例：

对ieee4潮流分析数据进行电力系统潮流计算，包含如下步骤：

s1：给出电压向量初始值，并给出收敛精度，为：ε＝10^-6；

s2：形成节点导纳矩阵，得：

s3：根据复数的节点功率方程按节点电压的极坐标形式分别对v和δ求偏导，得：

其中，表示节点复功率；表示节点复功率增量；diag表示取对应矢量的对角矩阵；表示电压向量；表示的共轭复数；表示节点导纳矩阵；表示的共轭复数；v表示电压幅值；表示电压相角；为实数向量；

s4：展开式(1)，得到复功率不平衡量表达式：

式中，δv为电压幅值增量，δδ为电压相角增量；

s5：把式(2)中的δv项先乘以diag(v)后，再除去diag(v)，得：

s6：令式(3)改为：

s7：从式(4)提取矩阵d、矩阵j的实部和虚部，得：

式中，re表示提取对应矩阵的实部；im表示提取对应矩阵的虚部；dr表示d的实部；jr表示j的实部；di表示d的虚部；ji表示j的虚部；

s8：将式(5)代入式(4)中，得：

s9：将式(6)整理为矩阵形式，系数矩阵即为雅可比：

式中，p表示有功功率；δp表示有功功率增量；q表示无功功率；δq表示无功功率增量；

按照上述步骤s3-s9，并处理pv节点电压恒定问题，形成雅可比矩阵，在第一次迭代过程中，雅可比矩阵为：

s10：进行迭代，达到收敛精度后输出计算结果，为：

v1＝0.9847e^-0.0087j，v2＝0.9648e^-0.1126j，v3＝1.1000e^0.1175j，v1＝1.0500e^0j。

运用本发明的方法进行上述电力系统潮流计算，其中的雅可比矩阵的形成过程得到了极大的简化，在实际应用中，能够极大提高大电网潮流计算效率，降低对电脑cpu、内存的占用。

综上所述，本发明在现有使用牛顿法求解潮流问题形成雅可比矩阵的过程中，从节点复功率方程出发，根据复数的节点功率方程，按节点电压的极坐标形式分别对v和δ求偏导；δv项先乘以v的对角矩阵后再除去；然后展开得到复功率不平衡量表达式；分别提取替换后矩阵的实部和虚部，并将其代入复功率不平衡量表达式；最后整理为矩阵形式，其系数矩阵即为雅可比矩阵。

以上所揭露的仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或变型，都应涵盖在本发明的保护范围之内。