CN106127237A - 基于谱聚类的含vsc‑hvdc交直流系统最优解列断面搜索方法 - Google Patents

基于谱聚类的含vsc‑hvdc交直流系统最优解列断面搜索方法 Download PDF

Info

Publication number
CN106127237A
CN106127237A CN201610436573.8A CN201610436573A CN106127237A CN 106127237 A CN106127237 A CN 106127237A CN 201610436573 A CN201610436573 A CN 201610436573A CN 106127237 A CN106127237 A CN 106127237A
Authority
CN
China
Prior art keywords
matrix
formula
neighboring rights
electromotor
notequal
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
Application number
CN201610436573.8A
Other languages
English (en)
Other versions
CN106127237B (zh
Inventor
李扬
张明泽
王振浩
李国庆
辛业春
王朝斌
王鹤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Original Assignee
Northeast Dianli University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeast Dianli University filed Critical Northeast Dianli University
Priority to CN201610436573.8A priority Critical patent/CN106127237B/zh
Publication of CN106127237A publication Critical patent/CN106127237A/zh
Application granted granted Critical
Publication of CN106127237B publication Critical patent/CN106127237B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2323Non-hierarchical techniques based on graph theory, e.g. minimum spanning trees [MST] or graph cuts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2321Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J5/00Circuit arrangements for transfer of electric power between ac networks and dc networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Discrete Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明是一种基于谱聚类的含VSC‑HVDC交直流系统最优解列断面搜索方法,其特点是,包括根据电力系统结构构建电力系统无向边权图G、根据潮流数据计算邻接权矩阵W和度矩阵A、判断系统是否因大扰动而失稳、根据发电机受扰后转子功角摇摆曲线将发电机组分为k个同调群,更新邻接权矩阵W和度矩阵A、根据VSC‑HVDC换流站落点,更新邻接权矩阵W和度矩阵A、计算规范化拉普拉斯矩阵LN、根据同调机群数k,通过基于k均值聚类算法对规范化拉普拉斯矩阵LN的前k最小特征值对应的特征向量进行聚类,确定最优解列断面等步骤。具有科学合理,适用性强,可靠性高等优点。

Description

基于谱聚类的含VSC-HVDC交直流系统最优解列断面搜索方法
技术领域
本发明属于电力系统暂态安全稳定与控制技术领域,尤其涉及一种基于谱聚类的含VSC-HVDC交直流系统最优解列断面搜索方法。
背景技术
从国内外近些年发生的大停电事故中可以得知,电力系统在运行和设计中尽管采取一系列措施提高系统稳定性,但仍有可能会遇到突发故障情况使系统丧失稳定,因此必须通过了解系统失去稳定后的现象并采取相应措施以减轻稳定破坏后所带来的伤害,迅速使系统恢复稳定运行。解列工作是限制电力系统遭受大扰动后崩溃的有效手段,通过将系统解列成几个相对独立的子系统避免大停电事故发生。解列工作的核心步骤是确定合理的解列断面,断面即连接各个独立子系统线路的集合,如果将失去稳定的系统在错误的断面处解列,不仅为解列后采用一系列控制措施恢复系统稳定增加难度,而且极有可能使系统进一步失稳直至崩溃,因此如何确定最优解列断面已经成为国内外电力系统方向研究的热点。
最优断面搜索问题可以视为约束组合优化问题,现有方法的目标函数可以分为两种主要类型,最小不平衡功率以及最小有功潮流冲击。最小不平衡功率方法最小化各个子系统内部不平衡功率,通过子系统内部切机、切负荷,使有功功率的输出与消耗近似相等。最小有功潮流冲击方法最小化解列后系统有功潮流的改变,由此也降低了有功潮流改变时的扰动、传输线路过载的危险,提高了岛屿暂态稳定性。同时,最优断面搜索必须满足大量的约束,如负载平衡、发电机一致性、输电线路的可用性、热限制、电压稳定、暂态稳定性等。然而在实际电力系统中每一条线路都是可能解,随着系统规模增大解空间呈几何指数O(2m)爆炸增长,属于非确定多项式(Non-deterministic Polynomial,NP)难题,满足所有以上约束求解方案将会变得十分困难,甚至无法确定可行解是否真实存在。但若只考虑以上约束条件的一部分,如负载平衡和发电机相干性,产生一组可行的候选解决方案,这组候选解决方案可以与其他约束相协调,最终找到满足所有约束的最优解决方案,这种近似考虑降低了解决问题的复杂性,可以应用于实际电力系统。
根据上述近似方法国内外已研究出很多最优断面搜索方案。通过启发式搜索方法可以克服这一难题,或通过仅解决该问题一个简化的网络模型,或选择原有的电力系统的子集进行求解,例如使用基于有序二元决策图的方法,电力系统模型被简化为含有少于大约40个节点。但化简网络模型缩小决策空间的同时,一些化简期间丢失的解决方案会可能比由该算法找到的最终解决方案更加合理。虽然启发式搜索方法通常具备灵活性高、计算速度快的优点,但因为这些方法趋向于局部收敛,而不是全局最小,因此不能保证最终解决方案是否是最优解。将求解最优断面的NP难题转化为其他相对容易求解的问题正是当前热门研究方向,谱聚类算法可以将实际NP难题转化为图论中最小切的问题,通过松弛法可以在求得在实数阈中的最优解,具有较高的计算效率并且不会丢失解集,但传统谱聚类算法如果仅考虑求解最小切而忽略切割后子图内部的密度分布,通常会出现某个子图内部只包含孤立节点的现象,不适用于应用到电力系统中解决实际问题,但已有研究表明若使用最小规范切准则可以避免节点单独成类现象发生。
随着电压源换流器型直流输电系统(Voltage Souce Converter Based HighVoltage Direct Current Transmission,VSC-HVDC)技术逐渐成熟,其相对于传统直流输电的优势已经逐渐显现。相比于传统直流输电,VSC-HVDC具备换流器无需无功补偿、不存在换相失败、可以为无源网络供电、同时独立调节有功功率和无功功率、连接异步系统等优点。交流系统解列成各个孤立子系统后,必然会出现某些子系统内部有功功率充裕但另一些子系统有功功率匮乏的情况,此时只能通过发电机调速器、低频减载等装置使子系统内部恢复稳定,而通过VSC-HVDC将两个子系统连接后,通过直流线路进行电力交换优化电源、负荷匹配程度,限制切机、切负荷,从而降低了大扰动对经济造成的损失。目前在国内VSC-HVDC属于快速发展建设阶段,福建省、浙江省等多个地区的电力系统都已建成并投产,因此关于含VSC-HVDC交直流系统最优解列断面搜索方法迫在眉睫。
发明内容
本发明的目的在于提供一种科学合理,适用性强,可靠性高,效果佳的基于谱聚类的含VSC-HVDC交直流系统最优解列断面搜索方法。
本发明的目的是这样实现的:一种基于谱聚类的含VSC-HVDC交直流系统最优解列断面搜索方法,其特征在于,它包括下列步骤:
步骤1:根据电力系统结构构建电力系统无向边权图G
根据电力系统结构构建电力系统无向边权图G,谱聚类算法是一种基于图论的聚类算法,图论中无向图模型G(V,VG,E,W)用来描述m台发电机,n条线路的系统,点集V={v1,...,vn},边集E中元素为eij(i,j=1,...,n)表示母线以及输电线,VG是V的子集包含那些有发电机直接相连的母线,矩阵W是一组边的权重,其中的元素为(1)式:
式中:wij是边eij的权值,边权值的大小代表两点相似程度的高低,从图论中可以得知,此时求无向图模型G(V,VG,E,W)切的代价c表示为(2)式:
c = 1 2 Σ c ( G s , G t ) = 1 2 Σ i ∈ G s , j ∈ G t w i j , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 2 )
由步骤1可知电力系统最优断面搜索问题可以转换为求无向图G的最小代价切的问题;
步骤2:根据潮流数据计算邻接权矩阵W和度矩阵A,若系统遭受大干扰后发电机组构成k个同调机群VG1,VG2,...,VGk,为了保持解列后各个独立子系统内部稳定,需要相对应解列成k个孤岛V1,V2,...,Vk。最小有功冲击为目标函数的表达式为(3)式:
P = m i n V 1 , V 2 , ... , V k ⋐ V ( 1 2 Σ i ∈ V s , j ∈ V t | P i j | ) , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 3 )
因此邻接权矩阵W由n个节点系统节点之间交换的有功功率来构造,同时必须考虑系统中的网损以产生用于谱聚类所需的对称无向图,为了确保矩阵W是对称的,矩阵W的元素定义为(IPijI+IPjiI)/2,Pij表示考虑网损后两个相邻的节点i,j之间传输的有功功率,表示为(4)式:
W i j = | P i j | + | P j i | 2 i f i ≠ j 0 i f i = j - - - ( 4 )
定义度矩阵A为对角阵,其只含有对角线元素,对角线元素为连接至该节点的所有边的权重的总和:
A i j = Σ j = 1 , j ≠ i n W i j , i = j 0 , i ≠ j - - - ( 5 )
通过这种方式定义,无向图的邻接权矩阵W和度矩阵A均为对称矩阵;
步骤3:判断系统是否因大扰动而失稳,若已经失稳进入步骤4,否则回到步骤1;
步骤4:根据发电机受扰后转子功角摇摆曲线将发电机组分为k个同调群,更新邻接权矩阵W和度矩阵A。为满足解列后同调机群必须处于同一子系统,而非同调机群务必处于不同子系统,需要增加成对约束条件来约束求解过程,因此需要将发电机分群信息作为约束监督信息对邻接权矩阵W和度矩阵A进行修改,发电机同调约束包括Must-Link约束与Cannot-Link约束,Must-Link约束表示任意两台同调发电机节点之间务必至少存在一条交流路径相连,Cannot-Link约束表示任意两台非同调发电机节点之间务必没有交流路径相连,因此邻接权矩阵W需按(6)式进行修改:
W ( V G i , V G j ) = ζ , ( V G i , V G j ) ∈ M u s t - L i n k 0 , ( V G i , V G j ) ∈ C a n n o t - L i n k , ( i , j = 1 , 2 , ... , k , i ≠ j ) - - - ( 6 )
式中:ζ为经验值,取邻接权矩阵W中最大元素的1.5倍,得到修改后的邻接权矩阵W,度矩阵A按公式(5)进行修改;
步骤5:根据VSC-HVDC换流站落点,更新邻接权矩阵W和度矩阵A。为迫使解列后各个子系统通过VSC-HVDC相连,确保解列后VSC-HVDC两端换流站不处于同一孤岛内部,直流落点约束属于Cannot-Link约束,表示两个换流站节点之间没有务必没有交流线路相连,邻接权矩阵W需按(7)式进行修改:
W(nVSC1,nVSC2)=0,(nVSC1,nVSC2)∈Cannot-link (7)
式中:nVSC1,nVSC2表示两个换流站节点;
步骤6:计算规范化拉普拉斯矩阵LN
LN=A-1(A-W) (8)
步骤7:基于Lanczos算法计算LNX=λX的前k个最小特征值对应的特征向量。根据最小规范切准则,求无向边权图G的最小规范切N-C可以表示为:
N - C = Σ c ( G s , G t ) v ( G s ) , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 9 )
式中v(Gs)表示Gs的节点的度之和,定义k个聚类指示向量βj=(β1j2j,...,βnj)T(j=1,2,...k),定义βij为:
令未规范化的拉普拉斯矩阵L=A-W,将k个聚类指示向量βj构成矩阵F,则求无向边权图G的最小规范切可以等效为:
m i n G 1 , G 2 , ... , G k T ( F T L F ) s . t . F T L F = E
式中T(FTLF)表示求矩阵FTLF的迹。由于求图划分准则的最优解属于NP难题,在多项式时间内无法快速求解,但通过将求解F的问题转化为考虑F的连续放松形式,根据Rayleigh-Ritz定理,F的解对应规范化拉普拉斯矩阵LN的前k最小特征值对应的特征向量;
步骤8:根据同调机群数k,通过基于k均值聚类算法对规范化拉普拉斯矩阵LN的前k最小特征值对应的特征向量进行聚类,确定最优解列断面;
本发明的一种基于谱聚类的含VSC-HVDC交直流系统最优解列断面搜索方法的优点体现在:使用谱聚类最小规范切(Normalized Cut,N-C)准则,能够避免发生节点单独成类现象;以最小有功潮流冲击为目标函数,降低了解列时刻有功潮流改变时的扰动、传输线路过载的危险,提高了岛屿暂态稳定性;根据发电机同调约束限制同调机组相连、非同调机组分离,确保解列工作可以平息振荡,解列后各个子系统内部保持同步;为进一步提高各个子系统的稳定性的同时切机切负荷尽可能少,加入直流落点约束使各个子系统能够通过VSC-HVDC相连。具有科学合理,适用性强,可靠性高等优点。
具体实施方式
下面结合具体实施方式对本发明作详细说明。
本发明的一种基于谱聚类的含VSC-HVDC交直流系统最优解列断面搜索方法,包括下列步骤:
步骤1:根据电力系统结构构建电力系统无向边权图G
根据电力系统结构构建电力系统无向边权图G,谱聚类算法是一种基于图论的聚类算法,图论中无向图模型G(V,VG,E,W)用来描述m台发电机,n条线路的系统,点集V={v1,...,vn},边集E中元素为eij(i,j=1,...,n)表示母线以及输电线,VG是V的子集包含那些有发电机直接相连的母线,矩阵W是一组边的权重,其中的元素为(1)式:
式中:wij是边eij的权值,边权值的大小代表两点相似程度的高低,从图论中可以得知,此时求无向图模型G(V,VG,E,W)切的代价c表示为(2)式:
c = 1 2 Σ c ( G s , G t ) = 1 2 Σ i ∈ G s , j ∈ G t w i j , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 2 )
由步骤1可知电力系统最优断面搜索问题可以转换为求无向图G的最小代价切的问题;
步骤2:根据潮流数据计算邻接权矩阵W和度矩阵A,若系统遭受大干扰后发电机组构成k个同调机群VG1,VG2,...,VGk,为了保持解列后各个独立子系统内部稳定,需要相对应解列成k个孤岛V1,V2,...,Vk。最小有功冲击为目标函数的表达式为(3)式:
P = m i n V 1 , V 2 , ... , V k ⋐ V ( 1 2 Σ i ∈ V s , j ∈ V t | P i j | ) , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 3 )
因此邻接权矩阵W由n个节点系统节点之间交换的有功功率来构造,同时必须考虑系统中的网损以产生用于谱聚类所需的对称无向图,为了确保矩阵W是对称的,矩阵W的元素定义为(IPijI+IPjiI)/2,Pij表示考虑网损后两个相邻的节点i,j之间传输的有功功率,表示为(4)式:
W i j = | P i j | + | P j i | 2 i f i ≠ j 0 i f i = j - - - ( 4 )
定义度矩阵A为对角阵,其只含有对角线元素,对角线元素为连接至该节点的所有边的权重的总和:
A i j = Σ j = 1 , j ≠ i n W i j , i = j 0 , i ≠ j - - - ( 5 )
通过这种方式定义,无向图的邻接权矩阵W和度矩阵A均为对称矩阵;
步骤3:判断系统是否因大扰动而失稳,若已经失稳进入步骤4,否则回到步骤1;
步骤4:根据发电机受扰后转子功角摇摆曲线将发电机组分为k个同调群,更新邻接权矩阵W和度矩阵A。为满足解列后同调机群必须处于同一子系统,而非同调机群务必处于不同子系统,需要增加成对约束条件来约束求解过程,因此需要将发电机分群信息作为约束监督信息对邻接权矩阵W和度矩阵A进行修改,发电机同调约束包括Must-Link约束与Cannot-Link约束,Must-Link约束表示任意两台同调发电机节点之间务必至少存在一条交流路径相连,Cannot-Link约束表示任意两台非同调发电机节点之间务必没有交流路径相连,因此邻接权矩阵W需按(6)式进行修改:
W ( V G i , V G j ) = ζ , ( V G i , V G j ) ∈ M u s t - L i n k 0 , ( V G i , V G j ) ∈ C a n n o t - L i n k , ( i , j = 1 , 2 , ... , k , i ≠ j ) - - - ( 6 )
式中:ζ为经验值,取邻接权矩阵W中最大元素的1.5倍,得到修改后的邻接权矩阵W,度矩阵A按公式(5)进行修改;
步骤5:根据VSC-HVDC换流站落点,更新邻接权矩阵W和度矩阵A。为迫使解列后各个子系统通过VSC-HVDC相连,确保解列后VSC-HVDC两端换流站不处于同一孤岛内部,直流落点约束属于Cannot-Link约束,表示两个换流站节点之间没有务必没有交流线路相连,邻接权矩阵W需按(7)式进行修改:
W(nVSC1,nVSC2)=0,(nVSC1,nVSC2)∈Cannot-link (7)
式中:nVSC1,nVSC2表示两个换流站节点;
步骤6:计算规范化拉普拉斯矩阵LN
LN=A-1(A-W) (8)
步骤7:基于Lanczos算法计算LNX=λX的前k个最小特征值对应的特征向量。根据最小规范切准则,求无向边权图G的最小规范切N-C可以表示为:
N - C = Σ c ( G s , G t ) v ( G s ) , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 9 )
式中v(Gs)表示Gs的节点的度之和,定义k个聚类指示向量βj=(β1j2j,...,βnj)T(j=1,2,...k),定义βij为:
令未规范化的拉普拉斯矩阵L=A-W,将k个聚类指示向量βj构成矩阵F,则求无向边权图G的最小规范切可以等效为:
m i n G 1 , G 2 , ... , G k T ( F T L F ) s . t . F T L F = E
式中T(FTLF)表示求矩阵FTLF的迹。由于求图划分准则的最优解属于NP难题,在多项式时间内无法快速求解,但通过将求解F的问题转化为考虑F的连续放松形式,根据Rayleigh-Ritz定理,F的解对应规范化拉普拉斯矩阵LN的前k最小特征值对应的特征向量;
步骤8:根据同调机群数k,通过基于k均值聚类算法对规范化拉普拉斯矩阵LN的前k最小特征值对应的特征向量进行聚类,确定最优解列断面;
本发明的保护范围并不局限于具体实施方式,本领域技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (1)

1.一种基于谱聚类的含VSC-HVDC交直流系统最优解列断面搜索方法,其特征在于,它包括下列步骤:
步骤1:根据电力系统结构构建电力系统无向边权图G
根据电力系统结构构建电力系统无向边权图G,谱聚类算法是一种基于图论的聚类算法,图论中无向图模型G(V,VG,E,W)用来描述m台发电机,n条线路的系统,点集V={v1,...,vn},边集E中元素为eij(i,j=1,...,n)表示母线以及输电线,VG是V的子集包含那些有发电机直接相连的母线,矩阵W是一组边的权重,其中的元素为(1)式:
W i j = w i j , ( i , j ) ∈ E 0 , ( i , j ) ∉ E - - - ( 1 )
式中:wij是边eij的权值,边权值的大小代表两点相似程度的高低,从图论中可以得知,此时求无向图模型G(V,VG,E,W)切的代价c表示为(2)式:
c = 1 2 Σ c ( G s , G t ) = 1 2 Σ i ∈ G s , j ∈ G t w i j , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 2 )
由步骤1可知电力系统最优断面搜索问题可以转换为求无向图G的最小代价切的问题;
步骤2:根据潮流数据计算邻接权矩阵W和度矩阵A,若系统遭受大干扰后发电机组构成k个同调机群VG1,VG2,...,VGk,为了保持解列后各个独立子系统内部稳定,需要相对应解列成k个孤岛V1,V2,...,Vk。最小有功冲击为目标函数的表达式为(3)式:
P = m i n V 1 , V 2 , ... , V k ⋐ V ( 1 2 Σ i ∈ V s , j ∈ V t | P i j | ) , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 3 )
因此邻接权矩阵W由n个节点系统节点之间交换的有功功率来构造,同时必须考虑系统中的网损以产生用于谱聚类所需的对称无向图,为了确保矩阵W是对称的,矩阵W的元素定义为(|Pij|+|Pji|)/2,Pij表示考虑网损后两个相邻的节点i,j之间传输的有功功率,表示为(4)式:
W i j = | P i j | + | P j i | 2 i f i ≠ j 0 i f i = j - - - ( 4 )
定义度矩阵A为对角阵,其只含有对角线元素,对角线元素为连接至该节点的所有边的权重的总和:
A i j = Σ j = 1 , j ≠ i n W i j , i = j 0 , i ≠ j - - - ( 5 )
通过这种方式定义,无向图的邻接权矩阵W和度矩阵A均为对称矩阵;
步骤3:判断系统是否因大扰动而失稳,若已经失稳进入步骤4,否则回到步骤1;
步骤4:根据发电机受扰后转子功角摇摆曲线将发电机组分为k个同调群,更新邻接权矩阵W和度矩阵A。为满足解列后同调机群必须处于同一子系统,而非同调机群务必处于不同子系统,需要增加成对约束条件来约束求解过程,因此需要将发电机分群信息作为约束监督信息对邻接权矩阵W和度矩阵A进行修改,发电机同调约束包括Must-Link约束与Cannot-Link约束,Must-Link约束表示任意两台同调发电机节点之间务必至少存在一条交流路径相连,Cannot-Link约束表示任意两台非同调发电机节点之间务必没有交流路径相连,因此邻接权矩阵W需按(6)式进行修改:
W ( V G i , V G j ) = ζ , ( V G i , V G j ) ∈ M u s t - L i n k 0 , ( V G i , V G j ) ∈ C a n n o t - L i n k , ( i , j = 1 , 2 , ... , k , i ≠ j ) - - - ( 6 )
式中:ζ为经验值,取邻接权矩阵W中最大元素的1.5倍,得到修改后的邻接权矩阵W,度矩阵A按公式(5)进行修改;
步骤5:根据VSC-HVDC换流站落点,更新邻接权矩阵W和度矩阵A。为迫使解列后各个子系统通过VSC-HVDC相连,确保解列后VSC-HVDC两端换流站不处于同一孤岛内部,直流落点约束属于Cannot-Link约束,表示两个换流站节点之间没有务必没有交流线路相连,邻接权矩阵W需按(7)式进行修改:
W(nVSC1,nVSC2)=0,(nVSC1,nVSC2)∈Cannot-link (7)
式中:nVSC1,nVSC2表示两个换流站节点;
步骤6:计算规范化拉普拉斯矩阵LN
LN=A-1(A-W) (8)
步骤7:基于Lanczos算法计算LNX=λX的前k个最小特征值对应的特征向量。根据最小规范切准则,求无向边权图G的最小规范切N-C可以表示为:
N - C = Σ c ( G s , G t ) v ( G s ) , ( s , t = 1 , 2 , ... , k , s ≠ t ) - - - ( 9 )
式中v(Gs)表示Gs的节点的度之和,定义k个聚类指示向量βj=(β1j2j,...,βnj)T(j=1,2,...k),定义βij为:
β i j = 1 v ( G j ) , i ∈ G j 0 , i ∉ G j , ( i = 1 , 2 , ... , n ; j = 1 , 2 , ... , k ) - - - ( 10 )
令未规范化的拉普拉斯矩阵L=A-W,将k个聚类指示向量βj构成矩阵F,则求无向边权图G的最小规范切可以等效为:
m i n G 1 , G 2 , ... , G k T ( F T L F ) s . t . F T L F = E
式中T(FTLF)表示求矩阵FTLF的迹。由于求图划分准则的最优解属于NP难题,在多项式时间内无法快速求解,但通过将求解F的问题转化为考虑F的连续放松形式,根据Rayleigh-Ritz定理,F的解对应规范化拉普拉斯矩阵LN的前k最小特征值对应的特征向量;
步骤8:根据同调机群数k,通过基于k均值聚类算法对规范化拉普拉斯矩阵LN的前k最小特征值对应的特征向量进行聚类,确定最优解列断面。
CN201610436573.8A 2016-06-17 2016-06-17 基于谱聚类的含vsc-hvdc交直流系统最优解列断面搜索方法 Expired - Fee Related CN106127237B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610436573.8A CN106127237B (zh) 2016-06-17 2016-06-17 基于谱聚类的含vsc-hvdc交直流系统最优解列断面搜索方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610436573.8A CN106127237B (zh) 2016-06-17 2016-06-17 基于谱聚类的含vsc-hvdc交直流系统最优解列断面搜索方法

Publications (2)

Publication Number Publication Date
CN106127237A true CN106127237A (zh) 2016-11-16
CN106127237B CN106127237B (zh) 2019-11-19

Family

ID=57469853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610436573.8A Expired - Fee Related CN106127237B (zh) 2016-06-17 2016-06-17 基于谱聚类的含vsc-hvdc交直流系统最优解列断面搜索方法

Country Status (1)

Country Link
CN (1) CN106127237B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330809A (zh) * 2017-04-12 2017-11-07 广东电网有限责任公司电力调度控制中心 基于半监督谱聚类算法的黑启动分区方法
CN107994569A (zh) * 2017-11-30 2018-05-04 天津大学 一种基于谱聚类分析的互联电网热稳定关键断面识别方法
CN109038635A (zh) * 2018-08-02 2018-12-18 南瑞集团有限公司 一种计及直流调制能力的解列断面优化选择方法及系统
CN109995023A (zh) * 2019-03-01 2019-07-09 上海电力学院 基于谱聚类区域划分的电力系统稳定性的判断方法和系统
CN110336266A (zh) * 2019-05-07 2019-10-15 山东大学 交直流混联受端电网的直流分群规划方法
CN111585280A (zh) * 2020-03-06 2020-08-25 山东大学 阻断多直流换相失败连锁的电网临时解列方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7750501B2 (en) * 2005-10-27 2010-07-06 Continental Automotive Systems Us, Inc. System and method of over voltage control for a power system
US8848401B2 (en) * 2010-10-27 2014-09-30 Alstom Technology Ltd. Modular multilevel power electronic converter having selectively definable circulation path
CN104077393A (zh) * 2014-07-01 2014-10-01 武汉大学 一种基于半监督谱聚类的最优解列断面搜索方法
CN104090985A (zh) * 2014-07-25 2014-10-08 武汉大学 一种基于电气距离的主动解列最优断面搜索方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7750501B2 (en) * 2005-10-27 2010-07-06 Continental Automotive Systems Us, Inc. System and method of over voltage control for a power system
US8848401B2 (en) * 2010-10-27 2014-09-30 Alstom Technology Ltd. Modular multilevel power electronic converter having selectively definable circulation path
CN104077393A (zh) * 2014-07-01 2014-10-01 武汉大学 一种基于半监督谱聚类的最优解列断面搜索方法
CN104090985A (zh) * 2014-07-25 2014-10-08 武汉大学 一种基于电气距离的主动解列最优断面搜索方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王乙斐: "基于Dijkstra算法的最优解列断面快速搜索方法", 《电力自动化设备》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330809A (zh) * 2017-04-12 2017-11-07 广东电网有限责任公司电力调度控制中心 基于半监督谱聚类算法的黑启动分区方法
CN107994569A (zh) * 2017-11-30 2018-05-04 天津大学 一种基于谱聚类分析的互联电网热稳定关键断面识别方法
CN109038635A (zh) * 2018-08-02 2018-12-18 南瑞集团有限公司 一种计及直流调制能力的解列断面优化选择方法及系统
CN109995023A (zh) * 2019-03-01 2019-07-09 上海电力学院 基于谱聚类区域划分的电力系统稳定性的判断方法和系统
CN110336266A (zh) * 2019-05-07 2019-10-15 山东大学 交直流混联受端电网的直流分群规划方法
CN110336266B (zh) * 2019-05-07 2020-09-04 山东大学 交直流混联受端电网的直流分群规划方法
CN111585280A (zh) * 2020-03-06 2020-08-25 山东大学 阻断多直流换相失败连锁的电网临时解列方法及系统

Also Published As

Publication number Publication date
CN106127237B (zh) 2019-11-19

Similar Documents

Publication Publication Date Title
CN106127237A (zh) 基于谱聚类的含vsc‑hvdc交直流系统最优解列断面搜索方法
CN105098762B (zh) 一种含分布式电源配电网的孤岛划分方法
CN109598377B (zh) 一种基于故障约束的交直流混合配电网鲁棒规划方法
Xia et al. Distributed control method for economic dispatch in islanded microgrids with renewable energy sources
Ustun et al. Implementation of Dijkstra's algorithm in a dynamic microgrid for relay hierarchy detection
CN102244394A (zh) 基于归一化谱聚类和约束谱聚类的两阶段主动解列方法
Linda et al. Optimal design of multi-machine power system stabilizer using robust ant colony optimization technique
Dou et al. Transient control for micro-grid with multiple distributed generations based on hybrid system theory
CN107069801B (zh) 一种基于最小用户停电损失的配电网孤岛划分方法
CN106849162B (zh) 大量可再生能源并网的主动配电网自适应鲁棒优化方法
Mishra et al. A scheduled intentional islanding method based on ranking of possible islanding zone
Zhang et al. SVC damping controller design based on novel modified fruit fly optimisation algorithm
Zhang et al. Hierarchical flexible operation approach on a VSC-MTDC interconnected hybrid grid with a high share of renewable power
CN109768543B (zh) 一种基于混合整数线性规划的弹性保底网架搜索建模方法
CN106712059B (zh) 基于凸优化理论的电力系统主动解列最优断面搜索方法
Ding et al. A new controlled islanding algorithm based on spectral clustering
Žarković et al. Reliability improvement of distribution system through distribution system planning: MILP vs. GA
Brahim et al. Load frequency control of hybrid power system using classical PID controller
CN116722536A (zh) 基于数据驱动的柔性配电网分布鲁棒扩展规划方法及系统
Matute et al. Transmission expansion planning considering the impact of distributed generation
Bento et al. A new pmu placement method for the assessment of the voltage stability margin
Zhang et al. The Voltage Stabilizing Control Strategy of Off‐Grid Microgrid Cluster Bus Based on Adaptive Genetic Fuzzy Double Closed‐Loop Control
Wang et al. A Novel Method for Islanding in Active Distribution Network Considering Distributed Generation.
Pourmirasghariyan et al. Application of max flow-min cut theory to find the best placement of electronic-based DC-PFCs for enhancing static security in MT-HVDC meshed grids
Shi et al. Distribution network post-disaster power supply recovery strategy based on path description

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Li Yang

Inventor after: Zhang Mingze

Inventor after: Wang Zhenhao

Inventor after: Li Guoqing

Inventor after: Xin Yechun

Inventor after: Wang Chaobin

Inventor after: Wang He

Inventor after: Liu Ruitong

Inventor before: Li Yang

Inventor before: Zhang Mingze

Inventor before: Wang Zhenhao

Inventor before: Li Guoqing

Inventor before: Xin Yechun

Inventor before: Wang Chaobin

Inventor before: Wang He

CB03 Change of inventor or designer information
TA01 Transfer of patent application right

Effective date of registration: 20170323

Address after: Jilin City, Jilin province Changchun ship 132012 Camp Road No. 169

Applicant after: Northeast Dianli University

Applicant after: Electric Power Research Institute of State Grid Liaoning Electric Power Co., Ltd.

Address before: Jilin City, Jilin province Changchun ship 132012 Camp Road No. 169

Applicant before: Northeast Dianli University

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20191119

Termination date: 20200617

CF01 Termination of patent right due to non-payment of annual fee