基于改进分裂h-k聚类方法的风电场机群划分方法

文档序号:10594338阅读:348来源:国知局
基于改进分裂 h-k 聚类方法的风电场机群划分方法
【专利摘要】一种基于改进分裂H?K聚类方法的风电场机群划分方法,包括:采用分裂层次聚类算法,将风机的参数划分成k个簇,其中k是风电场划分数目;在所有簇中寻找具有最大直径的簇C;找出簇C中与其它点平均相异度最大的一个点p,并将其放入一个新的簇中;在簇C中寻找目标点,该目标点到新的簇中的点的距离不大于到簇C中最近的点,并将目标点放入新的簇中;计算不同k值时的轮廓系数,根据轮廓系数确定最优聚类数;在最优聚类数的状态下,计算出聚类中心来代替k均值聚类算法中的随机聚类中心以进行聚类;计算簇C中的剩余点到各个聚类中心的距离并将各个剩余点分配进入最近的类别中,分配完成后计算新的聚类中心;然后以聚类结果为准建立新的风电场模型进行仿真模拟。
【专利说明】
基于改进分裂H-K聚类方法的风电场机群划分方法
技术领域
[0001] 本发明设及风力发电系统领域,更具体地说,本发明设及一种基于改进分裂H-K 化ierarchical K-means)聚类方法的风电场机群划分方法,其可应用于建立大型风电场等 值模型,运用轮廓系数并加 W合适的准则解决机群划分数目问题,利用本发明提出的聚类 方法完成风电场聚类,最终建立等值模型。
【背景技术】
[0002] 随着风电技术的不断发展和发电成本的不断降低,风力发电技术受到的关注程度 越来越高,但相对于传统能源来说,风电可W说是"劣质电源"因为风电是具有间歇性和随 机波动性的特殊电源,因此其输出功率也具有很强的不确定性。当大规模的风电场并网运 行后会对电网的稳定性和电能质量带来巨大的影响。目前在分析大规模风电场并网后产生 的影响,首先需要解决的问题是大规模风电场建模。早期研究人员是对风电场中的每台风 机进行详细建模,但是运种方法随着风电场规模的不断扩大,逐渐开始变得越来越困难,因 此研究建立一个相对简单并能够刻画风电场的模型变得很重要,即建立风电场对电网的外 特性模型。目前对于风电场动态模型主要使用简化的等值模型主要有两种方式即单机等值 模型和多机等值模型。相对来说多机等值模型的准确度更高。
[0003] 但就目前来说,等值模型的建立都存在相应的问题,对于机群划分数目没有明确 的标准,采用的都是经验方法居多。因此,需要提供一种符合实际的需求的划分方式和简明 的划分指标,确保建立模型的准确性。

【发明内容】

[0004] 本发明所要解决的技术问题是针对现有技术中存在上述缺陷,提供一种基于聚类 算法的大规模风电场建模,通过轮廓系数确定最优的聚类数目,并详细分析风机运行机理 选择合适的聚类指标,在保证模型精确的同时降低多机模型复杂程度。就目前来说关于轮 廓系数在风电场等值建模方面没有详细的描述,在相关方面没有深入探究。
[0005] 为了实现上述技术目的,根据本发明,提供了一种基于改进分裂H-K聚类方法的风 电场机群划分方法,包括:
[0006] 第一步骤:采用分裂层次聚类算法,将风机的参数划分成k个簇,其中k是风电场划 分数目;
[0007] 第二步骤:在所有簇中寻找具有最大直径的簇C用于数据的分裂;
[000引第=步骤:找出簇C中与其它点平均相异度最大的一个点P,并将其放入一个新的 簇中;
[0009] 第四步骤:在簇C中寻找目标点,该目标点到新的簇中的点的距离不大于到簇C中 最近的点,并将目标点放入新的簇中;
[0010] 第五步骤:计算不同k值时的轮廓系数,根据轮廓系数确定最优聚类数;
[0011] 第六步骤:在最优聚类数的状态下,计算出聚类中屯、来代替k-means聚类化均值聚 类)算法中的随机聚类中屯、W进行聚类;
[0012] 第屯步骤:计算簇C中的剩余点到各个聚类中屯、的距离并将各个剩余点分配进入 最近的类别中,分配完成后计算新的聚类中屯、;
[0013] 第八步骤:然后W聚类结果为准建立新的风电场模型进行仿真模拟。
[0014] 优选地,第五步骤中轮廓系数的计算公式为:
[0015]
[0016] 其中m是簇中的元素 Xi到本簇中其它元素之间的平均距离,nil是簇中的元素 Xi到 其它簇中元素平均距离的最小值。
[0017] 优选地,S(i)的值域在[-1,1]之间。
[001引优选地,k为不小于2的整数。
[0019] 优选地,在所有簇中,簇C中存在欧式距离最大的两个点。
[0020] 优选地,在簇C中没有分配到新的簇中的目标点时结束第四步骤并执行第五步骤。
[0021] 根据本发明的基于改进分裂H-K聚类方法的风电场机群划分方法能够适合大规模 风电场等值模型,其基于H-K聚类算法的风电场建模,是将k-means聚类方法和层次聚类方 法在数据处理上的优缺点结合在一起,先使用层次聚类算法获取初始信息,然后使用k- means聚类算法进一步完善聚类过程。在本发明中将H-K聚类方法和使用轮廓系数优化聚类 数目运两者结合到一起,明确划分机群的数目运样可W使建立的风电场模型更加符合实际 情况。
[0022] 本发明基于H-K聚类算法的风电场建模,给予机群划分确立一个明确的标准,W此 为基础界定大规模风电场进行多机等值建模时所需要的等值机型的数目,再次之后选择合 适的聚类指标,运用H-K聚类算法进行风电场聚类建模。确保建模的高效和准确。
【附图说明】
[0023] 结合附图,并通过参考下面的详细描述,将会更容易地对本发明有更完整的理解 并且更容易地理解其伴随的优点和特征,其中:
[0024] 图1示意性地示出了根据本发明优选实施例的基于改进分裂H-K聚类方法的风电 场机群划分方法的流程图。
[0025] 图2示意性地示出了根据本发明优选实施例的按照变量聚类的轮廓值。
[0026] 需要说明的是,附图用于说明本发明,而非限制本发明。注意,表示结构的附图可 能并非按比例绘制。并且,附图中,相同或者类似的元件标有相同或者类似的标号。
【具体实施方式】
[0027] 为了使本发明的内容更加清楚和易懂,下面结合具体实施例和附图对本发明的内 容进行详细描述。
[002引对整体方法进行详细描述,首先在MATLAB中搭建仿真模型,建立一个风电场的详 细模型。在不同风速状况下进行风场的仿真模拟,并保存风电场相关的输出曲线,在深入分 析风机的运行机理后选择合适的指标作为初始数据输入算法中进行详尽的计算,附图1为 本发明中算法部分的流程图,算法按照如下步骤进行处理:
[0029] 第一步骤SI:采用分裂层次聚类算法,将风机的参数划分成k个簇,其中k是风电场 划分数目。其中为了贴合实际并减少运算量,k值可W从2开始选择根据风场规模设置上限; 即,k为不小于2的整数。
[0030] 第二步骤S2:在所有簇中寻找具有最大直径的簇C(即,在所有簇中,簇C中存在欧 式距离最大的两个点)用于数据的分裂。
[0031] 第=步骤S3:找出簇C中与其它点平均相异度最大的一个点P,并将其放入一个新 的簇中。
[0032] 第四步骤S4:在簇C中寻找目标点,该目标点到新的簇中的点的距离不大于到簇C 中最近的点,并将目标点放入新的簇中。在簇C中没有分配到新的簇中的目标点时结束第四 步骤S4并执行第五步骤S5。
[0033] 第五步骤S5:计算不同k值时的轮廓系数,根据轮廓系数确定最优聚类数,其中轮 廓系数的计算公式为:
[0034]
[0035] 其中m是簇中的元素 Xi到本簇中其它元素之间的平均距离,nil是簇中的元素 Xi到 其它簇中元素平均距离的最小值,而且S(i)的值域在[-1,1]之间;
[0036] 轮廓系数越接近于1说明该元素在一个紧密独立的集体中;越接近于-1则说明元 素不属于该簇;若轮廓系数为0则表明样本集中不存在自然簇结构,而是均匀分布的结构。
[0037] 例如,在不同k值状态下,计算各台风机参数相对于本簇中屯、的轮廓值,并且可W 建立一个轮廓值曲线图,W轮廓值之和最大并且相对曲线均比较倾向1的k值为最优聚类 数。
[0038] 第六步骤S6:在最优聚类数的状态下,即确定的最优的风电场划分数目下,计算出 聚类中屯、来代替k-means聚类算法中的随机聚类中屯、W进行聚类。
[0039] 第屯步骤S7:计算簇C中的剩余点到各个聚类中屯、的距离并将各个剩余点分配进 入最近的类别中,分配完成后计算新的聚类中屯、。优选地,若第屯步骤S7计算的聚类中屯、的 发生变化则重新计算聚类中屯、再次进行剩余点分配直至收敛。
[0040] 第八步骤S8:然后W聚类结果为准建立新的风电场模型进行仿真模拟,并且优选 地还可W与之前的所保存的详细模型的输出曲线进行对比分析。
[0041] 下面结合具体事例,对本发明做进一步描述,W24台风机构成一个风电场进行分 析,风机初始数据设定见表1(数据均为标么值),定风速是指在12m/s是风机输出功率,在两 种风速模式运行,第一种为基本风和渐进风的叠加;第二种为波动风速W恒定风速和高斯 白噪声信号合成。选择聚类指标进行聚类,在不同聚类数下的轮廓值如图2所示。
[0042] 表1风机初始参数
[0043]
[0
[0045] 根据图2所示的计算结果显示,聚成3类时结果最优,然后设定等值风机模型共有 为討几等值模型,等值模型参数计算依照等值模式进行。
[0046] 需要说明的是,除非特别说明或者指出,否则说明书中的术语"第一"、"第二"、"第 立"等描述仅仅用于区分说明书中的各个组件、元素、步骤等,而不是用于表示各个组件、元 素、步骤之间的逻辑关系或者顺序关系等。
[0047] 可W理解的是,虽然本发明已W较佳实施例披露如上,然而上述实施例并非用W 限定本发明。对于任何熟悉本领域的技术人员而言,在不脱离本发明技术方案范围情况下, 都可利用上述掲示的技术内容对本发明技术方案作出许多可能的变动和修饰,或修改为等 同变化的等效实施例。因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对 W上实施例所做的任何简单修改、等同变化及修饰,均仍属于本发明技术方案保护的范围 内。
【主权项】
1. 一种基于改进分裂H-K聚类方法的风电场机群划分方法,其特征在于包括: 第一步骤:采用分裂层次聚类算法,将风机的参数划分成k个簇,其中k是风电场划分数 目; 第二步骤:在所有簇中寻找具有最大直径的簇C用于数据的分裂; 第三步骤:找出簇C中与其它点平均相异度最大的一个点p,并将其放入一个新的簇中; 第四步骤:在簇C中寻找目标点,该目标点到新的簇中的点的距离不大于到簇C中最近 的点,并将目标点放入新的簇中; 第五步骤:计算不同k值时的轮廓系数,根据轮廓系数确定最优聚类数; 第六步骤:在最优聚类数的状态下,计算出聚类中心来代替k均值聚类算法中的随机聚 类中心以进行聚类; 第七步骤:计算簇C中的剩余点到各个聚类中心的距离并将各个剩余点分配进入最近 的类别中,分配完成后计算新的聚类中心; 第八步骤:然后以聚类结果为准建立新的风电场模型进行仿真模拟。2. 根据权利要求1所述的基于改进分裂Η-K聚类方法的风电场机群划分方法,其特征在 于,第五步骤中轮廓系数的计算公式为:其中m是簇中的元素xi到本簇中其它兀素之间的肀均跑离,πη是簇中的元素Xi到其它簇 中元素平均距离的最小值。3. 根据权利要求2所述的基于改进分裂Η-K聚类方法的风电场机群划分方法,其特征在 于,S(i)的值域在[_1,1]之间。4. 根据权利要求1或2所述的基于改进分裂Η-K聚类方法的风电场机群划分方法,其特 征在于,k为不小于2的整数。5. 根据权利要求1或2所述的基于改进分裂Η-K聚类方法的风电场机群划分方法,其特 征在于,在所有簇中,簇C中存在欧式距离最大的两个点。6. 根据权利要求1或2所述的基于改进分裂Η-K聚类方法的风电场机群划分方法,其特 征在于,在簇C中没有分配到新的簇中的目标点时结束第四步骤并执行第五步骤。
【文档编号】G06Q50/06GK105956318SQ201610334655
【公开日】2016年9月21日
【申请日】2016年5月19日
【发明人】朱成亮, 刘三明, 王致杰, 殷建炜, 潘磊
【申请人】上海电机学院
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1