一种基于多参考站的海上测风塔数据代表年订正方法
技术领域
本发明属于风力发电计算技术领域,尤其涉及一种基于多参考站的海上测 风塔数据代表年订正方法。
背景技术
近年来我国海洋风能产业快速发展。研究和摸清海洋风能资源状况对于海 上风电场的选址、规划等工作具有重要价值。通常方法是收集海上测风塔数据 来分析和评估风能资源。但由于海上测风塔的数据时间一般仅有1-2年,而海 上风电场的设计运行期一般为25年,因此需收集周围参考站的长期数据,来对 海上测风塔数据进行代表年订正。
实际工作中可选择的参考站来源非常多,例如当地的气象站可作为参考站 可提供长期数据,另外再分析资料集CFSR、MERRA2、ERA5等也可作为参考站提 供长期数据,还有商业公司3TIER、Vortex、AWS等中尺度和小尺度资料集均可 作为参考站提供长期数据。这些不同来源的参考站各有互补,在实际工作中均 有大量应用。然而由于现有的技术是基于唯一参考站的长期数据来对海上测风 塔数据代表年订正,首先对唯一参考站与海上测风塔的同期数据进行相关性分 析,然而建立一元一次线性方程,进而求出订正量,最终完成代表年订正,具 体见《风电场风能资源评估方法》GB/T 18710-2002、《风电场工程风能资源测 量与评估技术规范》NB/T 31147-2018等标准规范。
由于现有的技术无法充分、合理利用不同来源的多个参考站,一方面造成 了数据资源的浪费,另一方面也会导致海上测风塔数据代表年订正的结果出现 较大误差。因此尽快提出基于多参考站的海上测风塔数据代表年订正方法的意 义重大。
发明内容
本发明的目的是提供一种基于多参考站的海上测风塔数据代表年订正方法, 以解决海上测风塔数据代表年订正的结果出现较大误差的问题。
本发明提供了一种基于多参考站的海上测风塔数据代表年订正方法,包括 如下步骤:
步骤一,获取海上测风塔单层高度的风速数据、风向数据,所述风速、风 向数据为逐10分钟的风速、风向数据,时间长度为至少1年;获取海上测风塔 半径100公里范围内的多个参考站单层高度的风速数据、风向数据,所述风速、 风向数据为逐1小时的风速、风向数据,时间长度为至少30年,且与海上测风 塔的同期时间在1年以上;
步骤二,根据步骤一获取的风速、风向数据,对海上测风塔的逐10分钟的 风速、风向数据,顺序按照每小时0分、10分、20分、30分、40分、50分进 行等权重加权平均,得到海上测风塔的逐1小时的风速、风向数据;
对每个参考站,依次计算参考站逐1小时的风速与海上测风塔的逐1小时 的风速的同期时间1年以上数据之间的相关系数;对于相关系数低于0.6情况, 则剔除对应参考站及数据;对于相关系数大于和等于0.6情况,则保留对应参考 站及数据;
对经上述步骤剔除后的多个参考站的逐1小时的风速、风向数据,采用下 式归一化:
其中,X
t是变量X在t小时的数据,
是X
t归一化后的数值,X
min、X
max、 X
median、X
mean分别是X
t中的最小值、最大值、中位值和平均值;
对海上测风塔逐1小时的风速Vt、风向Dt数据,以及经上述步骤归一化后 的多个参考站的逐1小时的风速、风向数据,建立数据库,并对其中同期时间1 年以上数据按时间段前0~50%的数据划分为训练集,将时间段后50%~100%的 数据划分为测试集;
步骤三,根据步骤二得到的训练集和测试集数据,对归一化后的多个参考 站的逐1小时的风速Vi,t、风向Di,t建立联合场矩阵Ui,t,其中i为参考站的编号; 采用经验正交函数对联合场矩阵Ui,t进行分解,提取出空间模态和对应时间主成 分Pj,t,其中j是空间模态编号,总数量为n的空间模态的合计方差贡献率至少 达到90%;
步骤四,根据步骤二得到的训练集的海上测风塔逐1小时的风速Vt、风向 Dt数据,和步骤三得到的训练集归一化后的多个参考站的数据经验正交函数分 解得到的主成分Pj,t数据,采用下式建立回归方程:
其中,αj、βj和χj分别是回归方程的系数,ε和δ分别是回归拟合残差;
步骤五,根据步骤三获得的测试集归一化后的多个参考站的数据经验正交 函数分解的得到的主成分Pj,t数据,利用步骤四建立的回归方程,进行代入计算, 从而得到多个参考站回归结果的海上测风塔逐1小时的风速Vt′、风向Dt′数据。 将其与步骤二获得的测试集的海上测风塔逐1小时的风速Vt、风向Dt数据进行 对比,利用下式计算出Vt′和Vt之间的均方根误差E:
其中,m是测试集数据的时间总长度;
若均方根误差E大于等于0.5m/s,则返回步骤三,对空间模态总数量n增 加1;然后执行步骤四,建立回归方程;最后执行步骤五,进行误差检验;若均 方根误差E小于0.5m/s,执行步骤六;
步骤六,根据步骤三得到的归一化后的多个参考站的逐1小时的风速、风向 数据,时间长度为至少30年;利用步骤四建立的回归方程,进行代入计算,从 而得到多个参考站回归的海上测风塔逐1小时的风速Vt′、风向Dt′数据,时间长度 为至少30年;最后按照自然年选取1年的风速、风向数据,要求该年的年平均风 速与时间长度为至少30年的风速Vt′平均值之间的绝对误差在所有自然年中为最 小,即为代表年订正结果。
借由上述方案,通过基于多参考站的海上测风塔数据代表年订正方法,可 以解决海上测风塔数据代表年订正的结果出现较大误差的问题。该方法具有通 用性,适用于各种来源的参考站及数据,适用于各种离岸距离的近海和远海的 海上测风塔及数据。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术 手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例详细说明 如后。
附图说明
图1是本发明基于多参考站的海上测风塔数据代表年订正方法的流程图;
图2是本发明一实施例中测试集2014年1月1日至2014年12月31日多 个参考站回归结果的海上测风塔逐1小时的风速Vt′数据与步骤二获得的测试 集2014年1月1日至2014年12月31日的海上测风塔逐1小时的风速Vt数据 的对比图;
图3是本发明一实施例中多个参考站回归的海上测风塔逐1小时的风速 Vt′在1986年1月1日至2015年12月31日的各自然年平均风速。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以 下实施例用于说明本发明,但不用来限制本发明的范围。
参图1所示,本实施例提供了一种基于多参考站的海上测风塔数据代表年订 正方法,包括如下步骤:
步骤一,数据采集。获取海上测风塔单层高度的风速数据、风向数据,所 述风速、风向数据为逐10分钟的风速、风向数据,时间长度为至少1年;获取 海上测风塔半径100公里范围内的多个参考站单层高度的风速数据、风向数据, 所述风速、风向数据为逐1小时的风速、风向数据,时间长度为至少30年,且 与海上测风塔的同期时间在1年以上;
步骤二,数据处理。根据步骤一获取的风速、风向数据,首先对海上测风 塔的逐10分钟的风速、风向数据,顺序按照每小时0分、10分、20分、30分、 40分、50分进行等权重加权平均,得到海上测风塔的逐1小时的风速、风向数 据;
然后对每个参考站,依次计算参考站逐1小时的风速与海上测风塔的逐1 小时的风速的同期时间1年以上数据之间的相关系数;对于相关系数低于0.6情 况,则剔除对应参考站及数据;对于相关系数大于和等于0.6情况,则保留对应 参考站及数据;
之后对经上述步骤剔除后的多个参考站的逐1小时的风速、风向数据,采 用下式归一化:
其中,X
t是变量X在t小时的数据,
是X
t归一化后的数值,X
min、X
max、 X
median、X
mean分别是X
t中的最小值、最大值、中位值和平均值;
最后对海上测风塔逐1小时的风速Vt、风向Dt数据,以及经上述步骤归一 化后的多个参考站的逐1小时的风速、风向数据,建立数据库,并对其中同期 时间1年以上数据按时间段前0~50%的数据划分为训练集,将时间段后 50%~100%的数据划分为测试集;
步骤三,经验正交函数分解。根据步骤二得到的训练集和测试集数据,对 归一化后的多个参考站的逐1小时的风速Vi,t、风向Di,t建立联合场矩阵Ui,t,其 中i为参考站的编号;采用经验正交函数对联合场矩阵Ui,t进行分解,提取出空 间模态和对应时间主成分Pj,t,其中j是空间模态编号,总数量为n的空间模态 的合计方差贡献率至少达到90%;
步骤四,建立回归方程。根据步骤二得到的训练集的海上测风塔逐1小时 的风速Vt、风向Dt数据,和步骤三得到的训练集归一化后的多个参考站的数据 经验正交函数分解得到的主成分Pj,t数据,采用下式建立回归方程:
其中,αj、βj和χj分别是回归方程的系数,ε和δ分别是回归拟合残差;
步骤五,误差检验。根据步骤三获得的测试集归一化后的多个参考站的数 据经验正交函数分解的得到的主成分Pj,t数据,利用步骤四建立的回归方程,进 行代入计算,从而得到多个参考站回归结果的海上测风塔逐1小时的风速Vt′、 风向Dt′数据。将其与步骤二获得的测试集的海上测风塔逐1小时的风速Vt、风 向Dt数据进行对比,利用下式计算出Vt′和Vt之间的均方根误差E:
其中,m是测试集数据的时间总长度;
若均方根误差E大于等于0.5m/s,则返回步骤三,对空间模态总数量n增 加1;然后执行步骤四,建立回归方程;最后执行步骤五,进行误差检验;若均 方根误差E小于0.5m/s,执行步骤六;
步骤六,完成代表年订正。根据步骤三得到的归一化后的多个参考站的逐1 小时的风速、风向数据,时间长度为至少30年;利用步骤四建立的回归方程, 进行代入计算,从而得到多个参考站回归的海上测风塔逐1小时的风速Vt′、风向 Dt′数据,时间长度为至少30年;最后按照自然年选取1年的风速、风向数据,要 求该年的年平均风速与时间长度为至少30年的风速Vt′平均值之间的绝对误差在 所有自然年中为最小,即为代表年订正结果。
通过基于多参考站的海上测风塔数据代表年订正方法,可以解决海上测风 塔数据代表年订正的结果出现较大误差的问题。该方法具有通用性,适用于各 种来源的参考站及数据,适用于各种离岸距离的近海和远海的海上测风塔及数 据。
下面通过具体实例对本发明作进一步说明。
某海上测风塔离岸距离5km、水深12m,塔高90m。
步骤一,数据采集。获取海上测风塔90m高度的风速数据、风向数据,时 间长度为2013年1月1日至2014年12月31日,满足至少1年。获取海上测 风塔半径100公里范围内的响水气象站的10m高度、CFSR的4个格点10m高 度、MERRA2的4个格点50m高度、ERA5的4个格点100m高度、3TIER的4 个格点80m高度、Vortex的4个格点100m高度、AWS的4个格点80m高度的共25个参考站风速数据、风向数据。时间长度为1986年1月1日至2015年12 月31日,满足至少30年。
步骤二,数据处理。获得的训练集为2013年1月1日至2013年12月31日 的数据,测试集为2014年1月1日至2014年12月31日的数据。
步骤三,经验正交函数分解。
步骤四,建立回归方程。
步骤五,误差检验。将获得的测试集2014年1月1日至2014年12月31日 多个参考站回归结果的海上测风塔逐1小时的风速Vt′数据,与步骤二获得的测 试集2014年1月1日至2014年12月31日的海上测风塔逐1小时的风速Vt数 据,对比结果见图2。可见Vt′和Vt两者十分接近,均方根误差E为0.18m/s,误 差很小、满足要求。
步骤六,完成代表年订正。获得了多个参考站回归的海上测风塔逐1小时 的风速Vt′、风向Dt′数据,时间长度为30年。其中风速Vt′在1986年1月1日至 2015年12月31日的各自然年平均风速见图3。可见1999年1月1日至1999 年12月31日的年平均风速与时间长度30年的风速Vt′平均值之间的绝对误差 0.04m/s,在所有自然年中为最小,即为代表年订正结果。
以上所述仅是本发明的优选实施方式,并不用于限制本发明,应当指出, 对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还 可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。