CN102999788B - 基于多层图分割的电力系统主动解列断面快速搜索方法 - Google Patents
基于多层图分割的电力系统主动解列断面快速搜索方法 Download PDFInfo
- Publication number
- CN102999788B CN102999788B CN201210442714.9A CN201210442714A CN102999788B CN 102999788 B CN102999788 B CN 102999788B CN 201210442714 A CN201210442714 A CN 201210442714A CN 102999788 B CN102999788 B CN 102999788B
- Authority
- CN
- China
- Prior art keywords
- node
- alligatoring
- power
- graph
- pqk
- 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
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明属于电力系统稳定与控制领域。为实现解列断面的快速搜索,本发明采取的技术方案是,基于多层图分割的电力系统主动解列断面快速搜索方法,包括如下步骤:第一阶段:对原始图G0中的点或边进行合并,得到简化图G1,重复该步骤直到最终简化图Gn只包含少量的节点;第二阶段:对简化图Gn进行初始分割,得到初始解列面;第三阶段:将简化图Gn依次地还原到原始图G0,并在每一步还原过程中通过优化手段改善分区结果。第一阶段按照粗化对系统最优解的解空间的影响程度将粗化过程分为等效粗化、基于潮流追踪的粗化和基于重权边匹配的粗化三个子步骤。本发明主要应用于电力系统稳定与控制。
Description
技术领域
本发明属于电力系统稳定与控制领域,具体讲,基于多层图分割的电力系统主动解列断面快速搜索方法。
背景技术
失步解列作为保证电力系统安全运行的重要措施,是保证整个电网不致完全崩溃的最后一道防线。传统的失步解列根据预想场景对解列装置进行离线整定,然后在线基于本地量测和判据进行解列操作,因此,该方法存在着整定复杂度高、多套装置相互配合困难以及实际条件与预想场景差异大等缺点[1]。而随着通讯和监控手段的发展,失步解列的研究更加倾向于具有集中和协调特征的在线主动解列控制[2]。
主动解列需要快速地确定系统最优的解列断面,使得解列后的系统孤岛能够保持同步运行且功率供求基本平衡。针对该问题已有诸多学者从不同角度提出了快速搜索方法,如[4-11]。于电力系统庞大的节点数,对原始系统进行简化是解列断面快速搜索算法的基础:[4][5][6]基于图的等效原理对系统冗余节点进行了合并;在此基础上,[6]和[7]分别基于电气距离计算和潮流追踪技术确定发电机和紧密负荷的依附关系从而可对系统进行大幅简化。[8]提出一种基于弱联接理论的决策空间预筛方法以有效降低预决策的空间的规模。[9]提出一种基于调度分区的电力系统解列割集搜索算法。在简化图上实现解列断面搜索:[4][8]利用有序二元决策图来求解满足机组同步和平衡约束的系统划分;[10]基于最大流最小割来求解该问题;[11]提出了一种基于含连通图约束的背包问题的图分割求解方法。
上述方法的应用往往受掣于系统的规模,而由于系统解列问题从本质上是图论中的平衡图分割的问题,因此图论中大型图划分方法的发展对本问题的研究具有重要的参考意义。
参考文献/著作:
[1]高鹏,王建全,甘德强,等.电力系统失步解列综述[J],电力系统自动化,2005,29(19):90-96.
[2]沈沉,吴佳耘,乔颖,等.电力系统主动解列控制方法的研究[J],中国电机工程学报,2006(13):1-6.
[3]GuangyueXu,VijayVittal,etal.ControlledIslandingDemonstrationsontheWECCSystem[J].IEEETransonPowerSystems,2011,26(1):334-343.
[4]KaiSun,ZhengDa-Zhong,andLuQiang.Splittingstrategiesforislandingoperationoflarge-scalepowersystemsusingOBDD-basedmethods[J].IEEETransonPowerSystems,2003,18(2):912-923.
[5]6uangyueXu,andVijayVittal.SlowCoherencyBasedCutsetDeterminationAlgorithmforLargePowerSystems[J].IEEETransonPowerSystems,2010,25(2):877-884.
[6]吴学娟,沈沉,向学军,等.主动解列策略求解过程中的网络化简.中国电机工程学报[J],2008(7):7-12.
[7]C.G.Wang,B.H.Zhang,etal.ANovelReal-TimeSearchingMethodforPowerSystemSplittingBoundary[J].IEEETransonPowerSystems,2010.,25(4):1902-1909.
[8]乔颖,沈沉,卢强.大电网解列决策空间筛选及快速搜索方法[J].中国电机工程学报,2008(22):23-28.
[9]刘源祺,刘玉田.基于调度分区的电力系统解列割集搜索算法[J].电力系统自动化,2008,32(11):20-24.
[10]XiaomingWang,andVijayVittal.Systemislandingusingminimalcutsetswithminimumnetflow[C].IEEEPESPowerSystemsConference&Exposition,379-384,2004.
[11]林济铿,王旭东,李胜文,等.基于含连通图约束的背包问题的图分割方法[J].中国电机工程学报,2012(10):134-141.
[12]RomualdBrillout.Amultilevelframeworkforbisectionheuristics.KarlsruheInstituteofTechnology,December,2009.
发明内容
本发明旨在克服现有技术的不足,实现解列断面的快速搜索,本发明采取的技术方案是,
基于多层图分割的电力系统主动解列断面快速搜索方法,包括如下步骤:
第一阶段:对原始图G0中的点或边进行合并,得到简化图G1,重复该步骤直到最终简化图Gn只包含少量的节点;
第二阶段:对简化图Gn进行初始分割,得到初始解列面;
第三阶段:将简化图Gn依次地还原到原始图G0,并在每一步还原过程中通过优化手段改善分区结果。
第一阶段按照粗化对系统最优解的解空间的影响程度将粗化过程分为等效粗化、基于潮流追踪的粗化和基于重权边匹配的粗化三个子步骤。
等效粗化:
a、对多回并行线路进行合并;
b、将节点度为1的冗余节点归并至邻接点;
c、将权重为0的冗余节点归并至任一邻接点;
重复步骤a、b、c直至系统中不再存在上述冗余支路或节点,进入下一步基于潮流追踪的粗化过程。
基于潮流追踪的粗化进一步细化为:用公式[Vjk,Pjk]=PFT_BFS(vi,Φ,k)表示过程:以节点vi为起点,在节点集合Φ中进行潮流追踪,在追踪深度为k的情况下,可追踪到的节点集合为Vik,对应追踪到的功率为Pik,其中i为节点编号,k为追踪深度;
1)将发电机根据同调性划分为n个集合,置i=1;
2)将第i组同调机群中的所有发电节点加入集合F={vFj},j=1:m其中m为F中节点的个数,vFj为发电机编号为j的节点;置k=1;
3)对集合{vFj},j=1:m中全部节点并行地执行过程[Vjk,Pjk]=PFT_BFS(vFj,Φ,k)(1),其中Pjk为追踪到的节点j的功率,式(1)意为通过广度优先的方法,以节点vFj为起点,在节点集合Φ中进行潮流追踪,在追踪深度为k的情况下,可追踪到的节点集合为Vjk,对应追踪到的功率为Pjk若所有过程均搜索结束,其中Pjk为节点编号为i的节点功率,跳至步骤5;否则,判断集合Vjk之间有无交集:若无,k=k+1,重复步骤3;若有,进入步骤4;
4)假设其中Vpk、Vqk为以vpk,vqk为追踪起点追踪到的节点集合,p、q为节点编号,为空集,令z=k,k为追踪深度,分别继续执行[Vpz,Ppz]=PFT_BFS(vFp,Vpk∪Vqk,z)和[Vqz,Pqz]=PFT_BFS(vFq,Vpk∪Vqk,z)直至追踪完毕,此时,z=zc;令Vpqk=Vpzc∪Vqzc,其中Vpzc、Vqzc为追踪起点分别为节点编号为p、q的节点,追踪深度为zc的节点集合,Vpqk为集合Vpzc与Vqzc的并集,并将Vpzc∩Vqzc中的节点重复追踪到的功率相累加,得到Vpqk对应的追踪功率Ppqk,Ppqk为集合Vpqk中相应节点的功率;进一步,原过程[Vpk,Ppk]=PFT_BFS(vFp,Φ,k)和[Vqk,Pqk]=PFT_BFS(vFq,Φ,k)将被[Vpqk,Ppqk]=PFT_BFS(Vpqk,Φ,k)统一取代,k=k+1;返回步骤3;
5)i=i+1,若i≤n,返回步骤2;否则所有同步机群的搜索已结束;
通过潮流追踪之后,负荷节点vi归属发电节点vj的隶属度可定义为sij=Pij/Pi,即由发电节点vj供给负荷节点vi的功率Pij占节点vi总负荷Pi的比例,通过设置隶属度的阈值可将负荷与对应的发电节点进行归并,从而围绕多个同步机群可形成多个广义节点。
基于重权边匹配的粗化:对节点通过重权边匹配(HeavyEdgeMatching,HEM)的方法加以合并。
第二阶段:对Gn进行初始分割,得到初始解列面,具体为,选择图生长法的一种改进形式贪婪图生长法(GreedyGraphGrowingPartitioning,GGGP)作为初始分区算法,即以N个节点为起点以BFS为基础生长出N片区域直至覆盖整个图,而在生长过程中优先加入使分区平衡度得到显著改善的节点。
针对任意系统的一个边界节点vi用集合Sm(m=1,...,Ns)表示与原分区相连的第m条通路所包含的所有节点,Ns为通路总数,其中Sm中的节点又可以分为两类,优化前与vi同属一个分区的记为Sm1,否则记为Sm2用集合Rn,n=1,...,NR,表示与交换分区相连的第n条通路所包含的所有节点,其中优化前与vi同属一个分区的记为Rm1,否则记为Rm2,则对于vi,同时考虑其保留在原分区和交换至其他分区的连通性约束条件,留在原分区的节点需要将通路上与其连接的节点一并进行约束,以此保证其连通性,同理,对于交换至其他分区的节点,也需将通路上与其连接的节点一并进行约束,按此求解可以保证优化解的拓扑连通性。
本发明的技术特点及效果:
本发明通过粗化、初始分区和还原优化三个过程实现了最优解列断面的确定,从而摆脱了传统方法在处理大规模系统的限制。基于潮流追踪的粗化手段可大幅提高粗化效率,针对该问题,通过基于同步机组群的优化追踪方法可将传统的潮流追踪效率提高70%以上。与前人的工作相比,本文通过在还原优化中考虑拓扑结构约束可在改善子分区不平衡度时保证各自的连通性,避免了解列断面搜索结束后的额外合并工作。
附图说明
图1为多层图分割理论示意图。
图2为新英格兰39节点系统局部示意图,描述于同步机组群的优化追踪算法。
图3为传统解列方法的局限性示例图。
图4为场景一下IEEE-118节点分区图。
图5为场景一下粗化-还原过程中节点数与不平衡度的变化示意图。
具体实施方式
运用图论中的平衡图分割的方法来寻找电网解列面对于本课题有重要的意义。平衡图的划分方法主要包括谱分析方法、几何图方法和多层图分割方法等。其中,谱分析方法计算量较大(需要计算特征值和特征向量),而基于几何图分割的方法需要坐标信息辅助,因此上述两种方法的使用受到极大的限制。本发明以多层图分割理论为框架提出了一种电力系统主动解列的断面搜索方法,该方法结合了图论和电力系统的自身特点,通过粗化、初始分区和还原优化三步实现解列断面的快速搜索。在粗化过程中,针对失步解列的特殊问题,提出一种基于同步机组群的潮流追踪方法,提高了粗化过程的效率。使用贪婪图生长算法完成初始分区,并在还原优化过程中提出一种考虑拓扑结构约束的0-1规划方法保证较高的划分质量和子系统的连通性。
本发明以多层图分割理论作为电力系统主动解列断面搜索方法的主框架。在粗化、初始分区和还原优化三个过程中通过结合电力系统的实际特点提出相应的算法改进以实现解列断面的快速搜索。
基于多层图分割理论的电力系统主动解列断面快速搜索方法的主要阶段:
第一阶段:对原始图G0中的点或边进行合并,得到简化图G1,可重复该步骤直到最终简化图Gn只包含少量的节点。
第二阶段:对Gn进行初始分割,得到初始解列面。
第三阶段:将Gn依次地还原到G0,并在每一步还原过程中通过优化手段改善分区结果。
关键技术一:系统的粗化。随着系统规模的增加,解列原始策略空间呈几何级数急剧增长,但实际可行解空间的规模却始终维持在较小的数量级上。因此在图的分割之前基于图论和电力系统运行特点对系统进行粗化处理是十分必要的。本发明按照粗化对系统最优解的解空间的影响程度将粗化过程分为等效粗化、基于潮流追踪的粗化和基于重权边匹配的粗化三步。
1.等效粗化:等效粗化基于等效图原理将冗余节点或支路进行归并以减小系统规模,粗化过程并不会缩减最优解的解空间。简要总结如下:
a、对多回并行线路进行合并;
b、将节点度为1的冗余节点归并至邻接点;
c、将权重为0的冗余节点归并至任一邻接点。
重复步骤a、b、c直至系统中不再存在上述冗余支路或节点,可进入下一步粗化。
2.基于潮流追踪的粗化:对于图的进一步粗化,很多文献倾向于这样一种逻辑:即紧急情况,机组应该保证就近重要负荷的供电,因此在解列计算前可将发电机和与其有紧密电气联系的负荷预先归并。例如分别利用计算电气距离和潮流追踪的方法获得发电机和负荷的依附关系。其中,由于潮流追踪可以明确获得发电节点与负荷的电气联系,可作为本发明对系统进行进一步粗化的基础。由于系统在解列后同调的机组仍需保持对就近负荷的供电,所以该粗化过程不会明显减小系统合理解的解空间。
传统的潮流追踪技术可以确定某一电源到某一负荷点的有功功率大小,但主动解列问题中的潮流追踪存在以下两个特点:
a.同调的机组需要划分到同一分区中,因此不必区分某一负荷是由同调机组中的哪一台供给的;
b.同调的机组在分布上往往很接近,因此与其有紧密电气联系的负荷集合可能较为相似,倘若单独遍历,重复率很高。
基于以上两点考虑,本发明提出了一种基于同步机组群的优化追踪算法,并通过图2对算法的基本思想进行描述:
图2取自新英格兰39节点系统的局部,本发明假设发电机节点v35和发电机节点v36为同调的信息已给出。传统追踪方法基于广度优先搜索(BreadthFirstSearch,BFS)分别以节点v35和节点v36为起点对其供电负荷进行追踪,其追踪路径分别如a和b所示:
a.v35→v22→{v25,v21}→{v24,v16}→{v15,v17}→v18
b.v36→v25→v24→v16→{v15,v17}→v18
追踪路径a表示从节点v35开始按照以上节点顺序进行追踪,最后追踪到节点v18,追踪路径b表示从节点v36开始按照以上节点顺序进行追踪,最后追踪到节点v18,从节点通过路径a与b的对比可以发现:在搜索进行到第二步时,两条路径均已搜索到公共点v25,此时若将{v35,v36,v22,v25,v21,v24}视为一个广义节点,则由节点v35和节点v36供给的功率将均通过该广义节点流出,以此为基础继续进行潮流追踪则可避免重复遍历剩余的负荷节点。
传统潮流追踪的算法在此不再赘述,仅通过式1描述潮流追踪的过程:
[Vik,Pik]=PFT_BFS(vi,Φ,k)(1)
式(1)意为通过广度优先的方法,以节点vi为起点,在节点集合Φ中进行潮流追踪,在追踪深度为k的情况下,可追踪到的节点集合为Vik,对应追踪到的功率为Pik。当Vik随着k的增加不再增大时,追踪过程结束。以此为基础,本发明改进算法的描述如下:
1.将发电机根据同调性划分为n个集合,置i=1;
2.将第i组同调机群中的所有发电节点加入集合F={vFj},j=1:m其中m为F中节点的个数,vFj为发电机编号为j的节点;置k=1;
3.对集合{vFj},j=1:m中全部节点并行地执行过程[Vjk,Pjk]=PFT_BFS(vFj,Φ,k)(1),其中Pjk为追踪到的节点j的功率,式(1)意为通过广度优先的方法,以节点vFj为起点,在节点集合Φ中进行潮流追踪,在追踪深度为k的情况下,可追踪到的节点集合为Vjk,对应追踪到的功率为Pjk若所有过程均搜索结束,其中Pjk为节点编号为i的节点功率,跳至步骤5;否则,判断集合Vjk之间有无交集:若无,k=k+1,重复步骤3;若有,进入步骤4;
4.假设其中Vpk、Vqk为以vpk,vqk为追踪起点追踪到的节点集合,p、q为节点编号,为空集,令z=k,k为追踪深度,分别继续执行[Vpz,Ppz]=PFT_BFS(vFp,Vpk∪Vqk,z)和[Vqz,Pqz]=PFT_BFS(vFq,Vpk∪Vqk,z)直至追踪完毕,此时,z=zc;令Vpqk=Vpzc∪Vqzc,其中Vpzc、Vqzc为追踪起点分别为节点编号为p、q的节点,追踪深度为zc的节点集合,Vpqk为集合Vpzc与Vqzc的并集,并将Vpzc∩Vqzc中的节点重复追踪到的功率相累加,得到Vpqk对应的追踪功率Ppqk,Ppqk为集合Vpqk中相应节点的功率;进一步,原过程[Vpk,Ppk]=PFT_BFS(vFp,Φ,k)和[Vqk,Pqk]=PFT_BFS(vFq,Φ,k)将被[Vpqk,Ppqk]=PFT_BFS(Vpqk,Φ,k)统一取代,k=k+1;返回步骤3;
5.i=i+1,若i≤n,返回步骤2;否则所有同步机群的搜索已结束;
通过潮流追踪之后,负荷节点vi归属发电节点vj的隶属度可定义为sij=Pij/Pi,即由发电节点vj供给负荷节点vi的功率Pij占vi总负荷Pi的比例。通过设置隶属度的阈值可将负荷与对应的发电节点进行归并,从而围绕多个同步机群可形成多个广义节点。
3.基于重权边匹配的粗化:通过潮流追踪进行粗化后,系统规模往往可以大幅减小,若系统仍具有一定的规模,则需要对节点通过重权边匹配(HeavyEdgeMatching,HEM)的方法加以合并。HEM方法的启发式逻辑为通过将权重较大的边的两端节点进行合并以降低图的总边权,从而使割集边的权重总和更小。该归并方法对电力系统解列断面优化也具有明确的意义,即潮流非常重的支路往往并非解列的优解,因此将其两端节点合并不会对最优解产生较大的影响。该步骤可重复进行直至最终简化图小于进行初始分区的规模设定。
关键技术二:粗化图的初始分区。图的初始划分可采用多种方法:目前主要包括谱分析法、几何分割法和图生长法。其中,谱分析法需要计算第二小特征值的特征向量,几何分割法需要坐标信息辅助。对本发明而言,前两种方法或费时或有较多限制,且无法保证子图的连通性。因此,本发明选择图生长法的一种改进形式贪婪图生长法(GreedyGraphGrowingPartitioning,GGGP)作为初始分区算法,即以N个节点为起点以BFS为基础生长出N片区域直至覆盖整个图,而在生长过程中优先加入使分区平衡度得到显著改善的节点。该方法适合于本文问题主要基于以下两点考虑:
1)基于BFS的GGGP可保证生长子图的连通性。
2)GGGP对初始点的选择比较敏感,但电力系统的特点为以发电节点为核心供一片负荷区域。因此,以粗化后含同步机群的广义节点作为初始点可回避初始点选择难的问题。
通过GGGP可将最终简化系统划分为具有较好平衡度的为N个连通子图。而分区的平衡度可通过还原优化过程进一步改善。
关键技术三:粗化图的还原。在还原过程中,通过解散合并节点来提高图的自由度可继续改善解列断面的优化解。多层图分割理论中主要通过Kernighan-Lin(KL)算法、边界KL算法和FM算法来完成还原优化。上述算法存在以下两点缺陷:1)每次只允许移动1个或着交换2个节点,易陷入局部最优解;2)无法保证节点移动后的连通性。下面通过图3加以说明:
以图3(a)为例,割集1将系统初始划分为G1和G2左右两个子图,其不平衡功率分别为-0.5和0.5。边界节点集合为{v3,v4}和{v7,v8,v9},若从中移动1个节点或两边交换一对节点,均无法改善不匹配度;但实际上,将节点v4和{v8,v9}交换即是最优解,如割集2所示。边界的优化尚可通过进一步的还原进行改善,与之相比,子图的连通性则是一种更强的约束。例如,图3(b)将图3(a)中节点v3和v4的权重加以交换。若仅考虑不平衡度,通过两边交换节点v3和v7可使G1、G2的不平衡度降至0(如割集3所示),但割集3将G2分割为非连通的两部分,这是解列断面搜索所不能接受的解。基于上述两点原因,本文提出一种考虑拓扑结构的0-1规划方法进行还原过程中的解列断面优化。
(1)当不考虑拓扑结构时,该优化问题可描述为:
满足以下条件:
i=1,2…,N(2b)
其中,N为分区个数,Ii为i分区内的不平衡度,w(vm)为节点编号为m的节点权重值,xm为bool变量,xm=1表示节点保留在原分区内,xm=0表示通过割集支路交换到相连接分区。集合Ai表示其他分区内可通过割集与分区i相连的点的集合;Bi和Ci分别表示分区i内与其他分区有、无直接连接的点的集合。因此,Ii的变化与Ai和Bi中节点的交换情况有关,如式(2b)所示。但为保证节点交换后分区本身的连通性,需进一步考虑网络的拓扑结构。
(2)当考虑拓扑结构时:
在问题2的求解过程中,连通性的破坏是由于边界点不加约束、随意移动导致的。因此,对于任意边界点,其移动结果必须保证两个条件:1)Bi中边界点的移动不会导致Ci连通性的破坏;2)边界点不论保留在原分区还是交换至其它分区,均需保持与目标分区的连通性,即至少存在一支可行通路与目标分区相连。对条件1而言,若边界点不属于Ci的点割则可自然满足,否则,需将该边界点锁定防止移动。条件2可通过在优化问题中补充连通性约束加以保证。
首先将分区i的边界点集合Bi分为Tia和Tib两类:Tia表示Bi中与Ci相连的子集;Tib表示Bi/Tia。对两类节点施加约束条件的情况有所不同:当节点留在原分区时,Tia中的节点与原分区一定是连通的,无需添加约束,而Tib中的节点与原分区的连通性则需视Tia中节点选择确定;而当节点交换至其它分区时,无论Tia或Tib中的节点均需考虑与目标分区的连通性补充约束条件。
约束条件的补充方法以图3加以说明,由于解列断面优化时仅涉及边界点间的拓扑关系,因此将其单独表示为图3(c),其中,节点v3、v4∈T1a,节点v7∈T2a,及节点v8、v9∈T2b。以其中拓扑关系最为复杂的v9为例加以说明:
1.若节点v9留在G2
通过网络拓扑可得节点v9与区域G2除边界点外直接相连的路径有两条:9→8→7和9→3→8→7。用变量表示1-xi,xi为bool变量,xi=1表示节点保留在原分区内,xi=0表示
通过割集支路交换到相连接分区。若v9和G2之间不存在任何通路可表示为式(3):
此时v9不存在保留在G2的可能性,x9=0;而至少存在一条通路可表示为式(4):
此时v9={0,1}。综上,保证v9留在G2且不会形成孤岛的约束如式(5)所示:
2.若v9交换至G1
v9与G1之间的通路为:9→3和9→8→7→4两条。与情况a的分析类似,为保证v9交换到G1后仍保持连通性的约束条件为式(6):
通过上述的分析,针对任意系统的一个边界节点vi用集合Sm(m=1,…,NS)表示与原分区相连的第m条通路所包含的所有节点,NS为通路总数。其中Sm中的节点又可以分为两类,优化前与vi同属一个分区的记为Sm1,否则记为Sm2用集合Rn(n=1,…,NR)表示与交换分区相连的第n条通路所包含的所有节点,其中优化前与vi同属一个分区的记为Rm1,否则记为Rm2。则对于vi,同时考虑其保留在原分区和交换至其他分区的连通性约束条件如式(7)所示:
上式中vp、vq为节点编号分别为p、q的节点,用变量表示1-xp,xp为bool变量,xp=1表示节点编号为p的节点保留在原分区内,xp=0表示通过割集支路交换到相连接分区,xq同理,通过将式(7)加入问题(3)求解即可保证优化解的拓扑连通性。以图3(c)为例,该优化问题共有变量5个,需添加约束7个。而通过对大量算例的研究表明,实际系统网架结构的稀疏性使得与割集关联的边界节点并不多,而且边界点中出现Tib类型节点的概率也比较小,因此,系统的优化变量和约束条件的个数可以保证在一个较小的规模内,并不会随着系统的增大而有明显增大。
下面通过IEEE-118节点算例(如图4所示)来给出本发明的一个最佳实施方式。
系统的总发电量为3800.36MW,总负荷量为3650MW,其中,线路损耗可按照潮流追踪的结果分摊到负荷上,从而将系统转换为无损图。根据文献[4]的发电机分组结果,本文将发电机分为{10,12,25,26,31}、{49,54,59,61,65,66,69,80}和{87,89,100,103,111}共3组。本算例设置为第1组的发电机对第2组和第3组的发电机失步。并对本算例进行解列断面搜索。
如表1所示,我们将原始118系统定义为图G0,通过逐步的粗化过程得到G4,对G4进行初始分区并将分区结果定义为G5。由图5可知,基于潮流追踪的粗化可大幅减少系统节点数目,提高粗化效率。本算例中,G3和G4图由两次基于HEM的粗化生成,以达到初始分区的规模要求。而在还原优化过程中,如表1和图5所示,初始分区G5的不平衡度为147.46,G5到G4的第一次优化将广义节点{v33,v35,v37,v39}和节点v24移动至左侧分区,不平衡度可降至28.07;由G4到G3的还原过程中,广义节点{v33,v35,v37,v39}解散为{v33,v35,v37}和v39,优化过程通过将v19和v39交换将不平衡度降至0.5659;而剩余的还原优化过程已不能进一步优化该结果,最终解列断面如图4所示。下图中{vi,vm,vn}-vk表示由节点vi,vm和vn组成的广义节点{vi,vm,vn}与节点vk相连接的边;
表1场景一的解列分区结果
Tab.1thepartitionresultofscenario1
本发明的测试环境使用计算机的配置为IntelDual2.0GHz,2G内存;使用C++搭建主平台并编制相关程序。在计算时间上:潮流追踪部分,传统方法平均耗时约为4ms,而本文优化方法约为1ms左右,而针对更大的系统,优化的效果更加明显(以IEEE300节点系统为例,未优化时为24ms,优化后对于三个同步机组群仅需5.5ms,两个同步机组群仅需3ms);其余粗化过程耗时小于1ms;还原优化过程G5→G4和G4→G3两步分别耗时<1ms和5ms。
Claims (4)
1.一种基于多层图分割的电力系统主动解列断面快速搜索方法,其特征是,包括如下步骤:
第一阶段:对原始图G0中的点或边进行合并,得到简化图G1,重复该步骤直到最终简化图Gn只包含少量的节点;
第二阶段:对简化图Gn进行初始分割,得到初始解列面;
第三阶段:将简化图Gn依次地还原到原始图G0,并在每一步还原过程中通过优化手段改善分区结果;
第一阶段按照粗化对系统最优解的解空间的影响程度将粗化过程分为等效粗化、基于潮流追踪的粗化和基于重权边匹配的粗化三个子步骤;等效粗化具体步骤为:
a、对多回并行线路进行合并;
b、将节点度为1的冗余节点归并至邻接点;
c、将权重为0的冗余节点归并至任一邻接点;
重复步骤a、b、c直至系统中不再存在上述冗余支路或节点,进入下一步基于潮流追踪的粗化过程。
2.如权利要求1所述的基于多层图分割的电力系统主动解列断面快速搜索方法,其特征是,基于潮流追踪的粗化进一步细化为:
用公式[Vjk,Pjk]=PFT_BFS(vi,Φ,k)表示过程:以节点vi为起点,在节点集合Φ中进行潮流追踪,在追踪深度为k的情况下,可追踪到的节点集合为Vik,对应追踪到的功率为Pik,其中i为节点编号,k为追踪深度;
1)将发电机根据同调性划分为n个集合,置i=1;
2)将第i组同调机群中的所有发电节点加入集合F={vFj},j=1:m其中m为F中节点的个数,vFj为发电机编号为j的节点;置k=1;
3)对集合{vFj},j=1:m中全部节点并行地执行过程:
[Vjk,Pjk]=PFT_BFS(vFj,Φ,k)(1)
其中Pjk为追踪到的节点j的功率,式(1)意为通过广度优先的方法,以节点vFj为起点,在节点集合Φ中进行潮流追踪,在追踪深度为k的情况下,可追踪到的节点集合为Vjk,对应追踪到的功率为Pjk若所有过程均搜索结束,其中Pjk为节点编号为i的节点功率,跳至步骤5;否则,判断集合Vjk之间有无交集:若无,k=k+1,重复步骤3;若有,进入步骤4;
4)假设其中Vpk、Vqk为以vpk,vqk为追踪起点追踪到的节点集合,p、q为节点编号,为空集,令z=k,k为追踪深度,分别继续执行[Vpz,Ppz]=PFT_BFS(vFp,Vpk∪Vqk,z)和[Vqz,Pqz]=PFT_BFS(vFq,Vpk∪Vqk,z)直至追踪完毕,此时,z=zc;令Vpqk=Vpzc∪Vqzc,其中Vpzc、Vqzc为追踪起点分别为节点编号为p、q的节点,追踪深度为zc的节点集合,Vpqk为集合Vpzc与Vqzc的并集,并将Vpzc∩Vqzc中的节点重复追踪到的功率相累加,得到Vpqk对应的追踪功率Ppqk,Ppqk为集合Vpqk中相应节点的功率;进一步,原过程[Vpk,Ppk]=PFT_BFS(vFp,Φ,k)和[Vqk,Pqk]=PFT_BFS(vFq,Φ,k)将被[Vpqk,Ppqk]=PFT_BFS(Vpqk,Φ,k)统一取代,k=k+1;返回步骤3;
5)i=i+1,若i≤n,返回步骤2;否则所有同步机群的搜索已结束;
通过潮流追踪之后,负荷节点vi归属发电节点vj的隶属度可定义为sij=Pij/Pi,即由发电节点vj供给负荷节点vi的功率Pij占节点vi总负荷Pi的比例,通过设置隶属度的阈值可将负荷与对应的发电节点进行归并,从而围绕多个同步机群可形成多个广义节点;
基于重权边匹配的粗化:对节点通过重权边匹配(HeavyEdgeMatching,HEM)的方法加以合并。
3.如权利要求1所述的基于多层图分割的电力系统主动解列断面快速搜索方法,其特征是,第二阶段:对Gn进行初始分割,得到初始解列面,具体为,选择图生长法的一种改进形式-贪婪图生长法(GreedyGraphGrowingPartitioning,GGGP)作为初始分区算法,即以N个节点为起点以BFS为基础生长出N片区域直至覆盖整个图,而在生长过程中优先加入使分区平衡度得到显著改善的节点。
4.如权利要求1所述的基于多层图分割的电力系统主动解列断面快速搜索方法,其特征是,针对任意系统的一个边界节点vi用集合Sm表示与原分区相连的第m条通路所包含的所有节点,m=1,…,NS,NS为通路总数,其中Sm中的节点又可以分为两类,优化前与vi同属一个分区的记为Sm1,否则记为Sm2用集合Rn,n=1,…,NR;用集合Rn表示与交换分区相连的第n条通路所包含的所有节点,n=1,…,NR,NR为与交换分区相连的通路总数,其中优化前与vi同属一个分区的记为Rm1,否则记为Rm2,则对于vi,同时考虑其保留在原分区和交换至其他分区的连通性约束条件,留在原分区的节点需要将通路上与其连接的节点一并进行约束,以此保证其连通性,同理,对于交换至其他分区的节点,也需将通路上与其连接的节点一并进行约束,按此求解可以保证优化解的拓扑连通性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210442714.9A CN102999788B (zh) | 2012-11-06 | 2012-11-06 | 基于多层图分割的电力系统主动解列断面快速搜索方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210442714.9A CN102999788B (zh) | 2012-11-06 | 2012-11-06 | 基于多层图分割的电力系统主动解列断面快速搜索方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102999788A CN102999788A (zh) | 2013-03-27 |
CN102999788B true CN102999788B (zh) | 2016-03-30 |
Family
ID=47928335
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210442714.9A Expired - Fee Related CN102999788B (zh) | 2012-11-06 | 2012-11-06 | 基于多层图分割的电力系统主动解列断面快速搜索方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102999788B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246933A (zh) * | 2013-04-02 | 2013-08-14 | 中国南方电网有限责任公司 | 基于虚拟阻抗矩阵和边界节点迁移的在线校核分区方法 |
CN104866611B (zh) * | 2015-06-04 | 2018-03-27 | 苗伟威 | 一种电力系统主动解列断面的实用化快速搜索方法 |
CN105515057B (zh) * | 2015-12-14 | 2018-05-29 | 国家电网公司 | 基于潮流追踪的主动解列策略空间缩减方法及系统 |
CN106383960B (zh) * | 2016-09-28 | 2019-07-02 | 天津大学 | 基于最小割集分析法的多电飞机电力系统可靠性分析方法 |
CN108281963B (zh) * | 2018-01-17 | 2019-08-27 | 浙江大学 | 一种适用于含多个柔性多状态开关的配电网分区方法 |
CN111784106B (zh) * | 2020-05-26 | 2022-07-01 | 国电南瑞南京控制系统有限公司 | 一种基于受扰轨迹的主动解列断面离线筛选方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101741611A (zh) * | 2009-12-03 | 2010-06-16 | 哈尔滨工业大学 | 基于MLkP/CR算法的无向图分割方法 |
-
2012
- 2012-11-06 CN CN201210442714.9A patent/CN102999788B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101741611A (zh) * | 2009-12-03 | 2010-06-16 | 哈尔滨工业大学 | 基于MLkP/CR算法的无向图分割方法 |
Non-Patent Citations (1)
Title |
---|
A fast and high quality multilevel scheme for partitioning irregular graphs;George Karipys and Vipin Kumar;《SIAM Journal on Scientific Computing》;19981231;第363-370页及图1、2 * |
Also Published As
Publication number | Publication date |
---|---|
CN102999788A (zh) | 2013-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102999788B (zh) | 基于多层图分割的电力系统主动解列断面快速搜索方法 | |
CN102280885B (zh) | 一种将电力系统解列成多个子系统的方法 | |
CN103532140B (zh) | 一种含分布式电源的配电网故障后供电恢复的方法及系统 | |
CN104077393B (zh) | 一种基于半监督谱聚类的最优解列断面搜索方法 | |
CN103093037B (zh) | 基于主从问题交替优化的电力系统解列断面搜索方法 | |
CN102799950B (zh) | 基于粒子群算法的船舶电网重构优化方法 | |
CN102708250A (zh) | 基于树形分层双向迭代的电力系统数字混合仿真方法 | |
CN102136721A (zh) | 一种计及多目标优化的供电恢复推理搜索方法 | |
CN103401241B (zh) | 一种基于分布式电源自组网策略的微电网供电恢复方法 | |
CN104466959A (zh) | 电力系统关键线路辨识方法和系统 | |
CN104362624A (zh) | 含分布式电源配电网故障的主网与孤岛恢复同步算法 | |
CN103746366A (zh) | 一种基于节点重要度的电力系统恢复方法 | |
Zhai et al. | Distributed optimal power flow for VSC-MTDC meshed AC/DC grids using ALADIN | |
CN108667027B (zh) | 含柔性直流的交直流系统潮流转移搜索与量化分析方法 | |
CN111799800A (zh) | 一种交直流混合配电网潮流计算方法 | |
CN108233367A (zh) | 含分布式电源配电网最优孤岛划分计算方法 | |
CN107846018A (zh) | 基于路径搜索算法的配网最大供电能力评估方法 | |
CN104934982B (zh) | 一种多直流馈入系统的直流分区方法 | |
CN110929403B (zh) | 一种基于环路的微电网拓扑规划方法 | |
CN109768543B (zh) | 一种基于混合整数线性规划的弹性保底网架搜索建模方法 | |
Li et al. | Distributed neurodynamic optimization for multi-energy management with time-varying external disturbances considering time-varying emission limitations and load demand in multi-microgrid | |
CN110120673A (zh) | 基于戴维南等值参数辨识的分布式输配协同无功优化方法及系统 | |
CN105870917A (zh) | 一种电力电网供给系统 | |
CN110752594B (zh) | 基于改进凝聚层次聚类法的配电网路径优化方法 | |
CN106712059A (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 | ||
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: 20160330 Termination date: 20201106 |