CN106532712A - 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 - Google Patents
含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 Download PDFInfo
- Publication number
- CN106532712A CN106532712A CN201611131677.4A CN201611131677A CN106532712A CN 106532712 A CN106532712 A CN 106532712A CN 201611131677 A CN201611131677 A CN 201611131677A CN 106532712 A CN106532712 A CN 106532712A
- Authority
- CN
- China
- Prior art keywords
- node
- delta
- power
- voltage
- impedance
- 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)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (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
本发明公开了一种含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法,在进行潮流计算之前,先对小阻抗支路进行处理,把此小阻抗支路变成两个支路串联的形式,其中一个支路的阻抗为z1=rij+j(xij+xc),另一个支路的阻抗为z2=‑jxc,这里xc为电力系统正常电抗均值xav。经过串联补偿把一个小阻抗变成两个阻抗较大的支路,可以提高潮流计算的收敛性。本发明能够计算任意阻抗的小阻抗支路,甚至阻抗为零的支路。本发明是对小阻抗支路数据进行处理,此处理过程是在潮流计算之前,对数据进行预处理,不需要改变潮流计算程序。因此本发明特别适合于对老的潮流计算程序进行改造。
Description
技术领域
本发明涉及含小阻抗支路电力系统的潮流计算方法,特别是一种电力系统的直角坐标牛顿法潮流计算方法。
背景技术
电力系统潮流计算是研究电力系统稳态运行的一项基本计算,它根据电力系统给定的运行条件和网络结构确定整个电力系统的运行状态。潮流计算也是电力系统其他分析的基础,如安全分析、暂态稳定分析等都要用到潮流计算。由于具有收敛可靠、计算速度较快及内存需求适中的优点,牛顿法成为当前潮流计算的主流方法。牛顿法分为极坐标和直角坐标两种形式,两种形式的牛顿法潮流计算都在电力系统中得到了广泛的应用。
在直角坐标牛顿法潮流计算中,节点i的电压采用直角坐标表示为:
对正常电力网络,直角坐标牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,直角坐标牛顿法潮流计算就可能发散。电力系统小阻抗支路可分为小阻抗线路和小阻抗变压器支路,在数学模型上线路可以看作变比为1:1的变压器,因此下面分析时仅以小阻抗变压器支路为例分析。小阻抗变压器模型见图1,变压器的非标准变比k位于节点i侧,阻抗位于标准变比侧。变压器阻抗zij=rij+jxij很小,导纳为
式中,yij、gij、bij分别为节点i和节点j之间小阻抗支路的导纳、电导和电纳;rij、xij分别为节点i和节点j之间小阻抗支路的电阻和电抗。
由于小阻抗支路li-j的阻抗很小,支路的电压降也很小,因此变压器两端节点的电压应满足:
如图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、xij分别为节点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给定的电压幅值。
平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差。
求各节点功率或电压偏差中绝对值最大的值,称为最大不平衡量,如果最大不平衡量的绝对值小于给定的收敛精度,转步骤F,否则执行步骤D。
D、形成雅可比矩阵J
雅可比矩阵J的元素(i≠j时)计算公式如下:
雅可比矩阵J的元素(i=j时)计算公式如下:
PQ节点按式(16)-(19)计算雅可比矩阵元素;PV节点按式(16)、(17)、(20)、(21)计算雅可比矩阵元素。
平衡节点不计算雅可比矩阵元素。
E、解修正方程及修正电压实部e、虚部f
潮流计算的基本方程(6)和(9)是非线性方程组,通常采用逐次线性化方法迭代求解。线性化得到的方程称为修正方程,用来求电压实部和虚部的修正量。
修正方程为:
式中,J为雅可比矩阵;ΔP和ΔQ分别为有功功率和无功功率偏差列向量;ΔV2为电压幅值偏差列向量;Δe和Δf分别为电压相量的实部和虚部修正量列向量;为有功功率偏差函数列向量对电压相量实部列向量转置的偏导矩阵,上标T为转置符号。
电压修正公式为:
式中,上标t表示第t次迭代。
F、输出节点及支路数据。
对正常电力网络,牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,牛顿法潮流计算就可能发散。而电力系统中小阻抗支路普遍存在,收敛性是电力系统潮流计算这类非线性问题的最重要指标,计算不收敛就无法得到方程的解。因此改善直角坐标牛顿法潮流计算针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。
中国专利ZL201611094297.8提出了一种随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法,该方法首次迭代时所有PQ节点与后续各次迭代采用不同的雅可比矩阵计算方法,首次迭代时所有PQ节点计算雅可比元素时采用由给定值Pis和Qis计算的ai和bi值,首次迭代时所有PV节点及后续各次迭代时所有节点计算雅可比元素时仍采用传统方法,有效解决了电力系统包含小阻抗支路的潮流计算的发散问题,但迭代次数仍然较多,需要进一步改进。
现在仍然在使用的一些直角坐标牛顿法潮流计算程序不具备处理小阻抗支路的能力,但对这些老的潮流计算程序进行改造往往比较麻烦。有些程序可能采用比较老的编程语言编写,不易修改,有的潮流计算版本甚至没有源文件,无法修改。
发明内容
为解决现有技术存在的上述问题,本发明要提出一种含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法,以改善直角坐标牛顿法潮流计算方法分析含有小阻抗支路电力系统的潮流计算的收敛性。
为了实现上述目的,本发明的技术方案如下:含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法,在进行潮流计算之前,先对小阻抗支路进行处理,把此小阻抗支路变成两个支路串联的形式,其中一个支路的阻抗为z1=rij+j(xij+xc),另一个支路的阻抗为z2=-jxc,这里xc为电力系统中电抗绝对值大于小电抗阈值的支路的电抗绝对值的平均值,称为电力系统正常电抗均值xav。经过串联补偿把一个小阻抗变成两个阻抗较大的支路,可以提高潮流计算的收敛性。方案包括以下步骤:
A、输入原始数据和初始化电压;
B、对小阻抗支路进行串联补偿;
用电力系统正常电抗均值对小阻抗支路进行补偿,得到两个阻抗较大的支路,其中一个支路的阻抗为z1=r+j(x+xc),另一个支路的阻抗为z2=-jxc,为电力系统正常电抗均值xav。
小阻抗支路串联补偿的方法,包括以下步骤:
B1、读入电力系统所有线路和变压器支路数据,设置小电阻阈值rmin和小电抗阈值xmin。
B2、计算电力系统正常电抗均值xav。
B3、设置支路计数初值m=1。
B4、设置新增加支路和节点计数初值p=0。
B5、取支路m的首末节点号i和j、电阻rij、电抗xij、变比k。
B6、判断是否满足rij≤rmin且xij≤xmin的条件,如果不满足转至步骤B12。
B7、令xc=xav。
B8、令p=p+1。
B9、增加节点号为n+p的节点,节点类型设置为PQ节点,节点电源有功功率和无功功率及负荷有功功率和无功功率都设为0。
B10、设支路m的末节点号为n+p、电抗为xij+xc,其它不变。
B11、增加支路n+p,令其首末节点号分别为n+p和j、电阻为0、电抗为-xc、变比为1.0。
B12、令m=m+1。
B13、判断m是否大于支路数l,如果m不大于l,则返回到步骤B5;否则,转至步骤C。
C、初始化电压;
D、形成节点导纳矩阵;
E、设置迭代计数t=0;
F、计算功率及电压偏差,求最大不平衡量ΔWmax;
G、判断最大不平衡量绝对值|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤K;否则,执行步骤H;
H、形成雅可比矩阵;
按式(10)-(21)计算雅可比矩阵元素,式(16)-(19)中的ai和bi按式(7)计算。
I、解修正方程及修正电压实部e、虚部f;
J、令t=t+1,返回步骤F进行下一次迭代;
K、输出节点及支路数据。
与现有技术相比,本发明具有以下有益效果:
1、本发明提出的对小阻抗支路进行串联补偿方法,把小阻抗支路变成两个阻抗较大的支路,有效改善了直角坐标牛顿法潮流计算在分析含有小阻抗支路系统时的收敛性。能够计算任意阻抗的小阻抗支路,甚至阻抗为零的支路。迭代次数比现有方法少,计算速度稍快。
2、由于本发明不仅能有效解决常规直角坐标牛顿法潮流方法分析含有小阻抗支路系统的收敛性问题,同时也能对正常系统进行潮流计算,因此没有不良影响。
3、本发明是对小阻抗支路数据进行处理,此处理过程是在潮流计算之前,对数据进行预处理,不需要改变潮流计算程序。因此本发明方法特别适合于对老的潮流计算程序进行改造,这些程序可能采用比较老的编程语言编写,不易修改,甚至有的潮流计算版本没有源文件,无法修改。
附图说明
本发明共有附图4张。其中:
图1是电力系统小阻抗变压器模型示意图。
图2是现有直角坐标牛顿法潮流计算的流程图。
图3是本发明直角坐标牛顿法潮流计算的流程图。
图4是本发明串联补偿的流程图。
具体实施方式
下面结合附图对本发明进行进一步地说明,按照图3-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。计算时各方法均采用了稀疏矩阵技术和节点优化编号技术。潮流计算的收敛精度ε为0.00001。作为对比,采用现有专利ZL201611094297.8所提出方法对该实际大型电网进行了潮流计算。
计算结果见表1,其中计算时间是在同一计算机环境的计算结果。
表1不同潮流方法的迭代结果
方法 | 现有专利方法 | 本发明方法 |
迭代次数 | 6次收敛 | 5次收敛 |
计算时间(ms) | 41.368 | 40.765 |
由表1可见,对于修改后的445节点实际电力系统算例,现有发明专利迭代6次收敛,本发明方法的迭代次数少1次。本发明经过串联补偿,增加了41个节点和41条支路,内存占用量增加不到10%。本发明的迭代次数比现有专利方法少1次,但由于本发明增加10%左右的节点,每次迭代的时间要增加10%左右,所以总时间比现有专利方法减少得不多。
为了验证本发明处理零阻抗支路的能力,把l118-125的阻抗改为0+j0,本发明仍然能正常计算,迭代次数为5次;现有专利方法因为零除计算失败。
本发明可以采用任何一种编程语言和编程环境实现,如C语言、C++、FORTRAN、Delphi等。开发环境可以采用Visual C++、Borland C++Builder、Visual FORTRAN等。
本发明不局限于本实施例,任何在本发明披露的技术范围内的等同构思或者改变,均列为本发明的保护范围。
Claims (1)
1.用于含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法,其特征在于:采用串联补偿方式对小阻抗支路进行补偿,把每个小阻抗支路变成两个正常阻抗的支路,然后进行直角坐标牛顿法潮流计算;具体包括以下步骤:
A、原始数据输入;
B、对小阻抗支路进行串联补偿;
用电力系统正常电抗均值对小阻抗支路进行补偿,得到两个阻抗较大的支路,其中一个支路的阻抗为z1=r+j(x+xc),另一个支路的阻抗为z2=-jxc,这里xc为电力系统正常电抗均值xav;
小阻抗支路串联补偿的方法,包括以下步骤:
B1、读入电力系统所有线路和变压器支路数据,设置小阻抗阈值电阻rmin和阈值电抗xmin;
B2、计算电力系统正常电抗均值xav;
B3、设置支路计数初值m=1;
B4、设置新增加支路和节点计数初值p=0;
B5、取支路m的首末节点号i和j、电阻r、电抗x、变比k;
B6、判断是否满足r≤rmin且x≤xmin的条件,如果不满足转至步骤B12;
B7、令xc=xav;
B8、令p=p+1;
B9、增加节点号为n+p的节点,节点类型设置为PQ节点,节点电源有功功率和无功功率及负荷有功功率和无功功率都设为0;
B10、设支路m的末节点号为n+p、电抗为x+xc,其它不变;
B11、增加支路n+p,令其首末节点号分别为n+p和j、电阻为0、电抗为-xc、变比为1.0;
B12、令m=m+1;
B13、判断m是否大于支路数l,如果m不大于l,则返回到步骤B5;否则,转至步骤C;
C、初始化电压;
根据电力系统节点的特点,潮流计算把电力系统节点分成3类:节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点;节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点;节点电压幅值和电压相角已知,节点有功功率和无功功率未知的节点称为平衡节点;
电压初始化采用平启动,即PV节点和平衡节点的电压实部取给定值,PQ节点的电压实部取1.0;所有电压的虚部都取0.0;单位采用标幺值;
D、形成节点导纳矩阵;
E、设置迭代计数t=0;
F、计算功率及电压偏差,求最大不平衡量ΔWmax;
PQ节点的功率偏差计算公式为:
式中,Pis、Qis分别为节点i给定的注入有功功率和无功功率,Pis为电源有功功率与负荷有功功率之差,Qis为电源无功功率与负荷无功功率之差;ei和fi分别为节点i的电压相量的实部和虚部;ai、bi分别为节点i的计算注入电流相量的实部和虚部,为
式中,n为电力系统的节点数;
PV节点的有功功率及电压偏差计算公式为:
式中,Vis为节点i给定的电压幅值;
平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差;
求各节点功率或电压偏差中绝对值最大的值,称为最大不平衡量ΔWmax;
G、判断最大不平衡量绝对值|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤K;否则,执行步骤H;
H、形成雅可比矩阵J
当i≠j时,雅可比矩阵J的元素计算公式如下:
当i=j时,雅可比矩阵J的元素计算公式如下:
PQ节点按式(10)-(13)计算雅可比矩阵元素;PV节点按式(10)、(11)、(14)、(15)计算雅可比矩阵元素;平衡节点不计算雅可比矩阵元素;
I、解修正方程及修正电压实部e、虚部f
潮流计算的基本方程(1)和(3)是非线性方程组,采用逐次线性化方法迭代求解;线性化得到的方程称为修正方程,用来求电压实部和虚部的修正量;
修正方程为:
式中,J为雅可比矩阵;ΔP和ΔQ分别为有功功率和无功功率偏差列向量;ΔV2为电压幅值偏差列向量;Δe和Δf分别为电压相量的实部和虚部修正量列向量;为有功功率偏差函数列向量对电压相量实部列向量转置的偏导矩阵,上标T为转置符号;
电压修正公式为:
式中,上标t表示第t次迭代;
J、令t=t+1,返回步骤F进行下一次迭代;
K、输出节点及支路数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611131677.4A CN106532712B (zh) | 2016-12-09 | 2016-12-09 | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611131677.4A CN106532712B (zh) | 2016-12-09 | 2016-12-09 | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106532712A true CN106532712A (zh) | 2017-03-22 |
CN106532712B CN106532712B (zh) | 2019-02-22 |
Family
ID=58342887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611131677.4A Expired - Fee Related CN106532712B (zh) | 2016-12-09 | 2016-12-09 | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106532712B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109586268A (zh) * | 2018-11-05 | 2019-04-05 | 南昌大学 | 一种基于支路电阻偏差的牛顿-拉夫逊法直流电网潮流计算方法 |
CN110265983A (zh) * | 2019-07-02 | 2019-09-20 | 广东电网有限责任公司 | 输电线路的保护方法、装置、设备及存储介质 |
CN110649624A (zh) * | 2019-11-09 | 2020-01-03 | 国网吉林省电力有限公司 | 一种电力系统潮流并行计算方法 |
CN111413890A (zh) * | 2020-04-01 | 2020-07-14 | 国网四川省电力公司技能培训中心 | 基于数字物理混合仿真的智能变电站一次平台仿真方法 |
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-12-09 CN CN201611131677.4A patent/CN106532712B/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 (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109586268A (zh) * | 2018-11-05 | 2019-04-05 | 南昌大学 | 一种基于支路电阻偏差的牛顿-拉夫逊法直流电网潮流计算方法 |
CN109586268B (zh) * | 2018-11-05 | 2022-02-11 | 南昌大学 | 一种基于支路电阻偏差的牛顿-拉夫逊法直流电网潮流计算方法 |
CN110265983A (zh) * | 2019-07-02 | 2019-09-20 | 广东电网有限责任公司 | 输电线路的保护方法、装置、设备及存储介质 |
CN110265983B (zh) * | 2019-07-02 | 2021-05-28 | 广东电网有限责任公司 | 输电线路的保护方法、装置、设备及存储介质 |
CN110649624A (zh) * | 2019-11-09 | 2020-01-03 | 国网吉林省电力有限公司 | 一种电力系统潮流并行计算方法 |
CN111413890A (zh) * | 2020-04-01 | 2020-07-14 | 国网四川省电力公司技能培训中心 | 基于数字物理混合仿真的智能变电站一次平台仿真方法 |
CN111413890B (zh) * | 2020-04-01 | 2023-07-11 | 国网四川省电力公司技能培训中心 | 基于数字物理混合仿真的智能变电站一次平台仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106532712B (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN104037764B (zh) | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 | |
CN104037763B (zh) | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 | |
CN106709243B (zh) | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 | |
CN106532712A (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
CN104899396B (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
CN109617080B (zh) | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 | |
CN107196306B (zh) | 基于Matlab稀疏矩阵的快速分解法潮流计算方法 | |
CN104600697A (zh) | 一种计及温度影响的拟直流最优潮流方法 | |
CN104636829A (zh) | 一种提高电力系统温度最优潮流计算效率的解耦算法 | |
CN106856327B (zh) | 一种含小阻抗支路电网串联补偿快速分解法潮流计算方法 | |
CN104022507B (zh) | 一种直角坐标牛顿法潮流计算方法 | |
CN101621200B (zh) | 一种直角坐标牛顿法潮流计算方法 | |
CN106229988A (zh) | 一种基于Matlab的极坐标牛顿法潮流计算方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN109494748B (zh) | 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法 | |
CN102780220A (zh) | 含有pv恒定型分布式电源的配电网潮流计算方法 | |
CN106712029B (zh) | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 | |
CN106529089B (zh) | 用于含小阻抗支路电网的补偿法快速分解法潮流计算方法 | |
CN106786605B (zh) | 一种含小阻抗电网串联补偿直角坐标牛顿法潮流计算方法 | |
CN104300572A (zh) | 一种配电网中分布式电源的配置方法 | |
CN107658880B (zh) | 基于关联矩阵运算的快速分解法系数矩阵计算方法 | |
CN104269872A (zh) | 一种变压器三相节点导纳矩阵奇异的处理方法 | |
CN106786604B (zh) | 一种含小阻抗电网的串联补偿极坐标牛顿法潮流计算方法 | |
CN103972886B (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: 20191209 |