CN106886564B - 一种基于空间聚类订正nwp风能图谱的方法及装置 - Google Patents

一种基于空间聚类订正nwp风能图谱的方法及装置 Download PDF

Info

Publication number
CN106886564B
CN106886564B CN201710002216.5A CN201710002216A CN106886564B CN 106886564 B CN106886564 B CN 106886564B CN 201710002216 A CN201710002216 A CN 201710002216A CN 106886564 B CN106886564 B CN 106886564B
Authority
CN
China
Prior art keywords
data
wind
clustering
nwp
correcting
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.)
Active
Application number
CN201710002216.5A
Other languages
English (en)
Other versions
CN106886564A (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.)
New Energy Polytron Technologies Inc
Original Assignee
New Energy Polytron Technologies Inc
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 New Energy Polytron Technologies Inc filed Critical New Energy Polytron Technologies Inc
Priority to CN201710002216.5A priority Critical patent/CN106886564B/zh
Publication of CN106886564A publication Critical patent/CN106886564A/zh
Application granted granted Critical
Publication of CN106886564B publication Critical patent/CN106886564B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/28—Databases characterised by their database models, e.g. relational or object models
    • G06F16/284—Relational databases
    • G06F16/285—Clustering or classification
    • 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/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Probability & Statistics with Applications (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Primary Health Care (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Tourism & Hospitality (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种基于空间聚类订正NWP风能图谱的方法级装置，方法包括对气象站和测风塔的实测数据进行分类以及处理，建立实测数据序列；计算NWP风能图谱，并对NWP网格点进行空间聚类，将风能图谱划分为不同的区域；利用建立的实测数据序列对NWP风能图谱里的每个格点的数据进行订正。本发明用实测数据对数值模拟风速做订正,提高风速的准确率；用K值空间聚类对区域进行划分，使格点风速能找到属性最接近的实测站点进行订正，提高了订正的可靠性和合理性；通过订正，提高了风资源图谱的精确度，为宏观选址提供可靠的依据。

Description

一种基于空间聚类订正NWP风能图谱的方法及装置
技术领域
本发明属于风电功率技术领域，尤其是涉及一种基于空间聚类订正NWP风能图谱的方法及装置。
背景技术
风力发电是我国重要的后续能源之一,大力发展风电对我国能源供应和环境保护极为重要。NWP模式通过大量计算得出的风能图谱为宏观选址提供依据，帮助投资者找到风资源丰富适合开发的区域。解决人工现场踏勘工作量大、耗时长、可靠性低的问题。
目前行业常规的风资源图谱大多是利用全球环流的再分析数据，采用中尺度大气模式,对所在区域进行模拟计算，获得分辨率1～3km范围的风资源图谱。数值模式是在一定的初始和边界条件下，通过数值求解大气运动方程，从而预知大气运动状态的技术。受初始边界条件及地形、粗糙度描述的准确性，次网格尺度物理过程的引入和数值方法自身的误差等因素的限制，数值模式的误差是客观存在的，并具有较大的不确定性。而且中尺度数值模式水平分辨率太低，无法体现出局地地形和粗糙度对风速分布的影响。由于这些原因，现在的风资源图谱存在较大的误差。
数值模式是在一定的初始和边界条件下，通过数值求解大气运动方程，从而预知大气运动状态的技术。受初始边界条件及地形、粗糙度描述的准确性，次网格尺度物理过程的引入和数值方法自身的误差等因素的限制，数值模式的误差是客观存在的，并具有较大的不确定性。而且中尺度数值模式水平分辨率太低，无法体现出局地地形和粗糙度对风速分布的影响。由于这些原因，现在的风资源图谱存在较大的误差。
用实测数据订正风资源风速分布时，研究人员通常找距离最近的气象站数据来进行订正。气象要素是典型的时空数据，在时间和空间上有很强的相关性。用最邻近的点订正只考虑了空间距离，忽略了周围环境差异造成的较近两点气象特征差异大、相关性低。
空间聚类是指将空间数据集中的对象分成由相似对象组成的类。同类中的对象间具有较高的相似度，而不同类中的对象间差异较大。通过空间聚类，将属性相似、空间邻近的点聚成一簇，划为一个区域，用这一区域里的气象站数据订正NWP风速序列。聚类后的这一区域不仅距离相近，气象属性相似度高，从理论上来说同一簇里的测风塔与格点之间风速相关性较高，能有效的订正NWP数值模拟风速序列。
发明内容
有鉴于此，本发明旨在提出一种基于空间聚类订正NWP风能图谱的方法，用实测数据对NWP数值模拟风速序列进行订正后，以减小数值模式产生的误差，极大的提高风能图谱的准确性，为风场前期选址提供依据。
为达到上述目的，本发明的技术方案是这样实现的：
一种基于空间聚类订正NWP风能图谱的方法，包括如下步骤：
(1)对气象站和测风塔的实测数据进行分类以及处理，建立实测数据序列；
(2)计算NWP风能图谱，并对NWP网格点进行空间聚类，将风能图谱划分为不同的区域；
(3)利用步骤(1)建立的实测数据序列对NWP风能图谱里的每个格点的数据进行订正。
进一步的，所述步骤(1)具体包括如下步骤：
(a1)对实测数据进行分类，分为气象站数据和测风塔数据；
(a2)对气象站自动观测数据进行检验，删除掉无效数据；
(a3)对测风塔数据进行完整性和合理性分析，检验出缺测数据和不合理数据，整理出一套连续一年完整的逐小时测风数据；
(a4)测风塔风速序列推长，利用风场附近长期气象站30年的观测数据，将风场测风塔数据推长，补齐历史30年风速序列；
(a5)循环a.1-a.4步骤直至所有实测数据处理完毕。
进一步的，所述步骤(a2)中无效数据包括不刷新死数，即连续6小时数值不变、以及越限数据，即风速值超过40m/s。
进一步的，所述步骤(2)具体包括如下步骤：
(b1)采用NCEP发布的水平分辨率2.5°×2.5°的历史30年全球再分析数据，利用WRF数值模式降尺度，得到目标区域的1km×1km的逐小时风资源数据；
(b2)提取每个格点的10m高度层的30年平均风速V(i,j)，其中i为经度，j为纬度；
(b3)将N个V(i,j)分为k个簇；
(b4)选择k个聚类中心a1,a2,a3……ak,用这k个点的平均风速建立k个空间聚类表V1,V2,V3……Vk；
(b5)按照最小距离法则逐个对样本V(i,j)进行分类,计算其与每个代表点的距离，将其指派到离聚类中心最近的组里；
(b6)用各聚类表计算聚类均值V，并用来作为各聚类新的代表点；
(b7)循环b4、b5直到聚类均值V不变或代表点未发生变化，则说明函数收敛，停止计算；
(b8)将每个聚类表里的元素边界连在一起，将要订正的图谱划分成为k个区域m1,m2,m3……mk；
(b9)变换k值，循环b4～b8，选取效果较好的一次k值进行空间聚类。
进一步的，所述步骤(3)具体包括如下步骤：
(c1)对风能图谱里的每个格点,找到所属区域范围内最近的三个实测站点n1,n2,n3；
(c2)根据站点风速序列与模拟的风速序列的相关关系，分别用这三个站点的风速序列对NWP模拟的风速序列进行订正，得到3组订正结果；
(c3)将3组订正结果根据站点与格点的距离进行距离加权，得到最终的订正结果。
相对于现有技术，本发明所述的一种基于空间聚类订正NWP风能图谱的方法具有以下优势：
(1)本发明用实测数据对数值模拟风速做订正,提高风速的准确率；
(2)本发明用K值空间聚类对区域进行划分,使格点风速能找到属性最接近的实测站点进行订正,提高了订正的可靠性和合理性，通过订正，能使风资源图谱的精确度提高，为宏观选址提供可靠的依据。
本发明的另一目的在于提出一种基于空间聚类订正NWP风能图谱的装置，以提高风能图谱的准确性，为风场前期选址提供依据。
为达到上述目的，本发明的技术方案是这样实现的：
一种基于空间聚类订正NWP风能图谱的装置，包括
用于对气象站和测风塔的实测数据进行分类以及处理，建立实测数据序列的实测数据序列建立装置；
用于计算NWP风能图谱，并对NWP网格点进行空间聚类，将风能图谱划分为不同的区域的空间聚类装置以及风能图谱划分装置；
用于利用实测数据序列对NWP风能图谱里的每个格点的数据进行订正的订正装置。
进一步的，所述实测数据序列建立装置包括
用于对实测数据进行分类，分为气象站数据和测风塔数据的数据分类装置；
用于对气象站自动观测数据进行检验，删除掉无效数据的数据检验装置；
用于对测风塔数据进行完整性和合理性分析，检验出缺测数据和不合理数据，整理出一套连续一年完整的逐小时测风数据的数据整理装置；
用于对测风塔风速序列推长，利用风场附近长期气象站30年的观测数据，将风场测风塔数据推长，补齐历史30年风速序列的数据补齐装置。
进一步的，所述空间聚类装置以及风能图谱划分装置包括
用于采用NCEP发布的水平分辨率2.5°×2.5°的历史30年全球再分析数据，利用WRF数值模式降尺度，得到目标区域的1km×1km的逐小时风资源数据的目标区域获得装置；
用于提取每个格点的10m高度层的30年平均风速V(i,j)，其中i为经度，j为纬度的格点数据提取装置；
用于将N个V(i,j)分为k个簇的划分装置；
用于选择k个聚类中心a1,a2,a3……ak,用这k个点的平均风速建立k个空间聚类表V1,V2,V3……Vk；的空间聚类表建立装置；
用于按照最小距离法则逐个对样本V(i,j)进行分类,计算其与每个代表点的距离，将其指派到离聚类中心最近的组里的分类装置；
用于用各聚类表计算聚类均值V，并用来作为各聚类新的代表点的聚类均值计算装置；
用于循环分类装置以及聚类均值计算装置直到聚类均值V不变或代表点未发生变化，则说明函数收敛，停止计算的聚类循环装置；
用于将每个聚类表里的元素边界连在一起，将要订正的图谱划分成为k个区域m1,m2,m3……mk的图谱划分装置；
用于变换k值，选取效果较好的一次k值进行空间聚类的聚类结果选取装置。
进一步的，所述订正装置包括
用于对风能图谱里的每个格点,找到所属区域范围内最近的三个实测站点n1,n2,n3的实测点查找装置；
用于根据站点风速序列与模拟的风速序列的相关关系，分别用这三个站点的风速序列对NWP模拟的风速序列进行订正，得到3组订正结果的风速序列订正装置；
用于将3组订正结果根据站点与格点的距离进行距离加权，得到最终的订正结果的最终结果订正装置。
所述一种基于空间聚类订正NWP风能图谱的装置与上述一种基于空间聚类订正NWP风能图谱的方法相对于现有技术所具有的优势相同，在此不再赘述。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例所述的一种基于空间聚类订正NWP风能图谱的方法示意图。
具体实施方式
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面将参考附图并结合实施例来详细说明本发明。
如图1所示，本发明是采用下述的技术方案实现的：
A、实测数据处理，所述方法包括：
a.1数据类型分类，分为气象站数据和测风塔数据；
a.2气象站自动观测数据检验，删除掉无效数据，无效数据包括：(1)不刷新死数(连续6小时数值不变)；(2)越限数据(风速值超过40m/s)；
a.3对无效数据进行修订，根据邻近气象站测风数据与该气象站测风数据的相关关系，补齐缺测和无效数据；
a.4根据GB/T18710-2002标准检查测风塔原始数据，对其进行完整性和合理性分析，检验出缺测数据和不合理数据，经过适当处理，整理出一套连续一年完整的逐小时测风数据；
a.5测风塔风速序列推长，利用风场附近长期气象站30年的观测数据，将风场测风塔数据推长，补齐历史30年风速序列；
循环a.1-a.5步骤直至所有实测数据处理完毕。
B、NWP网格点空间聚类：
b.1计算NWP风能图谱。采用NCEP发布的水平分辨率2.5°×2.5°的历史30年全球再分析数据，利用WRF(Weather Research and Forecast Model)数值模式降尺度，得到目标区域的1km×1km的逐小时风资源数据；
b.2提取每个格点的10m高度层的30年平均风速V(i,j)，其中i为经度，j为纬度；
b.3采用K均值(K-Means)空间聚类算法对所有格点进行聚类。将N个V(i,j)分为k个簇，k值需事先设定，k值大小应小于站点个数；
b.4初始化,选择k个聚类中心a1,a2,a3……ak,用这k个点的平均风速建立k个空间聚类表V1,V2,V3……Vk；
b.5按照最小距离法则逐个对样本V(i,j)进行分类,计算其与每个代表点的距离，将其指派到离聚类中心最近的组里；
b.6用各聚类表计算聚类均值V，并用来作为各聚类新的代表点(更新代表点)
b.7循环b5、b6,直到聚类均值V不变或代表点未发生变化，则说明函数收敛，停止计算；
b.8区域划分，将每个聚类表里的元素边界连在一起，将要订正的图谱划分成为k个区域m1,m2,m3……mk；
b.9变换k值，循环b4～b8，选取效果较好的一次k值进行空间聚类。C、区域内实测站点订正：
c.1对风能图谱里的每个格点,找到所属区域范围内最近的三个实测站点n1,n2,n3；
c.2根据站点风速序列与模拟的风速序列的相关关系，分别用这三个站点的风速序列对NWP模拟的风速序列进行订正，得到3组订正结果；
c.3将3组订正结果根据站点与格点的距离进行距离加权，得到最终的订正结果。
本发明还提出一种基于空间聚类订正NWP风能图谱的装置，包括
(1)用于对气象站和测风塔的实测数据进行分类以及处理，建立实测数据序列的实测数据序列建立装置；
所述实测数据序列建立装置包括
用于对实测数据进行分类，分为气象站数据和测风塔数据的数据分类装置；
用于对气象站自动观测数据进行检验，删除掉无效数据的数据检验装置；
用于对测风塔数据进行完整性和合理性分析，检验出缺测数据和不合理数据，整理出一套连续一年完整的逐小时测风数据的数据整理装置；
用于对测风塔风速序列推长，利用风场附近长期气象站30年的观测数据，将风场测风塔数据推长，补齐历史30年风速序列的数据补齐装置
(2)用于计算NWP风能图谱，并对NWP网格点进行空间聚类，将风能图谱划分为不同的区域的空间聚类装置以及风能图谱划分装置；
所述空间聚类装置以及风能图谱划分装置包括
用于采用NCEP发布的水平分辨率2.5°×2.5°的历史30年全球再分析数据，利用WRF数值模式降尺度，得到目标区域的1km×1km的逐小时风资源数据的目标区域获得装置；
用于提取每个格点的10m高度层的30年平均风速V(i,j)，其中i为经度，j为纬度的格点数据提取装置；
用于将N个V(i,j)分为k个簇的划分装置；
用于选择k个聚类中心a1,a2,a3……ak,用这k个点的平均风速建立k个空间聚类表V1,V2,V3……Vk；的空间聚类表建立装置；
用于按照最小距离法则逐个对样本V(i,j)进行分类,计算其与每个代表点的距离，将其指派到离聚类中心最近的组里的分类装置；
用于用各聚类表计算聚类均值V，并用来作为各聚类新的代表点的聚类均值计算装置；
用于循环分类装置以及聚类均值计算装置直到聚类均值V不变或代表点未发生变化，则说明函数收敛，停止计算的聚类循环装置；
用于将每个聚类表里的元素边界连在一起，将要订正的图谱划分成为k个区域m1,m2,m3……mk的图谱划分装置；
用于变换k值，选取效果较好的一次k值进行空间聚类的聚类结果选取装置。
(3)用于利用实测数据序列对NWP风能图谱里的每个格点的数据进行订正的订正装置，所述订正装置包括
用于对风能图谱里的每个格点,找到所属区域范围内最近的三个实测站点n1,n2,n3的实测点查找装置；
用于根据站点风速序列与模拟的风速序列的相关关系，分别用这三个站点的风速序列对NWP模拟的风速序列进行订正，得到3组订正结果的风速序列订正装置；
用于将3组订正结果根据站点与格点的距离进行距离加权，得到最终的订正结果的最终结果订正装置。
本方案通过K均值空间聚类给格点快速找到属性接近的实测站点。在方法上面，应用其他空间聚类方法也可以得结果。但大多只考虑空间目标的邻近性，忽视了空间目标的属性相似性，或者考虑了属性相似性忽视了目标的空间邻近性。相比较本方案兼顾空间位置和属性双重特性。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于空间聚类订正NWP风能图谱的方法，其特征在于：包括如下步骤：
(1)对气象站和测风塔的实测数据进行分类以及处理，建立实测数据序列；
(2)计算NWP风能图谱，并对NWP网格点进行空间聚类，将风能图谱划分为不同的区域；
(b1)采用NCEP发布的水平分辨率2.5°×2.5°的历史30年全球再分析数据，利用WRF数值模式降尺度，得到目标区域的1km×1km的逐小时风资源数据；
(b2)提取每个格点的10m高度层的30年平均风速V(i,j)，其中i为经度，j为纬度；
(b3)将N个V(i,j)分为k个簇；
(b4)选择k个聚类中心a1,a2,a3……ak,用这k个点的平均风速建立k个空间聚类表V1,V2,V3……Vk；
(b5)按照最小距离法则逐个对样本V(i,j)进行分类,计算其与每个代表点的距离，将其指派到离聚类中心最近的组里；
(b6)用各聚类表计算聚类均值V，并用来作为各聚类新的代表点；
(b7)循环b4、b5直到聚类均值V不变或代表点未发生变化，则说明函数收敛，停止计算；
(b8)将每个聚类表里的元素边界连在一起，将要订正的图谱划分成为k个区域m1,m2,m3……mk；
(b9)变换k值，循环b4～b8，选取误差最小的一次k值进行空间聚类；
(3)利用步骤(1)建立的实测数据序列对NWP风能图谱里的每个格点的数据进行订正。
2.根据权利要求1所述的一种基于空间聚类订正NWP风能图谱的方法，其特征在于：所述步骤(1)具体包括如下步骤：
(a1)对实测数据进行分类，分为气象站数据和测风塔数据；
(a2)对气象站自动观测数据进行检验，删除掉无效数据；
(a3)对测风塔数据进行完整性和合理性分析，检验出缺测数据和不合理数据，整理出一套连续一年完整的逐小时测风数据；
(a4)测风塔风速序列推长，利用风场附近长期气象站30年的观测数据，将风场测风塔数据推长，补齐历史30年风速序列；
(a5)循环a.1-a.4步骤直至所有实测数据处理完毕。
3.根据权利要求2所述的一种基于空间聚类订正NWP风能图谱的方法，其特征在于：所述步骤(a2)中无效数据包括不刷新死数，即连续6小时数值不变、以及越限数据，即风速值超过40m/s。
4.根据权利要求1所述的一种基于空间聚类订正NWP风能图谱的方法，其特征在于：所述步骤(3)具体包括如下步骤：
(c1)对风能图谱里的每个格点,找到所属区域范围内最近的三个实测站点n1,n2,n3；
(c2)根据站点风速序列与模拟的风速序列的相关关系，分别用这三个站点的风速序列对NWP模拟的风速序列进行订正，得到3组订正结果；
(c3)将3组订正结果根据站点与格点的距离进行距离加权，得到最终的订正结果。
5.一种基于空间聚类订正NWP风能图谱的装置，其特征在于：包括
用于对气象站和测风塔的实测数据进行分类以及处理，建立实测数据序列的实测数据序列建立装置；
用于计算NWP风能图谱，并对NWP网格点进行空间聚类，将风能图谱划分为不同的区域的空间聚类装置以及风能图谱划分装置；
用于利用实测数据序列对NWP风能图谱里的每个格点的数据进行订正的订正装置。
6.根据权利要求5所述的一种基于空间聚类订正NWP风能图谱的装置， 其特征在于：所述实测数据序列建立装置包括
用于对实测数据进行分类，分为气象站数据和测风塔数据的数据分类装置；
用于对气象站自动观测数据进行检验，删除掉无效数据的数据检验装置；
用于对测风塔数据进行完整性和合理性分析，检验出缺测数据和不合理数据，整理出一套连续一年完整的逐小时测风数据的数据整理装置；
用于对测风塔风速序列推长，利用风场附近长期气象站30年的观测数据，将风场测风塔数据推长，补齐历史30年风速序列的数据补齐装置。
7.根据权利要求5所述的一种基于空间聚类订正NWP风能图谱的装置， 其特征在于：所述空间聚类装置以及风能图谱划分装置包括
用于采用NCEP发布的水平分辨率2.5°×2.5°的历史30年全球再分析数据，利用WRF数值模式降尺度，得到目标区域的1km×1km的逐小时风资源数据的目标区域获得装置；
用于提取每个格点的10m高度层的30年平均风速V(i,j)，其中i为经度，j为纬度的格点数据提取装置；
用于将N个V(i,j)分为k个簇的划分装置；
用于选择k个聚类中心a1,a2,a3……ak,用这k个点的平均风速建立k个空间聚类表V1,V2,V3……Vk；的空间聚类表建立装置；
用于按照最小距离法则逐个对样本V(i,j)进行分类,计算其与每个代表点的距离，将其指派到离聚类中心最近的组里的分类装置；
用于用各聚类表计算聚类均值V，并用来作为各聚类新的代表点的聚类均值计算装置；
用于循环分类装置以及聚类均值计算装置直到聚类均值V不变或代表点未发生变化，则说明函数收敛，停止计算的聚类循环装置；
用于将每个聚类表里的元素边界连在一起，将要订正的图谱划分成为k个区域m1,m2,m3……mk的图谱划分装置；
用于变换k值，选取误差最小的一次k值进行空间聚类的聚类结果选取装置。
8.根据权利要求5所述的一种基于空间聚类订正NWP风能图谱的装置， 其特征在于：所述订正装置包括
用于对风能图谱里的每个格点,找到所属区域范围内最近的三个实测站点n1,n2,n3的实测点查找装置；
用于根据站点风速序列与模拟的风速序列的相关关系，分别用这三个站点的风速序列对NWP模拟的风速序列进行订正，得到3组订正结果的风速序列订正装置；
用于将3组订正结果根据站点与格点的距离进行距离加权，得到最终的订正结果的最终结果订正装置。
CN201710002216.5A 2017-01-03 2017-01-03 一种基于空间聚类订正nwp风能图谱的方法及装置 Active CN106886564B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710002216.5A CN106886564B (zh) 2017-01-03 2017-01-03 一种基于空间聚类订正nwp风能图谱的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710002216.5A CN106886564B (zh) 2017-01-03 2017-01-03 一种基于空间聚类订正nwp风能图谱的方法及装置

Publications (2)

Publication Number Publication Date
CN106886564A CN106886564A (zh) 2017-06-23
CN106886564B true CN106886564B (zh) 2020-02-14

Family

ID=59176651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710002216.5A Active CN106886564B (zh) 2017-01-03 2017-01-03 一种基于空间聚类订正nwp风能图谱的方法及装置

Country Status (1)

Country Link
CN (1) CN106886564B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107943747B (zh) * 2017-11-17 2021-01-08 南京航空航天大学 基于二维导热微分方程对多连通区域自动分解的方法
CN111325376A (zh) * 2018-12-14 2020-06-23 北京金风科创风电设备有限公司 风速预测方法及装置
CN110208251B (zh) * 2019-06-20 2022-01-21 安徽创谱仪器科技有限公司 等离子体发射光谱干扰校正方法
CN111680408A (zh) * 2020-05-26 2020-09-18 中国能源建设集团广东省电力设计研究院有限公司 一种用于海上风电的风资源图谱绘制方法及装置
CN113239318B (zh) * 2021-05-17 2022-12-20 中国气象局乌鲁木齐沙漠气象研究所 一种区域数值预报模式的土壤湿度初值订正方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102005760A (zh) * 2010-11-18 2011-04-06 西北电网有限公司 一种通用风电功率短期预报方法
CN102570453A (zh) * 2012-01-06 2012-07-11 甘肃省电力公司风电技术中心 基于多套数值天气预报源的短期风电功率预测方法及系统
CN103400230A (zh) * 2013-08-08 2013-11-20 上海电机学院 一种风电功率预测系统及方法
CN103514341A (zh) * 2012-06-14 2014-01-15 华锐风电科技(集团)股份有限公司 基于数值天气预报和计算流体动力学的风资源评估方法
CN104112167A (zh) * 2014-06-06 2014-10-22 国家电网公司 可发电风资源分布的获取方法
CN104112062A (zh) * 2014-06-05 2014-10-22 清华大学 基于插值方法的风资源分布的获取方法
CN104574303A (zh) * 2014-12-26 2015-04-29 河海大学 基于空间聚类的机载LiDAR点云地面滤波方法
US9037521B1 (en) * 2015-01-23 2015-05-19 Iteris, Inc. Modeling of time-variant threshability due to interactions between a crop in a field and atmospheric and soil conditions for prediction of daily opportunity windows for harvest operations using field-level diagnosis and prediction of weather conditions and observations and user input of harvest condition states

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102005760A (zh) * 2010-11-18 2011-04-06 西北电网有限公司 一种通用风电功率短期预报方法
CN102570453A (zh) * 2012-01-06 2012-07-11 甘肃省电力公司风电技术中心 基于多套数值天气预报源的短期风电功率预测方法及系统
CN103514341A (zh) * 2012-06-14 2014-01-15 华锐风电科技(集团)股份有限公司 基于数值天气预报和计算流体动力学的风资源评估方法
CN103400230A (zh) * 2013-08-08 2013-11-20 上海电机学院 一种风电功率预测系统及方法
CN104112062A (zh) * 2014-06-05 2014-10-22 清华大学 基于插值方法的风资源分布的获取方法
CN104112167A (zh) * 2014-06-06 2014-10-22 国家电网公司 可发电风资源分布的获取方法
CN104574303A (zh) * 2014-12-26 2015-04-29 河海大学 基于空间聚类的机载LiDAR点云地面滤波方法
US9037521B1 (en) * 2015-01-23 2015-05-19 Iteris, Inc. Modeling of time-variant threshability due to interactions between a crop in a field and atmospheric and soil conditions for prediction of daily opportunity windows for harvest operations using field-level diagnosis and prediction of weather conditions and observations and user input of harvest condition states

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"聚类分析在短期风电功率预测模型中的应用";高阳 等;《电器与能效管理技术》;20150228;第12-15页 *

Also Published As

Publication number Publication date
CN106886564A (zh) 2017-06-23

Similar Documents

Publication Publication Date Title
CN106886564B (zh) 一种基于空间聚类订正nwp风能图谱的方法及装置
Zhong et al. A city-scale estimation of rooftop solar photovoltaic potential based on deep learning
Sun et al. Estimating the spatial distribution of solar photovoltaic power generation potential on different types of rural rooftops using a deep learning network applied to satellite images
CN111427100B (zh) 一种台风中心定位方法、装置及台风路径生成方法
CN105095589B (zh) 一种山区电网风区分布图绘制方法
US20150039228A1 (en) Ultra-short-term forecasting method including real-time monitoring of the effect of upper and lower courses
Al-Yahyai et al. Nested ensemble NWP approach for wind energy assessment
CN103268572A (zh) 一种千万千瓦级大型风电基地测风网络的微观选址方法
Lou et al. CIE Standard Sky classification by accessible climatic indices
CN103955521B (zh) 一种风电场机群划分方法
Nouri et al. Moroccan wind farm potential feasibility. Case study
Amaya-Martínez et al. A statistical analysis of wind speed distribution models in the Aburrá Valley, Colombia
CN110163882A (zh) 一种基于激光点云数据的电力线悬挂点空间坐标提取方法
CN110134907B (zh) 一种降雨缺失数据填补方法、系统及电子设备
Hofierka et al. The spatial distribution of photovoltaic power plants in relation to solar resource potential: the case of the Czech Republic and Slovakia
CN111080003B (zh) 一种适应复杂山地风电场的测风塔规划选址方法
CN111666725A (zh) 一种适应非复杂地形风电场的测风塔规划选址方法及系统
Wang et al. New estimation method of wind power density with three‐parameter Weibull distribution: A case on Central Inner Mongolia suburbs
Weekes et al. Low‐cost wind resource assessment for small‐scale turbine installations using site pre‐screening and short‐term wind measurements
Chen et al. An estimation framework of regional rooftop photovoltaic potential based on satellite remote sensing images.
Özdemir et al. SPAN: An open-source plugin for photovoltaic potential estimation of individual roof segments using point cloud data
CN111815020A (zh) 一种基于太阳辐射气候特征识别的南墙辐射预测方法
CN112001090A (zh) 一种风场数值模拟方法
CN108415101B (zh) 秒级探空资料的稀疏化方法
CN103473476B (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
CB02 Change of applicant information

Address after: No.227, 2f, building 1, building materials City, Xisanqi, Haidian District, Beijing 100096

Applicant after: New energy Polytron Technologies Inc

Address before: 100081, No. 19, building 68, College Road, Haidian District, Beijing, four, 407

Applicant before: Beijing Guo Nengxin system control technique company limited

CB02 Change of applicant information
CB03 Change of inventor or designer information

Inventor after: Xiang Jie

Inventor after: Yong Zheng

Inventor after: Dong Fen

Inventor after: He Jiangfeng

Inventor before: Xiang Jie

Inventor before: Yong Zheng

Inventor before: Dong Fen

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant