复杂地形风电场风速与功率预测方法、系统、设备及存储介质与流程

本发明属于风力发电的能效评估
技术领域：
，具体涉及一种复杂地形风电场风速与功率预测方法、系统、设备及存储介质。
背景技术：
：随着近年来持续进行的风电大规模开发，陆上适合风电开发的风资源丰富、地形简单的区域越来越少，风资源相对丰富、地形复杂的区域成为当前和进后陆上风电开发的重点。复杂地形风电场中高低起伏的地形地貌，加剧了地表对大气边界层流动的撞击、环绕、分离、阻滞作用，很大程度上改变了不同高度层面的平均流速与湍流强度，其风况直接影响着各风电机组的功率输出和发电量。由于风能本身属于随机波动的不稳定能源，加之在复杂地形下，各机组所在位置的风况与风质量存在着一定程度的差异，大规模的风电并入电网系统，必将会对系统的稳定性带来新的挑战。故风电场在生产、运维的过程中，通常需要对风场区域的平均风速进行预测，进而根据机组容量对未来一段时间内风电场所能输出的功率大小进行预测，以便电力生产调度机构需要对未来数小时的风电输出功率有所了解，合理安排调度计划，减轻风电的出力波动对电网电压稳定与功率平衡的影响，同时实现风能的科学与最有利用，提升风电场的运行效益。目前，在役复杂地形风电场的风功率预测，主要是基于风电场所在区域中尺度气象信息中的风速与风向预报数据，利用物理模拟计算和统计分析方法，对装机位置附近的风速进行短期预测，进而结合风电机组的功率曲线计算出实际输出风电功率，预测出整个风电场的功率。然而，由于受到测风装置的安装数量以及维护成本等条件限制，复杂地形风电场通常只在1～2处代表性地理位置进行实地测风，模拟计算和统计分析得到的风电场中各风机点位处的风况信息难以得到验证与校核，加之风机风速和风电功率之间存在非线性关系，使得风速预测值微小误差将产生较大的功率值误差，从而给电网系统的安全性、稳定性和可控性造成一定程度的影响。技术实现要素：针对上述
背景技术：
中提到的难点问题，本发明的目的在于提供一种复杂地形风电场风速与功率预测方法、系统、设备及存储介质，该方法简化了传统物理模拟预测模型的复杂程度，避免了工时较长的在线数值模拟，节省了风电场运维系统中的计算资源，从数理统计的角度为复杂地形风电场基于风机点位的精细化风速、功率预测提供了一种新的解决途径。为了达到上述目的，本发明采用的技术方案是：一种复杂地形风电场风速与功率预测方法，包括以下步骤：获取目标风电场区域轮毂高度处的中尺度气象预报历史数据；由不同来流风向扇区下，风电场区域性中尺度气象历史风速与各机组机舱风速计测得的同期风速数据的相关关系，建立整场气象风速与各机组点位风速之间的拟合函数；根据短期气象预报风速和拟合函数预测出风电场离散机位处的精确风速，并结合各机组实际功率曲线得到功率预测。作为本发明的进一步改进，所述获取目标风电场区域轮毂高度处的中尺度气象预报历史数据具体是：获取目标风电场区域轮毂高度处的中尺度气象预报历史数据，并从中提取风电场运行投产以来至今的平均风速与风向的时间序列数据，建立气象历史风况数据集。作为本发明的进一步改进，对气象历史风况数据进行不同来流风向扇区下的筛选与归类得到每个来流扇区方位下的子数据集；根据获得的每个来流扇区方位下子数据集中的日期和时间序列，获取风电场中各风电机组的同期scada运行数据，并从中提取机舱风速计测量得到风速数据，建立相应的整场机组风速数据集；将每个来流扇区方位下气象子数据集与同期整场风电机组风速数据集合并，建立相应的整场气象预报与机组风速测量数据集；针对整场气象预报与机组风速测量数据集计算得到气象预报各来流风向下，整场气象预报风速与各风机点位风速之间的相关关系；根据各风机点位风速与整场气象预报风速之间的相关关系，建立基于气象预报风向与风速的各机组预测风速函数关系库。作为本发明的进一步改进，对气象历史风况数据进行不同来流风向扇区下的筛选与归类得到每个来流扇区方位下的子数据集，具体包括以下步骤：来流风向扇区依照经典风向玫瑰图，建立十六方位来流风向，每个方位所代表的风向扇区中心角度，及各角度所代表的来流扇区角度范围的相互关系；以气象历史风况数据集m中的气象预报轮毂高度处风向数据为基础，所划分的16方位来流风向进行数据筛选，获得每个来流扇区方位下的子数据集。作为本发明的进一步改进，整场气象预报风速与各风机点位风速之间的相关关系采用函数关系进行拟合，最终根据各拟合函数的均方差与确定系数的优劣程度，建立最优的风速拟合函数：di_vj＝fi(vp)式中，di_vj为各风力机(其中j取1至风机总数)点位处的风速，vp为整场气象预报风速。作为本发明的进一步改进，根据短期气象预报风速和拟合函数预测出风电场离散机位处的精确风速具体为：结合风电场运维系统接入的未来数小时实时气象预报风向与风速数据，从得到各机组预测风速函数关系，计算当前气象预报风况下，各机组点位未来数小时相应的预测风速。作为本发明的进一步改进，结合各机组实际功率曲线得到功率预测具体为：根据计算得到的预报风向下，各机组点位未来数小时相应的预测风速，结合各机组的风速-功率曲线表，采用风速内部差值法，得到各机组预测风速下对应的功率预测值，进而对各机组功率预测值进行加和，得到整场的功率输出预测值。一种复杂地形风电场风速与功率预测系统，包括：获取模块，用于获取目标风电场区域轮毂高度处的中尺度气象预报历史数据；计算模块，用于由不同来流风向扇区下，风电场区域性中尺度气象历史风速与各机组机舱风速计测得的同期风速数据的相关关系，建立整场气象风速与各机组点位风速之间的拟合函数；预测模块，用于根据短期气象预报风速和拟合函数预测出风电场离散机位处的精确风速，并结合各机组实际功率曲线得到功率预测。一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述复杂地形风电场风速与功率预测方法的步骤。一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述复杂地形风电场风速与功率预测方法的步骤。与现有技术方法相比，本发明的有益效果是：本发明提供的一种基于气象数据与机组实测数据拟合的复杂地形风电场风速与功率预测方法，该方法以风电场所在区域气象预报中的历史风速与风向数据为基础，结合各风电机组机舱风速计所测量的相应风机点位处的同期风速与风向数据，研究不同来流风向扇区下，区域气象预报风速与各个风机点位处风速的拟合函数关系，通过该函数关系预测未来一段时间内各机组位置的风速，进而根据功率曲线精确预测各机组以及整场的功率输出，以降低风电场生产运维中短期预测功率与实际功率存在的偏差。通过历史数据建立整场气象风速与各机组点位风速之间的拟合函数关系，进而根据不同风向下短期气象预报风速预测出风电场离散机位处的风速，并结合各机组实际功率曲线得到功率预测。与传统物理模拟方法相比，该方法简化了的预测模型的复杂程度，避免了工时较长的在线数值模拟，节省了风电场运维系统中的计算资源。本发明从数理统计的角度对复杂地形风电场的精细化风速与功率输出进行了预测，简化了传统物理模型的复杂程度，同时可减少短期功率预测与实际功率输出存在的偏差，为在役复杂地形风电场运维中的能效预测环节提供了有效的解决方法。进一步地，该方法有效地降低了采用物理模型进行计算时，由于复杂地形引起的湍流效应对流场风速模拟结果的不确定性问题，减少了风电场生产运维中短期预测功率与实际功率存在的偏差，从数理统计的角度为复杂地形风电场基于风机点位的精细化风速、功率预测提供了一种新的解决方法。附图说明图1为本发明所涉及的一种基于气象数据与机组实测数据拟合的复杂地形风电场风速与功率预测方法的流程图；图2为气象预报“东北(d3)”风向下，整场轮毂高度处的预报风速与第1台风机(wt1)轮毂高度处的风速相关关系分布与一阶多项式函数拟合曲线图；图3为气象预报“东北(d3)”风向下，整场轮毂高度处的预报风速与第10台风机(wt10)轮毂高度处的风速相关关系分布与一阶指数函数拟合曲线图；图4为本发明复杂地形风电场风速与功率预测系统结构示意图；图5为本发明电子设备结构示意图。具体实施方式为了使本发明的目的、技术方案和优点更加清楚明白，下面结合附图和具体实施例对本发明进一步详细说明：结合图1，本发明第一个目的在于提供一种基于气象数据与机组实测数据拟合的复杂地形风电场风速与功率预测方法，具体包括以下步骤：步骤一：获取目标风电场区域轮毂高度处的中尺度气象预报历史数据，并从中提取风电场运行投产以来至今的平均风速与风向的时间序列数据，建立气象历史风况数据集m，其数据结构如表1所示；表1日期时间气象预报轮毂高度处风向(°)气象预报轮毂高度处风速(m/s)2020-10-0100:00817.292020-10-0100:10897.852020-10-0100:20807.532020-10-0100:30826.87…………步骤二：对步骤1建立的气象历史风况数据集m进行不同来流风向扇区下的筛选与归类。2.1)所述来流风向扇区，可依照经典风向玫瑰图，选用十六方位来流风向，其中每个方位所代表的风向扇区中心角度di(其中i取1至16)，及各角度所代表的来流扇区角度范围如表2所示；表216方位来流风向对应的风向扇区中心角度与扇区角度范围2.2)以气象历史风况数据集m中的气象预报轮毂高度处风向数据为基础，按照2.1)中所划分的16方位来流风向进行数据筛选，获得每个来流扇区方位下的子数据集m_di(其中i取1至16)。例如，气象历史风况数据集m中，落于“东北”风扇区内的数据组成了相应的子数据集m_d3，其数据结构如表3所示。表3日期时间“东北”风向扇区(45°±11.25°)气象预报轮毂高度处风速(m/s)2020-10-0102:0055°4.672020-10-0102:1052°4.052020-10-0102:2048°4.082020-10-0102:3044°4.84…………步骤三：根据步骤二中获得的每个来流扇区方位下子数据集(m_di)中的日期和时间序列，获取风电场中各风电机组的同期scada运行数据，并从中提取机舱风速计测量得到风速数据，建立相应的整场机组风速数据集di_wts。例如，东北风向扇区下气象风况子数据集m_d3对应的同期整场风电机组风速数据集为d3-wts，其数据结构如表4所示：表4步骤四：将步骤二中获得的每个来流扇区方位下气象子数据集(m_di)与步骤三中获得的同期整场风电机组风速数据集(di_wts)合并，建立相应的整场气象预报与机组风速测量数据集m_di_wts。例如，m_d3_wts数据结构如表5所示：表5步骤五：针对步骤四建立的数据集m_di_wts，通过数理统计的方法，研究气象预报各来流风向di下，整场气象预报风速与各风机点位风速之间的相关关系，可采用多项式、对数、指数等函数关系进行拟合，最终根据各拟合函数的均方差(rootmeansquareerror，rmse)与确定系数(coefficientofdetermination，r2)的优劣程度，建立最优的风速拟合函数：di_vj＝fi(vp)式(1)式(1)中，di-vj为各风力机(其中j取1至风机总数)点位处的风速，vp为整场气象预报风速。例如，本实施例中，气象预报“东北(d3)”风向下轮毂高度处的预报风速与第1台风机(wt1)轮毂高度处的风速相关关系如图2所示，综合评价各拟合曲线的均方差与确定系数，最终采用的最优曲线拟合为一阶多项式函数：d3_v1＝f3(vp)＝0.8709·vp+0.587气象预报“东北(d3)”风向下轮毂高度处的预报风速与第10台风机(wt10)轮毂高度处的风速相关关系如图3所示，综合评价各拟合曲线的均方差与确定系数，最终采用的最优曲线拟合为一阶指数函数：步骤六：根据步骤五得到的各风机点位风速与整场气象预报风速之间的相关关系，建立基于气象预报风向与风速的各机组预测风速函数关系库。例如，本实施例中，基于气象预报的风向扇区，建立得到的各机组位置预测风速与气象预报风速的关系库结构如表6所示。表6步骤七：结合风电场运维系统接入的未来数小时实时气象预报风向与风速数据，从步骤六中得到的各机组预测风速函数关系，计算当前气象预报风况下，各机组点位未来数小时相应的预测风速。例如，若风电场当前接入的未来某时刻的气象预报风向为49.25°(属于东北风d3扇区)、预报风速为vp＝8.5m/s，则根据表6，采用d3_v1的计算公式，得到第1台机组轮毂高度处同期的预测风速为d3_v1＝0.8709·vp+0.587＝7.99m/s；采用d3_v10的计算公式，得到第10台机组轮毂高度处同期的预测风速为步骤八：根据步骤七中计算得到的预报风向下，各机组点位未来数小时相应的预测风速，结合各机组的风速-功率曲线表，采用风速内部差值法，得到各机组预测风速下对应的功率预测值，进而对各机组功率预测值进行加和，得到整场的功率输出预测值。例如，对于风电场当前接入的未来某时刻区域性气象预报数据：vp＝8.5m/s来流风向49.25°，通过拟合函数计算得到第1台机组轮毂高度处在该时刻的预测风速为d3_v1＝7.99m/s，根据该机型风速-功率曲线表查得：7.5m/s风速对应1212kw输出功率，8.0m/s风速对应1418kw输出功率，采用风速内部差值法，得到d3_v1对应的输出功率为：通过拟合函数计算得到第10台机组轮毂高度处在该时刻的预测风速为d3_v10＝8.01m/s，根据该机型风速-功率曲线表查得：8.0m/s风速对应1418kw输出功率，8.5m/s风速对应1604kw输出功率，采用风速内部差值法，得到d3_v10对应的输出功率为：如图4所示，本发明第二个目的是提供一种复杂地形风电场风速与功率预测系统，包括：获取模块，用于获取目标风电场区域轮毂高度处的中尺度气象预报历史数据；计算模块，用于由不同来流风向扇区下，风电场区域性中尺度气象历史风速与各机组机舱风速计测得的同期风速数据的相关关系，建立整场气象风速与各机组点位风速之间的拟合函数；预测模块，用于根据短期气象预报风速和拟合函数预测出风电场离散机位处的精确风速，并结合各机组实际功率曲线得到功率预测。如图5所示，本发明第三个目的是提供一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述复杂地形风电场风速与功率预测方法的步骤。本发明第四个目的是提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述复杂地形风电场风速与功率预测方法的步骤。本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。当前第1页12