CN103106344B - 一种建立电力系统聚类负荷模型的方法 - Google Patents

一种建立电力系统聚类负荷模型的方法 Download PDF

Info

Publication number
CN103106344B
CN103106344B CN201310035067.4A CN201310035067A CN103106344B CN 103106344 B CN103106344 B CN 103106344B CN 201310035067 A CN201310035067 A CN 201310035067A CN 103106344 B CN103106344 B CN 103106344B
Authority
CN
China
Prior art keywords
cluster
subregion
load
contribution degree
initial value
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
Application number
CN201310035067.4A
Other languages
English (en)
Other versions
CN103106344A (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201310035067.4A priority Critical patent/CN103106344B/zh
Publication of CN103106344A publication Critical patent/CN103106344A/zh
Application granted granted Critical
Publication of CN103106344B publication Critical patent/CN103106344B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种建立电力系统聚类负荷模型的方法:先将时序负荷曲线排序得到负荷持续时间曲线,然后根据负荷水平对充裕度指标的贡献度将负荷持续时间曲线分成高贡献度、中等贡献度、低贡献度三个分区,如果分区为高贡献度区,则对该分区采用层次聚类算法选取聚类中心初始值;如果分区为中等贡献度区,则对该分区采用均值-标准差法选取聚类中心初始值;如果分区为低贡献度区,则对该分区根据经验或者随机确定聚类中心初始值;通过定义改进效率指标,将改进效率作为收敛条件从而确定K-均值聚类算法中的聚类数。根据本发明方法得出的聚类负荷模型在用于电力系统充裕度评估时具有很高的计算精度和快速的收敛特性。

Description

一种建立电力系统聚类负荷模型的方法
技术领域
本发明属于电力系统可靠性评估技术领域,具体涉及一种建立电力系统聚类负荷模型的方法。
背景技术
电力系统充裕度评估能够为制定电力系统的发展规划和系统的运行调度提供参考依据。根据系统状态抽取方法的不同,电力系统充裕度评估算法可分为解析法、状态抽样法和状态持续时间抽样法。其中状态持续时间抽样法直接采用年时序负荷曲线作为负荷模型,而解析法和状态抽样法则采用多级负荷模型,多级负荷模型的准确性将直接影响到可靠性评估结果的准确性。
在形成电力系统充裕度评估中的聚类负荷模型时,为了节约聚类时间和提高聚类模型精度,可以先将电力系统的年时序负荷按照从大到小的顺序排列,然后对排序后的负荷进行聚类从而获取多级负荷模型。通常采用K-均值聚类技术来建立与排序后的负荷曲线对应的多级水平负荷模型,然而K-均值聚类算法存在聚类中心初始值和聚类数难以选择的问题。针对聚类中心初始值难以选择的问题,基于层次聚类的K-均值聚类改进算法通过先应用层次聚类算法得到一个初始的划分,并计算每个类内对象的均值从而得到了初始聚类中心;基于均值-标准差的K-均值聚类改进算法则通过利用均值、标准差来选取聚类中心初始值。研究结果表明,如果单独采用其中的某种改进算法来建立电力系统年时序负荷的聚类模型,所得的聚类负荷模型在用于电力系统充裕度评估时模型精度不够高。采用基于层次聚类的K-均值聚类算法虽然可以对高水平部分的负荷进行精确的聚类,但同时也对低水平部分的负荷进行了细致却几乎无益的聚类,而且实际电力负荷模型通常采用一年的8760个小时的负荷水平为原始数据,采用基于层次聚类的K-均值聚类算法对8760个小时的负荷进行聚类计算将耗费大量的计算时间。而基于均值-标准差的K-均值聚类算法在建立初始聚类中心时对分散负荷点的处理比较粗糙,无法对相对稀疏的高水平负荷或者低水平负荷精确的聚类,而这其中高负荷水平部分的粗糙聚类将会使得计算结果出现较大的误差。针对聚类数的选择问题,2011年第39卷第10期《电力系统保护与控制》中的“基于二分法的聚类负荷模型及其在电力系统可靠性评估中的应用”一文提出利用二分法来快速确定聚类数,将显著性水平作为聚类方案好坏的评价标准,没有考虑到聚类数对电力系统充裕度评估计算速度的影响,因此得出的聚类数太多,不利于快速计算得出充裕度指标。
因此,需要一种新的建立电力系统聚类负荷模型的方法以解决上述问题。
发明内容
发明目的:本发明针对现有技术中电力系统聚类负荷模型建立方法在用于电力系统充裕度评估时的缺陷,提供一种计算精度高、收敛速度快的建立电力系统聚类负荷模型的方法。
技术方案:为解决上述技术问题,本发明的建立电力系统聚类负荷模型的方法采用如下技术方案:
一种建立电力系统聚类负荷模型的方法,包括以下步骤:
1)、输入电力系统时序负荷数据;
2)、对所述时序负荷数据进行降序排列得到负荷持续时间曲线;
3)、将步骤2)得到的所述负荷持续时间曲线按照负荷水平对可靠性影响的贡献划分成不同的分区,分区的数目为n;
4)、设定第i个分区的聚类数的初始值和改进效率阈值ηis,其中,i=1,2,…,n;
5)、将第i个分区的聚类数加倍,其中,i=1,2,…,n;
6)、对不同的分区选择相应的方法来选取初始聚类中心初始值,并根据步骤5)所得的聚类数进行聚类;
7)、计算第i个分区的改进效率值ηi,其中,ηi值由下式计算得到:
η i = | x i - x j | x j / | t i - t j | t j
其中,xi、xj为聚类数增加前和增加后的电力系统可靠性指标,ti、tj为聚类数增加前和增加后的指标计算时间;
8)、重复步骤5)-7),直到ηi≤ηis,其中,ηi为第i个分区的改进效率值,ηis为第i个分区的改进效率阈值,i=1,2,…,n;
9)、重复步骤5)-8),直到所有的分区均完成聚类,输出聚类负荷模型。
更进一步的,步骤3)中所述分区的数目n=3,分别为高贡献度区、中等贡献度区和低贡献度区。
更进一步的,所述高贡献度区的负荷为0.8~1,所述中等贡献度区的负荷为0.5~0.8,所述低贡献度区的负荷为0.3~0.5。
更进一步的,步骤4)中所述高贡献度区的聚类数初始值为1,改进效率阈值为η1s=0.12;所述中等贡献度区的聚类数初始值为1,改进效率阈值为η2s=0.15;所述低贡献度区的聚类数初始值为1,改进效率阈值为η3s=0.20。
更进一步的,步骤6)中,如果分区为高贡献度区,则对该分区采用层次聚类算法选取聚类中心初始值;如果分区为中等贡献度区,则对该分区采用均值-标准差法选取聚类中心初始值;如果分区为低贡献度区,则对该分区根据经验或者随机确定聚类中心初始值。
更进一步的,步骤1)中所述时序负荷数据为年时序负荷数据。
有益效果:本发明的建立电力系统聚类负荷模型的方法具有以下效果:
1、充分利用了基于层次聚类的K-均值聚类改进算法聚类精度高的特点,实现了对高水平负荷的细致划分从而保证负荷模型在电力系统充裕度评估中具有较高的计算精度。
2、充分利用了基于均值-标准差的K-均值聚类改进算法聚类效率高的特点,使得聚类负荷模型在具有较高的计算精度的同时提高了聚类效率。
3、克服了采用二分法时所确定的聚类数过多而导致计算时间成本过大的问题;
4、本发明提出的负荷聚类数的确定方法能够在满足一定的计算精度的同时节省时间成本,实现了计算精度和计算速度的综合权衡。
附图说明
图1是本发明的建立电力系统聚类负荷模型的方法的流程示意图。
具体实施方式
根据下述实施例,可以更好地理解本发明。然而,本领域的技术人员容易理解,实施例所描述的内容仅用于说明本发明,而不应当也不会限制权利要求书中所详细描述的本发明。
实施例1
请参阅图1所示,针对IEEE RTS79系统给出的年时序负荷数据,建立电力系统聚类负荷模型的方法的具体步骤如下:
1)输入IEEE RTS79系统给出的年时序负荷数据;
2)对步骤1)中的年时序负荷进行降序排列得到负荷持续时间曲线;
3)将步骤2)中的负荷持续时间曲线分成3个分区,即分区的数量n=3,其中第1个分区为负荷在0.8~1之间,此为高贡献度区;第2个分区为负荷在0.5~0.8之间,此为中等贡献度区;第3个分区为负荷在0.3~0.5之间,此为低贡献度区;
4)3个分区的改进效率阈值分别取为:高贡献度区的改进效率阈值η1s=0.12,中等贡献度区的改进效率阈值η2s=0.15,低贡献度区的改进效率阈值η3s=0.20;各个分区的的聚类数的初始值均取为1;
5)将第i个(i=1,2,3)分区的聚类数加倍;
6)对第1个分区采用层次聚类算法选取聚类中心初始值;对第2个分区采用均值-标准差法选取聚类中心初始值;对第3个分区根据经验指定聚类中心初始值。
7)计算第i个(i=1,2,3)分区的改进效率值ηi
η i = | x i - x j | x j / | t i - t j | t j
其中,xi、xj为聚类数增加前和增加后的电力系统可靠性指标(如电力不足时间期望),ti、tj为聚类数增加前和增加后的指标计算时间。
8)重复步骤5)-7),直到ηi≤ηis,其中,ηi为第i个分区的改进效率值,ηis为第i个分区的改进效率阈值,i=1,2,3;
9)重复步骤5)-8),直到3个分区均完成聚类,输出聚类负荷模型。
计算结果
采用本发明给出的步骤将IEEE RTS79系统提供的年时序负荷进行聚类,得到的聚类结果是:第一个分区和第二个分区在聚类数为8时,满足收敛条件;第三个分区在聚类数为2时,满足收敛条件;总的聚类数为18。最终形成的负荷聚类结果如下表1所示。
表1本发明方法所得的聚类负荷模型
为了更好的说明本发明方法得出的聚类负荷模型在相同的计算时间下具有更高的计算精度,分别采用基于层次聚类的K-均值聚类算法和基于均值-标准差的K-均值聚类算法在聚类数同样为18的情况下进行了聚类,聚类结果分别如表2和表3所示:
表2基于层次聚类的K-均值聚类算法的负荷模型
表3基于均值-标准差的K-均值聚类算法的负荷模型
分别利用本发明方法得出的聚类负荷模型和单独使用基于层次聚类的K-均值聚类算法或基于均值-标准差的K-均值聚类算法得出的负荷模型,利用状态抽样法(最大抽样次数设置为100000次)计算了系统峰荷为2850MW时可靠性指标,计算结果如表4所示。
表4采用不同聚类负荷模型的发电充裕度计算结果
文献R.Billinton,Li Wenyuan.Reliability Assessment of Electric Power Systems UsingMonte Carlo Methods.New York and London:Plenum Press,1994.中给出了系统峰荷为2850MW时,利用状态持续时间抽样法和年时序负荷模型对IEEE RTS79系统的结果:LOLE=9.3716小时/年,EENS=1197.4448MWh/年。将此发电系统充裕度指标计算结果同表4中采用三种聚类负荷模型时得出的结果相比较可知,本发明方法得出的聚类负荷模型具有更高的计算精度和更快的收敛速度。
本发明的建立电力系统聚类负荷模型的方法,克服了单独采用基于层次聚类的K-均值聚类算法和基于均值-标准差的K-均值聚类算法获取电力系统可靠性评估中的聚类负荷模型时的不足,既保证了对高水平负荷的细致划分从而改善了负荷模型的聚类精度,又避免了对低水平负荷的细致划分从而节省了计算时间;克服了采用二分法时所确定的聚类数过多而导致计算时间成本过大的问题;本发明给出的负荷聚类数的确定方法能够在满足一定的计算精度的同时节省时间成本,实现了计算精度和计算速度的综合权衡。

Claims (5)

1.一种建立电力系统聚类负荷模型的方法,其特征在于,包括以下步骤:
1)、输入电力系统的时序负荷数据;
2)、对年时序负荷数据进行降序排列得到负荷持续时间曲线;
3)、将步骤2)得到的所述负荷持续时间曲线按照负荷水平对可靠性影响的贡献划分成不同的分区,分区的数目为n;
4)、设定第i个分区的聚类数的初始值和改进效率阈值ηis,其中,i=1,2,…,n;
5)、将第i个分区的聚类数加倍,其中,i=1,2,…,n;
6)、对不同的分区选择相应的方法来选取初始聚类中心初始值,并根据步骤5)所得的聚类数进行聚类;
7)、计算第i个分区的改进效率值ηi,其中,ηi值由下式计算得到:
其中,xi、xj为聚类数增加前和增加后的电力系统可靠性指标,ti、tj为聚类数增加前和增加后的指标计算时间;
8)、重复步骤5)-7),直到ηi≤ηis,其中,ηi为第i个分区的改进效率值,ηis为第i个分区的改进效率阈值,i=1,2,…,n;
9)、重复步骤5)-8),直到所有的分区均完成聚类,输出聚类负荷模型。
2.如权利要求1所述的建立电力系统聚类负荷模型的方法,其特征在于,步骤3)中所述分区的数目n=3,分别为高贡献度区、中等贡献度区和低贡献度区。
3.如权利要求2所述的建立电力系统聚类负荷模型的方法,其特征在于,步骤4)中所述高贡献度区的聚类数初始值为1,改进效率阈值为η1s=0.12;所述中等贡献度区的聚类数初始值为1,改进效率阈值为η2s=0.15;所述低贡献度区的聚类数初始值为1,改进效率阈值为η3s=0.20。
4.如权利要求2所述的建立电力系统聚类负荷模型的方法,其特征在于,步骤6)中,如果分区为高贡献度区,则对该分区采用层次聚类算法选取聚类中心初始值;如果分区为中等贡献度区,则对该分区采用均值-标准差法选取聚类中心初始值;如果分区为低贡献度区,则对该分区根据经验或者随机确定聚类中心初始值。
5.如权利要求1所述的建立电力系统聚类负荷模型的方法,其特征在于,步骤1)中所述时序负荷数据为年时序负荷数据。
CN201310035067.4A 2013-01-30 2013-01-30 一种建立电力系统聚类负荷模型的方法 Expired - Fee Related CN103106344B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310035067.4A CN103106344B (zh) 2013-01-30 2013-01-30 一种建立电力系统聚类负荷模型的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310035067.4A CN103106344B (zh) 2013-01-30 2013-01-30 一种建立电力系统聚类负荷模型的方法

Publications (2)

Publication Number Publication Date
CN103106344A CN103106344A (zh) 2013-05-15
CN103106344B true CN103106344B (zh) 2015-09-16

Family

ID=48314198

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310035067.4A Expired - Fee Related CN103106344B (zh) 2013-01-30 2013-01-30 一种建立电力系统聚类负荷模型的方法

Country Status (1)

Country Link
CN (1) CN103106344B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104951654A (zh) * 2015-06-05 2015-09-30 华南理工大学 基于控制变量抽样的大规模风电场可靠性评估方法
CN106022578A (zh) * 2016-05-12 2016-10-12 国网天津市电力公司 基于数据高维化和K-means聚类的居民用电峰谷平时段划分方法
CN106055918B (zh) * 2016-07-26 2020-02-21 天津大学 一种电力系统负荷数据的辨识及修正方法
CN106446967A (zh) * 2016-10-28 2017-02-22 国网福建省电力有限公司 一种新型电力系统负荷曲线聚类方法
CN107918779A (zh) * 2017-08-02 2018-04-17 北京国电通网络技术有限公司 一种构建多元负荷聚类模型方法及系统
CN108429257B (zh) * 2018-03-26 2021-01-19 东北电力大学 基于负荷点聚类分区并计及地理因素的配电网网架规划方法
CN110994614A (zh) * 2019-12-26 2020-04-10 南方电网科学研究院有限责任公司 一种直流输电系统的可靠性评估方法、系统及设备
CN111640149A (zh) * 2020-06-01 2020-09-08 中国人民解放军63653部队 一种评估点云建模误差的分析方法
CN113064677B (zh) * 2021-02-23 2023-06-20 青岛海尔科技有限公司 应用运行方法及装置、存储介质及电子装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101789598A (zh) * 2010-03-05 2010-07-28 湖北省电力试验研究院 一种电力系统负荷的建模方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102779223B (zh) * 2011-05-13 2016-03-30 富士通株式会社 短期电力负荷预测的方法及装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101789598A (zh) * 2010-03-05 2010-07-28 湖北省电力试验研究院 一种电力系统负荷的建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于二分法的聚类负荷模型及其在电力系统可靠性评估中的应用;鲁国起等;《电力系统保护与控制》;20110516;第39卷(第10期);第31-36页 *
面向运行规划可靠性的综合聚类负荷模型;李小燕等;《电力系统自动化》;20100425;第34卷(第8期);第56-60页 *

Also Published As

Publication number Publication date
CN103106344A (zh) 2013-05-15

Similar Documents

Publication Publication Date Title
CN103106344B (zh) 一种建立电力系统聚类负荷模型的方法
CN103020423B (zh) 基于copula函数获取风电场出力相关特性的方法
CN102945223B (zh) 一种构建多个风电场出力联合概率分布函数的方法
CN103424654B (zh) 一种敏感设备电压暂降敏感度的评估方法
CN103872395B (zh) 一种锂离子动力电池梯次利用的配组方法
CN104079003B (zh) 含光伏电源配电网的概率潮流计算方法
CN110599006B (zh) 基于场景分析的配电网运行风险评估方法
CN104332665B (zh) 一种电池配组方法
CN103106314B (zh) 太阳能光伏电源输出功率的时序概率建模方法
CN103579700A (zh) 一种锂离子电池分选配组方法
CN105488538A (zh) 一种基于改进k均值聚类算法的天空图像云空辨识方法
CN109921426A (zh) 基于cv-kde的风电并网系统概率潮流计算方法
CN104504455A (zh) 一种梯级蓄能控制下水电站群长期优化调度方法
CN104537233B (zh) 一种基于核密度估计的配电网伪量测生成方法
CN106385039A (zh) 用以平抑光伏功率波动的混合储能系统设计方法
CN106056305B (zh) 一种基于状态聚类的发电系统可靠性快速评估方法
CN102800016B (zh) 基于信用等级与违约损失率匹配的信用评级调整方法
CN103606111A (zh) 一种综合电压合格率的评估方法
CN106094972B (zh) 一种基于函数模型的光伏发电系统最大功率点跟踪方法
CN108694475A (zh) 基于混合模型的短时间尺度光伏电池发电量预测方法
CN107392439A (zh) 一种电力系统可靠性评估的组合抽样方法
CN105095658A (zh) 水文流量波动情势识别方法及系统
CN103326396B (zh) 一种测定电力系统风电运行出力上限值的方法
CN112993984B (zh) 一种分布式光伏电站在有电气连接情况下的大扰动聚合方法
CN104715148B (zh) 一种基于Markov过程的可再生能源电力建模方法

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: 20150916

Termination date: 20200130