CN109802392B - 大规模配电网潮流计算方法及装置 - Google Patents
大规模配电网潮流计算方法及装置 Download PDFInfo
- Publication number
- CN109802392B CN109802392B CN201910141531.5A CN201910141531A CN109802392B CN 109802392 B CN109802392 B CN 109802392B CN 201910141531 A CN201910141531 A CN 201910141531A CN 109802392 B CN109802392 B CN 109802392B
- Authority
- CN
- China
- Prior art keywords
- current
- node
- load
- constant
- sub
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 71
- 239000011159 matrix material Substances 0.000 claims abstract description 108
- 238000002347 injection Methods 0.000 claims abstract description 42
- 239000007924 injection Substances 0.000 claims abstract description 42
- 238000004891 communication Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000000354 decomposition reaction Methods 0.000 claims description 12
- 239000000243 solution Substances 0.000 description 13
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 239000002699 waste material Substances 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明实施例提供一种大规模配电网潮流计算方法及装置,所述方法包括:将目标配电网络分解成若干个子网,获取用于求解子网间联络支路电流值的方程;计算各子网ZIP负荷模型的恒阻抗、恒电流和恒功率参数;根据恒阻抗参数计算各子网的节点导纳矩阵,并对节点导纳矩阵进行改进;根据当前迭代的节点电压,以及恒电流和恒功率参数,计算当前迭代的节点注入电流;根据改进后的节点导纳矩阵、当前迭代的节点注入电流和用于求解子网间联络支路电流的方程,获得当前迭代的联络支路电流;若前后两次迭代的联络支路电流间的差值大于收敛判据,则根据当前迭代的联络支路电流获得下一次迭代的节点电压相量。本发明实施例兼顾了潮流并行计算的效率与灵活性。
Description
技术领域
本发明实施例涉及电力系统分析与计算技术领域,更具体地,涉及一种大规模配电网潮流计算方法及装置。
背景技术
随着电力系统自动化水平不断提高,配电系统作为电力网络的重要组成部分,在电力系统中越来越受到重视,尤其是在优化配电网络中对配电网潮流进行计算。随着分布式电源、电动汽车等设备大量接入,配电网从传统的无源网络逐步变为有源的主动配电网,结构更为复杂,传统的潮流计算速度已无法满足大规模复杂配电网实时计算的要求。因此,研究配电网潮流并行算法,提高潮流并行计算效率已成为研究大规模配电系统的关键。
早期的潮流并行算法研究主要包括分块法、多重因子化法、稀疏矢量法和逆矩阵法。这些方法均是以线性方程组的前代回代和因子分解过程为切入点实现并行的,由于受进程数限定,存在当导纳矩阵维数不能被进程数整除时无法实现并行的问题,并行计算效率低。
J.R.Marti等人提出一种多区域戴维南等值法,利用网络中存在的集中参数线路或元件将网络分割为多个子网,先求出子网之间联络支路上流过的电流,再将联络支路电流纳入各子网的节点电压方程最后完成整个网络的求解。但采用该方法进行网络分割时,必须在集中参数元件上进行,且在并行计算时要求各子网采用相同的计算理论,限制了求解算法的灵活性。
发明内容
为了解决现有技术存在的并行计算效率低下、分网求解灵活性差的问题,本发明实施例提供一种大规模配电网潮流计算方法及装置。
第一方面,本发明实施例提供一种大规模配电网潮流计算方法,包括:
通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;
在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算各子网ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;
根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;
对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;
将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入所述用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;
计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,
若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量。
第二方面,本发明实施例提供一种大规模配电网潮流计算装置,包括:
分解模块,用于通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;
负荷模型求解模块,用于在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算各子网ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;
导纳矩阵获取模块,用于根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;
节点注入电流计算模块,用于对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;
联络支路电流计算模块,用于将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入所述用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;
收敛判断模块,用于计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,
若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量。
第三方面,本发明实施例提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如第一方面所提供的大规模配电网潮流计算方法的步骤。
第四方面,本发明实施例提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如第一方面所提供的大规模配电网潮流计算方法的步骤。
本发明实施例提供的大规模配电网潮流计算方法及装置,在节点处进行网络分解,有效保证各子网的计算规模均衡,避免了子网间互相等待浪费计算效率,提高了潮流并行计算效率;网络分解不受网络元件及计算方式限制,在方法上更具备灵活性;将负荷建模为ZIP模型,使得潮流计算结果更符合配电网实际运行情况,提高了潮流计算的收敛性和精度;提出了一种改进的节点导纳矩阵,使得配电网节点导纳矩阵可以直接求逆,从而进一步地提升了并行计算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的大规模配电网潮流计算方法的流程示意图;
图2为本发明实施例提供的配电网络分解示意图;
图3为本发明实施例提供的子网间联络示意图;
图4为本发明实施例提供的大规模配电网潮流计算装置的结构示意图;
图5为本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了便于理解,首先对本发明实施例中涉及到的名词进行解释。
节点:若干元件的汇合点,节点的特征是每一个瞬时电压相等,节点和系统短路发生点基本一致,因有可能发生短路而需要被列为短路点进行计算的点均被看作一个节点,比如各条母线、发电机出口、变压器高低压侧等。
节点的类型包括:PQ节点,比如变电所母线、负荷节点、浮游节点;PV节点,比如有一定无功功率储备的发电机母线或变电所、带有补偿设备的母线;平衡节点,也叫参考节点,一般该节点只有一个。
如图1所示,为本发明实施例提供的大规模配电网潮流计算方法的流程示意图,包括:
步骤100,通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;
可以理解的是,针对目标配电网络,选取多个分裂节点将所述目标配电网络分解成若干个子网。目标配电网络分解后,可以进一步获取各子网的线路拓扑参数和各子网中的负荷对应的负荷参数。
分裂节点的选取原则:子网规模的均衡性,即要保证分解后的各子网的计算时间和计算规模是相对均衡的;并行计算的效率,即进行分解后的并行计算效率要比分解前的计算效率有所提升;还可以考虑各种约束条件,计算出合适的分裂节点。本发明实施例不对选取分裂节点的具体方法作限制。
联络支路电流值反映了子网之间的信息交互,是完成潮流并行计算的重要参数。因此,将目标配电网进行分解后,还需获取求解联络支路电流值的方程。
根据基尔霍夫电流定律,可以列出各子网的支路电流方程。
根据同一分裂节点在各子网中电压相等的关系,可以列出各子网在分裂节点处的节点电压方程。
然后联立所述支路电流方程与节点电压方程,从而获得用于求解子网间联络支路电流值的方程。
下面以配电网络被分解为两个子网为例说明获得求解子网间联络支路电流值的方程的具体过程。设选取分裂节点为ε,将网络分解为a、b两部分,如图2所示,为配电网络分解示意图。
图3为子网间联络示意图,分裂节点ε分到子网a、b中,子网a、b仅通过分裂节点间的联络支路电流iε进行信息交互,假设电流方向如图所示。
由基尔霍夫电流定律可得子网a与子网b的支路电流方程为:
式(1)、(2)中,Yta、Ytb分别为子网a与子网b的节点导纳矩阵; 分别为子网a与子网b的k次迭代节点电压相量;分别为子网a与子网b的k次迭代节点注入电流相量;Lab为子网a中各节点与联络支路电流的连接关系,Lba为子网b中各节点与联络支路电流的连接关系,Lab、Lba均为元素非0即1的列向量,在分裂节点处元素为1。
由同一分裂节点在各子网中电压相等的关系,可得子网a与子网b在分裂点处的节点电压方程:
联立式(1)~(3),可得:
值得说明的是,当子网的个数大于2时,仍可以直接依据基尔霍夫电流定律列出各子网的支路电流方程,根据同一分裂节点在各子网中电压相等的关系列出各子网在分裂节点处的节点电压方程,只是方程的形式更复杂,方程数更多。
步骤200、在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;
在传统输配电系统潮流计算中负荷多处理为恒功率模型,然而随着越来越多的电压敏感性负荷接入系统中,单一的恒功率负荷模型不能真实地反映负荷的性质,从而无法准确地反映电网的潮流分布状况,导致求解结果的精度受到影响,收敛性差。
考虑负荷功率与其上电压的负荷电压静特性的潮流计算结果更符合配电网实际运行情况,尤其当配电网直接与负荷相连时。负荷电压静特性模型有ZIP模型和幂函数模型等形式。为了方便求解节点导纳矩阵与节点注入电流相量,本发明实施例采用ZIP负荷模型。
具体地,建立ZIP负荷模型为:
式(9)~(11)中:为负荷复功率相量;PL、QL分别为负荷实际有功功率相量、负荷实际无功功率相量,PN、QN分别为负荷额定有功功率相量、负荷额定无功功率相量;U为负荷节点电压相量,UN为负荷节点额定电压相量;a1+b1+c1=1,a2+b2+c2=1。
将式(10)~(11)代入式(9)得:
将式(12)简写为:
进一步地,对比式(12)和(13),可得:
获取各子网的负荷参数PN、QN,UN,将各子网的负荷参数PN、QN,UN,a1,b1,c1,a2,b2,c2代入式(14)~(16),可以得到各子网ZIP负荷模型的恒阻抗参数ZL、恒电流参数IL和恒功率参数
步骤300、根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;
具体地,本发明实施例中的节点导纳矩阵为考虑负荷恒阻抗部分后的节点导纳矩阵。
考虑负荷恒阻抗部分的自阻抗矩阵为:
ZT=Zt+ZL (17)
式(17)中,ZT为考虑负荷恒阻抗部分的自阻抗矩阵;Zt为未考虑负荷恒阻抗部分的自阻抗矩阵;ZL为负荷恒阻抗部分矩阵,即前述步骤中的各子网ZIP负荷模型的恒阻抗参数。
首先,根据现有技术可以计算获得未考虑负荷恒阻抗部分的自阻抗矩阵Zt。
然后,根据所述未考虑负荷恒阻抗部分的自阻抗矩阵Zt与所述ZIP负荷模型的恒阻抗参数ZL利用式(17)计算获得考虑负荷恒阻抗部分后的自阻抗矩阵ZT。
根据上述方法计算出来的节点导纳矩阵为奇异矩阵,而从用于求解子网间联络支路电流值的式(4)可以得知,计算联络支路电流值时会用到节点导纳矩阵的逆,而对奇异矩阵直接求逆是没有意义的,因此,还需对所获得的节点导纳矩阵进行改进,以使得改进后的节点导纳矩阵为非奇异矩阵,从而提高潮流计算效率。
具体地,设一系统共有n节点,则该系统的电力网络方程可描述为:
式(5)中,Yii为节点导纳矩阵的对角元,为节点i的自导纳,其值等于与该节点i直接相连的所有支路的导纳总和;Yij为节点导纳矩阵的非对角元(i≠j),为互导纳,其值等于连接节点i、j支路导纳总和的负数,且Yij=Yji。
式(5)可简写为:
IB=YBUB (6)
式(6)中,IB为节点注入电流列向量,规定流入网络为正;YB为电力网络的节点导纳矩阵;UB为节点电压列向量。
现假设一配电网共有n个节点,其中节点1为平衡节点,节点2~n为PQ节点,设平衡节点电压固定为Us(相角为0)。式(4)需对节点导纳矩阵求逆,而配电网节点导纳矩阵YB为奇异矩阵,直接求逆没有实际意义,因此需对其进行改进,将式(5)改写为:
式(7)中,修改节点导纳矩阵第一行元素,为满足等式IB=YBUB,节点电流相量与节点电压相量也作相应的修改。改进节点导纳矩阵如式(8)所示,式(8)可直接求逆,且不影响子网间联络支路电流值的求解。
步骤400、对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;
具体地,根据现有技术可得仅考虑负荷恒功率部分的节点注入电流相量的计算公式为:
那么,基于ZIP负荷模型的节点注入电流相量为:
可以理解的是,首先根据现有技术计算出仅考虑负荷恒功率部分的节点注入电流相量,然后根据所述仅考虑负荷恒功率部分的节点注入电流相量与所述ZIP负荷模型的恒电流参数,利用式(19)可得当前迭代的节点注入电流相量Ik。
步骤500、将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;
具体地,在前述步骤中已经通过联立支路电流方程与节点电压方程获得了用于求解子网间联络支路电流值的方程,将改进后的节点导纳矩阵和当前迭代的节点注入电流相量作为已知量,代入所述用于求解子网间联络支路电流值的方程,即可获得当前迭代的联络支路电流值。
步骤600、计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,
若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量。
具体地,比较前后两次迭代计算获得的联络支路电流值之间的差值Δi,若差值Δi小于等于预设收敛判据δ,则说明计算收敛了,可以结束迭代过程,并将所述当前迭代的联络支路电流值作为需求的联络支路电流值。
若所述差值Δi大于预设收敛判据,则说明计算结果还不够精确,需继续进行迭代,将所述当前迭代的联络支路电流值回代到各子网的支路电流方程中,获得下一次迭代的节点电压相量,然后开始下一次迭代。
然后,开始下一次迭代,即重复步骤400-500,根据所述下一次迭代的节点电压相量Uk+1,以及所述ZIP负荷模型的恒电流参数和恒功率参数,利用公式(18)~(19)计算获得下一次迭代的节点注入电流相量Ik+1,然后,将改进后的节点导纳矩阵Yt和下一次迭代的节点注入电流相量Ik+1代入用于求解子网间联络支路电流值的方程,获得下一次迭代的联络支路电流值
值得说明的是,第一次迭代的节点电压相量根据如下步骤获得:
获取各子网的线路拓扑参数,根据所述各子网的线路拓扑参数在各子网内部进行潮流计算,获得各子网的节点电压相量初始值,将所述各子网的节点电压相量初始值作为第一次迭代的节点电压相量。
本发明实施例提供的大规模配电网潮流计算方法,在节点处进行网络分解,有效保证各子网的计算规模均衡,避免了子网间互相等待浪费计算效率,提高了潮流并行计算效率;网络分解不受网络元件及计算方式限制,在方法上更具备灵活性;将负荷建模为ZIP模型,使得潮流计算结果更符合配电网实际运行情况,提高了潮流计算的收敛性和精度;提出了一种改进的节点导纳矩阵,使得配电网节点导纳矩阵可以直接求逆,从而进一步地提升了并行计算效率。
如图4所示,为本发明实施例提供的大规模配电网潮流计算装置的结构示意图,包括:分解模块410、负荷模型求解模块420、导纳矩阵获取模块430、节点注入电流计算模块440、联络支路电流计算模块450和收敛判断模块460,其中,
分解模块410,用于通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;
负荷模型求解模块420,用于在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算各子网ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;
导纳矩阵获取模块430,用于根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;
节点注入电流计算模块440,用于对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;
联络支路电流计算模块450,用于将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入所述用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;
收敛判断模块460,用于计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,
若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量。
该装置用于实现前述实施例中的大规模配电网潮流计算方法。因此,在前述实施例中所述的大规模配电网潮流计算方法的描述和定义,可以用于本发明实施例中各个执行模块的理解,在此不再赘述。
本发明实施例提供的大规模配电网潮流计算装置,在节点处进行网络分解,有效保证各子网的计算规模均衡,避免了子网间互相等待浪费计算效率,提高了潮流并行计算效率;网络分解不受网络元件及计算方式限制,在方法上更具备灵活性;将负荷建模为ZIP模型,使得潮流计算结果更符合配电网实际运行情况,提高了潮流计算的收敛性和精度;提出了一种改进的节点导纳矩阵,使得配电网节点导纳矩阵可以直接求逆,从而进一步地提升了并行计算效率。
图5为本发明实施例提供的电子设备的实体结构示意图,如图5所示,该电子设备可以包括:处理器(processor)510、通信接口(Communications Interface)520、存储器(memory)530和通信总线540,其中,处理器510,通信接口520,存储器530通过通信总线540完成相互间的通信。处理器510可以调用存储在存储器530上并可在处理器510上运行的计算机程序,以执行上述各方法实施例所提供的大规模配电网潮流计算方法,例如包括:通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算各子网ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入所述用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量。
此外,上述的存储器530中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本发明实施例还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现上述各方法实施例提供的大规模配电网潮流计算方法,例如包括:通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算各子网ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入所述用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.一种大规模配电网潮流计算方法,其特征在于,包括:
通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;
在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算各子网ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;
根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;
对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;
将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入所述用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;
计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,
若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量;
其中,所述根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵的步骤,具体为:
计算未考虑负荷恒阻抗部分的自阻抗矩阵,根据所述未考虑负荷恒阻抗部分的自阻抗矩阵与所述ZIP负荷模型的恒阻抗参数,计算获得考虑负荷恒阻抗部分后的自阻抗矩阵;
根据所述考虑负荷恒阻抗部分后的自阻抗矩阵获得考虑负荷恒阻抗部分后的节点导纳矩阵;
其中,所述对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵的步骤,具体为:
其中,Yii为节点导纳矩阵的对角元,为节点i的自导纳,其值等于与该节点i直接相连的所有支路的导纳总和;Yij为节点导纳矩阵的非对角元(i≠j),为互导纳,其值等于连接节点i、j支路导纳总和的负数,且Yij=Yji;
其中,所述根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量的步骤,具体为:
根据当前迭代的节点电压相量和所述ZIP负荷模型的恒功率参数计算仅考虑负荷恒功率部分的节点注入电流相量;
根据所述仅考虑负荷恒功率部分的节点注入电流相量与所述ZIP负荷模型的恒电流参数计算获得当前迭代的节点注入电流相量;
其中,所述当前迭代的节点注入电流相量的计算公式为:
2.根据权利要求1所述的方法,其特征在于,所述获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程的步骤,具体为:
根据基尔霍夫电流定律获得各子网的支路电流方程;
根据同一分裂节点在各子网中电压相等的关系,获得各子网在分裂节点处的节点电压方程;
联立所述支路电流方程与节点电压方程,获得用于求解子网间联络支路电流值的方程。
3.根据权利要求1所述的方法,其特征在于,还包括:
获取各子网的线路拓扑参数,根据所述各子网的线路拓扑参数在各子网内部进行潮流计算,获得各子网的节点电压相量初始值,将所述各子网的节点电压相量初始值作为第一次迭代的节点电压相量。
4.一种大规模配电网潮流计算装置,其特征在于,包括:
分解模块,用于通过选取多个分裂节点将目标配电网络分解成若干个子网,获取各子网的支路电流方程、各子网在分裂节点处的节点电压方程以及用于求解子网间联络支路电流值的方程;
负荷模型求解模块,用于在各子网内建立ZIP负荷模型,获取各子网的负荷参数,并根据所述负荷参数计算各子网ZIP负荷模型的恒阻抗参数、恒电流参数和恒功率参数;
导纳矩阵获取模块,用于根据所述ZIP负荷模型的恒阻抗参数计算获得各子网的节点导纳矩阵,并对计算获得的节点导纳矩阵进行改进,以使得改进后的所述节点导纳矩阵为非奇异矩阵;
节点注入电流计算模块,用于对于任一次迭代的节点电压相量,根据当前迭代的节点电压相量,以及所述ZIP负荷模型的恒电流参数和恒功率参数,计算获得当前迭代的节点注入电流相量;
联络支路电流计算模块,用于将改进后的节点导纳矩阵和所述当前迭代的节点注入电流相量代入所述用于求解子网间联络支路电流值的方程,获得当前迭代的联络支路电流值;
收敛判断模块,用于计算当前迭代的联络支路电流值与上一次迭代的联络支路电流值之间的差值,若所述差值小于等于预设收敛判据,则停止迭代,并将当前迭代的联络支路电流值作为需求的联络支路电流值;或者,
若所述差值大于预设收敛判据,则将所述当前迭代的联络支路电流值回代到所述各子网的支路电流方程中,获得下一次迭代的节点电压相量;
其中,所述导纳矩阵获取模块具体用于:
计算未考虑负荷恒阻抗部分的自阻抗矩阵,根据所述未考虑负荷恒阻抗部分的自阻抗矩阵与所述ZIP负荷模型的恒阻抗参数,计算获得考虑负荷恒阻抗部分后的自阻抗矩阵;
根据所述考虑负荷恒阻抗部分后的自阻抗矩阵获得考虑负荷恒阻抗部分后的节点导纳矩阵;
其中,Yii为节点导纳矩阵的对角元,为节点i的自导纳,其值等于与该节点i直接相连的所有支路的导纳总和;Yij为节点导纳矩阵的非对角元(i≠j),为互导纳,其值等于连接节点i、j支路导纳总和的负数,且Yij=Yji;
其中,所述节点注入电流计算模块具体用于:
根据当前迭代的节点电压相量和所述ZIP负荷模型的恒功率参数计算仅考虑负荷恒功率部分的节点注入电流相量;
根据所述仅考虑负荷恒功率部分的节点注入电流相量与所述ZIP负荷模型的恒电流参数计算获得当前迭代的节点注入电流相量;
其中,所述当前迭代的节点注入电流相量的计算公式为:
5.一种电子设备,其特征在于,包括:
至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中:
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至3任一所述的方法。
6.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至3任一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910141531.5A CN109802392B (zh) | 2019-02-26 | 2019-02-26 | 大规模配电网潮流计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910141531.5A CN109802392B (zh) | 2019-02-26 | 2019-02-26 | 大规模配电网潮流计算方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109802392A CN109802392A (zh) | 2019-05-24 |
CN109802392B true CN109802392B (zh) | 2020-10-13 |
Family
ID=66561345
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910141531.5A Active CN109802392B (zh) | 2019-02-26 | 2019-02-26 | 大规模配电网潮流计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109802392B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111523800B (zh) * | 2020-04-22 | 2023-03-31 | 中车株洲电力机车研究所有限公司 | 一种地铁潮流计算中节点电导矩阵的快速计算方法 |
CN111463782B (zh) * | 2020-04-29 | 2021-12-17 | 南方电网科学研究院有限责任公司 | 一种电压敏感负荷模型及参数辨识方法 |
CN111614099B (zh) * | 2020-05-22 | 2023-03-28 | 国家电网有限公司 | 一种电力系统潮流仿真计算方法及系统 |
CN112202177B (zh) * | 2020-09-29 | 2022-03-04 | 国网河南省电力公司洛阳供电公司 | 一种计算配电网支路系统潮流的方法 |
CN112290604B (zh) * | 2020-10-15 | 2023-06-30 | 珠海博威电气股份有限公司 | 计及负荷特性的配电网协调优化方法和装置、存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100373733C (zh) * | 2004-01-06 | 2008-03-05 | 中国电力科学研究院 | 电力系统潮流分网并行计算方法 |
CN100576223C (zh) * | 2007-12-12 | 2009-12-30 | 国电南瑞科技股份有限公司 | 基于计算模型拼接的互联系统分布式潮流计算方法 |
JP5743195B2 (ja) * | 2011-04-15 | 2015-07-01 | 富士電機株式会社 | ループ状系統用潮流計算方法及び潮流計算装置 |
CN103683284B (zh) * | 2013-12-26 | 2016-08-17 | 国家电网公司 | 一种含分布式电源的配电网三相不平衡潮流计算方法 |
CN108306299B (zh) * | 2018-01-29 | 2022-11-11 | 中国电力科学研究院有限公司 | 一种配电网异步迭代分布式三相潮流计算方法及系统 |
-
2019
- 2019-02-26 CN CN201910141531.5A patent/CN109802392B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109802392A (zh) | 2019-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109802392B (zh) | 大规模配电网潮流计算方法及装置 | |
Wang et al. | Explicit conditions on existence and uniqueness of load-flow solutions in distribution networks | |
JP5815382B2 (ja) | 平衡配電系統のための電力潮流解析 | |
JP6025575B2 (ja) | 三相電力潮流解析のためのシステムおよび方法 | |
Ju et al. | An extension of FBS three-phase power flow for handling PV nodes in active distribution networks | |
CN103066593B (zh) | 含多类型分布式电源的弱环配电网三相潮流计算方法 | |
Lee | Voltage stability assessment using equivalent nodal analysis | |
CN111625914A (zh) | 一种短路电流计算方法及系统 | |
Jabr et al. | A Fortescue approach for real-time short circuit computation in multiphase distribution networks | |
CN107342586B (zh) | 一种用于确定配电网静态等值的方法及装置 | |
CN102866315B (zh) | 含逆变型分布式电源配电网的对称故障分析方法 | |
CN111682546B (zh) | 基于灵敏度分析的直流潮流改进算法 | |
JP5933427B2 (ja) | 非接地配電系統の三相電力潮流解析のための方法およびシステム | |
Fernandes et al. | Symmetrical components based state estimator for power distribution systems | |
CN112234598B (zh) | 一种电磁暂态仿真初始化方法 | |
CN111725812B (zh) | 大规模配电系统的潮流计算方法、装置、设备和存储介质 | |
Abdel-Akher et al. | Development of unbalanced three-phase distribution power flow analysis using sequence and phase components | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
Lourenco et al. | Steady-state solution for power networks modeled at bus section level | |
CN111834996B (zh) | 电网线损计算方法和装置 | |
Husain et al. | Load flow analysis of radial and mesh distribution system using ZIP model | |
Fahmy et al. | Analytical computation of power grids’ sensitivity coefficients with voltage-dependent injections | |
CN105610156B (zh) | 一种多线路并发合环方法 | |
Vemprala et al. | Open-source poly-phase distribution system power flow analysis tool (DxFlow) | |
CN108134392B (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 |