CN109617080A - 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 - Google Patents
基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 Download PDFInfo
- Publication number
- CN109617080A CN109617080A CN201811626642.7A CN201811626642A CN109617080A CN 109617080 A CN109617080 A CN 109617080A CN 201811626642 A CN201811626642 A CN 201811626642A CN 109617080 A CN109617080 A CN 109617080A
- Authority
- CN
- China
- Prior art keywords
- node
- formula
- voltage
- power
- jacobian matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
Abstract
本发明公开了一种基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法,通过首次迭代采用与以后各次迭代过程不同的雅可比矩阵计算方法,解决了直角坐标牛顿法潮流计算在求解含有小阻抗支路电力系统时的收敛性问题。首次迭代时,PQ节点采用节点注入功率给定值计算出的电流相量实部和虚部计算雅可比元素,而PV节点采用节点注入有功功率给定值和负荷无功给定值计算出的电流相量实部和虚部计算雅可比元素,不使用电源无功给定值,避免了PV节点电源无功给定值不确定给潮流计算收敛性的不利影响。采用常规直角坐标牛顿法潮流计算不收敛时,本发明能可靠收敛,且比现有专利技术迭代次数少。本发明也能对正常电力系统进行潮流计算,无不良影响。
Description
技术领域
本发明涉及一种电力系统的直角坐标牛顿法潮流计算方法,特别适合含小阻抗支路系统的直角坐标牛顿法潮流计算。
背景技术
电力系统潮流计算是研究电力系统稳态运行的一项基本计算,它根据电力系统给定的运行条件和网络结构确定整个电力系统的运行状态。潮流计算也是电力系统其他分析的基础,如安全分析、暂态稳定分析等都要用到潮流计算。由于具有收敛可靠、计算速度较快及内存需求适中的优点,牛顿法成为当前潮流计算的主流算法。牛顿法分为极坐标和直角坐标两种形式,其中直角坐标牛顿法潮流计算不需要三角函数计算,计算量相对小一些。
在直角坐标牛顿法潮流计算中,节点i的电压采用直角坐标表示为:
对正常电力网络,直角坐标牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,直角坐标牛顿法潮流计算就可能发散。电力系统小阻抗支路可分为小阻抗线路和小阻抗变压器支路,在数学模型上线路可以看作变比为1:1的变压器,因此下面分析时仅以小阻抗变压器支路为例分析。小阻抗变压器模型见图1,变压器的非标准变比k位于节点i侧,阻抗位于标准变比侧。变压器阻抗zij=rij+jxij很小,导纳为
式中,yij、gjj、bjj分别为节点i和节点j之间小阻抗支路的导纳、电导和电纳;rij、xjj分别为节点i和节点j之间小阻抗支路的电阻和电抗。
由于小阻抗支路lij的阻抗很小,支路的电压降也很小,因此变压器两端节点的电压应满足:
如图2所示,现有直角坐标牛顿法潮流计算方法,主要包括以下步骤:
A、输入原始数据和初始化电压
根据电力系统节点的特点,潮流计算把电力系统节点分成3类:节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点;节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点;节点电压幅值和电压相角已知,节点有功功率和无功功率未知的节点称为平衡节点。
电压初始化采用平启动,即PV节点和平衡节点的电压实部取给定值,PQ节点的电压实部取1.0;所有电压的虚部都取0.0。这里单位采用标幺值。
B、形成节点导纳矩阵
设节点i和节点j原来的自电导与自电纳分别为Gi0、Bi0、Gj0、Bj0,在它们之间增加一条小阻抗支路后的自导纳和互导纳分别为:
式中,Yii、Yjj分别为节点i和节点j的自导纳;Yij为节点i和节点j之间的互导纳;rij、xjj分别为节点i和节点j之间小阻抗支路的电阻和电抗;k为节点i和节点j之间小阻抗支路的变比(如果是输电线支路,变比为1);
C、计算功率及电压偏差
PQ节点的功率偏差计算公式为:
式中,n为电力系统的节点数,ei、fi分别为节点i的电压相量的实部和虚部,Gim、Bim分别为节点导纳矩阵元素Yim的实部和虚部,Pis、Qis分别为节点i给定的注入有功功率和无功功率,为
式中,PiG为节点i的电源有功功率,PiL为节点i的负荷有功功率,QiG为节点i的电源无功功率,QiL为节点i的负荷无功功率。
PV节点的有功功率及电压偏差计算公式为:
式中,Vis为节点i给定的电压幅值。
平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差。
平衡节点的电源有功功率和无功功率以及PV节点的电源无功功率采用式(9)计算。
计算各节点功率或电压偏差的最大值,称为最大不平衡量,如果最大不平衡量的绝对值小于给定的收敛精度,转步骤F,否则执行步骤D。
D、形成雅可比矩阵J
雅可比矩阵J的元素(i≠j时)计算公式如下:
PQ节点按式(10)-(13)计算雅可比矩阵元素;PV节点按式(10)、(11)、(14)、(15)计算雅可比矩阵元素;平衡节点不计算雅可比矩阵元素。
雅可比矩阵J的元素(i=j时)计算公式如下:
式中,ai、ci分别为节点i的计算注入电流相量的实部和虚部,为
PQ节点按式(16)-(19)计算雅可比矩阵元素;PV节点按式(16)、(17)、(20)、(21)计算雅可比矩阵元素;平衡节点不计算雅可比矩阵元素。
E、解修正方程及修正电压实部e、虚部f
修正方程为:
式中,ΔW为修正方程不平衡量列向量,ΔP为有功功率不平衡量列向量,ΔQ为无功功率不平衡量列向量,ΔV2为电压平方不平衡量列向量,J为雅可比矩阵,Δe为电压相量实部列向量,Δf为电压相量虚部列向量。
电压修正公式为:
式中,上标(t)表示第t次迭代。
F、输出节点及支路数据。
对正常电力网络,牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,牛顿法潮流计算就可能发散。而电力系统中小阻抗支路普遍存在,收敛性是电力系统潮流计算这类非线性问题的最重要指标,计算不收敛就无法得到方程的解。因此改善直角坐标牛顿法潮流计算针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。
中国专利ZL201410299531.5披露了一种修改常规直角坐标牛顿法潮流计算雅可比矩阵的方法,该方法依据是:潮流计算收敛时,式(6)中ΔPi、ΔQi都趋近于0,则ai和ci等于由节点注入功率给定值Pis和Qis计算出的ais和cis,即
中国专利ZL201410299531.5方法计算雅可比元素时采用由给定值Pis和Qis计算的ai和ci值,改善了潮流计算的收敛性,有效解决了含有电阻为0的小阻抗支路电力系统潮流计算的发散问题。但当小阻抗支路的电阻不为0时,该方法迭代次数增加,收敛性变差,甚至不收敛。
中国专利ZL201410315785.1提出了一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法,该方法首次迭代和后续各次迭代采用不同的雅可比矩阵计算方法,首次迭代计算雅可比元素时采用由给定值Pis和Qis计算的ai和ci值,后续各次迭代计算雅可比元素时仍采用传统方法,有效解决了含电阻不为0的小阻抗支路电力系统潮流计算的发散问题,但该方法在电阻不为0的小阻抗支路较多时,迭代次数增加,收敛性变差。
发明内容
为解决现有技术存在的上述问题,本发明要提出基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法,该方法可以改善直角坐标牛顿法潮流计算分析含有电阻不为0的小阻抗支路电力系统的收敛性。
为了实现上述目的,本发明根据PV节点的电源无功功率不是给定值的特点,提出了基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法来改善潮流计算收敛性。本发明的首次迭代时PQ节点采用由给定值Pis和Qis计算得到的ai和ci值形成雅可比矩阵元素,首次迭代时PV节点则采用给定值Pis和QiL计算得到的ai和ci值形成雅可比矩阵元素;后续各次迭代时所有节点都采用传统方法计算雅可比矩阵元素。
本发明的技术方案如下:基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法,包括以下步骤:
A、输入原始数据和初始化电压;
B、形成节点导纳矩阵;
C、设置迭代计数t=0;
D、计算功率及电压偏差,求最大不平衡量ΔWmax;
E、判断最大不平衡量|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤I;否则,执行步骤F;
F、形成雅可比矩阵;
除首次迭代外,雅可比矩阵计算方法仍采用传统方法。首次迭代的雅可比矩阵计算方法根据节点类型采用不同方法。对于PQ节点,因潮流计算收敛时,计算功率Pi和Qi与给定值Pis和Qis相等,故采用式(25)计算注入电流相量的实部和虚部效果较好;由于PV节点的电源无功功率QiG不是给定的,原始数据输入时该值输入的是任意值,因此Qis与潮流计算收敛时的计算无功功率Qi相差很大,且误差不可预期,采用式(25)计算注入电流相量的实部和虚部效果可能较差,可采用式(26)计算。
式中,Pis为节点i给定的注入有功功率,QiL为节点i的负荷无功功率。
形成雅可比矩阵元素的具体步骤如下:
F1、按式(10)-(15)计算i≠j时的雅可比矩阵元素;
F2、令i=1;
F3、判断是否满足t=0,如果不满足转步骤F7;
F4、判断是否满足节点i是PQ节点的条件,如果不满足转步骤F5;如果满足,则按式(25)计算节点i的注入电流相量的实部ai和虚部ci,然后转步骤F8;
F5、判断是否满足节点i是PV节点的条件,如果不满足转步骤F9;
F6、按式(26)计算节点i的注入电流相量的实部ai和虚部ci,转步骤F8;
F7、按式(22)计算节点i的注入电流相量的实部ai和虚部ci;
F8、按式(16)-(21)计算i=j时的雅可比矩阵元素;
F9、令i=i+1;
F10、判断i是否大于节点数n,如果i不大于n转步骤F3;否则转步骤G;
G、解修正方程及修正各节点电压实部和虚部;
H、令t=t+1,返回步骤D进行下一次迭代;
I、输出节点及支路数据。
下面分两种情况分析首次迭代的情况。
(1)小阻抗支路两端节点都是PQ节点
首次迭代时,与小阻抗支路有关的修正方程为:
式中,Ai、Aj、Bi、Bj为与Δek、Δfk相关的项(k=1,…,n且k≠i,j);Pi0、Pj0、Qi0、Qj0为除小阻抗支路lij外节点的计算功率。
式(27)-(30)中考虑到首次迭代时,电压为电压初值,即电压初值实部为1.0,虚部为0.0。得:
|gij|和|bij|较大,与之相比,其它项较小,式(31)-(34)忽略较小量,得
-(gij/k2)Δei+(gij/k)Δej+(bij/k2)Δfi-(bij/k)Δfj≈-gij/k2+gij/k (35)
-gijΔej+(gij/k)Δei+bijΔfj-(bij/k)Δfi≈-gij+gij/k (36)
(bij/k2)Δei-(bij/k)Δej+(gij/k2)Δfi-(gij/k)Δfj≈bij/k2-bij/k (37)
bijΔej-(bij/k)Δei+gijΔfj-(gij/k)Δfi≈bij-bij/k (38)
式(35)乘以bij与式(37)乘以gij相加,得
式(39)中由于得
Δfi≈kΔfj (40)
由于初值则电压虚部修正后满足公式(2)。
式(37)乘以bij,再与式(35)乘以gij相减,得
式(41)中由于得
Δei/k2-Δej/k≈1/k2-1/k (42)
式(42)整理,得
(1-Δei)≈k(1-Δej) (43)
式(43)中,考虑电压实部初值首次迭代后电压实部为
式(44)满足公式(2)。
式(31)乘以k再加式(32),得
式(33)乘以k再加式(34),得
这样式(31)-(34)经过变换得到式(40)、(44)、(45)、(46),而式(40)、(44)、(45)、(46)已经不存在小阻抗了,且满足小阻抗支路两端电压关系式(2)。由于小阻抗的影响已经不存在了,因此首次迭代时端点都是PQ节点的小阻抗不会引起潮流计算的不收敛。
(2)小阻抗支路两端节点分别为PV节点和PQ节点
设支路首端节点i为PV节点,末端节点j为PQ节点。首次迭代时,与小阻抗支路有关的修正方程为
式(47)-(50)中考虑到首次迭代时,电压为电压初值,即节点i电压初值实部为Vis,虚部为0.0,节点j电压初值实部为1.0,虚部为0.0,得:
-2VisΔei=0 (53)
式(53)代入式(51)、式(52)、式(54),得
式(56)-(57)忽略较小量,得
-gijΔej+bijΔfj-(bij/k)Δfi≈-gij+gijVis/k (58)
bijΔej+gijΔfj-(gij/k)Δfi≈bij-bijVis/k (59)
式(58)乘以bij与式(59)乘以gij相加,得
式(60)中由于得
Δfi≈kΔfj (61)
由于初值则电压虚部修正后满足公式(2)。
式(59)乘以bij与式(58)乘以gij相减,得
式(62)中由于得
1-Δej≈Vis/k (63)
由式(63),再考虑电压实部初值以及修正量首次迭代后电压实部为
式(64)满足公式(2)。
式(56)乘以(Vis/k)再加式(55),得
式(61)和式(63)代入式(57),得
(cjs+Bj0)Δej+(-ajs+Gj0)Δfj+Bj≈Qjs+Bj0-Qj0 (66)
这样式(53)、式(55)~(57)经过变换得到式(61)、(64)、(65)、(66),而式(61)、(64)、(65)、(66)已经不存在小阻抗了,且满足小阻抗支路两端电压关系式(2)。由于小阻抗的影响已经不存在了,因此首次迭代时端点分别为PV节点和PQ节点的小阻抗不会引起潮流计算的不收敛。
下面分两种情况分析第2次迭代的情况。
(1)小阻抗支路两端节点都是PQ节点
第2次迭代时,与小阻抗支路有关的修正方程为:
把式(22)代入到式(67)-(70),得:
考虑到首次迭代以后,小阻抗支路两端节点电压已经满足 把此电压关系代入式(71)-(74),得:
式(75)-(78)忽略较小量,得
-(gijej+bijfj)Δei/k+(gijej+bijfj)Δej+(bijej-gijfj)Δfi/k+(gijfj-bijej)Δfj≈0 (79)
-(gijej+bijfj)Δej+(gijej+bijfj)Δei/k+(bijej-gijfj)Δfj+(gijfj-bijej)Δfi/k≈0 (80)
(bijej-gijfj)Δei/k+(gijfj-bijej)Δej+(gijej+bijfj)Δfi/k-(gijej+bijfj)Δfj≈0 (81)
(bijej-gijfj)Δej+(gijfj-bijej)Δei/k+(gijej+bijfj)Δfj-(gijej+bijfj)Δfi/k≈0 (82)
式(80)乘以bij与式(82)乘以gij相加,得
式(83)中由于得
-fjΔej+fjΔei/k+ejΔfj-ejΔfi/k≈0 (84)
式(82)乘以bij,再与式(80)乘以gij相减,得
式(85)中由于得
ejΔej-ejΔei/k+fjΔfj-fjΔfi/k≈0 (86)
式(84)乘以ej与式(86)乘以fj相加,得
式(87)中由于得
Δfi≈kΔfj (88)
由于首次迭代后有则修正后满足公式(2)。
式(88)代入式(84),得
Δei≈kΔej (89)
由于首次迭代后有则修正后满足公式(2)。
式(75)加式(76),得
式(77)加式(78),得
这样式(75)~(78)经过变换得到式(88)~(91),而式(88)~(91)已经不存在小阻抗了,且满足小阻抗支路两端电压关系式(2)。由于小阻抗的影响已经不存在了,因此第2次迭代时端点都是PQ节点的小阻抗不会引起潮流计算的不收敛。
(2)小阻抗支路两端节点分别为PV节点和PQ节点
设支路首端节点i为PV节点,末端节点j为PQ节点。第2次迭代时,与小阻抗支路有关的修正方程除第3个方程改为式(92)外,其它3个方程仍为式(75)、式(76)和式(78)。
除了如何得出关于无功偏差的方程式(91)外,其它推导过程同前。
式(88)和式(89)代入式(78),得
这样式(75)、式(76)、式(92)、式(78)经过变换得到式(88)、(89)、(92)、(93),而式(88)、(89)、(92)、(93)已经不存在小阻抗了,且满足小阻抗支路两端电压关系式(2)。由于小阻抗的影响已经不存在了,因此第2次迭代时端点分别为PV节点和PQ节点的小阻抗不会引起潮流计算的不收敛。
同理可证第2次后各次迭代时小阻抗不会引起潮流计算的不收敛。
综上所述本发明方法能够使得潮流计算收敛,但对PV节点的不同处理可能影响迭代次数。
由此可见,本发明解决了直角坐标牛顿法潮流计算在求解含有小阻抗支路系统时的收敛性问题。采用现有常规直角坐标牛顿法潮流计算不收敛时,本算法能够可靠收敛,且比现有专利技术迭代次数少。
与现有技术相比,本发明具有以下有益效果:
1、本发明通过首次迭代采用与以后各次迭代过程不同的雅可比矩阵计算方法,解决了直角坐标牛顿法潮流计算在求解含有小阻抗支路电力系统时的收敛性问题。在首次迭代过程中,PQ节点采用节点注入功率给定值计算出的电流相量的实部和虚部计算雅可比元素,而PV节点仅采用节点注入有功功率给定值和负荷无功功率给定值计算出的电流相量的实部和虚部计算雅可比元素,不使用电源无功功率给定值,避免了PV节点的电源无功功率给定值不确定给潮流计算收敛性的不利影响。采用常规直角坐标牛顿法潮流计算不收敛时,本发明能够可靠收敛,且比现有专利技术迭代次数少。
2、由于本发明不仅能有效解决常规直角坐标牛顿法潮流计算分析含有小阻抗支路电力系统的收敛性问题,同时也能对正常电力系统进行潮流计算,没有不良影响。
附图说明
本发明共有附图4张。其中:
图1是电力系统小阻抗变压器模型示意图。
图2是直角坐标牛顿法潮流计算的流程图。
图3是现有发明直角坐标牛顿法潮流计算的流程图。
图4是本发明直角坐标牛顿法潮流计算雅可比矩阵形成方法的流程图。
具体实施方式
下面结合附图对本发明进行进一步地说明。根据图1所示的小阻抗变压器模型,采用图2所示的直角坐标牛顿法潮流计算的流程图和图4所示的直角坐标牛顿法潮流计算雅可比矩阵形成方法的流程图,对一个实际大型电网进行了潮流计算。该实际大型电网有445个节点,含有大量的小阻抗支路。其中,x≤0.001的小阻抗支路有49条,x≤0.0001的小阻抗支路有41条,x≤0.00001的小阻抗支路有22条。其中阻抗值最小的是节点118和节点125之间的小阻抗支路l118-125为x=0.00000001,变比k=0.9565,k位于节点118侧。潮流计算的收敛精度为0.00001。为了验证本发明计算含电阻不为0的小阻抗支路电力系统的收敛性,把小阻抗支路l118-125、l60-122及l287-310的电阻改为r=0.0001。
作为对比,同时采用常规直角坐标牛顿法潮流算法及现有专利算法(专利号为ZL201410315785.1,如图3所示)对该实际大型电网进行了潮流计算,迭代次数见表1。
表1不同潮流方法的迭代结果
方法 | 常规算法 | ZL201410315785.1算法 | 本发明算法 |
迭代结果 | 不收敛 | 7次收敛 | 4次收敛 |
由表1可见,对于修改后的445节点实际电力系统算例,常规直角坐标牛顿法潮流算法不收敛,本发明算法和专利ZL201410315785.1算法都能够收敛,但本发明算法的迭代次数少3次。
不同潮流计算方法各次迭代最大不平衡量见表2。单位为标幺值。
表2不同潮流方法各次迭代最大不平衡量
迭代序号 | 常规算法 | ZL201410315785.1算法 | 本发明算法 |
0 | -4754.570367135 | -4754.570367135 | -4754.570367135 |
1 | -3451593.823720038 | -11.138394991 | -0.507179900 |
2 | -886651.468310079 | -6.163450054 | -0.068516394 |
3 | -222023.112200678 | -1.441071252 | -0.001325088 |
4 | -55754.415245002 | -0.106199006 | -0.000000792 |
5 | -13972.568194423 | -0.006353455 | |
6 | -6386.835620506 | -0.000141863 | |
7 | -6585.38761914 | -0.000000062 | |
8 | -378994.776907351 | ||
9 | -98508.025841226 | ||
10 | -37917.863557986 |
由表2可知,3种方法首次迭代前最大不平衡量相同且很大。首次迭代后,现有专利方法最大不平衡量明显减少,迭代7次收敛;本发明方法最大不平衡量减少速度更快,迭代4次收敛;而常规方法的最大不平衡量则变大,最终发散。
本发明可以采用任何一种编程语言和编程环境实现,如C语言、C++、FORTRAN、Delphi等。开发环境可以采用Visual C++、Borland C++Builder、Visual FORTRAN等。
本发明不局限于本实施例,任何在本发明披露的技术范围内的等同构思或者改变,均列为本发明的保护范围。
Claims (1)
1.基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法,其特征在于:包括以下步骤:
A、输入原始数据和初始化电压;
根据电力系统节点的特点,潮流计算把电力系统节点分成3类:节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点;节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点;节点电压幅值和电压相角已知,节点有功功率和无功功率未知的节点称为平衡节点;
初始化电压采用平启动,即PV节点和平衡节点的电压实部取给定值,PQ节点的电压实部取1.0;所有电压的虚部都取0.0;这里单位采用标幺值;
B、形成节点导纳矩阵;
设节点i和节点j原来的自电导与自电纳分别为Gi0、Bi0、Gj0、Bj0,在它们之间增加一条小阻抗支路后的自导纳和互导纳分别为:
式中,Yii、Yjj分别为节点i和节点j的自导纳;Yij为节点i和节点j之间的互导纳;rij、xjj分别为节点i和节点j之间小阻抗支路的电阻和电抗;k为节点i和节点j之间小阻抗支路的变比,如果是输电线支路,则变比k为1;
C、设置迭代计数t=0;
D、计算功率及电压偏差,求最大不平衡量ΔWmax;
PQ节点的功率偏差计算公式为:
式中,Pis、Qis分别为节点i给定的注入有功功率和无功功率;ei、fi分别为节点i的电压相量的实部和虚部;ai、ci分别为节点i的计算注入电流相量的实部和虚部,具体表达式为:
式中,n为电力系统的节点数;Gim、Bim分别为节点导纳矩阵元素Yim的实部和虚部;
PV节点的有功功率及电压偏差计算公式为:
式中,Vis为节点i给定的电压幅值;
平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差;
E、判断最大不平衡量|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤I;否则,执行步骤F;
其特征在于:还包括以下步骤:
F、形成雅可比矩阵J;
F1、按式(7)-(12)计算i≠j时的雅可比矩阵元素;
当i≠j时,雅可比矩阵J的元素计算公式如下:
F2、令i=1;
F3、判断是否满足t=0,如果不满足转步骤F7;
F4、判断是否满足节点i是PQ节点的条件,如果不满足转步骤F5;如果满足,则按式(13)计算节点i的注入电流相量的实部ai和虚部ci,然后转步骤F8;
F5、判断是否满足节点i是PV节点的条件,如果不满足转步骤F9;
F6、按式(14)计算节点i的注入电流相量的实部ai和虚部ci,转步骤F8;
式中,Pis为节点i给定的注入有功功率,QiL为节点i的负荷无功功率;
F7、按式(5)计算节点i的注入电流相量的实部ai和虚部ci;
F8、按式(15)-(20)计算i=j时的雅可比矩阵元素;
F9、令i=i+1;
F10、判断i是否大于节点数n,如果i不大于n转步骤F3;否则转步骤G;
G、解修正方程及修正电压实部e、虚部f;
修正方程为:
式中,ΔP为有功功率偏差列向量;ΔQ为无功功率偏差列向量;ΔV2为电压偏差列向量;Δe为电压相量实部修正量列向量;Δf为电压相量虚部修正量列向量;J为雅可比矩阵;
电压修正公式为:
式中,上标(t)表示第t次迭代;
H、令t=t+1,返回步骤D进行下一次迭代;
I、输出节点及支路数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811626642.7A CN109617080B (zh) | 2018-12-28 | 2018-12-28 | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811626642.7A CN109617080B (zh) | 2018-12-28 | 2018-12-28 | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109617080A true CN109617080A (zh) | 2019-04-12 |
CN109617080B CN109617080B (zh) | 2021-10-26 |
Family
ID=66011924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811626642.7A Active CN109617080B (zh) | 2018-12-28 | 2018-12-28 | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109617080B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110601202A (zh) * | 2019-08-28 | 2019-12-20 | 潘协印 | 一种电力系统潮流计算方法 |
CN111900738A (zh) * | 2020-08-14 | 2020-11-06 | 广东电网有限责任公司广州供电局 | 基于补偿算法的三相不平衡潮流计算方法 |
CN112260288A (zh) * | 2020-10-22 | 2021-01-22 | 云南电网有限责任公司 | 一种电力系统中节点电压的调整方法及装置 |
CN113297813A (zh) * | 2021-05-21 | 2021-08-24 | 华南理工大学 | 基于改进三阶牛顿迭代法的电-气能源系统潮流计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027642A1 (en) * | 2005-07-15 | 2007-02-01 | Chang Gung University | Method for Calculating Power Flow Solution of a Power Transmission Network that Includes Interline Power Flow Controller (IPFC) |
CN104037764A (zh) * | 2014-07-03 | 2014-09-10 | 大连海事大学 | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
CN106356859A (zh) * | 2016-09-29 | 2017-01-25 | 大连海事大学 | 一种基于Matlab的直角坐标牛顿法潮流计算方法 |
CN106532711A (zh) * | 2016-11-30 | 2017-03-22 | 大连海事大学 | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 |
CN106712029A (zh) * | 2016-12-09 | 2017-05-24 | 大连海事大学 | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 |
CN109066693A (zh) * | 2018-09-04 | 2018-12-21 | 河南理工大学 | 一种具有高阶收敛速率的柔性系统潮流计算方法 |
-
2018
- 2018-12-28 CN CN201811626642.7A patent/CN109617080B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027642A1 (en) * | 2005-07-15 | 2007-02-01 | Chang Gung University | Method for Calculating Power Flow Solution of a Power Transmission Network that Includes Interline Power Flow Controller (IPFC) |
CN104037764A (zh) * | 2014-07-03 | 2014-09-10 | 大连海事大学 | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
CN106356859A (zh) * | 2016-09-29 | 2017-01-25 | 大连海事大学 | 一种基于Matlab的直角坐标牛顿法潮流计算方法 |
CN106532711A (zh) * | 2016-11-30 | 2017-03-22 | 大连海事大学 | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 |
CN106712029A (zh) * | 2016-12-09 | 2017-05-24 | 大连海事大学 | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 |
CN109066693A (zh) * | 2018-09-04 | 2018-12-21 | 河南理工大学 | 一种具有高阶收敛速率的柔性系统潮流计算方法 |
Non-Patent Citations (4)
Title |
---|
YUBIN YAO, BENHUA GU: "A Newton Power Flow for Study Purpose", 《2012 INTERNATIONAL CONFERENCE ON FUTURE ENERGY, ENVIRONMENT, AND MATERIALS》 * |
YUBIN YAO, JUQING SHAN: "A Improved Newton Power Flow in Rectangular Form for Systems with Small Impedance Branches", 《2009 2ND INTERNATIONAL CONFERENCE ON POWER ELECTRONICS AND INTELLIGENT TRANSPORTATION SYSTEM (PEITS)》 * |
姚玉斌, 刘莉: "基于快速分解法的连续潮流法", 《哈尔滨工业大学学报》 * |
覃智君: "大规模交直流系统潮流计算的实用化模型", 《中国电机工程学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110601202A (zh) * | 2019-08-28 | 2019-12-20 | 潘协印 | 一种电力系统潮流计算方法 |
CN111900738A (zh) * | 2020-08-14 | 2020-11-06 | 广东电网有限责任公司广州供电局 | 基于补偿算法的三相不平衡潮流计算方法 |
CN111900738B (zh) * | 2020-08-14 | 2022-02-11 | 广东电网有限责任公司广州供电局 | 基于补偿算法的三相不平衡潮流计算方法 |
CN112260288A (zh) * | 2020-10-22 | 2021-01-22 | 云南电网有限责任公司 | 一种电力系统中节点电压的调整方法及装置 |
CN112260288B (zh) * | 2020-10-22 | 2023-05-23 | 云南电网有限责任公司 | 一种电力系统中节点电压的调整方法及装置 |
CN113297813A (zh) * | 2021-05-21 | 2021-08-24 | 华南理工大学 | 基于改进三阶牛顿迭代法的电-气能源系统潮流计算方法 |
CN113297813B (zh) * | 2021-05-21 | 2022-07-05 | 华南理工大学 | 基于改进三阶牛顿迭代法的电-气能源系统潮流计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109617080B (zh) | 2021-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN104037764B (zh) | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 | |
CN109617080A (zh) | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 | |
US8326594B2 (en) | Power flow analysis for balanced power distribution systems | |
CN104037763B (zh) | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 | |
CN106709243B (zh) | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 | |
CN109494748A (zh) | 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法 | |
CN106532712B (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
CN103956735B (zh) | 一种分布式发电系统的谐波潮流分析方法 | |
CN104022507B (zh) | 一种直角坐标牛顿法潮流计算方法 | |
CN104113061A (zh) | 一种含分布式电源的配电网三相潮流计算方法 | |
CN104636829A (zh) | 一种提高电力系统温度最优潮流计算效率的解耦算法 | |
CN106856327B (zh) | 一种含小阻抗支路电网串联补偿快速分解法潮流计算方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN105046588A (zh) | 一种基于网损迭代的改进直流动态最优潮流的计算方法 | |
CN102780220A (zh) | 含有pv恒定型分布式电源的配电网潮流计算方法 | |
CN106712029B (zh) | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 | |
CN107994567B (zh) | 一种广义快速分解状态估计方法 | |
CN111049144A (zh) | 基于定变化率矩阵的有源配电网潮流计算pv节点处理方法 | |
CN111355236A (zh) | 一种计及中性点电压变量的配电网三相潮流计算方法 | |
CN106529089B (zh) | 用于含小阻抗支路电网的补偿法快速分解法潮流计算方法 | |
CN108183487B (zh) | 一种基于线性数学模型的配电网潮流快速分析方法 | |
CN106786605B (zh) | 一种含小阻抗电网串联补偿直角坐标牛顿法潮流计算方法 | |
CN104218577A (zh) | 一种基于节点电压的主动配电网三相潮流计算方法 | |
CN107658880B (zh) | 基于关联矩阵运算的快速分解法系数矩阵计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220119 Address after: 116026 No. 1, Ling Hai Road, Liaoning, Dalian Patentee after: Dalian Maritime University Patentee after: Shanghai Shipping Research Institute Co., Ltd Address before: 116026 1 Linghai Road, hi tech park, Dalian, Liaoning Patentee before: Dalian Maritime University |