法律状态公告日
法律状态信息
法律状态
2020-03-27
未缴年费专利权终止 IPC(主分类):G06F19/00 授权公告日:20160120 终止日期:20190408 申请日:20130408
专利权的终止
2016-01-20
授权
授权
2013-08-14
实质审查的生效 IPC(主分类):G06F19/00 申请日:20130408
实质审查的生效
2013-07-17
公开
公开
技术领域
本发明涉及一种测风数据处理方法,尤其是涉及一种基于风速相关性的风电场测风仪风速缺失数据插补方法。
背景技术
由于当今世界严重的环境污染及传统能源的稀缺,清洁能源的开发和研究已经显得越来越重要。作为清洁能源之一的风能资源,目前在诸多国家都处于快速发展的阶段。一方面,为了获得更多的风能资源,风能资源评估是建造风电场的第一个关键步骤,在具有开发潜能的风电场中,完备的历史风速数据对于评估长期的平均风速及风功率密度极其重要;另一方面,对基于历史数据的统计学的风速预测模型来说,缺失风速或是不合理的风速显然会对建模产生误导的作用,也就是说可靠的风速数据集对提高风速预测的精度有很大的帮助。在实际环境中,由于某些不可控的原因,例如设备失灵或数据传输错误,风速数据并不总是可获得并且是准确的。缺失的数据应当得到插补,而不合理的数据应当被剔除并插补。
国标GB/T19710-2002是数据检验的主要参考,国标中规定了测风数据集中必须包含完整性与合理性检验两个部分。国标中建议使用备用数据来补全缺失数据。但在实际应用中,对于一个测风塔来说,在一个高度安装多个测风仪是不现实的。因为一个测风塔多个测风仪同时出现缺失数据的可能性很小,因为当某个高度出现缺失数据时,利用不同高度风速的关系来插补数据是当前插补风速数据最主要的办法。
最常用的方法是利用Hellman曲线,即将风速随高度的变化大致看成幂率分布:
>
式中的α为风切变系数,通常令α为1/7,故Hellman幂率函数通常也被称为1/7幂函数。然而这一方法存在很大问题,因为风切变系数并不是一个定值。特别是当风出现反层现象时(较小高度的风速大于较大高度的风速),这一方法插值会出现较大的误差。王有禄等人提出了一种利用一个风场多个测风塔的线性回归模型用于切变系数的确定。F.
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种实时性好、插值精度高的基于风速相关性的风电场测风仪风速缺失数据插补方法。
本发明的目的可以通过以下技术方案来实现:
一种基于风速相关性的风电场测风仪风速缺失数据插补方法,该方法根据垂直高度上不同高度风速的关系,对缺失数据进行插补,具体包括以下步骤:
1)导入目标风电场各高度的历史风速数据;
2)根据历史风速数据获得该风电场每个高度风速的范围、每小时风速的变化范围及两个高度间风速差值的变化范围,并保存;
3)导入实时测量数据,根据步骤2)中各范围判断当前测量数据是否存在不合理的风速数据,若是,则剔除不合理风速数据,并执行步骤5),若否,则执行步骤4);
4)根据时间顺序判断是否存在数据缺失,若是,则执行步骤5),若否,则结束;
5)利用相关阶模型对缺失数据进行插补。
所述的不合理风速数据是指该风速数据满足以下任一个条件:
a)超出步骤2)中每个高度风速的范围;
b)超出步骤2)中每小时风速的变化范围;
c)超出步骤2)中两个高度间风速差值的变化范围。
所述的相关阶模型为:
>
>
>
>
高度h0上时间段k1到k2的数据为缺失数据块,v(h0,k)为高度h0上第k个时间的数据,且k∈[k1,k2],在k1-1到k1+δ1以及k2+1到k2+δ2时间段数据是完整的,
>
>
m(·)、ρ(·)、σ(·)分别为向量的均值、标准差及两个向量的相关系数,δj为参数,j∈{1,2},M为可利用的参考高度的个数。
与现有技术相比,本发明具有以下优点:
1.实时性好,因为相关阶模型中的参数并不是固定的,而是与其它高度的风速信息有关。这就与传统的幂率函数插值的方法不同,即便出现反层现象,相关阶模型也可得到有较好的效果。
2.插值的精度高,实验表明,当缺失风速的个数小于9时,相关阶模型的绝对误差百分比(MAPE)会比1/7幂率分布的方法提高8%左右,均方根误差(RMSE)提高0.6m/s。
附图说明
图1为本发明的流程示意图;
图2为缺失数据的示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,一种基于风速相关性的风电场测风仪风速缺失数据插补方法,该方法根据垂直高度上不同高度风速的关系,对缺失数据进行插补,不同高度的风速大小不一致主要是因为地表的粗糙度对其影响不一致找到的,较大高度的风受地表影响较小,而较小高度的风受地表影响较大。该方法具体包括以下步骤:
1)导入目标风电场各高度的历史风速数据;
2)根据历史风速数据获得该风电场每个高度风速的范围、每小时风速的变化范围及两个高度间风速差值的变化范围,并保存;
3)导入实时测量数据,根据步骤2)中各范围判断当前测量数据是否存在不合理的风速数据,若是,则剔除不合理风速数据,并执行步骤5),若否,则执行步骤4);
4)根据时间顺序判断是否存在数据缺失,若是,则执行步骤5),若否,则结束;
5)利用相关阶模型对缺失数据进行插补。
所述的不合理风速数据是指该风速数据满足以下任一个条件:
a)超出步骤2)中每个高度风速的范围;
b)超出步骤2)中每小时风速的变化范围;
c)超出步骤2)中两个高度间风速差值的变化范围。
所述的相关阶模型为:
>
>
>
>
高度h0上时间段k1到k2的数据为缺失数据块,v(h0,k)为高度h0上第k个时间的数据,即为求得的插补的数据,且k∈[k1,k2],为自然数据,在k1-1到k1+δ1以及k2+1到k2+δ2时间段数据是完整的,如图2所示,
>
>
m(·)、ρ(·)、σ(·)分别为向量的均值、标准差及两个向量的相关系数,δj为参数,j∈{1,2},M为可利用的参考高度的个数。
参数δj的确定方法为:两个高度间风速的相关系数
对于一个测风塔来说,通常测风仪可以安装在70m、60m、50m、30m、10m这五个高度上。假设一个高度的风速出现缺失,可以使用4个其他高度的风速信息进行插补,为测风数据集中必须包含完整性与合理性提供了基础。
机译: 风向/风速监视装置,风向/风速监视方法,风向/风速监视程序,记录程序的记录介质,风向/风速监视系统以及风向/风速测量仪
机译: 在高风速条件下运行风电场的方法,涉及通过发电技术来运行风电场
机译: 用于校正风向和风速的系统,用于提示风向和风速的系统,用于校正风向和风速的方法,用于指示风向和风速的方法以及程序