CN107944682A - 基于Matlab矩阵运算的潮流计算导纳矩阵计算方法 - Google Patents
基于Matlab矩阵运算的潮流计算导纳矩阵计算方法 Download PDFInfo
- Publication number
- CN107944682A CN107944682A CN201711136624.6A CN201711136624A CN107944682A CN 107944682 A CN107944682 A CN 107944682A CN 201711136624 A CN201711136624 A CN 201711136624A CN 107944682 A CN107944682 A CN 107944682A
- Authority
- CN
- China
- Prior art keywords
- array
- admittance
- matrix
- branch
- node
- 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
- 239000011159 matrix material Substances 0.000 title claims abstract description 116
- 238000004364 calculation method Methods 0.000 title claims abstract description 49
- 238000000205 computational method Methods 0.000 title claims abstract description 8
- 238000003491 array Methods 0.000 claims abstract description 3
- 230000005611 electricity Effects 0.000 claims description 22
- 230000005540 biological transmission Effects 0.000 claims description 20
- 230000017105 transposition Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 10
- 238000000034 method Methods 0.000 description 25
- 230000006870 function Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Pure & Applied Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Entrepreneurship & Innovation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Educational Administration (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Software Systems (AREA)
- Development Economics (AREA)
- Computing Systems (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于Matlab矩阵运算的潮流计算导纳矩阵计算方法,完全采用矩阵运算计算导纳矩阵,包括以下步骤:读支路原始数据数组;形成有向支路的首末节点号数组I0和J0;计算支路π形等值电路串联支路导纳数组YB;形成部分互导纳矩阵;计算有向支路首节点和末节点的自导纳增量数组;形成自导纳对角矩阵Y0D;形成节点导纳矩阵Y;节点导纳矩阵Y的自导纳追加无功补偿导纳值。本发明在Matlab平台实现,便于科研人员使用Matlab提供的各种工具和函数对计算结果进行测试和分析。本发明提出的导纳矩阵计算完全采用矩阵运算和复数运算,减少了程序代码,简化了编程,使得程序更清晰;使用矩阵运算也大大提高了计算速度。
Description
技术领域
本发明涉及一种电力系统潮流计算的导纳矩阵计算方法,特别是一种适合研究目的使用的基于Matlab的潮流计算的导纳矩阵计算方法。
背景技术
电力系统潮流计算是研究电力系统稳态运行的一项基本计算,它根据给定的运行条件和网络结构确定整个网络的运行状态。潮流计算也是电力系统其他分析的基础,如安全分析、暂态稳定分析等都要用到潮流计算。潮流计算是电力系统分析的基本分析工具,科研人员经常以潮流计算为基础进行进一步地研究。实用的商业软件采用C语言等高级编程语言编写且使用稀疏矩阵技术和节点优化编号等高级技术。这些技术虽然能大幅度提高潮流计算的速度、降低内存占用量,但编程非常麻烦且难以修改和维护,不易增加新的功能,因而不适合科研人员用于研究目的使用。
Matlab软件以矩阵为最基本的数据单位,可以方便地处理各种矩阵和向量运算,也可以很方便自然地处理复数类型,其指令表达式与数学中常用的形式很接近,还有大量常见实用的函数,给编程带来很大便利。Matlab软件简单易用、代码短小易操作,易于编程和调试,计算功能强大,同时还具有非常强大的可视化图形处理和交互式功能,为科学研究以及工程应用提供了一种高效的编程工具,目前已经成为许多科学领域的基本工具和首选平台,在各种科学和工程计算领域得到了广泛的应用。为了适应越来越多的科研人员需要在Matlab平台上以潮流计算为基础进行进一步地研究的需求,迫切需要一种基于Matlab软件的易于编程、修改和调试的潮流计算方法。
目前常用的牛顿法和快速分解法潮流计算方法都是以节点电压法为基础的,都需要求取节点导纳矩阵。
节点导纳矩阵为:
式中,Yik为节点导纳矩阵元素,当下标i≠k时,为节点i和节点k之间的互导纳,当下标i=k时,为节点i的自导纳;n为节点数。
可以采用追加支路法计算导纳矩阵元素,即依次扫描所有支路,每扫描一条支路在原有导纳矩阵元素基础上增加一个导纳增量。由于输电线路和变压器支路都属于支路,通常把输电线路和变压器支路都作为支路数据统一输入。作为区分,变压器支路的非标准变比侧的节点号加个负号。
输电线路采用如图2所示的π形等值电路,设第m条支路的串联支路导纳为
ym=1/zm=1/(rm+jxm) (2)
式中,rm、xm和zm分别为输电线路等值电路的电阻、电抗和阻抗。
追加第m条输电线路时,导纳矩阵元素的计算公式如下:
式中,下标i、j分别表示支路的首节点号im和末节点号jm去掉负号后的节点号,bm为输电线路等值电路的对地电纳,符号“←”表示右端计算结果赋值给左端变量。
变压器支路采用如图3所示的理想变压器串联一等值阻抗的等值电路表示,根据变比及等值阻抗处于位置不同,分为4种情况,图3(a)和图3(b)为等值阻抗位于标准变比侧(即1侧),图3(c)和图3(d)为等值阻抗位于非标准变比侧(即km侧)。为了降低程序设计的复杂程度,通常把图3(c)和图3(d)的变比转换成(1/km):1形式,从而把图3(c)和图3(d)所示的等值电路变成图3(a)和图3(b)所示的等值电路。图4(a)和图4(b)分别为图3(a)和图3(b)所示电路的π形等值电路,其中ym=1/(rm+jxm)。
如图3(a)所示,当变比km位于首节点i侧时,根据其π形等值电路图4(a)得到第m条变压器支路的导纳矩阵元素的计算公式如下:
式中,km为变压器支路的变比。
如图3(b)所示,当变比km位于末节点j侧时,根据其π形等值电路图4(b)得到第m条变压器支路的导纳矩阵元素的计算公式如下:
如图1-5所示,现有潮流计算方法,主要包括以下步骤:
A、原始数据输入和电压初始化;
原始数据包括线路和变压器支路数据、节点注入有功功率和无功功率、节点电压幅值、节点无功补偿数据,以及收敛精度、最大迭代次数;
输电线路和变压器支路都作为支路数据统一输入,作为区分,变压器支路的非标准变比侧的节点号加个负号;
根据电力系统节点的特点,潮流计算把电力系统节点分成3类:节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点;节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点;节点电压幅值和电压相角已知,节点有功功率和无功功率未知的节点称为平衡节点。
电压初始化采用平启动,即PV节点和平衡节点的电压幅值取给定值,PQ节点的电压幅值取1.0;所有电压的相角都取0.0。这里相角单位为弧度,其他量单位采用标幺值。
B、形成节点导纳矩阵;
形成节点导纳矩阵的步骤如下:
B1、设置支路计数m=1;
B2、取支路m的首节点号im、末节点号jm,并令i=|im|、j=|jm|;
B3、取支路m的电阻rm、电抗xm,并令ym=1/(rm+jxm);
B4、判断支路m的首节点号im、末节点号jm是否都大于0,如果不满足转至步骤B7;
B5、取支路m的对地电纳bm;
B6、按式(2)和式(3)计算输电线路两端节点对应的节点导纳矩阵元素;
B7、判断支路m的首节点号im是否小于0,如果不满足转至步骤B10;
B8、取变压器支路m的变比km;
B9、按式(2)和式(4)计算变压器支路两端节点对应的节点导纳矩阵元素;
B10、判断支路m的末节点号jm是否小于0,如果不满足转至步骤B13;
B11、取变压器支路m的变比km;
B12、按式(2)和式(5)计算变压器支路两端节点对应的节点导纳矩阵元素;
B13、令m=m+1。
B14、判断m是否大于支路数l,如果m不大于l,则返回到步骤B2;否则,转步骤C。
C、潮流计算迭代主程序;
根据采用的潮流计算的方法不同,可以采用极坐标牛顿法、直角坐标牛顿法、快速分解法进行潮流计算。
D、计算平衡节点的有功功率和无功功率及PV节点的无功功率;
平衡节点的有功功率和无功功率及PV节点的无功功率未知,需要计算求出。
E、计算各支路有功功率和无功功率;
F、输出计算结果,结束。
直接采用上述原理实现的潮流计算软件计算速度较慢,商业使用的潮流计算软件采用稀疏矩阵技术和节点优化编号技术,比较复杂,不适合科研人员以此为基础进一步进行科学研究。因此,中国专利CN201710557623.2、CN201710557642.5和CN201710557622.8分别提出基于Matlab的极坐标牛顿法潮流计算方法、直角坐标牛顿法潮流计算方法和快速分解法潮流计算方法,可以充分利用Matlab特有的擅长矩阵运算和复数运算的特点,并采用Matlab的稀疏矩阵技术和方程求解算法,设计出了简洁又有较快计算速度的潮流计算方法,为以潮流计算为基础进行进一步研究的科研人员提供了3种易于修改和维护的潮流计算方法;中国专利CN201710942153.1在上述专利的基础上提出了基于关联矩阵和矩阵运算的导纳矩阵计算方法,进一步提高了潮流计算的速度。但中国专利CN201710942153.1方法计算导纳矩阵时所使用的关联矩阵未实现矩阵运算,导纳矩阵计算速度相对较慢,仍有待进一步提高计算速度。
发明内容
为解决现有技术存在的上述问题,本发明要提出基于Matlab矩阵运算的导纳矩阵计算方法,充分利用Matlab特有的擅长矩阵运算的特点,实现提高潮流计算的计算速度的目的。
为了实现上述目的,本发明的技术方案如下:基于Matlab矩阵运算的潮流计算导纳矩阵计算方法,完全采用矩阵运算计算导纳矩阵。下面推导完全采用矩阵运算计算导纳矩阵的公式。
变压器支路的π形等值电路不是对称的,追加支路法计算导纳矩阵元素时,支路的首末节点的自导纳增量不同。为了计算方便,采用有向支路概念,先规定电力网络中支路的方向,输电线路的方向为首节点指向末节点,变压器支路的方向为变压器非标准变比km侧的节点指向标准变比1侧的节点,图3(a)所示的变压器支路的方向为首节点指向末节点,图3(b)所示的变压器支路的方向为末节点指向首节点。
根据定义,有向支路的首末节点号数组为:
I0=(I<0)·*abs(I)+(I>0)·*abs(J) (6)
J0=(I<0)·*J+(I>0)·*I (7)
式中,I、J分别为支路原来的首、末节点号数组,其中变压器支路的非标准变比km侧的节点号为负数,I<0表示I数组元素分别与0比较,成立为1,不成立为0,结果仍为一个数组,I0、J0分别为有向支路首、末节点号数组,abs为Matlab求数组元素绝对值函数,“.*”表示两数组对应元素相乘。
第m条输电线路等值电路的串联支路导纳为
ym=1/zm=1/(rm+jxm) (8)
式中,rm、xm和zm分别为输电线路等值电路的电阻、电抗和阻抗。
由图4(a)和图4(b)所示的变压器π形等值电路可知,第m条变压器支路π形等值电路的串联支路导纳为
y′m=1/zm/km=1/(rm+jxm)/km (9)
式中,rm、xm、zm和km分别为变压器支路等值电路的电阻、电抗、阻抗和变比。
输电线路可以看作变比km=1的变压器支路,则输电线路和变压器支路π形等值电路的串联支路导纳可以统一写成式(9)的形式。
用矩阵表示支路π形等值电路的串联支路导纳数组YB为
YB=1·/(R+jX)·/K (10)
式中,R为支路电阻数组,X为支路电抗数组,K为支路变比数组,输电线路的变比为1,“·/”表示两数组对应元素相除。
电力网络中的导纳矩阵的部分互导纳矩阵为
Y1=sparse(I0,J0,-YB,n,n) (11)
式中,sparse为Matlab形成稀疏矩阵函数,其参数分别为矩阵的行号数组、列号数组、元素值数组、行数、列数,I0、J0分别为有向支路首节点号数组、末节点号数组,n为节点数。
式(11)形成的部分互导纳矩阵是由一半互导纳元素构成的矩阵,另外一半互导纳元素可以通过矩阵转置得到。利用sparse函数形成稀疏矩阵时,如果两个元素的行列号相同,这两个元素会合并。
电力网络中导纳矩阵的自导纳元素由输电线路和变压器的π形等值电路串联支路和对地支路计算。
由图2可知,输电线路等值电路两侧节点的对地导纳相同,两侧节点自导纳增量为
ys=ye=y′m+jbm/2 (12)
式中,ys为有向支路首节点对地电纳,ye为有向支路末节点对地电纳。
由图4(a)和4(b)可知,变压器支路π形等值电路非标准变比km侧节点的自导纳增量为
变压器支路π形等值电路标准变比1侧节点的自导纳增量为
ye=ym=y′mkm (14)
综合式(12)~式(14),并写成矩阵形式得到有向支路首节点的自导纳增量数组YS、末节点的自导纳增量数组YE为
YS=YB·/K+jB/2 (15)
YE=YB·*K+jB/2 (16)
式中,B为对地电纳数组,变压器支路B值为0。
利用Matlab的形成稀疏矩阵函数sparse由YS和YE生成导纳矩阵的自导纳元素构成的对角矩阵为:
Y0D=sparse(I0,I0,YS,n,n)+sparse(J0,J0,YE,n,n) (17)
式中,Y0D为YS和YE生成的导纳矩阵的对角矩阵。
使用式(17)生成矩阵时,连在相同节点的自导纳增量会合并,得到该节点总的自导纳。
电力网络支路的π形等值电路串联支路的形成节点互导纳矩阵加上各节点自导纳形成的对角矩阵得到由支路数据计算出的节点导纳矩阵,计算公式如下:
Y=Y1+Y1T+Y0D (18)
式中,上标T表示矩阵的转置。
本发明完全采用矩阵运算形成导纳矩阵,包括以下步骤:
B1、读支路首节点号数组I、末节点号数组J、电阻数组R、电抗数组X、对地电纳数组B、变压器变比数组K;
所述的首节点号数组I、末节点号数组J、电阻数组R、电抗数组X、对地电纳数组B、变压器变比数组K分别按顺序存放所有支路的首节点号im、末节点号jm、电阻rm、电抗xm、对地电纳bm、变压器变比km,变压器支路对地电纳bm为0,变压器非标准变比km侧的节点号加个负号,输电线路变比km设为1,其中下标m为支路序号;
B2、根据有向支路方向的定义由式(6)和式(7)形成有向支路的首末节点号数组I0和J0;
B3、用式(10)计算支路π形等值电路串联支路导纳数组YB;
B4、由串联支路形成部分互导纳矩阵;
用式(11)计算由支路π形等值电路串联支路形成的部分互导纳矩阵Y1;
B5、计算有向支路首节点和末节点的自导纳增量数组;
分别用式(15)和式(16)计算有向支路首节点和末节点的自导纳增量数组YS和YE;
B6、用式(17)由自导纳增量数组YS和YE形成自导纳对角矩阵Y0D;
B7、形成节点导纳矩阵Y;
由支路的π形等值电路串联支路形成的部分互导纳矩阵Y1和自导纳对角矩阵Y0D按式(18)计算节点导纳矩阵Y。
B8、节点导纳矩阵Y的自导纳追加无功补偿导纳值。
与现有技术相比,本发明具有以下有益效果:
1、本发明提出的方法在Matlab平台实现,便于科研人员使用Matlab提供的各种工具和函数对计算结果进行测试和分析。
2、本发明提出的导纳矩阵计算完全采用矩阵运算和复数运算,减少了程序代码,简化了编程,使得程序更加清晰;使用矩阵运算也大大提高了计算速度。
附图说明
本发明共有附图6张。其中:
图1是现有潮流计算的流程图。
图2是输电线路的等值电路图。
图3是变压器支路的等值电路图。
图4是变压器支路的π形等值电路图。
图5是现有潮流计算形成导纳矩阵的流程图。
图6是本发明形成导纳矩阵的流程图。
具体实施方式
下面结合附图对本发明进行进一步地说明,按照图1和图6所示流程对一个10428节点实际系统算例进行了计算,该算例有10428个节点、10436条支路。
采用本发明方法和两项已有专利方法对10428节点实际系统算例进行了计算,潮流计算采用极坐标牛顿法,计算时相角单位为弧度,其他量采用标幺值,收敛精度为0.00001。3种潮流计算方法分别为:
方法1:中国专利CN201710557623.2方法,导纳矩阵计算采用循环结构及复数运算;
方法2:中国专利CN201710942153.1方法,导纳矩阵计算采用复数运算和矩阵运算,其中关联矩阵采用循环结构;
方法3:本发明方法,完全采用矩阵运算和复数运算计算导纳矩阵。
3种方法的潮流计算和导纳矩阵计算的计算时间见表1,潮流计算的计算时间不包括数据读入和输出的时间。
表1 3种极坐标牛顿法潮流计算计算时间比较
潮流计算方法 | 潮流计算计算时间(s) | 导纳矩阵计算时 |
方法1 | 2.6085 | 间(s)2.0626 |
方法2 | 0.9413 | 0.3987 |
方法3 | 0.5506 | 0.0091 |
从表1可见,中国专利CN201710557623.2方法计算导纳矩阵时间较长,占潮流计算计算时间的79.1%,占潮流计算大部分计算时间;中国专利CN201710942153.1方法采用关联矩阵和矩阵运算技术计算导纳矩阵明显提高了导纳矩阵的计算速度,导纳矩阵计算时间占潮流计算计算时间的42.4%;本发明完全采用矩阵运算计算导纳矩阵进一步提高了导纳矩阵的计算速度,导纳矩阵计算时间仅占潮流计算计算时间的1.65%。本发明导纳矩阵计算时间为专利CN201710942153.1方法的1/44。
本发明可以在任何版本的MATLAB编程语言实现,但建议使用较新版本的MATLAB语言。
本发明不局限于本实施例,任何在本发明披露的技术范围内的等同构思或者改变,均列为本发明的保护范围。
Claims (1)
1.基于Matlab矩阵运算的潮流计算导纳矩阵计算方法,其特征在于:包括以下步骤:
B1、读支路首节点号数组I、末节点号数组J、电阻数组R、电抗数组X、对地电纳数组B、变压器变比数组K;
所述的首节点号数组I、末节点号数组J、电阻数组R、电抗数组X、对地电纳数组B、变压器变比数组K分别按顺序存放所有支路的首节点号im、末节点号jm、电阻rm、电抗xm、对地电纳bm、变压器变比km,变压器支路对地电纳bm为0,变压器非标准变比km侧的节点号加个负号,输电线路变比km设为1,其中下标m为支路序号;
B2、形成有向支路的首末节点号数组I0和J0;
规定输电线路的方向为首节点指向末节点,变压器支路的方向为变压器非标准变比km侧的节点指向标准变比1侧的节点;
有向支路的首节点号数组I0和末节点号数组J0分别为:
I0=(I<0)·*abs(I)+(I>0)·*abs(J) (1)
J0=(I<0)·*J+(I>0)·*I (2)
式中,I、J分别为支路原来的首、末节点号数组,其中变压器支路的非标准变比km侧的节点号为负数,I<0表示I数组元素分别与0比较,成立为1,不成立为0,结果仍为一个数组,I0、J0分别为有向支路首、末节点号数组,abs为Matlab求数组元素绝对值函数,“.*”表示两数组对应元素相乘;
B3、计算支路π形等值电路串联支路导纳数组;
支路π形等值电路的串联支路导纳数组YB为:
YB=1·/(R+jX)·/K (3)
式中,R为支路电阻数组,X为支路电抗数组,K为支路变比数组,输电线路的变比为1,“·/”表示两数组对应元素相除;
B4、由串联支路形成部分互导纳矩阵;
由串联支路形成的部分互导纳矩阵为:
Y1=sparse(I0,J0,-YB,n,n) (4)
式中,sparse为Matlab形成稀疏矩阵函数,其参数分别为矩阵的行号数组、列号数组、元素值数组、行数、列数,I0、J0分别为有向支路首节点号数组、末节点号数组,n为节点数;
B5、计算有向支路首节点和末节点的自导纳增量数组;
有向支路首节点的自导纳增量数组YS、末节点的自导纳增量数组YE分别为:
YS=YB·/K+jB/2 (5)
YE=YB·*K+jB/2 (6)
式中,B为对地电纳数组,变压器支路B值为0;
B6、由自导纳增量数组YS和YE形成自导纳对角矩阵Y0D;
由YS和YE生成的自导纳对角矩阵Y0D为:
Y0D=sparse(I0,I0,YS,n,n)+sparse(J0,J0,YE,n,n) (7)
B7、形成节点导纳矩阵Y;
由支路形成节点导纳矩阵Y为:
Y=Y1+Y1T+Y0D (8)
式中,上标T表示矩阵的转置;
B8、节点导纳矩阵Y的自导纳追加无功补偿导纳值,结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711136624.6A CN107944682B (zh) | 2017-11-16 | 2017-11-16 | 基于Matlab矩阵运算的潮流计算导纳矩阵计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711136624.6A CN107944682B (zh) | 2017-11-16 | 2017-11-16 | 基于Matlab矩阵运算的潮流计算导纳矩阵计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107944682A true CN107944682A (zh) | 2018-04-20 |
CN107944682B CN107944682B (zh) | 2020-04-07 |
Family
ID=61931409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711136624.6A Expired - Fee Related CN107944682B (zh) | 2017-11-16 | 2017-11-16 | 基于Matlab矩阵运算的潮流计算导纳矩阵计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107944682B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108964063A (zh) * | 2018-08-24 | 2018-12-07 | 南京工程学院 | 一种电力系统潮流计算中非标准变比变压器的建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289905A1 (en) * | 2012-04-23 | 2013-10-31 | Institute Of Nuclear Energy Research Atomic Energy Council, Executive Yuan | Distribution power flow analysis system and method |
CN103441496A (zh) * | 2013-09-04 | 2013-12-11 | 博爱县电业公司 | 基于matlab的电力系统潮流计算方法 |
CN106602570A (zh) * | 2017-01-25 | 2017-04-26 | 大连海事大学 | 一种基于Matlab的快速分解法潮流计算方法 |
CN107194131A (zh) * | 2017-07-10 | 2017-09-22 | 大连海事大学 | 基于Matlab稀疏矩阵的极坐标牛顿法潮流计算方法 |
-
2017
- 2017-11-16 CN CN201711136624.6A patent/CN107944682B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289905A1 (en) * | 2012-04-23 | 2013-10-31 | Institute Of Nuclear Energy Research Atomic Energy Council, Executive Yuan | Distribution power flow analysis system and method |
CN103441496A (zh) * | 2013-09-04 | 2013-12-11 | 博爱县电业公司 | 基于matlab的电力系统潮流计算方法 |
CN106602570A (zh) * | 2017-01-25 | 2017-04-26 | 大连海事大学 | 一种基于Matlab的快速分解法潮流计算方法 |
CN107194131A (zh) * | 2017-07-10 | 2017-09-22 | 大连海事大学 | 基于Matlab稀疏矩阵的极坐标牛顿法潮流计算方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108964063A (zh) * | 2018-08-24 | 2018-12-07 | 南京工程学院 | 一种电力系统潮流计算中非标准变比变压器的建模方法 |
CN108964063B (zh) * | 2018-08-24 | 2020-03-31 | 南京工程学院 | 一种电力系统潮流计算中非标准变比变压器的建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107944682B (zh) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106356859B (zh) | 一种基于Matlab的直角坐标牛顿法潮流计算方法 | |
Ghatak et al. | A fast and efficient load flow technique for unbalanced distribution system | |
CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN106602570B (zh) | 一种基于Matlab的快速分解法潮流计算方法 | |
CN107196306B (zh) | 基于Matlab稀疏矩阵的快速分解法潮流计算方法 | |
CN104037763B (zh) | 一种适合含小阻抗支路系统的快速分解法潮流计算方法 | |
CN106229988B (zh) | 一种基于Matlab的极坐标牛顿法潮流计算方法 | |
CN104899396B (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
Seifi et al. | Power flow study and comparison of FACTS: Series (SSSC), Shunt (STATCOM), and Shunt-Series (UPFC) | |
CN106709243B (zh) | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 | |
CN105162109B (zh) | 基于灵敏度分析的直流潮流控制器优化配置方法 | |
CN107194131B (zh) | 基于Matlab稀疏矩阵的极坐标牛顿法潮流计算方法 | |
Jereminov et al. | An equivalent circuit formulation for three-phase power flow analysis of distribution systems | |
Semlyen | S-domain methodology for assessing the small signal stability of complex systems in nonsinusoidal steady state | |
CN106886617A (zh) | 一种含多vsc的多速率电磁暂态分网方法 | |
CN108599167A (zh) | 一种辐射状配电网的线性化潮流计算方法 | |
CN107944682A (zh) | 基于Matlab矩阵运算的潮流计算导纳矩阵计算方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN107181260B (zh) | 基于Matlab稀疏矩阵直角坐标牛顿法潮流计算方法 | |
CN107658880B (zh) | 基于关联矩阵运算的快速分解法系数矩阵计算方法 | |
CN107665184A (zh) | 基于关联矩阵运算的潮流计算导纳矩阵计算方法 | |
CN109659943A (zh) | 一种电力系统潮流计算的导纳矩阵计算方法 | |
CN107834562B (zh) | 基于Matlab矩阵运算的快速分解法系数矩阵计算法 | |
Osaloni et al. | Integrated AC/DC systems power flow solution using Newton-Raphson and Broyden approaches | |
CN107704686A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200407 Termination date: 20201116 |
|
CF01 | Termination of patent right due to non-payment of annual fee |