CN106532711A - 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 - Google Patents
随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 Download PDFInfo
- Publication number
- CN106532711A CN106532711A CN201611094297.8A CN201611094297A CN106532711A CN 106532711 A CN106532711 A CN 106532711A CN 201611094297 A CN201611094297 A CN 201611094297A CN 106532711 A CN106532711 A CN 106532711A
- Authority
- CN
- China
- Prior art keywords
- node
- delta
- voltage
- power
- formula
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法,首次迭代时PQ节点采用由给定值Pis和Qis计算的ai和bi值计算雅可比矩阵元素,PV节点以及后续各次迭代时所有节点则采用传统方法计算雅可比矩阵元素。本发明通过在首次迭代过程中PQ节点采用与以后各次迭代过程不同的雅可比矩阵计算方法,解决了直角坐标牛顿法潮流计算在分析含有小阻抗支路电力系统时的收敛性问题。采用常规直角坐标牛顿法潮流计算不收敛时,本发明能可靠收敛,且比现有专利技术迭代次数少。由于本发明不仅能有效解决常规直角坐标牛顿法潮流计算分析含有小阻抗支路电力系统的收敛性问题,同时也能对正常电力系统进行潮流计算,没有不良影响。
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节点的功率偏差计算公式为:
式中,Pis、Qis分别为节点i给定的注入有功功率和无功功率,Pis为电源有功功率与负荷有功功率之差,Qis为电源无功功率与负荷无功功率之差;ai、bi分别为节点i的计算注入电流相量的实部和虚部,为
式中,n为电力系统的节点数。
潮流计算收敛时,式(6)中ΔPi、ΔQi都趋近于0,因此ai和bi等于由给定值Pis和Qis计算出的ais和bis
PV节点的有功功率及电压偏差计算公式为:
式中,Vis为节点i给定的电压幅值。
平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差。
平衡节点的电源有功功率和无功功率以及PV节点的电源无功功率采用式(6)计算。
计算各节点功率或电压偏差的最大值,称为最大不平衡量,如果最大不平衡量的绝对值小于给定的收敛精度,转步骤F,否则执行步骤D。
D、形成雅可比矩阵J
雅可比矩阵J的元素(i≠j时)计算公式如下:
PQ节点按式(10)-(13)计算雅可比矩阵元素;PV节点按式(10)、(11)、(14)、和(15)计算雅可比矩阵元素;平衡节点不计算雅可比矩阵元素。
雅可比矩阵J的元素(i=j时)计算公式如下:
PQ节点按式(16)-(19)计算雅可比矩阵元素;PV节点按式(16)、(17)、(20)、(21)计算雅可比矩阵元素;平衡节点不计算雅可比矩阵元素。
E、解修正方程及修正电压实部e、虚部f
修正方程为:
式中,J为雅可比矩阵。
电压修正公式为:
式中,上标(t)表示第t次迭代。
F、输出节点及支路数据。
对正常电力网络,牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,牛顿法潮流计算就可能发散。而电力系统中小阻抗支路普遍存在,收敛性是电力系统潮流计算这类非线性问题的最重要指标,计算不收敛就无法得到方程的解。因此改善直角坐标牛顿法潮流计算针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。
中国专利ZL201410299531.5披露了一种通过修改常规直角坐标牛顿法潮流计算雅可比矩阵的方法,该方法计算雅可比元素时采用由给定值Pis和Qis计算的ai和bi值,改善了潮流计算的收敛性,有效解决了含有电阻为0的小阻抗支路电力系统潮流计算的发散问题。但当小阻抗支路的电阻不为0时,该方法迭代次数增加,收敛性变差,甚至不收敛。
中国专利ZL201410315785.1提出了一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法,该方法首次迭代和后续各次迭代采用不同的雅可比矩阵计算方法,首次迭代计算雅可比元素时采用由给定值Pis和Qis计算的ai和bi值,后续各次迭代计算雅可比元素时仍采用传统方法,有效解决了含电阻不为0的小阻抗支路电力系统潮流计算的发散问题,但该方法在电阻不为0的小阻抗支路较多时,迭代次数增加,收敛性变差。
发明内容
为解决现有技术存在的上述问题,本发明要提出一种随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法,该方法可以改善其分析含有电阻不为0的小阻抗支路电力系统的收敛性。
为了实现上述目的,本发明根据PV节点的电源无功功率不是给定值的特点,提出了一种直角坐标牛顿法潮流计算方法来改善潮流计算收敛性。本发明的首次迭代时PQ节点采用由给定值Pis和Qis计算的ai和bi值计算雅可比矩阵元素,首次迭代时PV节点以及后续各次迭代时所有节点则采用传统方法计算雅可比矩阵元素。本发明的技术方案如下:随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法,包括以下步骤:
A、输入原始数据和初始化电压;
B、形成节点导纳矩阵;
C、设置迭代计数t=0;
D、计算功率及电压偏差,求最大不平衡量ΔWmax;
E、判断最大不平衡量|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤I;否则,执行步骤F;
F、形成雅可比矩阵;
除首次迭代外,雅可比矩阵计算方法仍采用传统方法。首次迭代的雅可比矩阵计算方法根据节点类型采用不同方法。对于PQ节点,因潮流计算收敛时,计算功率Pi和Qi与给定值Pis和Qis相等,故采用式(8)计算注入电流相量的实部和虚部效果较好;由于PV节点的电源无功功率不是给定的,原始数据输入时该值输入是任意值,因此Qis与潮流计算收敛时的计算无功功率Qi相差很大,采用式(8)计算注入电流相量的实部和虚部效果可能较差,应该采用式(7)计算。
形成雅可比矩阵元素的具体步骤如下:
F1、按式(10)-(15)计算i≠j时的雅可比矩阵元素;
F2、令i=1;
F3、判断是否同时满足t=0和节点i是PQ节点的条件,如果不满足转步骤F4;如果满足,则按式(8)计算节点i的注入电流相量的实部ai和虚部bi,然后转步骤F5;
F4、按式(7)计算节点i的注入电流相量的实部ai和虚部bi;
F5、按式(16)-(21)计算i=j时的雅可比矩阵元素;
F6、令i=i+1;
F7、判断i是否大于节点数n,如果i不大于n转步骤F3;否则转步骤G;
G、解修正方程及修正电压实部e、虚部f;
H、令t=t+1,返回步骤D进行下一次迭代;
I、输出节点及支路数据。
与现有技术相比,本发明具有以下有益效果:
1、本发明通过在首次迭代过程中PQ节点采用与以后各次迭代过程不同的雅可比矩阵计算方法,解决了直角坐标牛顿法潮流计算在分析含有小阻抗支路电力系统时的收敛性问题。采用常规直角坐标牛顿法潮流计算不收敛时,本算法能够可靠收敛,且比现有专利技术迭代次数少。
2、由于本发明不仅能有效解决常规直角坐标牛顿法潮流计算分析含有小阻抗支路电力系统的收敛性问题,同时也能对正常电力系统进行潮流计算,没有不良影响。
附图说明
本发明共有附图4张。其中:
图1是电力系统小阻抗变压器模型示意图。
图2是直角坐标牛顿法潮流计算的流程图。
图3是现有发明直角坐标牛顿法潮流计算的流程图。
图4是本发明直角坐标牛顿法潮流计算的流程图。
具体实施方式
下面结合附图对本发明进行进一步地说明。根据图1所示的小阻抗变压器模型,采用图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)对该实际大型电网进行了潮流计算,迭代次数见表1。
表1 不同潮流方法的迭代结果
方法 | 常规算法 | ZL201410315785.1算法 | 本发明算法 |
迭代结果 | 不收敛 | 7次收敛 | 6次收敛 |
由表1可见,对于修改后的445节点实际电力系统算例,常规直角坐标牛顿法潮流算法不收敛,本发明算法和专利ZL201410315785.1算法都能够收敛,但本发明算法的迭代次数少1次。
不同潮流计算方法各次迭代最大不平衡量见表2。单位为标幺值。
表2 不同潮流方法各次迭代最大不平衡量
迭代序号 | 常规算法 | ZL201410315785.1算法 | 本发明算法 |
0 | -4754.570367135 | -4754.570367135 | -4754.570367135 |
1 | -3451593.823720038 | -11.138394991 | -3.264368583 |
2 | -886651.468310079 | -6.163450054 | -0.715148045 |
3 | -222023.112200678 | -1.441071252 | -0.076847277 |
4 | -55754.415245002 | -0.106199006 | -0.002294590 |
5 | -13972.568194423 | -0.006353455 | -0.000017499 |
6 | -6386.835620506 | -0.000141863 | -0.000000001 |
7 | -6585.38761914 | -0.000000062 | |
8 | -378994.776907351 | ||
9 | -98508.025841226 | ||
10 | -37917.863557986 |
由表2可知,3种方法首次迭代前最大不平衡量相同且很大。首次迭代后,现有专利方法最大不平衡量明显减少,迭代7次收敛;本专利方法最大不平衡量减少速度更快,迭代6次收敛;而常规方法的最大不平衡量则变大,最终发散。
PV节点的电源无功功率输入值与潮流计算的计算值及差值见表3,差值为输入值与计算值的差。单位为标幺值。
表3 PV节点的电源无功功率输入值与计算值及差值
由表3可见,PV节点的电源无功功率输入值与潮流计算的计算值有较大的差值,有的差值达到1.6以上,用此输入值计算出的注入电流相量的实部和虚部与潮流计算收敛时的值必定有很大差距,影响潮流计算的收敛性。
本算法可以采用任何一种编程语言和编程环境实现,如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给定的注入有功功率和无功功率,Pis为电源有功功率与负荷有功功率之差,Qis为电源无功功率与负荷无功功率之差;ei、fi分别为节点i的电压相量的实部和虚部;ai、bi分别为节点i的计算注入电流相量的实部和虚部,具体表达式为:
式中,n为电力系统的节点数;Gim、Bim分别为节点i和节点m之间的互导纳的实部和虚部,如果m=i时,Gii、Bii分别为节点i的自导纳的实部和虚部;
PV节点的有功功率及电压偏差计算公式为:
式中,Vis为节点i给定的电压幅值;
平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差;
E、判断最大不平衡量|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤I;否则,执行步骤F;
其特征在于:还包括以下步骤:
F、形成雅可比矩阵J;
F1、按式(7)-(12)计算i≠j时的雅可比矩阵元素;
当i≠j时,雅可比矩阵J的元素计算公式如下:
F2、令i=1;
F3、判断是否同时满足t=0和节点i是PQ节点的条件,如果不满足则转步骤F4;如果满足,则按式(13)计算节点i的注入电流相量的实部ai和虚部bi,然后转步骤F5;
F4、按式(5)计算节点i的注入电流相量的实部ai和虚部bi;
F5、按式(14)-(19)计算i=j时的雅可比矩阵元素;
F6、令i=i+1;
F7、判断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 |
---|---|---|---|
CN201611094297.8A CN106532711B (zh) | 2016-11-30 | 2016-11-30 | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611094297.8A CN106532711B (zh) | 2016-11-30 | 2016-11-30 | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106532711A true CN106532711A (zh) | 2017-03-22 |
CN106532711B CN106532711B (zh) | 2019-02-22 |
Family
ID=58354640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611094297.8A Expired - Fee Related CN106532711B (zh) | 2016-11-30 | 2016-11-30 | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106532711B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107957974A (zh) * | 2017-11-29 | 2018-04-24 | 广西大学 | 一种提取极坐标潮流方程雅可比矩阵的方法 |
CN109494748A (zh) * | 2018-12-28 | 2019-03-19 | 大连海事大学 | 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法 |
CN109617080A (zh) * | 2018-12-28 | 2019-04-12 | 大连海事大学 | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 |
CN109980648A (zh) * | 2018-09-07 | 2019-07-05 | 山东大学 | 交直流混联潮流计算的方法及其装置、存储介质及终端 |
CN110707701A (zh) * | 2019-09-19 | 2020-01-17 | 东北电力大学 | 一种电力系统潮流计算方法 |
CN111046336A (zh) * | 2019-12-27 | 2020-04-21 | 大连海事大学 | 一种雅可比矩阵变化的极坐标牛顿法潮流计算方法 |
CN112600200A (zh) * | 2020-12-08 | 2021-04-02 | 深圳供电局有限公司 | 电网节点的谐波潮流计算方法、装置、设备及存储介质 |
CN113297813A (zh) * | 2021-05-21 | 2021-08-24 | 华南理工大学 | 基于改进三阶牛顿迭代法的电-气能源系统潮流计算方法 |
CN113313369A (zh) * | 2021-05-21 | 2021-08-27 | 华南理工大学 | 基于改进雅可比矩阵的热网潮流计算方法 |
CN117154735A (zh) * | 2023-08-30 | 2023-12-01 | 国网江苏省电力有限公司扬州供电分公司 | 一种适用于新能源发电台区的潮流计算方法、设备及介质 |
CN111030122B (zh) * | 2019-12-27 | 2023-12-22 | 大连海事大学 | 首次迭代雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104022507A (zh) * | 2014-06-26 | 2014-09-03 | 大连海事大学 | 一种直角坐标牛顿法潮流计算方法 |
CN104037763A (zh) * | 2014-07-02 | 2014-09-10 | 大连海事大学 | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 |
CN104037764A (zh) * | 2014-07-03 | 2014-09-10 | 大连海事大学 | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
-
2016
- 2016-11-30 CN CN201611094297.8A patent/CN106532711B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104022507A (zh) * | 2014-06-26 | 2014-09-03 | 大连海事大学 | 一种直角坐标牛顿法潮流计算方法 |
CN104037763A (zh) * | 2014-07-02 | 2014-09-10 | 大连海事大学 | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 |
CN104037764A (zh) * | 2014-07-03 | 2014-09-10 | 大连海事大学 | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107957974A (zh) * | 2017-11-29 | 2018-04-24 | 广西大学 | 一种提取极坐标潮流方程雅可比矩阵的方法 |
CN107957974B (zh) * | 2017-11-29 | 2020-12-18 | 广西大学 | 一种提取极坐标潮流方程雅可比矩阵的方法 |
CN109980648A (zh) * | 2018-09-07 | 2019-07-05 | 山东大学 | 交直流混联潮流计算的方法及其装置、存储介质及终端 |
CN109980648B (zh) * | 2018-09-07 | 2023-04-07 | 山东大学 | 交直流混联潮流计算的方法及其装置、存储介质及终端 |
CN109617080A (zh) * | 2018-12-28 | 2019-04-12 | 大连海事大学 | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 |
CN109494748A (zh) * | 2018-12-28 | 2019-03-19 | 大连海事大学 | 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法 |
CN109494748B (zh) * | 2018-12-28 | 2021-10-26 | 大连海事大学 | 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法 |
CN110707701A (zh) * | 2019-09-19 | 2020-01-17 | 东北电力大学 | 一种电力系统潮流计算方法 |
CN110707701B (zh) * | 2019-09-19 | 2022-01-28 | 东北电力大学 | 一种电力系统潮流计算方法 |
CN111046336A (zh) * | 2019-12-27 | 2020-04-21 | 大连海事大学 | 一种雅可比矩阵变化的极坐标牛顿法潮流计算方法 |
CN111030122B (zh) * | 2019-12-27 | 2023-12-22 | 大连海事大学 | 首次迭代雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
CN112600200A (zh) * | 2020-12-08 | 2021-04-02 | 深圳供电局有限公司 | 电网节点的谐波潮流计算方法、装置、设备及存储介质 |
CN113297813A (zh) * | 2021-05-21 | 2021-08-24 | 华南理工大学 | 基于改进三阶牛顿迭代法的电-气能源系统潮流计算方法 |
CN113313369A (zh) * | 2021-05-21 | 2021-08-27 | 华南理工大学 | 基于改进雅可比矩阵的热网潮流计算方法 |
CN117154735A (zh) * | 2023-08-30 | 2023-12-01 | 国网江苏省电力有限公司扬州供电分公司 | 一种适用于新能源发电台区的潮流计算方法、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106532711B (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN104037764B (zh) | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 | |
CN104037763B (zh) | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 | |
CN106709243B (zh) | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 | |
CN104899396B (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
JP2012139090A (ja) | 平衡配電系統のための電力潮流解析 | |
CN109617080B (zh) | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 | |
CN104022507B (zh) | 一种直角坐标牛顿法潮流计算方法 | |
CN106532712B (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
CN103632046A (zh) | 一种电网潮流计算方法 | |
CN101621200B (zh) | 一种直角坐标牛顿法潮流计算方法 | |
CN106856327B (zh) | 一种含小阻抗支路电网串联补偿快速分解法潮流计算方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN109494748A (zh) | 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法 | |
Zhao et al. | Holomorphic embedding power flow for AC/DC hybrid power systems using Bauer's eta algorithm | |
CN106712029B (zh) | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 | |
CN106529089B (zh) | 用于含小阻抗支路电网的补偿法快速分解法潮流计算方法 | |
CN106786605B (zh) | 一种含小阻抗电网串联补偿直角坐标牛顿法潮流计算方法 | |
CN106786604B (zh) | 一种含小阻抗电网的串联补偿极坐标牛顿法潮流计算方法 | |
CN107658880B (zh) | 基于关联矩阵运算的快速分解法系数矩阵计算方法 | |
CN103904649A (zh) | 一种适用于含不同阶次模型的发电机聚合方法 | |
CN107846021B (zh) | 一种广义快速分解潮流方法 | |
CN106159946A (zh) | 变电站内无功设备对变压器绕组无功功率影响的计算方法 | |
Matthiss et al. | Fast Probabilistic Load Flow for Non-Radial Distribution Grids | |
CN111030122B (zh) | 首次迭代雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190222 Termination date: 20191130 |