CN101141064A - 通过交换边界节点状态和网损信息完成分布式潮流分析的方法 - Google Patents
通过交换边界节点状态和网损信息完成分布式潮流分析的方法 Download PDFInfo
- Publication number
- CN101141064A CN101141064A CNA2007101217928A CN200710121792A CN101141064A CN 101141064 A CN101141064 A CN 101141064A CN A2007101217928 A CNA2007101217928 A CN A2007101217928A CN 200710121792 A CN200710121792 A CN 200710121792A CN 101141064 A CN101141064 A CN 101141064A
- Authority
- CN
- China
- Prior art keywords
- boundary
- node
- regional
- network
- power grid
- 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
- 238000000034 method Methods 0.000 title claims abstract description 72
- 230000008569 process Effects 0.000 claims abstract description 32
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 186
- 239000013598 vector Substances 0.000 claims description 53
- 238000005206 flow analysis Methods 0.000 claims description 21
- 238000002347 injection Methods 0.000 claims description 19
- 239000007924 injection Substances 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 12
- 230000011218 segmentation Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 7
- 239000000243 solution Substances 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 4
- 238000002940 Newton-Raphson method Methods 0.000 claims description 2
- 230000006855 networking Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 15
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 3
- 239000003990 capacitor Substances 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 229910052708 sodium Inorganic materials 0.000 description 3
- 239000011734 sodium Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- VWDWKYIASSYTQR-UHFFFAOYSA-N sodium nitrate Chemical compound [Na+].[O-][N+]([O-])=O VWDWKYIASSYTQR-UHFFFAOYSA-N 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- Y04S10/54—
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种通过交换边界节点状态和网损信息完成分布式潮流分析的方法,包括:根据电力***实际调度管理方式构建分布式潮流分析***;从电力***拓扑连接关系出发,根据实际***的运行状况,采用带有边界区域的切分方法对互联电网进行划分,明确潮流计算服务器和协调计算服务器的计算对象和数据来源;区域潮流计算服务器和上级协调计算服务器计算用数据和参数初始化,协调计算服务器调用区域潮流计算服务器共同完成全网一体化潮流分解协调求解过程;服务器间调用流程,区域电网潮流计算服务器将计算结果反馈给协调计算服务器,以使协调计算服务器获得收敛的全网一体化潮流计算结果。
Description
技术领域
本发明涉及一种完成分布式潮流分析的方法,特别涉及一种通过交换边界节点状态和网损信息完成分布式潮流分析的方法,属于电力***分布式仿真技术领域。
背景技术
随着电力***互联不断加强,电网构成更加复杂,运行难度大大增加。发生局部故障时,如果处理不当,就会发展为危及整个***的重大事故。为了提高互联电力***的安全稳定性,对整个电网进行无简化的一体化仿真分析显得越来越有必要的。
我国电网由多级调度中心协调管理和控制,采用“分级管理,分层控制,分布处理”的管理体系。如图1所示,级别较低的区域电网调度中心,如省调、地调等,只负责维持所辖区域内的功率平衡,管理和维护***运行状态和参数;而级别较高的调度中心,如网调、国调等,则要负责区域间功率交换控制,协调区域电网调度活动。由于我国电力***具有广域分布、参数海量、模型复杂的特点,高级别调度中心无法直接获取所辖电网状态和参数,需要由低级别的调度中心层层转发或上报所需数据。这使得全国电网参数和实时状态的采集、同步、整合所需时间较长,且维护难度较大。同时,随着电力市场化进程的不断深入,在市场竞争中,各调度机构更加重视对自身数据和信息的保护,加厚了各个调度机构的信息壁垒,进一步阻碍了***内数据和信息的共享。因此,传统的集中式的全网一体化仿真计算难以在线实现。
采用分布式计算技术的电网分析计算方法可以在保持参与计算的各方自身数据和计算资源的独立性的同时获得全网一体化仿真分析结果,具有快速高效、信息安全、应用灵活和易于扩展的特点。分布式潮流计算的重要功能是可以在保留各调度中心计算独立性的同时获得全网统一的仿真结果,因而该技术有望成为解决大规模互联***一体化仿真的有效手段。目前,国内外对互联电网分布式潮流及其相关领域的研究,主要集中在以下几个方面:
1.合理的电力网络切分方法;
2.并行潮流求解算法向分布式环境中的移植和测试;
3.潮流分解协调求解过程的等值网络计算和***不平衡功率分配;
可以看出,现有的互联电网一体化潮流分解协调计算方法的研究,主要集中在如何根据实际电力***区域情况进行潮流计算的分解和对已有并行潮流算法改进方面。这些研究大多缺少对广域网络环境中的分布式计算环境所具有的高通信延时、数据和计算资源异构等特性的考虑,也没有建立能够统一处理全网潮流平衡和不平衡功率分配的分解协调计算模型,因而在实际电力***中应用能力不足。
发明内容
本发明的目的在于,提供一种通过交换边界节点状态和网损信息完成分布式潮流分析的方法,以联络线电流为协调变量,计算过程中只需要传递少量边界信息,可以灵活的处理全网不平衡功率分配的问题。
本发明是采用以下技术手段实现的:
一种通过交换边界节点状态和网损信息完成分布式潮流分析的方法,包括:
步骤1;根据电力***实际调度管理方式构建分布式潮流分析***;所述的电力***由区域调度中心和上级调度中心协调一致、信息互联;所述的分布式潮流计算***由区域调度中心内的潮流计算服务器和上级调度中心内的协调计算服务器通过网络连接构成;其中;潮流计算服务器,负责其所辖电网的潮流计算;而上级调度中心内的协调计算服务器,负责协调各区域电网潮流计算过程;
步骤2;从电力***拓扑连接关系出发,根据实际***的运行状况,采用带有边界区域的切分方法对***进行切分,按照实际网络连接情况,对互联电网进行划分,明确潮流计算服务器和协调计算服务器的计算对象和数据来源;
所述的***定义为(S0);包括两个区域***(S1)和(S2),通过联络线l相互连接;其中(B1)(B2)分别代表(S1)(S2)区域的边界节点,(Sln 1)(Sln 2)分别为(S1)(S2)中除了边界节点外其他节点组成的网络;将边界节点B1和B2分别***,虚拟出边界节点和并将联络线l与其两端的虚拟边界节点和单独看作一个SB,则区域起到连接其他所有区域作用,定义为边界区域;
所述区域S1和S2属于区域电网调度中心管辖范围,区域SB属于上级调度中心管辖范围;在这种切分方式下,全网潮流的收敛条件为:区域S1和S2中潮流计算达到收敛,它们内部节点功率均达到平衡,且边界节点状态满足
其中uB和θB是区域电网内潮流计算收敛结果中边界节点B1和B2的电压幅值和相角向量;和是边界区域SB内潮流计算收敛时边界节点B1和B2的电压幅值和相角向量;是区域电网潮流计算收敛结果中边界节点B1和B2上来自联络线的注入电流相量组成的向量,以流入为正方向;是边界区域SB中潮流计算收敛时从边界节点B1和B2向联络线上注入的电流相量组成的向量,以流出为正方向。
步骤3;区域潮流计算服务器和上级协调计算服务器计算用数据和参数初始化,包括以下步骤:
步骤3.1:区域电网调度中心内的潮流计算服务器初始化所辖电网的网络参数,包括所有参与潮流计算的网络节点的类型、负荷和发电机出力设定值、发电机机端节点电压幅值以及松弛节点相角等;
步骤3.2:区域电网调度中心内的潮流计算服务器初始化边界条件信息,包括边界节点的编号、名称等、边界条件的组成,即边界节点的注入电流相量组成的向量和区域电网松弛节点相角;
步骤3.3:区域电网调度中心内的潮流计算服务器初始化潮流计算的参数,包括所使用的迭代求解方法、最大迭代次数、节点电压上下限、用以判断潮流收敛的最小功率偏差;
步骤3.4:上级调度中心内的协调计算服务器初始化边界区域网络参数,包括边界节点的名称、联络线阻抗;
步骤3.5:上级调度中心内的协调计算服务器设定全网一体化潮流协调求解控制参数;最大迭代次数、各区域负担全网有功网损的比例系数,以及边界节点功率平衡精度要求,即判定边界节点功率平衡所需的边界节点上不平衡的电流向量最大幅值。
步骤4:协调计算服务器调用区域潮流计算服务器共同完成全网一体化潮流分解协调求解过程;服务器间调用流程所示,包括以下基本步骤:
步骤4.1:上级调度中心内的协调计算服务器按照步骤3.4、3.5初始化本地基础网络数据和计算参数,区域调度中心内区域潮流计算服务器按照步骤3.1~3.3初始化本地基础网络数据和计算参数;
步骤4.2:协调计算服务器设定初始边界条件,即区域电网边界节点上来自联络线的注入电流相量组成的向量IB及各区域电网松弛节点相角设定值组成的向量Θ,并按照步骤2中网络化分关系,将相应的注入电流和相角设定值发送给区域电网潮流计算服务器;
步骤4.3:区域电网潮流计算服务器启动本区域潮流计算过程,根据所得到的边界条件和内部节点参数求解区域潮流方程,从收敛的潮流结果中得到偏差边界节点的电压VB i和区域网损信息Ploss i,i=1,2,并将信息发送给协调计算服务器;
步骤4.4:协调计算服务器根据得到的区域潮流计算结果信息以及初始化过程中设定的参数,计算边界节点的电流不平衡量ΔI,和区域电网负担全网有功网损的不平衡量ΔP;
步骤4.5;若‖[ΔI,ΔP]T‖2≤ξ则判断全网潮流收敛,计算过程结束,其中‖·‖2表示向量的二范数,ξ为小的正常数;否则根据[ΔI,ΔP]T计算边界条件修正量[ΔIB,ΔΘ]T后更新边界条件设定值,返回步骤4.2;
步骤5:区域电网潮流计算服务器将计算结果反馈给协调计算服务器,以使协调计算服务器获得收敛的全网一体化潮流计算结果。
前述的根据计算边界节点的电流不平衡量和区域电网有功网损的不平衡量,判断全网一体化潮流是否一致收敛;需要同时满足以下基本条件:
A:边界节点功率达到平衡,即满足方程;
B:全网有功网损在各区域电网间合理分配,满足以下方程,
其中∏={ηj}是区域电网S1和S2负担全网有功网损的比例向量,i=1,2, Pi loss是从区域电网潮流计算结果中统计出的有功网损; 是由各区域电网有功网损相加得到的全网有功网损;ΔPloss i是区域电网负担全网有功网损的不平衡量,定义为区域电网有功网损不平衡量。
前述的边界条件包括:区域电网的边界母线上来自联络线的注入电流相量组成的向量IB以及区域电网松弛节点相角设定值Θ。
前述的区域潮流计算服务器根据给定的边界条件计算区域潮流;其中所述的区域电网S1:
设电网由n个节点和若干条支路组成,潮流计算模型可表述为非线性方程组形式,即
其中 i=1,2,...,n为节点电压相量,g为节点功率平衡方程,Pi SP和Qi SP是节点i的给定有功和无功功率,Gij+jBij是节点i和节点j之间的导纳,j∈i表示与i节点相联的所有节点,包括i节点自身。
前述的S1区域所有节点由内部节点Sln 1和边界节点B1两部分,对应的以 代表内部节点的电压相量组成的向量,以 代表边界节点电压相量组成的向量;对应地,可将方程(2)描述的节点功率平衡方程组也可以划分为两类,内部节点Sln 1的功率平衡方程组gln 1和边界节点的功率平衡方程组gB 1;如区域电网计算潮流时不选择边界节点作为松弛节点,在完成给定内部节点的中发电机节点的有功出力、节点电压幅值,负荷节点的有功和无功负荷,以及松弛节点的电压幅值和相角情况下,内部节点功率平衡方程gln 1与内部节点和边界节点电压相关,边界节点功率平衡方程gB 1与内部节点电压、边界节点电压及其来自联络线的注入电流相量相关,其表述为如下形式,
其中 是由边界节点B1上来自联络线的注入电流相量组成的向量,表示节点k上来自联络线的注入电流相量的共轭,其他变量意义同式(3)。
前述的完成步骤3.1中的基础数据初始化过程后,再给定IB 1和区域电网S1内松弛节点电压相角设定值θref 1即可采用普通的Newton-Raphson方法求解方程(4)式,得到区域电网S1内的潮流。
前述的步骤4.3中协调计算服务器需要根据边界条件、联络线参数和区域电网潮流计算结果来计算边界节点的电流不平衡量ΔI,和区域电网负担全网有功网损的不平衡量ΔP,即[ΔI,ΔP]T。具体计算流程如图6所示,包括以下步骤:
a:协调计算服务器设定边界条件IB和Θ,发送至区域潮流计算服务器;
b:区域潮流计算服务器按照步骤(4)式计算本区域电网潮流,从潮流结果中提取VB和Ploss信息,其中VB是边界节点电压相量组成的向量,Ploss是区域电网有功网损信息,并将这些信息发送到协调计算服务器;
c:协调计算服务器根据边界节点的电压VB,按照(5)式计算边界节点来自联络线的注入电流理论值 以流出边界节点为正方向,
d协调计算服务器根据设定的区域电网负担全网有功网损的比例向量∏和各区域电网潮流结果中统计出的有功网损信息Ploss,根据(2)式计算区域电网有功网损不平衡量组成的向量
前述的协调计算服务器采用改进的Jacobian-free Newton-GMRES(m)方法计算边界条件的修正量,并迭代直至‖[ΔI,ΔP]T‖2≤ξ,具体步骤如下:
a.令k=-1,令xB=[Re(IB),Im(IB),Θ]T为边界条件组成的向量,其中Re(IB)和Im(IB)分别表示边界节点来自联络线的注入电流相量实部和虚部组成的向量;给定初始边界条件xB 0,选择任意非奇异矩阵M0为预处理矩阵;
b.k=k+1,重复步骤a~i,直至‖[ΔI,ΔP]T‖2≤ξ满足或是k>max NIter,max NIter为Newton迭代次数限制,结束;
c.令 为边界条件为xB k时通过步骤a~e计算得到的边界节点电流不平衡量和区域有功网损不平衡量,为了计算需要将ΔI的实部和虚部分开;
d. l=1,ρ=β=‖r0‖2,v1=r0,errtolG=ε‖r0‖‖2>0,1>ε>0;
e.若ρ>errtolG并且l<maxGIter,maxGIter是GMRES迭代次数限制,则l=l+1,zl=Mkvl, vl+1=ΔGl/w,w为一个小的正常数,一般取10-5>w>0;
f.修正预处理矩阵
h.
i.修正预处理矩阵 返回步骤b;
前述的步骤f和i为内层和外层预处理矩阵修正。
本发明与现有技术相比,具有以下明显的优势和有益效果:
本发明通过交换边界节点状态和网损信息完成分布式潮流分析的方法,针对我国电力***的分布式监控技术及电网互联的现状,提出了通过交换边界节点状态和网损信息完成分布式潮流分析的方法,该方法以联络线电流为协调变量,计算过程中只需要传递少量边界信息,可以灵活的处理全网不平衡功率分配的问题,计算过程中各区域电网可以使用独立的相角参考节点,具有模型通用,数据接口简化、收敛性能高、协调求解过程所需通信次数较少等优点,适合在广域网络中的异构计算环境中应用。
附图说明
图1为现有技术多级调度中心连接示意图;
图2为分布式一体化潮流分析基本流程示意图;
图3为电力***分布式潮流分析***示意图;
图4为带边界区域的***切分方式示意图;
图5为全网一体化潮流分解协调求解过程示意图;
图6为协调计算服务器计算[ΔI,ΔP]T过程示意图;
图7为二个区域互联协调示意图;
图8为带边界区域的***切分方式示意图;
图9为IEEE39***的切分示意图。
具体实施方式
下面结合附图对本发明的具体实施例加以说明:
请参阅图2所示,为分布式一体化潮流分析基本流程的示意图;
为一种借助Jaobian-free Newton GMRES(m)迭代算法,在保留各个调度中心计算独立性的前提下,只需要通过交换边界节点状态和区域电网有功网损信息即可完成全网潮流的协调求解的方法。可以完成满足以下基本要求的全网一体化潮流计算:
a、全网所有节点都满足功率平衡;
b、全网有功网损在各个区域电网间合理分配,即各个区域电网按照一定的比例分担的全网有功网损;
c、各区域电网潮流计算过程中可以采用独立的电压相角参考节点;
请参阅图3所示,为全网一体化潮流分析连接示意图;其中,区域调度中心内的潮流计算服务器,负责其所辖电网的潮流计算;而上级调度中心内的协调计算服务器,负责协调各区域电网潮流计算过程,即通过改变区域电网潮流计算所需的边界条件来影响区域电网潮流计算结果,直至所有区域电网的潮流计算同时达到和全网集中求解收敛时一样的潮流结果。
请参阅图4所示,为带边界区域的***切分方式示意图;
从图中可以看出,***S0包括两个区域***S1和S2,它们通过联络线l相互连接。其中B1代表S1区域的边界节点,Sln 1为S1中除了边界节点外其他节点组成的网络。S2的情况类似。若将边界节点B1和B2分别***开来,虚拟出边界节点和并将联络线l与其两端的虚拟边界节点和单独看作一个SB,则区域起到连接其他所有区域作用,可称为边界区域。区域S1和S2属于区域电网调度中心管辖范围,而区域SB属于上级调度中心管辖范围。在这种切分方式下,全网潮流的收敛条件为:区域S1和S2中潮流计算达到收敛,它们内部节点功率均达到平衡,且边界节点状态满足
其中uB和θB是区域电网内潮流计算收敛结果中边界节点B1和B2的电压幅值和相角向量;和是边界区域SB内潮流计算收敛时边界节点B1和B2的电压幅值和相角向量;是区域电网潮流计算收敛结果中边界节点B1和B2上来自联络线的注入电流相量组成的向量,以流入为正方向;是边界区域SB中潮流计算收敛时从边界节点B1和B2向联络线上注入的电流相量组成的向量,以流出为正方向。
请参阅图5所示,全网一体化潮流分解协调求解过程示意图;可以看出,包括以下步骤;
a、上级调度中心内的协调计算服务器按照步骤3.4、3.5初始化本地基础网络数据和计算参数,区域调度中心内区域潮流计算服务器按照步骤3.1~3.3初始化本地基础网络数据和计算参数;
b、协调计算服务器设定初始边界条件,即区域电网边界节点上来自联络线的注入电流相量组成的向量IB及各区域电网松弛节点相角设定值组成的向量Θ,并按照步骤2中网络化分关系,将相应的注入电流和相角设定值发送(网络通信)给区域电网潮流计算服务器;
c、区域电网潮流计算服务器启动本区域潮流计算过程,根据所得到的边界条件和内部节点参数求解区域潮流方程,从收敛的潮流结果中得到偏差边界节点的电压VB i和区域网损信息Ploss i,i=1,2,并将这些信息发送(网络通信)给协调计算服务器;
d、协调计算服务器根据得到的区域潮流计算结果信息以及初始化过程中设定的参数,计算边界节点的电流不平衡量ΔI,和区域电网负担全网有功网损的不平衡量ΔP;
e、若‖[ΔI,ΔP]T‖2≤ξ则判断全网潮流收敛,计算过程结束,其中‖·‖2表示向量的二范数,ξ为小的正常数;否则根据[ΔI,ΔP]T计算边界条件修正量[ΔIB,ΔΘ]T后更新边界条件设定值,返回步骤(4.2)。
请参阅图6所示,为协调计算服务器计算[ΔI,ΔP]T过程示意图;
a、协调计算服务器设定边界条件IB和Θ,发送至区域潮流计算服务器;
b、区域潮流计算服务器按照步骤4.3.1~4.3.3计算本区域电网潮流,从潮流结果中提取VB和Ploss信息,其中VB是边界节点电压相量组成的向量,Ploss是区域电网有功网损信息,并将这些信息发送到协调计算服务器;
c、协调计算服务器根据边界节点的电压VB,按照(5)式计算边界节点来自联络线的注入电流理论值 以流出边界节点为正方向,
d、协调计算服务器根据设定的区域电网负担全网有功网损的比例向量∏和各区域电网潮流结果中统计出的有功网损信息Ploss,根据(2)式计算区域电网有功网损不平衡量组成的向量
在具体实施例中,本发明按照以下3个阶段实施
1、阶段1:从电力***拓扑连接关系出发,根据实际***的运行状况,采用带有边界区域的切分方法对***进行切分。
电力***具有区域运营、分布式管理的特点,区域电网之间耦合较弱,通过少数的联络线连接,一般的互联电网可以简化成的两区域互联***形式。如图7所示,***S0包括两个区域***S1和S2,它们通过联络线l相互连接。其中B1代表S1区域的边界节点,Sln 1为S1中除了边界节点外其他节点组成的网络。S2的情况类似。若将边界节点B1和B2分别***开来,虚拟出边界节点和并将联络线l与其两端的虚拟边界节点和单独看作一个SB,则区域起到连接其他所有区域作用,可称为边界区域。区域S1和S2属于区域电网调度中心管辖范围,而区域SB属于上级调度中心管辖范围。
请参阅图8所示,为带边界区域的***切分方式示意图;
请参阅图9所示,为IEEE39***的切分示意图。
以IEEE39节点***为例,将其分为三区域电网区域及一个边界区域,具体切分方式见图4。支路9~8、3~4、14~15、16~17为联络线,节点3、4、8、9、14、15、16、17为边界节点,它们共同边界区域;三个区域电网区域包含节点数情况如表1所示,各区域的节点和支路参数如表2~7所示
表1IEEE39节点***区域情况
区域 | 1 | 2 | 3 |
节点数 | 15 | 12 | 12 |
表2区域电网区域1节点参数
节点号 | 电压幅值 | 电压相角 | 有功出力 | 无功出力 | 有功负荷 | 无功负荷 | 并联电导 | 并联电钠 | 节点类型 |
1 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
2 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
3 | 1 | 0 | 0 | 0 | 3.22 | 0.024 | 0 | 0.1107 | PQ |
9 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0.1902 | PQ |
17 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0.0671 | PQ |
18 | 1 | 0 | 0 | 0 | 1.58 | 0.3 | 0 | 0 | PQ |
25 | 1 | 0 | 0 | 0 | 2.24 | 0.472 | 0 | 0 | PQ |
26 | 1 | 0 | 0 | 0 | 1.39 | 0.17 | 0 | 0 | PQ |
27 | 1 | 0 | 0 | 0 | 2.81 | 0.755 | 0 | 0 | PQ |
28 | 1 | 0 | 0 | 0 | 2.06 | 0.276 | 0 | 0 | PQ |
29 | 1 | 0 | 0 | 0 | 2.835 | 0.269 | 0 | 0 | PQ |
30 | 1.047 | 0 | 2.5 | 0 | 0 | 0 | 0 | 0 | PV |
37 | 1.027 | 0 | 5.4 | 0 | 0 | 0 | 0 | 0 | PV |
38 | 1.026 | 0 | 8.3 | 0 | 0 | 0 | 0 | 0 | PV |
39 | 1.03 | 0 | 10 | 0 | 11.04 | 2.5 | 0 | 0 | Vθ |
表3区域电网区域1支路参数
首节点 | 尾节点 | 电阻 | 电抗 | 充电电容 | 变比 |
2 | 1 | 0.0035 | 0.0411 | 0.6987 | 1 |
39 | 1 | 0.001 | 0.025 | 0.75 | 1 |
3 | 2 | 0.0013 | 0.0151 | 0.2572 | 1 |
25 | 2 | 0.007 | 0.0086 | 0.146 | 1 |
18 | 3 | 0.0011 | 0.0133 | 0.2138 | 1 |
39 | 9 | 0.001 | 0.025 | 1.2 | 1 |
18 | 17 | 0.0007 | 0.0082 | 0.1319 | 1 |
27 | 17 | 0.0013 | 0.0173 | 0.3216 | 1 |
26 | 25 | 0.0032 | 0.0323 | 0.513 | 1 |
27 | 26 | 0.0014 | 0.0147 | 0.2396 | 1 |
28 | 26 | 0.0043 | 0.0474 | 0.7802 | 1 |
29 | 26 | 0.0057 | 0.0625 | 1.029 | 1 |
29 | 28 | 0.0014 | 0.0151 | 0.249 | 1 |
2 | 30 | 0 | 0.0181 | 0 | 1.025 |
25 | 37 | 0.0006 | 0.0232 | 0 | 1.025 |
29 | 38 | 0.0008 | 0.0156 | 0 | 1.025 |
表4区域电网区域2节点参数
节点号 | 电压幅值 | 电压相角 | 有功出力 | 无功出力 | 有功负荷 | 无功负荷 | 并联电导 | 并联电钠 | 节点类型 |
4 | 1 | 0 | 0 | 0 | 5 | 1.84 | 0 | 0.1107 | PQ |
5 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
6 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
7 | 1 | 0 | 0 | 0 | 2.338 | 0.84 | 0 | 0 | PQ |
8 | 1 | 0 | 0 | 0 | 5.22 | 1.76 | 0 | 0.1902 | PQ |
10 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
11 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
12 | 1 | 0 | 0 | 0 | 0.085 | 0.88 | 0 | 0 | PQ |
13 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
14 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0.183 | PQ |
31 | 0.982 | 0 | 5.72 | 0 | 0.092 | 0.046 | 0 | 0 | Vθ |
32 | 0.983 | 0 | 6.5 | 0 | 0 | 0 | 0 | 0 | PV |
表5区域电网区域2支路参数
首节点 | 尾节点 | 电阻 | 电抗 | 充电电容 | 变比 |
5 | 4 | 0.0008 | 0.0128 | 0.1342 | 1 |
14 | 4 | 0.0008 | 0.0129 | 0.1382 | 1 |
6 | 5 | 0.0002 | 0.0026 | 0.0434 | 1 |
8 | 5 | 0.0008 | 0.0112 | 0.1476 | 1 |
7 | 6 | 0.0006 | 0.0092 | 0.113 | 1 |
11 | 6 | 0.0007 | 0.0082 | 0.1389 | 1 |
8 | 7 | 0.0004 | 0.0046 | 0.078 | 1 |
11 | 10 | 0.0004 | 0.0043 | 0.0729 | 1 |
13 | 10 | 0.0004 | 0.0043 | 0.0729 | 1 |
14 | 13 | 0.0009 | 0.0101 | 0.1723 | 1 |
6 | 31 | 0 | 0.025 | 0 | 1.07 |
12 | 11 | 0.0016 | 0.0435 | 0 | 1.006 |
12 | 13 | 0.0016 | 0.0435 | 0 | 1.006 |
10 | 32 | 0 | 0.02 | 0 | 1.07 |
表6区域电网区域3节点参数
节点号 | 电压幅值 | 电压相角 | 有功出力 | 无功出力 | 有功负荷 | 无功负荷 | 并联电导 | 并联电钠 | 节点类型 |
15 | 1 | 0 | 0 | 0 | 3.2 | 1.53 | 0 | 0.183 | PQ |
16 | 1 | 0 | 0 | 0 | 3.29 | 0.323 | 0 | 0.0671 | PQ |
19 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
20 | 1 | 0 | 0 | 0 | 6.8 | 1.03 | 0 | 0 | PQ |
21 | 1 | 0 | 0 | 0 | 2.74 | 1.15 | 0 | 0 | PQ |
22 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | PQ |
23 | 1 | 0 | 0 | 0 | 2.475 | 0.846 | 0 | 0 | PQ |
24 | 1 | 0 | 0 | 0 | 3.086 | -0.922 | 0 | 0 | PQ |
33 | 0.997 | 0 | 6.32 | 0 | 0 | 0 | 0 | 0 | PV |
34 | 1.012 | 0 | 5.08 | 0 | 0 | 0 | 0 | 0 | PV |
35 | 1.049 | 0 | 6.5 | 0 | 0 | 0 | 0 | 0 | Vθ |
36 | 1.063 | 0 | 5.6 | 0 | 0 | 0 | 0 | 0 | PV |
表7区域电网区域3支路参数
首节点 | 尾节点 | 电阻 | 电抗 | 充电电容 | 变比 |
16 | 15 | 0.0009 | 0.0094 | 0.171 | 1 |
19 | 16 | 0.0016 | 0.0195 | 0.304 | 1 |
21 | 16 | 0.0008 | 0.0135 | 0.2548 | 1 |
24 | 16 | 0.0003 | 0.0059 | 0.068 | 1 |
19 | 20 | 0.0007 | 0.0138 | 0 | 1.06 |
22 | 21 | 0.0008 | 0.014 | 0.2565 | 1 |
23 | 22 | 0.0006 | 0.0096 | 0.1846 | 1 |
24 | 23 | 0.0022 | 0.035 | 0.361 | 1 |
19 | 33 | 0.0007 | 0.0142 | 0 | 1.07 |
20 | 34 | 0.0009 | 0.018 | 0 | 1.009 |
22 | 35 | 0 | 0.0143 | 0 | 1.025 |
23 | 36 | 0.0005 | 0.0272 | 0 | 1 |
基于上述互联电力***切分方式,由区域电网潮流计算服务器负责计算区域电网潮流,由上级调度中心内协调计算服务器负责更新并向各区域电网发送边界条件,以及判断全网一体化潮流是否收敛。
2.阶段2:全网一体化潮流分解协调计算初始化过程。该过程包括区域电网潮流计算服务器初始化和上级调度中心内协调计算服务器初始化两部分,以下分别介绍。
(2.1)区域电网潮流计算服务器按照以下步骤完成初始化;
(2.1.1)读入本区域潮流计算的节点和支路参数,例如表2和3中所示参数;
(2.1.2)初始化边界条件信息,包括边界节点编号和名称,以及边界条件与边界节点之间的对应关系等;
(2.1.3)初始化潮流计算控制参数,包括使用的迭代求解方法、最大迭代次数、节点电压上下限、收敛判据等;
(2.2)上级调度中心内协调计算服务器按照以下步骤完成初始化;
(2.2.1)读入边界区域网络参数,包括边界节点名称、联络线阻抗等;
(2.2.2)初始化区域电网承担全网有功网损比例参数∏;
(2.2.3)设定采用JFNG(m)求解边界协调方程的控制参数,包括最大Newton迭代次数max NIter、最大GMRES迭代次数maxGIter、Newton迭代收敛精度ξ(边界节点电流平衡和区域负担全网有功网损平衡的精度要求)、GMRES迭代收敛精度ε、有限差分步长w等。
3.阶段3:采用改进JFNG(m)方法求解边界协调方程,实现全网潮流一体化求解。
具体求解过程如下:
(3.1)令k=-1,令xB=[Re(IB),Im(IB),Θ]T为边界条件组成的向量,其中Re(IB)和Im(IB)分别表示边界节点来自联络线的注入电流相量实部和虚部组成的向量;给定初始边界条件xB 0,选择任意非奇异矩阵M0为预处理矩阵,例如单位阵;
(3.2)k=k+1,重复步骤(3.2)~(3.9), 或是k>max NIter,结束;
(3.3)进入GMRES(m)迭代求解修正方程: 其中G′(xB k)表示边界协调方程在xB k的一阶倒数,即边界协调方程的Jacobian矩阵,G(xB k)为边界条件为xB k时边界协调函数的值,为了计算需要将ΔI的实部和虚部分开,令
(3.4) l=1,ρ=β=‖r0‖2,v1=r0,errtolG=ε‖r0‖2>0,1>ε>0;
(3.5)若ρ>errtolG并且l<max GIter,则l=l+1,计算zl=Mkvl, vl+1=ΔGl/w,w为一个小的正常数,一般取10-5>w>0;
(3.6)修正预处理矩阵
(3.8)
(3.9)修正预处理矩阵 返回步骤(3.2);
其中边界协调函数G(xB)的求值过程按照以下步骤完成:
(3.10)上级调度中心内的协调计算服务器设定边界条件,即区域电网边界节点上来自联络线的注入电流相量组成的向量IB及各区域电网松弛节点相角设定值组成的向量Θ,并按照地域从属关系,将相应的注入电流和相角设定值发送给区域电网潮流计算服务器;
(3.11)区域电网潮流计算服务器启动本区域潮流计算过程,根据所得到的边界条件和内部节点参数求解S1和S2潮流方程(4),从收敛的潮流结果中得到偏差边界节点的电压VB和区域网损信息Ploss,并将这些信息发送给协调计算服务器;
(3.12)协调计算服务器根据得到的区域潮流计算结果信息以及初始化过程中设定的参数,按照式(5)和(6)计算边界节点的电流不平衡量ΔI,并按照式
(2)计算区域电网负担全网有功网损的不平衡量ΔP;
对本发明的计算机仿真结果如下:
采用IEEE39节点***进行仿真。网络参数及区域情况如表1~7所示。算例按照以下条件进行:
(1)各区域电网潮流计算服务器参数为:
(1.1)各区域电网潮流从平启动开始,即PQ节点电压幅值和相角分别为1和0,PV节点电压相角设为0,松弛(Vθ)节点相角为0;
(1.2)采用Newton-Raphson迭代方法求解潮流方程(3);
(1.3)最大迭代次数为30,节点电压幅值上下限为1.5和0.5,潮流收敛精度10-7;
(2)协调侧计算服务器参数为:
(2.1)设定区域电网负担全网有功网损的比例向量为∏=[0,1,0],即规定全网所有的网损均由区域电网区域2内的松弛节点承担,这样的设置和IEEE39节点标准数据的设置是一致的;
(2.2)JFNG(m)算法的控制参数为,ξ=10-4,ε=0.1,max NIter=30,maxGlter=30,w=10-6。
按照以上条件,进行对IEEE39节点***一体化潮流分解协调计算。
表8分布式潮流计算测试结果
*** | 集中NR法迭代次数 | 改进JFNG(m)法 | 节点电压幅值最大偏差 | 节点电压相角最大偏差 | ||
Newton迭代 | GMRES迭代 | 求值次数 | ||||
IEEE39 | 3 | 4 | 0+7+4+4 | 19 | 6.37×10-7 | 3.94×10-6 |
从表8可看出:采用JFNG(m)算法求解边界协调方程,其外层迭代次数与传统的串行NR法相当,具有较高的收敛性。自适应预处理技术的使用可以显著提高内层迭代的收敛性,从而减少协调方程求解的次数,这将大大提高算法在分布式计算环境中的通信次数,从而使得算法适用于高延时的广域网络通信环境。本发明给出的全网一体化潮流分解协调求解算法计算精度较高,能够满足实用计算需要。
最后应说明的是:以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此,尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围当中。
Claims (9)
1.一种通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于包括:
步骤1;根据电力***实际调度管理方式构建分布式潮流分析***;
所述的电力***由区域调度中心和上级调度中心协调一致、信息互联;
所述的分布式潮流计算***由区域调度中心内的潮流计算服务器和上级调度中心内的协调计算服务器通过网络连接构成;其中:潮流计算服务器,负责其所辖电网的潮流计算;而上级调度中心内的协调计算服务器,负责协调各区域电网潮流计算过程;
步骤2;从电力***拓扑连接关系出发,根据实际***的运行状况,采用带有边界区域的切分方法对***进行切分,按照实际网络连接情况,对互联电网进行划分,明确潮流计算服务器和协调计算服务器的计算对象和数据来源;
所述的***定义为(S0);包括两个区域***(S1)和(S2),通过联络线l相互连接;其中(B1)(B2)分别代表(S1)(S2)区域的边界节点,(Sln 1)(Sln 2)分别为(S1)(S2)中除了边界节点外其他节点组成的网络;
所述区域S1和S2属于区域电网调度中心管辖范围,区域SB属于上级调度中心管辖范围;在这种切分方式下,全网潮流的收敛条件为:区域S1和S2中潮流计算达到收敛,它们内部节点功率均达到平衡,且边界节点状态满足
其中uB和θB是区域电网内潮流计算收敛结果中边界节点B1和B2的电压幅值和相角向量;和是边界区域SB内潮流计算收敛时边界节点B1和B2的电压幅值和相角向量;是区域电网潮流计算收敛结果中边界节点B1和B2上来自联络线的注入电流相量组成的向量,以流入为正方向;是边界区域SB中潮流计算收敛时从边界节点B1和B2向联络线上注入的电流相量组成的向量,以流出为正方向。
步骤3;区域潮流计算服务器和上级协调计算服务器计算用数据和参数初始化,包括以下步骤:
步骤3.1:区域电网调度中心内的潮流计算服务器初始化所辖电网的网络参数,包括所有参与潮流计算的网络节点的类型、负荷和发电机出力设定值、发电机机端节点电压幅值以及松弛节点相角等;
步骤3.2:区域电网调度中心内的潮流计算服务器初始化边界条件信息,包括边界节点的编号、名称等、边界条件的组成,即边界节点的注入电流相量组成的向量和区域电网松弛节点相角;
步骤3.3:区域电网调度中心内的潮流计算服务器初始化潮流计算的参数,包括所使用的迭代求解方法、最大迭代次数、节点电压上下限、用以判断潮流收敛的最小功率偏差;
步骤3.4:上级调度中心内的协调计算服务器初始化边界区域网络参数,包括边界节点的名称、联络线阻抗;
步骤3.5:上级调度中心内的协调计算服务器设定全网一体化潮流协调求解控制参数;最大迭代次数、各区域负担全网有功网损的比例系数,以及边界节点功率平衡精度要求,即判定边界节点功率平衡所需的边界节点上不平衡的电流向量最大幅值。
步骤4:协调计算服务器调用区域潮流计算服务器共同完成全网一体化潮流分解协调求解过程;服务器间调用流程,包括以下基本步骤:
步骤4.1:上级调度中心内的协调计算服务器按照步骤3.4、3.5初始化本地基础网络数据和计算参数,区域调度中心内区域潮流计算服务器按照步骤3.1~3.3初始化本地基础网络数据和计算参数;
步骤4.2:协调计算服务器设定初始边界条件,即区域电网边界节点上来自联络线的注入电流相量组成的向量IB及各区域电网松弛节点相角设定值组成的向量Θ,并按照步骤2中网络化分关系,将相应的注入电流和相角设定值发送给区域电网潮流计算服务器;
步骤4.3:区域电网潮流计算服务器启动本区域潮流计算过程,根据所得到的边界条件和内部节点参数求解区域潮流方程,从收敛的潮流结果中得到偏差边界节点的电压VB i和区域网损信息Ploss i,i=1,2,并将信息发送给协调计算服务器;
步骤4.4:协调计算服务器根据得到的区域潮流计算结果信息以及初始化过程中设定的参数,计算边界节点的电流不平衡量ΔI,和区域电网负担全网有功网损的不平衡量ΔP:
步骤4.5:若‖[ΔI,ΔP]T‖2≤ξ则判断全网潮流收敛,计算过程结束,其中‖.‖2表示向量的二范数,ξ为小的正常数;否则根据[ΔI,ΔP]T计算边界条件修正量[ΔIB,ΔΘ]T后更新边界条件设定值,返回步骤4.2;
步骤5:区域电网潮流计算服务器将计算结果反馈给协调计算服务器,以使协调计算服务器获得收敛的全网一体化潮流计算结果。
2.根据权利要求1所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:根据计算边界节点的电流不平衡量和区域电网有功网损的不平衡量,判断全网一体化潮流是否一致收敛;需要同时满足以下基本条件:
A:边界节点功率达到平衡,即满足方程;
B:全网有功网损在各区域电网间合理分配,满足以下方程,
其中∏={ηi}是区域电网S1和S2负担全网有功网损的比例向量,i=1,2, Pi loss是从区域电网潮流计算结果中统计出的有功网损; 是由各区域电网有功网损相加得到的全网有功网损;ΔPloss i是区域电网负担全网有功网损的不平衡量,定义为区域电网有功网损不平衡量。
3.根据权利要求1所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:所述的边界条件包括:区域电网的边界母线上来自联络线的注入电流相量组成的向量IB以及区域电网松弛节点相角设定值Θ。
4.根据权利要求1所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:所述的区域潮流计算服务器根据给定的边界条件计算区域潮流;其中所述的区域电网S1:
设电网由n个节点和若干条支路组成,潮流计算模型可表述为非线性方程组形式,即
其中 i=1,2,…,n为节点电压相量,g为节点功率平衡方程,Pi SP和Qi SP是节点i的给定有功和无功功率,Gij+jBij是节点i和节点j之间的导纳,j∈i表示与i节点相联的所有节点,包括i节点自身。
5.根据权利要求1所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:S1区域所有节点由内部节点Sln 1和边界节点B1两部分,对应的以 代表内部节点的电压相量组成的向量,以 代表边界节点电压相量组成的向量;对应地,可将方程(2)描述的节点功率平衡方程组也可以划分为两类,内部节点Sln 1的功率平衡方程组gln 1和边界节点的功率平衡方程组gB 1;如区域电网计算潮流时不选择边界节点作为松弛节点,在完成给定内部节点的中发电机节点的有功出力、节点电压幅值,负荷节点的有功和无功负荷,以及松弛节点的电压幅值和相角情况下,内部节点功率平衡方程gln 1内部节点和边界节点电压相关,边界节点功率平衡方程gB 1与内部节点电压、边界节点电压及其来自联络线的注入电流相量相关,其表述为如下形式,
6.根据权利要求1所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:当完成步骤3.1中的基础数据初始化过程后,再给定IB 1和区域电网S1内松弛节点电压相角设定值θref 1即可采用普通的Newton-Raphson方法求解方程(4)式,得到区域电网S1内的潮流。
7.根据权利要求1所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:步骤4.3中协调计算服务器需要根据边界条件、联络线参数和区域电网潮流计算结果来计算边界节点的电流不平衡量ΔI,和区域电网负担全网有功网损的不平衡量ΔP;包括以下步骤:
a:协调计算服务器设定边界条件IB和Θ,发送至区域潮流计算服务器;
b:区域潮流计算服务器按照(4)式计算本区域电网潮流,从潮流结果中提取VB和Ploss。信息,其中VB是边界节点电压相量组成的向量,Ploss是区域电网有功网损信息,并将这些信息发送到协调计算服务器;
c:协调计算服务器根据边界节点的电压VB,按照(5)式计算边界节点来自联络线的注入电流理论值 以流出边界节点为正方向,
d协调计算服务器根据设定的区域电网负担全网有功网损的比例向量∏和各区域电网潮流结果中统计出的有功网损信息Ploss,根据(2)式计算区域电网有功网损不平衡量组成的向量
8.根据权利要求1所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:所述的协调计算服务器采用改进的Jacobian-free Newton-GMRES(m)方法计算边界条件的修正量,并迭代直至‖[ΔI,ΔP]T‖2≤ξ,具体步骤如下:
a.令k=-1,令xB=[Re(IB),Im(IB),Θ]T为边界条件组成的向量,其中Re(IB)和Im(IB)分别表示边界节点来自联络线的注入电流相量实部和虚部组成的向量;给定初始边界条件xB 0,选择任意非奇异矩阵M0为预处理矩阵;
b.k=k+1,重复步骤a~i,直至‖[ΔI,ΔP]T‖2≤ξ满足或是k>max NIter,max NIter为Newton迭代次数限制,结束;
c.令 为边界条件为xB k时通过步骤a~e计算得到的边界节点电流不平衡量和区域有功网损不平衡量,为了计算需要将ΔI的实部和虚部分开;
d. l=1,ρ=β=‖r0‖2,v1=r0,errtolG=ε‖r0‖2>0,1>ε>0;
e.若ρ>errtolG并且l<max GIter,max GIter是GMRES迭代次数限制,则l=l+1,zl=Mkvl, vl+1=ΔGl/w,w为一个小的正常数,一般取10-5>w>0;
f.修正预处理矩阵
h.
i.修正预处理矩阵 返回步骤b;
9.根据权利要求8所述的通过交换边界节点状态和网损信息完成分布式潮流分析的方法,其特征在于:其中步骤f和i为内层和外层预处理矩阵修正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101217928A CN101141064B (zh) | 2007-09-14 | 2007-09-14 | 通过交换边界节点状态和网损信息完成分布式潮流分析的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101217928A CN101141064B (zh) | 2007-09-14 | 2007-09-14 | 通过交换边界节点状态和网损信息完成分布式潮流分析的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101141064A true CN101141064A (zh) | 2008-03-12 |
CN101141064B CN101141064B (zh) | 2010-12-15 |
Family
ID=39192864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101217928A Active CN101141064B (zh) | 2007-09-14 | 2007-09-14 | 通过交换边界节点状态和网损信息完成分布式潮流分析的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101141064B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101937490A (zh) * | 2010-08-27 | 2011-01-05 | 河北省电力研究院 | 一种规划态电网网损分析计算方法 |
CN101540504B (zh) * | 2009-04-30 | 2012-01-11 | 东北大学 | 一种基于变步长神经网络的潮流分析装置及方法 |
CN102385654A (zh) * | 2011-10-25 | 2012-03-21 | 中国电力科学研究院 | 一种基于vsc并网光伏、储能电站的机电暂态建模方法 |
CN101777764B (zh) * | 2010-01-05 | 2013-01-09 | 华北电力大学 | 基于区域控制偏差和机组发电偏差的调度反馈控制*** |
CN103198205A (zh) * | 2013-02-25 | 2013-07-10 | 中国电力科学研究院 | 适应于大规模数据潮流计算的初值计算方法 |
CN103353921A (zh) * | 2013-06-18 | 2013-10-16 | 国家电网公司 | 一种基于启发式自动网络分区电力***并行潮流确定方法 |
CN103530736A (zh) * | 2013-10-23 | 2014-01-22 | 华北电力大学 | 基于web-svg的分布式潮流计算***图形平台的构建方法 |
CN104680430A (zh) * | 2015-03-19 | 2015-06-03 | 国网上海市电力公司 | 基于传输节点容量模型的多电压等级电网可靠性评估方法 |
CN104716646A (zh) * | 2015-03-31 | 2015-06-17 | 湖南大学 | 一种基于注入电流的节点耦合度分析方法 |
CN105406459A (zh) * | 2014-09-10 | 2016-03-16 | 三菱电机株式会社 | 利用基于一致性分布式处理估计电网中的最优潮流的方法 |
CN105576648A (zh) * | 2015-11-23 | 2016-05-11 | 中国电力科学研究院 | 一种基于gpu-cpu异构计算平台的静态安全分析双层并行方法 |
CN105811420A (zh) * | 2016-05-05 | 2016-07-27 | 国网江苏省电力公司南京供电公司 | 一种分解协调式主配网一体化潮流计算方法 |
CN106130005A (zh) * | 2016-05-13 | 2016-11-16 | 国家电网公司 | 基于地理分区的220kV电网分区边界节点处理方法 |
CN106201985A (zh) * | 2016-07-07 | 2016-12-07 | 三峡大学 | 一种基于pq法的分布式并行潮流计算***开发方法 |
CN107134769A (zh) * | 2016-02-26 | 2017-09-05 | 中国电力科学研究院 | 一种电网模型边界规划方法 |
CN107482624A (zh) * | 2017-08-11 | 2017-12-15 | 国网江苏省电力公司经济技术研究院 | 一种基于upfc的电网损耗控制策略 |
CN107994583A (zh) * | 2017-11-14 | 2018-05-04 | 华北电力大学 | 直流联络线在协调层建模的交直流***分布式潮流算法 |
CN108400601A (zh) * | 2018-01-25 | 2018-08-14 | 国电南瑞科技股份有限公司 | 一种特高压直流近区电网国分省协同无功备用计算方法 |
CN108647830A (zh) * | 2018-05-16 | 2018-10-12 | 清华大学 | 面向跨区跨省电力交易的省内安全约束边界等值生成方法 |
CN108847672A (zh) * | 2018-07-02 | 2018-11-20 | 国网山西省电力公司电力科学研究院 | 基于分布式连续潮流的最小电压稳定裕的获取方法和*** |
CN109193664A (zh) * | 2018-09-12 | 2019-01-11 | 国电南瑞科技股份有限公司 | 含dg并网的配电网潮流计算方法 |
CN111262249A (zh) * | 2019-12-20 | 2020-06-09 | 国网新疆电力有限公司 | 一种电网潮流仿真方法和*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100373733C (zh) * | 2004-01-06 | 2008-03-05 | 中国电力科学研究院 | 电力***潮流分网并行计算方法 |
CN100481668C (zh) * | 2006-09-22 | 2009-04-22 | 天津大学 | 大规模电力***网络分割潮流并行计算方法 |
-
2007
- 2007-09-14 CN CN2007101217928A patent/CN101141064B/zh active Active
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101540504B (zh) * | 2009-04-30 | 2012-01-11 | 东北大学 | 一种基于变步长神经网络的潮流分析装置及方法 |
CN101777764B (zh) * | 2010-01-05 | 2013-01-09 | 华北电力大学 | 基于区域控制偏差和机组发电偏差的调度反馈控制*** |
CN101937490A (zh) * | 2010-08-27 | 2011-01-05 | 河北省电力研究院 | 一种规划态电网网损分析计算方法 |
CN102385654B (zh) * | 2011-10-25 | 2015-09-23 | 中国电力科学研究院 | 一种基于vsc并网光伏、储能电站的机电暂态建模方法 |
CN102385654A (zh) * | 2011-10-25 | 2012-03-21 | 中国电力科学研究院 | 一种基于vsc并网光伏、储能电站的机电暂态建模方法 |
CN103198205A (zh) * | 2013-02-25 | 2013-07-10 | 中国电力科学研究院 | 适应于大规模数据潮流计算的初值计算方法 |
CN103198205B (zh) * | 2013-02-25 | 2017-02-08 | 中国电力科学研究院 | 适应于大规模数据潮流计算的初值计算方法 |
CN103353921A (zh) * | 2013-06-18 | 2013-10-16 | 国家电网公司 | 一种基于启发式自动网络分区电力***并行潮流确定方法 |
CN103530736A (zh) * | 2013-10-23 | 2014-01-22 | 华北电力大学 | 基于web-svg的分布式潮流计算***图形平台的构建方法 |
CN105406459A (zh) * | 2014-09-10 | 2016-03-16 | 三菱电机株式会社 | 利用基于一致性分布式处理估计电网中的最优潮流的方法 |
CN105406459B (zh) * | 2014-09-10 | 2018-03-16 | 三菱电机株式会社 | 利用基于一致性分布式处理估计电网中的最优潮流的方法 |
CN104680430A (zh) * | 2015-03-19 | 2015-06-03 | 国网上海市电力公司 | 基于传输节点容量模型的多电压等级电网可靠性评估方法 |
CN104680430B (zh) * | 2015-03-19 | 2018-04-27 | 国网上海市电力公司 | 基于传输节点容量模型的多电压等级电网可靠性评估方法 |
CN104716646A (zh) * | 2015-03-31 | 2015-06-17 | 湖南大学 | 一种基于注入电流的节点耦合度分析方法 |
CN104716646B (zh) * | 2015-03-31 | 2016-02-24 | 湖南大学 | 一种基于注入电流的节点耦合度分析方法 |
CN105576648A (zh) * | 2015-11-23 | 2016-05-11 | 中国电力科学研究院 | 一种基于gpu-cpu异构计算平台的静态安全分析双层并行方法 |
CN107134769A (zh) * | 2016-02-26 | 2017-09-05 | 中国电力科学研究院 | 一种电网模型边界规划方法 |
CN107134769B (zh) * | 2016-02-26 | 2020-08-18 | 中国电力科学研究院 | 一种电网模型边界规划方法 |
CN105811420A (zh) * | 2016-05-05 | 2016-07-27 | 国网江苏省电力公司南京供电公司 | 一种分解协调式主配网一体化潮流计算方法 |
CN106130005B (zh) * | 2016-05-13 | 2021-09-28 | 国家电网公司 | 基于地理分区的220kV电网分区边界节点处理方法 |
CN106130005A (zh) * | 2016-05-13 | 2016-11-16 | 国家电网公司 | 基于地理分区的220kV电网分区边界节点处理方法 |
CN106201985A (zh) * | 2016-07-07 | 2016-12-07 | 三峡大学 | 一种基于pq法的分布式并行潮流计算***开发方法 |
CN106201985B (zh) * | 2016-07-07 | 2019-09-24 | 三峡大学 | 一种基于pq法的分布式并行潮流计算***开发方法 |
CN107482624B (zh) * | 2017-08-11 | 2019-09-13 | 国网江苏省电力公司经济技术研究院 | 一种基于upfc的电网损耗控制策略 |
CN107482624A (zh) * | 2017-08-11 | 2017-12-15 | 国网江苏省电力公司经济技术研究院 | 一种基于upfc的电网损耗控制策略 |
CN107994583B (zh) * | 2017-11-14 | 2020-05-08 | 华北电力大学 | 直流联络线在协调层建模的交直流***分布式潮流算法 |
CN107994583A (zh) * | 2017-11-14 | 2018-05-04 | 华北电力大学 | 直流联络线在协调层建模的交直流***分布式潮流算法 |
CN108400601A (zh) * | 2018-01-25 | 2018-08-14 | 国电南瑞科技股份有限公司 | 一种特高压直流近区电网国分省协同无功备用计算方法 |
CN108647830A (zh) * | 2018-05-16 | 2018-10-12 | 清华大学 | 面向跨区跨省电力交易的省内安全约束边界等值生成方法 |
CN108647830B (zh) * | 2018-05-16 | 2021-10-15 | 清华大学 | 面向跨区跨省电力交易的省内安全约束边界等值生成方法 |
CN108847672A (zh) * | 2018-07-02 | 2018-11-20 | 国网山西省电力公司电力科学研究院 | 基于分布式连续潮流的最小电压稳定裕的获取方法和*** |
CN109193664A (zh) * | 2018-09-12 | 2019-01-11 | 国电南瑞科技股份有限公司 | 含dg并网的配电网潮流计算方法 |
WO2020052104A1 (zh) * | 2018-09-12 | 2020-03-19 | 国电南瑞科技股份有限公司 | 含dg并网的配电网潮流确定方法及计算机存储介质 |
CN111262249A (zh) * | 2019-12-20 | 2020-06-09 | 国网新疆电力有限公司 | 一种电网潮流仿真方法和*** |
CN111262249B (zh) * | 2019-12-20 | 2024-01-30 | 国网新疆电力有限公司 | 一种电网潮流仿真方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN101141064B (zh) | 2010-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101141064A (zh) | 通过交换边界节点状态和网损信息完成分布式潮流分析的方法 | |
CN106549384B (zh) | 一种含upfc电力***的通用潮流计算方法 | |
CN106208102B (zh) | 一种基于辅助问题原理的主动配电网分布式无功优化方法 | |
CN106250640B (zh) | 一种适用于地区电网的分层动态等值方法 | |
CN103094900B (zh) | 考虑相序混合法的分布式发电配电网三相潮流计算方法 | |
CN111428351B (zh) | 基于前推回代法的电-热综合能源***潮流计算方法 | |
CN109802392B (zh) | 大规模配电网潮流计算方法及装置 | |
Zhang et al. | A distributed multi-control-center dynamic power flow algorithm based on asynchronous iteration scheme | |
CN108388688A (zh) | 一种主动配电网可靠性评估方法 | |
CN105184418B (zh) | 基于灵敏度和约束等值的互联电网最优潮流计算方法 | |
CN104113061A (zh) | 一种含分布式电源的配电网三相潮流计算方法 | |
CN102522744B (zh) | 一种过网网损分摊方法 | |
CN110120673A (zh) | 基于戴维南等值参数辨识的分布式输配协同无功优化方法及*** | |
CN109768543B (zh) | 一种基于混合整数线性规划的弹性保底网架搜索建模方法 | |
CN108390390A (zh) | 一种含优先消纳风电能源路由器的三相不平衡主动配电网随机模糊潮流方法 | |
CN112234598B (zh) | 一种电磁暂态仿真初始化方法 | |
CN106991229B (zh) | 一种针对复杂拓扑的风电场等值建模方法 | |
CN108599974A (zh) | 一种基于图论连通性的微电网分布式通讯拓扑设计方法 | |
CN111834996B (zh) | 电网线损计算方法和装置 | |
CN105787604A (zh) | 一种输配协同的配网态势快速感知方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN105305502B (zh) | 基于协调一致性的配电网络分布式电源控制方法及*** | |
CN107465195B (zh) | 一种基于微电网结合潮流计算的最优潮流双层迭代方法 | |
Zhao et al. | An asynchronous iterative distributed power flow algorithm based on the boundary bus injection power | |
CN109687431A (zh) | 一种考虑新能源随机性的主动配电网概率等值建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |