一种基于风速日波动特征的超短期风速修订方法、系统

1.本发明属于风电功率预测风速修订领域，尤其涉及一种基于风速日波动特征的超短期风速修订方法、系统。
背景技术：
：：2.现阶段风能成为了替代化石能源在传统能源结构中重要地位的首要选择。然而，风能具有间歇性、波动性与随机性，这三项特征给风电大规模安全并网带来了巨大的挑战。因此，风电功率预测是风电大规模运用的前提保障，风电功率预测的准确度直接影响了电网的稳定运营与有效调度。3.风电功率预测在过去的几十年中受到了研究界的大量关注，相关研究可以根据研究内容的重点分为两个方面，一方面是根据现有的风速预测数据实现预测风速转化为预测风电功率；另一方面是根据风速预测与实际风速间的误差，对风速预测值进行修订，从而实现风电功率预测精度的提升。4.目前风电功率预测的风速预测数据均由nwp（numericalweatherprediction,数值天气预报）提供，但受限于其物理参数化方案不完善、分辨率较低、地形不够精确等因素的影响，导致数值天气预报的预报结果误差较大、预报结果具有不确定性，不能满足风电场对风速预报的要求，因而有必要对nwp提供的预测风速进行修订，进一步提高超短期风速预测的准确度。因此，本发明的研究重点集中于第二个方面，使用机器学习算法对nwp风速预测值进行修订。技术实现要素：5.本发明所要解决的技术问题：提供一种基于风速日波动特征的超短期风速修订方法，将过往风速日数据、对应波动情况与nwp风速预测数据作为风速订正模型输入，并结合经过bo（bayesianoptimization,贝叶斯优化算法）的lstm（longshort-termmemorynetworks,长短期记忆神经网络）算法实现对nwp风速预测值的修订，能有效提高风速预测值的准确度。6.本发明为解决以上技术问题而采用以下技术方案：本发明提出的一种基于风速日波动特征的超短期风速修订方法，包括以下步骤：s1、收集风场scada系统记录的历史实测风速数据与nwp模型输出的对应风速预测数据。7.s2、利用风速波动特征统计方法和k-means聚类分析算法，构建vsda模型（volatilitystableday,风速日波动特征提取模型）。8.s3、将历史实测风速数据输入vsda模型，进行风速日分段、提取风速日数据段中风速波动特征，并根据风速日波动特征对风速段进行标记，为风速修订提供风速间的相关性与风速波动特征。9.s4、考虑到风速的随机性与间歇性使得风速数据非线性极强，将经过vsda模型生成的风速日数据与对应的波动标记和nwp风速预测数据作为风速订正模型输入，构建lstm神经网络模型，并采用贝叶斯优化算法对lstm神经网络的各个层数、神经元与训练步骤进行优化，建立bo-lstm网络。10.进一步的，步骤s3中，对风速段进行标记的具体步骤为：s301、以采样窗宽为24小时（96个数据点）、采样间隔为4小时（16个数据点）为采样标准，利用vsda模型对历史实测风速数据进行分段，得到大量风速日数据段。11.s302、将分段所得风速日数据段视作为一种时间序列数据，利用vsda模型提取风速脉冲强度和相似熵两类波动特征，描述过往风速时间序列数据的波动情况。12.s303、在一定高度上，用风速脉动值的均方根表示风速脉冲强度，反映了该时段风速波动的振幅情况，具体公式为：；其中，为脉冲强度，u为瞬时风速，为平均风速，为风速脉动值，n为采集风速的样本数。13.s304、相似熵是一种全新的时间序列复杂性测度方法，该方法把皮尔逊相关系数作为衡量时间序列内部相关性的手段。不同于样本熵中采用欧氏距离作为评判指标，相似熵采用的皮尔逊相关系数的取值在-1与1之间，因而可更加直观表现时间序列内部真实存在的相关性。相似熵在vsda模型中表征风速日数据的时序复杂度，反映风速的无序变化程度，提取相似熵的具体过程如下：s3041、构建一组维数为m的向量：；其中，x(i)为第i时刻的风速，i=1,2,…,n-m+1；n为风速段总时刻长度；m为维数，通常取4。14.定义向量x(i)与x(j)之间的相关系数为：；其中，；s3042、计算每个i值下的数目，即为，具体公式为：。15.其中，为相关性阈值，取值范围为0.4~0.6；num()为计数函数；i=1,2,…,n-m+1。16.s3043、对取平均值，记为，具体公式为：。17.s3044、将维数m增加1，重复步骤s3041‑ꢀs3043，计算的平均值，具体公式为：。18.s3045、计算时间序列相似熵，具体公式为：。19.s305、使用k-means聚类分析算法，根据各个数据段的两类波动特征进行聚类处理，得出风速段区中的平稳日与波动日的分类标准，将平稳日与波动日分别标记为0与1，以此作为风速波动标记，用于反映风速波动情况。20.进一步的，还包括步骤s5：利用建立的bo-lstm网络对预测数据进行修订，得到修订后的风速预测数据。21.进一步的，本发明还提出了一种基于风速日波动特征的超短期风速修订系统，包括：数据收集模块，用于收集风场历史实测风速数据与对应数值天气预报风速预测数据。22.风速日波动特征标记模块，用于利用风速波动特征统计方法和k-means聚类分析算法，构建vsda模型；将历史实测风速数据输入vsda模型，进行风速日分段、提取风速日数据段中风速波动特征，并对风速段进行标记。23.lstm模型优化模块，用于将经过vsda模型生成的风速日数据、对应的波动标记与数值天气预报风速预测数据作为风速订正模型输入，构建lstm模型，并采用贝叶斯优化算法对lstm模型的各个层数、神经元与训练步骤进行优化。24.进一步的，数据收集模块中，被配置以执行以下动作：风场历史实测风速数据来自scada系统，天气预报风速预测数据由nwp模型输出。25.进一步的，风速日波动特征标记模块中，被配置以执行以下动作：步骤1、以采样窗宽为24小时、采样间隔为4小时作为采样标准，利用vsda模型对历史实测风速数据进行分段，得到大量风速日数据段。26.步骤2、将分段所得风速日数据作为一种时间序列数据，利用vsda模型提取风速脉冲强度和相似熵两类波动特征。27.步骤3、在一定高度上，用风速脉动值的均方根表示风速脉冲强度，具体公式为：；其中，为脉冲强度，u为瞬时风速，为平均风速，为风速脉动值，n为采集风速的样本数。28.步骤4、提取相似熵的具体过程如下：步骤401、构建一组维数为m的向量：；其中，x(i)为第i时刻的风速，i=1,2,…,n-m+1，n为风速段总时刻长度。29.定义向量x(i)与x(j)之间的相关系数为：。30.其中，；步骤402、计算每个i值下的数目，即为，具体公式为：；其中，为相关性阈值，num()为计数函数，i=1,2,…,n-m+1。31.步骤403、对取平均值，记为，具体公式为：。32.步骤404、将维数m增加1，重复步骤s3041‑ꢀs3043，计算的平均值，具体公式为：。33.步骤405、计算时间序列相似熵，具体公式为：。34.步骤5、使用k-means聚类分析算法，根据各个数据段的两类波动特征进行聚类处理，得出风速段区中的平稳日与波动日的分类标准，将平稳日与波动日分别标记为0与1，以此作为风速波动标记，用于反映风速波动情况。35.进一步的，本发明还提出一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现前文所述的基于风速日波动特征的超短期风速修订方法的步骤。36.进一步的，本发明还提出一种计算机可读的存储介质，所述计算机可读的存储介质存储有计算机程序，所述计算机程序被处理器运行时执行前文所述的基于风速日波动特征的超短期风速修订方法。37.本发明采用以上技术方案，与现有技术相比，其显著技术效果如下：本发明提出的基于风速日波动特征的超短期风速修订方法（vsda-bo-lstm），通过vsda模型有效提取风速的波动特征与历史相关性并将其输出结果，即风速波动标记与风速日数据作为风速修订模型的部分输入；利用贝叶斯优化算法优化lstm神经网络提高修订结果的可靠度。vsda模型与贝叶斯优化算法均可有效地提高风速预测值的准确度，使得直接输入风电功率预测模型的风速预测值可靠性大幅提升，增强了风电功率预测的准确度。附图说明38.图1是本发明的整体流程图。39.图2是本发明实施例中经过vsda模型处理后的数据分布情况。40.图3是vsda-bo-lstm模型的2020年4月10-12日修订效果图。41.图4是vsda-bo-lstm模型的2020年7月10-12日修订效果图。42.图5是vsda-bo-lstm模型的2020年10月10-12日修订效果图。43.图6是vsda-bo-lstm模型的2021年1月10-12日修订效果图。44.图7是nwp风速预测值与vsda-bo-lstm模型、vsda-lstm模型以及bo-lstm模型对四组数据的风速修订值与实际风速值的绝对误差分布图。具体实施方式45.下面结合附图和实例对本发明的技术方案作进一步的说明。46.为实现上述目的，本发明提出了一种基于风速日波动特征的超短期风速修订方法（vsda-bo-lstm），如图1所示，包括以下步骤：s1、收集风场scada系统记录的历史实测风速数据与nwp模型输出的对应风速预测数据。47.s2、利用风速波动特征统计方法和k-means聚类分析算法，构建vsda模型。48.s3、将历史实测风速数据输入vsda模型，进行风速日分段、提取风速日数据段中风速波动特征，并对风速日波动特征进行标记，具体步骤为：s301、以采样窗宽为24小时（96个数据点）、采样间隔为4小时（16个数据点）为采样标准，利用vsda模型对历史实测风速数据进行分段，得到大量风速日数据段。49.s302、将分段所得风速日数据作为一种时间序列数据，利用vsda模型提取风速脉冲强度和相似熵两类波动特征，描述过往风速时间序列数据的波动情况。50.s303、在一定高度上，用风速脉动值的均方根表示风速脉冲强度，反映了该时段风速波动的振幅情况，具体公式为：；其中，为脉冲强度，u为瞬时风速，为平均风速，为风速脉动值，n为采集风速的样本数。51.s304、相似熵在vsda模型中表征风速日数据的时序复杂度，反映风速的无序变化程度，提取相似熵的具体过程如下：s3041、构建一组维数为m的向量：；其中，x(i)为第i时刻的风速，i=1,2,…,n-m+1，n为风速段总时刻长度。52.定义向量x(i)与x(j)之间的相关系数为：。53.其中，；s3042、计算每个i值下的数目，即为，具体公式为：。54.其中，为相关性阈值，取值范围为0.4~0.6；num()为计数函数；i=1,2,…,n-m+1。55.s3043、对取平均值，记为，具体公式为：。56.s3044、将维数m增加1，重复步骤s3041-s3043计算得到的平均值，具体公式为：。57.s3045、计算时间序列相似熵，具体公式为：。58.s305、使用k-means聚类算法，根据各个数据段的两类波动特征进行聚类处理，得出风速段区中的平稳日与波动日的分类标准，将平稳日与波动日分别标记为0与1，以此作为风速波动标记，用于反映风速波动情况。59.s5、利用步骤s4建立的优化lstm模型对预测数据进行修订，得到修订后的风速预测数据。60.本实施例中，数据来源于江苏如东某风电场2020年4月至2021年1月近十个月的实测数据与nwp预测数据，经过vsda模型处理后共计有529条稳定日数据与502条波动日数据，其分布情况如图2所示。为了检验本专利所提出模型的鲁棒性与各个模块的有效性，进行以下实验，实例验证具体为：一、验证本发明提出的vsda-bo-lstm超短期风速修订模型的鲁棒性为了验证本发明提出的vsda-bo-lstm超短期风速修订模型的鲁棒性，本实验选取各个季节的数据即4月、7月、10月、次年1月中随机三天的风速预测数据进行修订，将该模型修订后结果与真实值及nwp风速预测值进行对比分析。61.图3~图6为四季风速修订前后对比结果折线图，从图中可以看出，nwp风速预测数据仅反映出了实际风速的变化趋势，却无法反映出实际风速的时序复杂性与波动性，难以体现风速本身具有的随机性、间歇性与波动性的特点。而本发明提出的vsda-bo-lstm超短期风速修订模型在一年四季均有较为优秀的修订效果，具体表现在vsda-bo-lstm模型在实际工程运用中，保留了nwp风速预测的趋势特征的同时大幅提高了预测风速的准确性，使得修订值准确率远高于nwp风速预测结果，更加贴合实际的真实风速值。62.二、验证本发明提出的vsda-bo-lstm超短期风速修订模型中各模块的有效性为了验证所提模型的各个模块的有效性，分别将四组数据的vsda-lstm模型、bo-lstm模型、nwp风速预测值以及vsda-bo-lstm模型的四季风速值修订误差统计对比，同样采用了相对平均绝对误差（rmae）、相对均方根误差（rrmse）和平均绝对误差（mae）的三种误差评价指标来评价每组数据不同模型的预测结果精度。63.图7为四组数据经过各个模型修订后风速值与实际风速值的绝对误差分布图。由图可知vsda-bo-lstm风速修订模型误差均值远低于其他两种模型；同时相比于其他算法，vsda-bo-lstm的误差分布集中情况最好，25%-75%误差分布区间起止点最小且区间长度最短。分析可得，vsda-bo-lstm模型有效地调整了风速修订的误差分布，使得原本呈均匀分布的nwp风速预测值的误差经过修订后分布更为集中。64.表1为nwp模式、vsda-bo-lstm模型、vsda-lstm模型以及bo-lstm模型在四组试验数据下的平均绝对误差（mae）、相对均方根误差（rrmse）以及相对平均绝对误差（rmae）结果对比。65.表1nwp模式及三种模型风速预测修订误差66.从表1中可以看出，经过本发明提出的vsda-bo-lstm超短期风速修订模型修订后的风速数据对比未经修订前数据在准确度方面有巨大的提升，该算法将nwp模式风速数据误差降低超过50%以上；尤其在秋冬两季本文所提出的方法在修订风速预测值的能力上更加有效；同时，与vsda-lstm和bo-lstm模型相比，vsda-bo-lstm模型在mae、rmae与rrmse三项误差评价指标上数值最小。结果也证明了本发明所提出的基于风速日波动特征的超短期风速修订方法（vsda-bo-lstm）对nwp风速预测值的修订效果显著，说明了有效提取风速波动特征的vsda模型对nwp风速预测值的修订效果有着重要提升，使得修订后的风速值更加贴近真实风速，使得风速预报的精确度得到显著提高，证明了本发明所提模型的优越性。67.本发明实施例还提出一种基于风速日波动特征的超短期风速修订系统包括数据收集模块、风速日波动特征标记模块、lstm模型优化模块及可在处理器上运行的计算机程序。需要说明的是，上述系统中的各个模块对应本发明实施例所提供的方法的具体步骤，具备执行方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节，可参见本发明实施例所提供的方法。68.本发明实施例还提出一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序。需要说明的是，上述系统中的各个模块对应本发明实施例所提供的方法的具体步骤，具备执行方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节，可参见本发明实施例所提供的方法。69.本发明实施例还提出一种计算机可读的存储介质，所述计算机可读的存储介质存储有计算机程序。需要说明的是，上述系统中的各个模块对应本发明实施例所提供的方法的具体步骤，具备执行方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节，可参见本发明实施例所提供的方法。70.以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。当前第1页12当前第1页12
技术特征：
1.一种基于风速日波动特征的超短期风速修订方法，其特征在于，包括以下步骤：s1、收集风场历史实测风速数据与对应数值天气预报风速预测数据；s2、利用风速波动特征统计方法和k-means聚类分析算法，构建vsda模型；s3、将历史实测风速数据输入vsda模型，进行风速日分段、提取风速日数据段中风速波动特征，并对风速段进行标记；s4、将经过vsda模型生成的风速日数据、对应的波动标记与数值天气预报风速预测数据作为风速订正模型输入，构建lstm模型，并采用贝叶斯优化算法对lstm模型的各个层数、神经元与训练步骤进行优化。2.根据权利要求1所述的基于风速日波动特征的超短期风速修订方法，其特征在于，步骤s1中，风场历史实测风速数据来自scada系统，天气预报风速预测数据由nwp模型输出。3.根据权利要求1所述的基于风速日波动特征的超短期风速修订方法，其特征在于，步骤s3中，对风速波动特征进行标记的具体步骤为：s301、以采样窗宽为24小时、采样间隔为4小时作为采样标准，利用vsda模型对历史实测风速数据进行分段，得到大量风速日数据段；s302、将分段所得风速日数据作为一种时间序列数据，利用vsda模型提取风速脉冲强度和相似熵两类波动特征；s303、在一定高度上，风速脉动值的均方根表示风速脉冲强度，具体公式为：；其中，为脉冲强度，u为瞬时风速，为平均风速，为风速脉动值，n为采集风速的样本数；s304、提取相似熵的具体过程如下：s3041、构建一组维数为m的向量：；其中，x(i)为第i时刻的风速，i=1,2,
…
,n-m+1，n为风速段总时刻长度；定义向量x(i)与 x(j)之间的相关系数为：；其中，；s3042、计算每个i值下的数目，即为，具体公式为：；其中，为相关性阈值，num( )为计数函数, i=1,2,
…
, n-m+1；s3043、对取平均值，记为，具体公式为：
；s3044、将维数m增加1，重复步骤s3041-s3043，计算的平均值，具体公式为：；s3045、计算时间序列相似熵，具体公式为：；s305、使用k-means聚类分析算法，根据各个数据段的两类波动特征进行聚类处理，得出风速段区中的平稳日与波动日的分类标准，将平稳日与波动日分别标记为0与1，以此作为风速波动标记，用于反映风速波动情况。4.根据权利要求1所述的基于风速日波动特征的超短期风速修订方法，其特征在于，还包括步骤s5：利用步骤s4建立的优化lstm模型对预测数据进行修订，得到修订后的风速预测数据。5.一种基于风速日波动特征的超短期风速修订系统，其特征在于，包括：数据收集模块，用于收集风场历史实测风速数据与对应数值天气预报风速预测数据；风速日波动特征标记模块，用于利用风速波动特征统计方法和k-means聚类分析算法，构建vsda模型；将历史实测风速数据输入vsda模型，进行风速日分段、提取风速日数据段中风速波动特征，并对风速段进行标记；lstm模型优化模块，用于将经过vsda模型生成的风速日数据、对应的波动标记与数值天气预报风速预测数据作为风速订正模型输入，构建lstm模型，并采用贝叶斯优化算法对lstm模型的各个层数、神经元与训练步骤进行优化。6.根据权利要求5所述的基于风速日波动特征的超短期风速修订系统，其特征在于，数据收集模块中，被配置以执行以下动作：风场历史实测风速数据来自scada系统，天气预报风速预测数据由nwp模型输出。7.根据权利要求5所述的基于风速日波动特征的超短期风速修订系统，其特征在于，风速日波动特征标记模块中，被配置以执行以下动作：步骤1、以采样窗宽为24小时、采样间隔为4小时作为采样标准，利用vsda模型对历史实测风速数据进行分段，得到大量风速日数据段；步骤2、将分段所得风速日数据作为一种时间序列数据，利用vsda模型提取风速脉冲强度和相似熵两类波动特征；步骤3、在一定高度上，用风速脉动值的均方根表示风速脉冲强度，具体公式为：；其中，为脉冲强度，u为瞬时风速，为平均风速，为风速脉动值，n为采集风速的样本数；
步骤4、提取相似熵的具体过程如下：步骤401、构建一组维数为m的向量：；其中，x(i)为第i时刻的风速，i=1,2,
…
,n-m+1，n为风速段总时刻长度；定义向量x(i)与 x(j)之间的相关系数为：；其中，；步骤402、计算每个i值下的数目，即为，具体公式为：；其中，为相关性阈值，num( )为计数函数, i=1,2,
…
, n-m+1；步骤403、对取平均值，记为，具体公式为：；步骤404、将维数m增加1，重复步骤s3041
‑ꢀ
s3043，计算的平均值，具体公式为：；步骤405、计算时间序列相似熵，具体公式为：；步骤5、使用k-means聚类分析算法，根据各个数据段的两类波动特征进行聚类处理，得出风速段区中的平稳日与波动日的分类标准，将平稳日与波动日分别标记为0与1，以此作为风速波动标记，用于反映风速波动情况。8.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4中任一项所述方法的步骤。9.一种计算机可读的存储介质，所述计算机可读的存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器运行时执行所述权利要求1至4中任一项所述的方法。

技术总结
本发明公开了一种基于风速日波动特征的超短期风速修订方法、系统，该方法为：通过构建VSDA模型对历史风速数据进行风速日分段、波动特征提取以及分类，为风速修订算法提供风速日数据与风速波动标记，高度概括了风速间的相关性与风速波动特征；而后将过往风速日数据、对应波动标记与数值天气预报风速预测数据作为模型输入，结合经过贝叶斯优化的长短期记忆神经网络算法实现对数值天气预报风速预测值修订。本发明提出的修订方法可有效地提高风速预测值的准确度，使得风速预测值可靠性大幅提升，增强了风电功率预测的准确度，具有一定的实用价值。实用价值。实用价值。


技术研发人员：熊雄 邹瑞麟 王翼虎 陈一鸣 杨旺
受保护的技术使用者：南京信息工程大学
技术研发日：2023.06.05
技术公布日：2023/7/7