CN106709243B - 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 - Google Patents
含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 Download PDFInfo
- Publication number
- CN106709243B CN106709243B CN201611129772.0A CN201611129772A CN106709243B CN 106709243 B CN106709243 B CN 106709243B CN 201611129772 A CN201611129772 A CN 201611129772A CN 106709243 B CN106709243 B CN 106709243B
- Authority
- CN
- China
- Prior art keywords
- node
- voltage
- small impedance
- power
- branch
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种小阻抗支路PQ端点变雅可比矩阵的牛顿法潮流计算方法,包括以下步骤:根据支路电阻和电抗的大小确定两端节点所连支路类型T;小阻抗支路的PQ端点与其他节点采用不同的方法计算雅可比矩阵;解修正方程。本发明通过在迭代过程中小阻抗支路的PQ端点采用与其他节点不同的雅可比矩阵计算方法,解决了直角坐标牛顿法潮流计算在分析含有小阻抗支路电力系统时的收敛性问题。采用常规直角坐标牛顿法潮流计算不收敛时,本发明能够可靠收敛,且比现有专利技术迭代次数少。由于本发明不仅能有效解决常规直角坐标牛顿法潮流计算分析含有小阻抗支路电力系统的收敛性问题,同时也能对正常电力系统进行潮流计算,没有不良影响。
Description
技术领域
本发明涉及含小阻抗支路电力系统的潮流计算方法,特别是一种电力系统的极坐标牛顿法潮流计算方法。
背景技术
电力系统潮流计算是研究电力系统稳态运行的一项基本计算,它根据电力系统给定的运行条件和网络结构确定整个电力系统的运行状态。潮流计算也是电力系统其他分析的基础,如安全分析、暂态稳定分析等都要用到潮流计算。由于具有收敛可靠、计算速度较快及内存需求适中的优点,牛顿法成为当前潮流计算的主流方法。牛顿法分为极坐标和直角坐标两种形式,两种形式的牛顿法潮流计算都在电力系统中得到了广泛的应用。
在极坐标牛顿法潮流计算中,节点i的电压采用极坐标表示为:
对正常电力网络,极坐标牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,极坐标牛顿法潮流计算就可能发散。电力系统小阻抗支路可分为小阻抗线路和小阻抗变压器支路,在数学模型上线路可以看作变比为1:1的变压器,因此下面分析时仅以小阻抗变压器支路为例分析。小阻抗变压器模型见图1,变压器的非标准变比k位于节点i侧,阻抗位于标准变比侧。变压器阻抗zij=rij+jxij很小,导纳为
式中,yij、gij、bij分别为节点i和节点j之间小阻抗支路的导纳、电导和电纳;rij、xij分别为节点i和节点j之间小阻抗支路的电阻和电抗。
由于小阻抗支路li-j的阻抗很小,支路的电压降也很小,因此变压器两端节点的电压应满足:
式中,Ui、θi分别为节点i的电压幅值和相角;Uj、θj分别为节点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、计算功率偏差
节点功率公式为:
式中,Pi、Qi分别为节点i的节点有功功率和无功功率;Ui、Uk分别为节点i和节点k的节点电压幅值;θik=θi-θk,θi和θk分别为节点i和节点k的节点电压相角;Gik、Bik分别为节点导纳矩阵元素Yik的实部和虚部;n为电力系统的节点数。
节点功率偏差计算公式为:
式中,ΔPi、ΔQi分别为节点i的节点有功功率偏差和无功功率偏差;Pis、Qis分别为节点i给定的节点注入有功功率和注入无功功率;m为PQ节点数。
PV节点的无功功率不是给定值,不需要计算无功功率偏差;平衡节点不参与迭代计算,不需要计算有功功率偏差和无功功率偏差。
求各节点功率偏差中绝对值最大的值,称为最大不平衡量,如果最大不平衡量的绝对值小于给定的收敛精度,转步骤F,否则执行步骤D。
D、形成雅可比矩阵J
雅可比矩阵J的元素(i≠j时)计算公式如下:
雅可比矩阵J的元素(i=j时)计算公式如下:
E、解修正方程及修正电压幅值U、相角θ
潮流计算的基本方程(7)是非线性方程组,通常采用逐次线性化方法迭代求解。线性化得到的方程称为修正方程,用来求电压幅值和相角的修正量。
修正方程为:
式中,J为雅可比矩阵;H、N、M、L为雅可比矩阵的分块子矩阵;ΔP和ΔQ分别为有功功率和无功功率偏差列向量;Δθ为电压相角修正量列向量;ΔU/U为电压幅值修正量与电压幅值之比的列向量。
电压修正公式为:
式中,上标t表示第t次迭代。
F、输出节点及支路数据。
对正常电力网络,牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,牛顿法潮流计算就可能发散。而电力系统中小阻抗支路普遍存在,收敛性是电力系统潮流计算这类非线性问题的最重要指标,计算不收敛就无法得到方程的解。因此改善极坐标牛顿法潮流计算针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。
发明者在早期研究中提出了一种变雅可比牛顿法潮流计算方法,该方法采用给定值Pis和Qis代替计算值Pi和Qi计算雅可比矩阵的Nii和Lii元素,而雅可比矩阵的Hii和Mii元素仍由Pi和Qi计算。变雅可比牛顿法潮流计算方法有效解决了极坐标牛顿法潮流计算分析包含小阻抗支路的电力系统时的发散问题。
现在仍然在使用的一些极坐标牛顿法潮流计算程序不具备处理小阻抗支路的能力,但对这些老的潮流计算程序进行改造往往比较麻烦。有些程序可能采用比较老的编程语言编写,不易修改,有的潮流计算版本甚至没有源文件,无法修改。
发明内容
为解决现有技术存在的上述问题,本发明要提出一种含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法,以改善极坐标牛顿法潮流计算方法分析含有小阻抗支路电力系统的潮流计算的收敛性。
为了实现上述目的,本发明的技术方案如下:含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法,在进行潮流计算之前,先对小阻抗支路进行处理,把此小阻抗支路变成两个支路串联的形式,其中一个支路的阻抗为z1=rij+j(xij+xc),另一个支路的阻抗为z2=-jxc,这里xc为电力系统中电抗绝对值大于小电抗阈值的支路的电抗绝对值的平均值,称为电力系统正常电抗均值xav。经过串联补偿把一个小阻抗变成两个阻抗较大的支路,可以提高潮流计算的收敛性。方案包括以下步骤:
A、输入原始数据;
B、对小阻抗支路进行串联补偿;
用电力系统正常电抗均值对小阻抗支路进行补偿,得到两个阻抗较大的支路,其中一个支路的阻抗为z1=rij+j(xij+xc),另一个支路的阻抗为z2=-jxc,这里xc为电力系统正常电抗均值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、形成雅可比矩阵;
按式(8)-(15)计算雅可比矩阵元素。
I、解修正方程及修正电压幅值U、相角θ;
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。作为对比,采用发明者早期提出的方法对该实际大型电网进行了潮流计算。
计算结果见表1,其中计算时间是在同一计算机环境的计算结果。
表1不同潮流方法的迭代结果
方法 | 对比方法 | 本发明方法 |
迭代次数 | 5次收敛 | 5次收敛 |
计算时间(ms) | 33.153 | 38.987 |
由表1可见,对于修改后的445节点实际电力系统算例,对比方法和本发明方法的迭代次数都是5次。本发明方法经过串联补偿,增加了41个节点和41条支路,内存占用量增加不到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和q、电阻r、电抗x、变比α;
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和q、电阻为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;
节点功率公式为:
式中,Pi、Qi分别为节点i的节点有功功率和无功功率;Ui、Uk分别为节点i和节点k的节点电压幅值;θik=θi-θk,θi和θk分别为节点i和节点k的节点电压相角;Gik、Bik分别为节点导纳矩阵元素Yik的实部和虚部;n为电力系统的节点数;
节点功率偏差计算公式为:
式中,ΔPi、ΔQi分别为节点i的节点有功功率偏差和无功功率偏差;Pis、Qis分别为节点i给定的节点注入有功功率和注入无功功率;m为PQ节点数;
PV节点的无功功率不是给定值,不需要计算无功功率偏差;平衡节点不参与迭代计算,不需要计算有功功率偏差和无功功率偏差;
求各节点功率偏差中绝对值最大的值,称为最大不平衡量ΔWmax;
G、判断最大不平衡量绝对值|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤K;否则,执行步骤H;
H、形成雅可比矩阵J
当i≠q时,雅可比矩阵J的元素计算公式如下:
当i=q时,雅可比矩阵J的元素计算公式如下:
I、解修正方程及修正电压幅值U、相角θ
潮流计算的基本方程(2)是非线性方程组,采用逐次线性化方法迭代求解;线性化得到的方程称为修正方程,用来求电压幅值和相角的修正量;
修正方程为:
式中,J为雅可比矩阵;H、N、M、L为雅可比矩阵的分块子矩阵;ΔP和ΔQ分别为有功功率和无功功率偏差列向量;Δθ为电压相角修正量列向量;ΔU/U为电压幅值修正量与电压幅值之比的列向量;
电压修正公式为:
式中,上标t表示第t次迭代;
J、令t=t+1,返回步骤F进行下一次迭代;
K、输出节点及支路数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611129772.0A CN106709243B (zh) | 2016-12-09 | 2016-12-09 | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611129772.0A CN106709243B (zh) | 2016-12-09 | 2016-12-09 | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106709243A CN106709243A (zh) | 2017-05-24 |
CN106709243B true CN106709243B (zh) | 2019-04-05 |
Family
ID=58936589
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611129772.0A Expired - Fee Related CN106709243B (zh) | 2016-12-09 | 2016-12-09 | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106709243B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108964047A (zh) * | 2018-08-07 | 2018-12-07 | 河海大学 | 一种配电网状态估计方法 |
CN110046450B (zh) * | 2019-04-24 | 2023-04-07 | 深圳供电局有限公司 | 一种适合含超导电缆电网的牛顿法潮流计算的初值选取方法 |
CN111046336A (zh) * | 2019-12-27 | 2020-04-21 | 大连海事大学 | 一种雅可比矩阵变化的极坐标牛顿法潮流计算方法 |
CN111030122B (zh) * | 2019-12-27 | 2023-12-22 | 大连海事大学 | 首次迭代雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
CN111049146B (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 | 大连海事大学 | 一种直角坐标牛顿法潮流计算方法 |
CN104037764A (zh) * | 2014-07-03 | 2014-09-10 | 大连海事大学 | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
CN104037763A (zh) * | 2014-07-02 | 2014-09-10 | 大连海事大学 | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 |
-
2016
- 2016-12-09 CN CN201611129772.0A patent/CN106709243B/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 | 大连海事大学 | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 |
Non-Patent Citations (4)
Title |
---|
"小阻抗支路对潮流计算收敛性的影响和处理";韩平等;《电力系统保护与控制》;20090916;第37卷(第18期);第17-19、31页 |
"求解含小阻抗支路配电网潮流的牛顿法";初壮等;《电力系统及其自动化学报》;20160930;第28卷(第9期);第36-41页 |
"求解含有小阻抗支路系统潮流的一种新方法";姚玉斌等;《哈尔滨工业大学学报》;20010830;第33卷(第4期);第525-529页 |
"牛顿类潮流计算方法的收敛性分析";孙秋野等;《中国电机工程学报》;20140505;第34卷(第13期);第2196-2200页 |
Also Published As
Publication number | Publication date |
---|---|
CN106709243A (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106709243B (zh) | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 | |
CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN104037764B (zh) | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 | |
CN106356859B (zh) | 一种基于Matlab的直角坐标牛顿法潮流计算方法 | |
JP6025575B2 (ja) | 三相電力潮流解析のためのシステムおよび方法 | |
CN106532712B (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
CN104037763B (zh) | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 | |
CN109617080B (zh) | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 | |
CN104899396B (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
CN106856327B (zh) | 一种含小阻抗支路电网串联补偿快速分解法潮流计算方法 | |
CN107196306B (zh) | 基于Matlab稀疏矩阵的快速分解法潮流计算方法 | |
CN106229988B (zh) | 一种基于Matlab的极坐标牛顿法潮流计算方法 | |
CN104600697A (zh) | 一种计及温度影响的拟直流最优潮流方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN104022507A (zh) | 一种直角坐标牛顿法潮流计算方法 | |
CN109494748A (zh) | 基于节点类型和修正的雅可比矩阵的牛顿法潮流计算方法 | |
CN106529089B (zh) | 用于含小阻抗支路电网的补偿法快速分解法潮流计算方法 | |
CN107846022B (zh) | 基于ilutp预处理并行迭代法的大规模配电网潮流分析方法 | |
CN106712029B (zh) | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 | |
CN106786605B (zh) | 一种含小阻抗电网串联补偿直角坐标牛顿法潮流计算方法 | |
CN108183487B (zh) | 一种基于线性数学模型的配电网潮流快速分析方法 | |
CN106786604B (zh) | 一种含小阻抗电网的串联补偿极坐标牛顿法潮流计算方法 | |
CN107658880B (zh) | 基于关联矩阵运算的快速分解法系数矩阵计算方法 | |
CN107181260A (zh) | 基于Matlab稀疏矩阵直角坐标牛顿法潮流计算方法 | |
CN111049146B (zh) | 首次迭代雅可比矩阵改变的极坐标牛顿法潮流计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
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: 20190405 Termination date: 20191209 |