CN107994569A - 一种基于谱聚类分析的互联电网热稳定关键断面识别方法 - Google Patents
一种基于谱聚类分析的互联电网热稳定关键断面识别方法 Download PDFInfo
- Publication number
- CN107994569A CN107994569A CN201711248087.4A CN201711248087A CN107994569A CN 107994569 A CN107994569 A CN 107994569A CN 201711248087 A CN201711248087 A CN 201711248087A CN 107994569 A CN107994569 A CN 107994569A
- Authority
- CN
- China
- Prior art keywords
- branch
- power grid
- thermostabilization
- matrix
- key sections
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2323—Non-hierarchical techniques based on graph theory, e.g. minimum spanning trees [MST] or graph cuts
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Discrete Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于谱聚类分析的互联电网热稳定关键断面识别方法,包括:步骤一、给定电网的网络连接方式,各支路的参数,各条支路的热稳定传输限值,发电机出力和负荷水平,进行潮流计算;步骤二、根据步骤一所得潮流计算结果,计算各条支路的热稳定裕度;步骤三、根据电网结构,构建无向加权图G(V,E),设置图中各边的权值,用W表示,形成G(V,E)连接矩阵A和度矩阵D,计算拉普拉斯矩阵L;步骤四、计算L的特征值λ与特征向量rV,对L的所有特征值进行升序排序,提取出所有特征值中第2小的特征值λ2所对应的特征向量rV2;步骤五、根据电网的接线关系和rV2筛选出电网的热稳定关键断面的构成支路,将电网所有节点分为两个区域。
Description
技术领域
本发明属于电力系统分析方法,尤其涉及一种基于谱聚类分析的互联电网热稳定关键断面识别方法。
背景技术
随着区域互联的发展,现代电网的规模越来越大,节点和支路数据显著增加。随之而来的是,电网的潮流分布更加复杂,热稳定安全成为限制互联电网安全稳定运行的重要约束。由于电网的网络连接关系复杂,支路众多,在实际运行中,若要求电力系统调度人员同时监视所有线路的负载情况,则意味着巨大的工作量。考虑到实际电网中,热稳定安全问题往往出现在一些薄弱线路构成的关键断面上,而这些关键断面通常由少数几条支路构成,因此,对关键断面进行监控是电网运行中处理热稳定约束的重要手段之一。所谓的热稳定关键断面应具备以下一些特性:该断面可将整个网络分割为两个互不联通的子网,并且该断面的热稳定传输裕度较低,容易出现热稳越限问题。
发明内容
针对互联电网热稳定安全问题日益突出的挑战,本发明提出了一种基于谱聚类分析的互联电网热稳定关键断面识别方法。
为了解决上述技术问题,本发明一种基于谱聚类分析的互联电网热稳定关键断面识别方法,包括以下步骤:
步骤一、给定电网的网络连接方式,各支路的参数,各条支路的热稳定传输限值,发电机出力和负荷水平,进行潮流计算,获得电网中所有的支路的视在功率Sl;其中,各支路的参数包括支路的电阻、电抗、对地导纳和最大允许传输潮流
步骤二、根据步骤一所得潮流计算结果,计算各条支路的热稳定裕度;
步骤三、根据电网结构,构建无向加权图G(V,E),设置无向加权图G(V,E)中各边的权值,对应的权重矩阵用W表示,形成无向加权图G(V,E)的连接矩阵A和度矩阵D,计算拉普拉斯矩阵L;
步骤四、计算拉普拉斯矩阵L的特征值λ与特征向量rV,对拉普拉斯矩阵L的所有特征值λ进行升序排序,提取出所有特征值中第2小的特征值λ2所对应的特征向量rV2;
步骤五、根据电网的接线关系和特征向量rV2筛选出电网的热稳定关键断面的构成支路,根据rV2向量各元素的符号,将电网所有节点分为两个区域。
进一步讲,本发明基于谱聚类分析的互联电网热稳定关键断面识别方法步骤二的具体内容包括:
根据步骤一潮流计算得到的电网中所有的支路的视在功率Sl,其中,l表示第l支路的编号,并且支路l的两端节点分别用节点i和节点j表示,第l条支路的热稳定裕度ηl根据式(1)计算:
步骤三的具体内容包括:
步骤1)将电网中的所有节点构成无向加权图G的顶点的集合,用V表示,将电网中所有支路构成无向加权图G的边的集合,用E表示,对于存在并列双回或多回线路的情况,用一条边表示,即节点i和节点j之间存在两条或两条以上的支路,则在无向加权图G中,仅用一条边表示;
步骤2)计算无向加权图G中各边的权重矩阵W,其中,W矩阵的第i行第j列元素用Wij表示,计算公式如下:
步骤3)建立无向加权图G(V,E)的连接矩阵A、度矩阵D和拉普拉斯矩阵L;连接矩阵A的第i行第j列元素用Aij表示;度矩阵D的第i行第j列元素用Dij表示;拉普拉斯矩阵L的第i行第j列元素用Lij表示;具体计算公式如下:
Lij=Dij-Aij(5)
步骤五中,特征向量rV2的第i个元素用ci表示,步骤五的具体内容包括:
步骤1)对于第i个节点,若ci大于0,则节点i属于分区1,若若ci小于等于0,则节点i属于分区2;
步骤2)对于支路l,按照式(6)计算指标κl,其中i和j表示支路两端的节点编号,若κl小于0,则支路l属于电网热稳定关键断面,若κl大于等于0,则支路l不属于电网热稳定关键断面;所有满足κl<0的支路构成电网热稳定关键断面,如式(7)所示,
κl=cicj (6)
CritCut={lκl<0,l=1,2,...,nB} (7)
与现有技术相比,本发明的有益效果是:本发明方法可以考虑网络的拓扑结构和潮流负载水平,确定对电网热稳定安全具有重要影响的关键断面,为电网调度运行人员有效监控电网提供重要的辅助决策依据。
附图说明
图1是本发明用到示例的接线图;
图2是本发明对IEEE24系统识别出的热稳定关键断面。
具体实施方式
下面结合附图和具体实施算例对本发明技术方案作进一步详细描述。
本发明一种基于谱聚类分析的互联电网热稳定关键断面识别方法,对其实施流程详细说明如下:
步骤一:采用IEEE24节点测试系统验证本发明所提方法的有效性与正确性。该算例包括24个节点,38条线路,电网的接线如图1所示。给定电网的网络连接方式,各支路的参数(包括电阻、电抗、对地导纳和最大允许传输潮流),各条支路的热稳定传输限值,发电机出力和负荷水平,对该电网进行潮流计算,获得电网中所有的支路的视在功率。
步骤二:根据步骤一潮流计算得到的电网中所有的支路的视在功率Sl,其中,l表示第l支路的编号,并且支路l的两端节点分别用节点i和节点j表示,第l条支路的热稳定裕度ηl根据式(1)计算:
计算得到各支路的热稳定裕度如表1所示。
表1
步骤三:根据电网的接线结构,构建无向加权图G(V,E),图G的定点集合V共包含24个节点,由于该电网中有4组并列双回线,在将并列双回线合并为单一边后,图G的边集合E共包含34条边。根据支路热稳定裕度计算得到图G的权值矩阵W,连接矩阵A,度矩阵D和拉普拉斯矩阵L。具体过程是:
步骤1)将电网中的所有节点构成无向加权图G的顶点的集合,用V表示,将电网中所有支路构成无向加权图G的边的集合,用E表示,对于存在并列双回或多回线路的情况,用一条边表示,即节点i和节点j之间存在两条或两条以上的支路,则在无向加权图G中,仅用一条边表示;
步骤2)计算无向加权图G中各边的权重矩阵W,其中,W矩阵的第i行第j列元素用Wij表示,计算公式如下:
步骤3)建立无向加权图G(V,E)的连接矩阵A、度矩阵D和拉普拉斯矩阵L;连接矩阵A的第i行第j列元素用Aij表示;度矩阵D的第i行第j列元素用Dij表示;拉普拉斯矩阵L的第i行第j列元素用Lij表示;具体计算公式如下:
Lij=Dij-Aij(5)
步骤四、计算拉普拉斯矩阵L的特征值λ与特征向量rV,对拉普拉斯矩阵L的所有特征值λ进行升序排序,提取出所有特征值中第2小的特征值λ2所对应的特征向量rV2;所有特征值如表2所示,第2小的特征值为15.53,对应的特征向量rV2=[0.00544 -0.003690.02461 -0.007667 -0.013482 -0.011640 -0.92977 -0.095275 -0.01031 -0.023068 -0.003792 0.0055456 0.00941754-0.001939 0.1170937 0.101365 0.1449070.1424190.0725008 0.0598679 0.138161 0.1467022 0.042637 0.089969]。
表2
步骤五、根据电网的接线关系和特征向量rV2筛选出电网的热稳定关键断面的构成支路,根据rV2向量各元素的符号,将电网所有节点分为两个区域。具体过程是:
特征向量rV2的第i个元素用ci表示,步骤五的具体内容包括:
步骤1)把所有节点分为两个区:对于第i个节点,若ci大于0,则节点i属于分区1,若若ci小于等于0,则节点i属于分区2;
步骤2)确定步骤1)中两个区间的联络线,这些联络线就构成了热稳关键断面,对于支路l,按照式(6)计算指标κl,其中i和j表示支路两端的节点编号,若κl小于0,则支路l属于电网热稳定关键断面,若κl大于等于0,则支路l不属于电网热稳定关键断面;所有满足κl<0的支路构成电网热稳定关键断面,如式(7)所示,
κl=cicj (6)
CritCut={lκl<0,l=1,2,...,nB} (7)
根据电网结构和rV2计算得到各支路的κl指标如表3所示,表3中加粗字体表示的支路L1-3、L3-9、L11-14、L12-23和L13-23构成了该电网的热稳定关键断面。如图2所示,该断面将整个电网分割为互补联通的两个子网,并且其热稳定裕度偏低,需在运行中给予重点关注。
表3
尽管上面结合附图对本发明进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨的情况下,还可以做出很多变形,这些均属于本发明的保护之内。
Claims (4)
1.一种基于谱聚类分析的互联电网热稳定关键断面识别方法,其特征在于,包括以下步骤:
步骤一、给定电网的网络连接方式,各支路的参数,各条支路的热稳定传输限值,发电机出力和负荷水平,进行潮流计算,获得电网中所有的支路的视在功率Sl;其中,各支路的参数包括支路的电阻、电抗、对地导纳和最大允许传输潮流Sl M,
步骤二、根据步骤一所得潮流计算结果,计算各条支路的热稳定裕度;
步骤三、根据电网结构,构建无向加权图G(V,E),设置无向加权图G(V,E)中各边的权值,对应的权重矩阵用W表示,形成无向加权图G(V,E)的连接矩阵A和度矩阵D,计算拉普拉斯矩阵L;
步骤四、计算拉普拉斯矩阵L的特征值λ与特征向量rV,对拉普拉斯矩阵L的所有特征值λ进行升序排序,提取出所有特征值中第2小的特征值λ2所对应的特征向量rV2;
步骤五、根据电网的接线关系和特征向量rV2筛选出电网的热稳定关键断面的构成支路,根据rV2向量各元素的符号,将电网所有节点分为两个区域。
2.根据权利要求1所述一种基于谱聚类分析的互联电网热稳定关键断面识别方法,其中,步骤二的具体内容包括:
根据步骤一潮流计算得到的电网中所有的支路的视在功率Sl,其中,l表示第l支路的编号,并且支路l的两端节点分别用节点i和节点j表示,第l条支路的热稳定裕度ηl根据式(1)计算:
。
3.根据权利要求1所述一种基于谱聚类分析的互联电网热稳定关键断面识别方法,其中,步骤三的具体内容包括:
步骤1)将电网中的所有节点构成无向加权图G的顶点的集合,用V表示,将电网中所有支路构成无向加权图G的边的集合,用E表示,对于存在并列双回或多回线路的情况,用一条边表示,即节点i和节点j之间存在两条或两条以上的支路,则在无向加权图G中,仅用一条边表示;
步骤2)计算无向加权图G中各边的权重矩阵W,其中,W矩阵的第i行第j列元素用Wij表示,计算公式如下:
步骤3)建立无向加权图G(V,E)的连接矩阵A、度矩阵D和拉普拉斯矩阵L;连接矩阵A的第i行第j列元素用Aij表示;度矩阵D的第i行第j列元素用Dij表示;拉普拉斯矩阵L的第i行第j列元素用Lij表示;具体计算公式如下:
Lij=Dij-Aij (5)。
4.根据权利要求1所述一种基于谱聚类分析的互联电网热稳定关键断面识别方法,其中,步骤五中,特征向量rV2的第i个元素用ci表示,步骤五的具体内容包括:
步骤1)对于第i个节点,若ci大于0,则节点i属于分区1,若若ci小于等于0,则节点i属于分区2;
步骤2)对于支路l,按照式(6)计算指标κl,其中i和j表示支路两端的节点编号,若κl小于0,则支路l属于电网热稳定关键断面,若κl大于等于0,则支路l不属于电网热稳定关键断面;所有满足κl<0的支路构成电网热稳定关键断面,如式(7)所示,
κl=cicj (6)
CritCut={l|κl<0,l=1,2,...,nB} (7)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711248087.4A CN107994569A (zh) | 2017-11-30 | 2017-11-30 | 一种基于谱聚类分析的互联电网热稳定关键断面识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711248087.4A CN107994569A (zh) | 2017-11-30 | 2017-11-30 | 一种基于谱聚类分析的互联电网热稳定关键断面识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107994569A true CN107994569A (zh) | 2018-05-04 |
Family
ID=62035089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711248087.4A Pending CN107994569A (zh) | 2017-11-30 | 2017-11-30 | 一种基于谱聚类分析的互联电网热稳定关键断面识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107994569A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109285089A (zh) * | 2018-10-25 | 2019-01-29 | 天津大学 | 一种电力系统热稳定安全关键机组的筛选方法 |
CN113204854A (zh) * | 2021-06-03 | 2021-08-03 | 广西师范大学 | 一种基于发电机节点及网络加权拓扑的电网分区方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104077393A (zh) * | 2014-07-01 | 2014-10-01 | 武汉大学 | 一种基于半监督谱聚类的最优解列断面搜索方法 |
CN106127237A (zh) * | 2016-06-17 | 2016-11-16 | 东北电力大学 | 基于谱聚类的含vsc‑hvdc交直流系统最优解列断面搜索方法 |
-
2017
- 2017-11-30 CN CN201711248087.4A patent/CN107994569A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104077393A (zh) * | 2014-07-01 | 2014-10-01 | 武汉大学 | 一种基于半监督谱聚类的最优解列断面搜索方法 |
CN106127237A (zh) * | 2016-06-17 | 2016-11-16 | 东北电力大学 | 基于谱聚类的含vsc‑hvdc交直流系统最优解列断面搜索方法 |
Non-Patent Citations (2)
Title |
---|
杨文武 等: ""基于潮流转移的输电断面快速搜索"", 《四川电力技术》 * |
赵金利 等: ""电力系统电压稳定分区和关键断面的确定"", 《电力系统自动化》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109285089A (zh) * | 2018-10-25 | 2019-01-29 | 天津大学 | 一种电力系统热稳定安全关键机组的筛选方法 |
CN113204854A (zh) * | 2021-06-03 | 2021-08-03 | 广西师范大学 | 一种基于发电机节点及网络加权拓扑的电网分区方法 |
CN113204854B (zh) * | 2021-06-03 | 2022-11-29 | 广西师范大学 | 一种基于发电机节点及网络加权拓扑的电网分区方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107482626A (zh) | 一种区域电网关键节点识别方法 | |
CN105912835B (zh) | 一种计及风电场接入的电力线路脆弱性辨识模型建立方法 | |
CN107565601A (zh) | 一种分布式光伏电站集群的动态等值建模方法 | |
CN105656031B (zh) | 基于高斯混合分布特征的含风电电力系统安全风险评估方法 | |
CN112149280B (zh) | 含svg的新能源的多馈入系统电网强度获得方法 | |
CN107846014A (zh) | 一种基于随机特征和级数计算的电网概率脆弱性评估方法 | |
CN107947151A (zh) | 一种电力系统关键节点的识别方法 | |
CN109217296A (zh) | 一种基于加权潮流熵和电压稳定的电网关键节点辨识方法 | |
CN109034467A (zh) | 一种电力系统抗灾型骨干网架多目标优化方法 | |
CN107578137A (zh) | 一种考虑最小奇异值灵敏度熵的电网脆弱性评估方法 | |
CN107622360A (zh) | 一种综合考虑主客观因素的关键线路识别方法 | |
CN106096806A (zh) | 一种基于主观与客观权重集成的多馈入直流落点方案评估方法 | |
CN107994569A (zh) | 一种基于谱聚类分析的互联电网热稳定关键断面识别方法 | |
CN102709955B (zh) | 一种基于多断面潮流控制的方法 | |
CN113258607A (zh) | 基于不同渗透下的光伏发电系统动态离散等值模型建立方法 | |
CN106786602B (zh) | 一种配电网潮流计算方法 | |
CN109768543B (zh) | 一种基于混合整数线性规划的弹性保底网架搜索建模方法 | |
CN108899896B (zh) | 一种基于改进benders分解法的供电能力评估方法 | |
CN112653136B (zh) | 一种电力电子多馈入电力系统关键线路识别方法及系统 | |
CN110783913A (zh) | 基于群的计及预想事故集最佳电网拓扑在线优化方法 | |
CN110120673A (zh) | 基于戴维南等值参数辨识的分布式输配协同无功优化方法及系统 | |
CN111404163A (zh) | 一种电磁环网开环方法 | |
CN112103950B (zh) | 一种基于改进gn分裂算法的电网分区方法 | |
CN109617079A (zh) | 一种直流配电系统潮流解的存在性与稳定性分析方法 | |
CN112052543B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180504 |