CN113283035B - 双参数风力机机舱传递函数构建方法、系统、设备及存储介质 - Google Patents

双参数风力机机舱传递函数构建方法、系统、设备及存储介质 Download PDF

Info

Publication number
CN113283035B
CN113283035B CN202110656398.4A CN202110656398A CN113283035B CN 113283035 B CN113283035 B CN 113283035B CN 202110656398 A CN202110656398 A CN 202110656398A CN 113283035 B CN113283035 B CN 113283035B
Authority
CN
China
Prior art keywords
wind speed
wind
wind turbine
transfer function
data
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
CN202110656398.4A
Other languages
English (en)
Other versions
CN113283035A (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.)
Huaneng Hami Wind Power Generation Co ltd
Huaneng Xinjiang Energy Development Co ltd
North China Electric Power University
Xian Thermal Power Research Institute Co Ltd
Original Assignee
Huaneng Hami Wind Power Generation Co ltd
Huaneng Xinjiang Energy Development Co ltd
North China Electric Power University
Xian Thermal Power Research Institute Co Ltd
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 Huaneng Hami Wind Power Generation Co ltd, Huaneng Xinjiang Energy Development Co ltd, North China Electric Power University, Xian Thermal Power Research Institute Co Ltd filed Critical Huaneng Hami Wind Power Generation Co ltd
Priority to CN202110656398.4A priority Critical patent/CN113283035B/zh
Publication of CN113283035A publication Critical patent/CN113283035A/zh
Application granted granted Critical
Publication of CN113283035B publication Critical patent/CN113283035B/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
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/06—Wind turbines or wind farms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Wind Motors (AREA)

Abstract

本发明提供一种双参数风力机机舱传递函数构建方法、系统、设备及存储介质。方法包括，通过选择标杆风电机组，测量自由来流风速，获取同时段机舱风速和风电机组输出功率数据，形成数据集；对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理，进而基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数机舱传递函数关系。采用风力机机舱传递函数，从模型机理上增加了传递函数的考虑因素，解决了平坦地形条件下风力机运行控制和性能评价所需的自由来流风速参数缺失或估算准确性偏低的问题。同时减少了风电机组后评估时对遥感测风装置的依赖，为风力机运行控制优化和能效性能评估提供了一种实用的解决方案。

Description

双参数风力机机舱传递函数构建方法、系统、设备及存储介质
技术领域
本发明属于风力发电技术领域，具体涉及一种双参数风力机机舱传递函数构建方法、系统、设备及存储介质。
背景技术
风能工程中，自由来流是指风力机前方未受到叶轮旋转扰动的自然空气流动，自由来流风速是风力机主控系统运行、功率曲线等性能参数评估的重要依据。然而，风力机在运行期间通常是通过位于叶轮后机舱尾部的机舱风速计进行机组点位处风速的测量，进而在风力机监控及数据采集系统(SCADA)中记录。实践表明，受到机舱和叶片外形、叶轮转速、叶片桨距角、以及当地地形地貌等因素的影响，机舱风速计测得的机舱风速与自由来流风速存在着一定程度的误差，无法真实反映真实自由来流的情况。
对于没有自由来流风速测量装置的风力机，采用机舱传递函数(NTF，NacelleTransfer Function)是基于机舱风速反向推导自由来流风速的一种常用方法。机舱传递函数是指风力机机舱风速计测量风速(简称为机舱风速)与自由来流风速之间的函数关系。在确定的运行控制策略的条件下，风力机叶轮转速和桨距角等运行参数也是自由来流风速的函数，因此对于地形起伏较为平坦的风电场区域，通常假设机舱风速只是自由来流风速的单值函数。但在风力机的实际运行中，很可能出现相同的自由来流风速对应不同风力机运行状态的情况，此时理论上机舱风速是自由来流风速的多值函数，反之自由来流风速也是机舱风速的多值函数，需要考虑更多影响因素来确定自由来流风速与机舱风速的函数映射关系。
现有的机舱传递函数通常是基于自由来流风速与机舱风速的线性函数或多项式函数关系进行构造，其缺点是考虑的影响因素单一，特别是对于自由来流风速与风力机出力不是单值映射关系的风力机，根据机舱风速计算的自由来流风速偏差较大。
发明内容
本发明的目的在于提供双参数风力机机舱传递函数构建方法、系统、设备及存储介质，本发明方法解决了风力机运行控制和性能评价所需的自由来流风速参数缺失或估算准确性偏低的问题，提高了自由来流风速计算的准确性。
为了达到上述目的，本发明采用的技术方案是：
一种双参数风力机机舱传递函数构建方法，包括以下步骤：
通过选择标杆风电机组，测量自由来流风速，获取同时段机舱风速和风电机组输出功率数据，形成数据集；
对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理，进而基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数机舱传递函数关系。
作为本发明的进一步改进，选择标杆风电机组是从风电场区域主来流风向下的第一迎风排机组中，选取风力机。
作为本发明的进一步改进，测量自由来流风速是在所述标杆风电机组处安装测风装置，并进行校准与自由来流测风；
所述测风装置包括测风塔、地面激光测风雷达和机舱激光测风雷达；安装测风装置时，应对数据采集时间戳进行校准，使其与标杆风电机组SCADA系统中的时间戳保持一致，并提取自由来流风速的一定时间内的平均值。
作为本发明的进一步改进，形成数据集具体包括：
将测风装置测风期内每一个时间戳对应的叶轮前自由来流风速数据、机舱风速、机组输出功率进行汇总，建立同一时间下不同空间位置测风与相应功率输出的数据集；
对同一时间戳下的每条风速数据进行有效检测，如果记录中有缺失数据或超范围数据，则将该时间戳对应的整个数据序列视为无效数据；对无效数据记录进行剔除后形成的数据集作为标杆机组的自由来流风速-机舱风速-功率数据集M。
作为本发明的进一步改进，对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理具体包括：
将得到的自由来流风速-机舱风速-功率数据集M中，剔除风电机组输出功率为0的数据序列，得到功率有效风速数据集M_P；
将得到的功率有效风速数据集M_P进行修正，即保留风电机组输出功率大于a％额定功率值的数据序列，a为常数，形成修正的功率有效风速数据集M_P_R；
作为本发明的进一步改进，基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数函数关系具体包括：
将数据集中数据进行归一化处理，无量纲参数x为机组SCADA系统中机舱风速与最大机舱风速的比值，无量纲参数y为机组实际输出功率与额定功率的比值，无量纲参数z为自由来流风速与最大自由来流风速的比值；
将机舱风速、机组输出功率的无量纲参数(x、y)作为自变量，以自由来流风速的无量纲参数(z)为因变量，采用二元五次多项式函数，对三者形成的三维空间关系进行拟合，得到双参数机舱传递函数方程：
z＝fs(x,y)＝p00+p10x+p01y+p20x2+p11xy+p02y2+p30x3+p21x2y+p12xy2+p03y3+p40x4+p31x3y+p22x2y2+p13xy3+p04y4+p50x5+p41x4y+p32x3y2+p23x2y3+p14xy4+p05y5
其中，p00、p10、p01、p20、p11、p02、p30、p21、p12、p03、p40、p31、p22、p13、p04、p50、p41、p32、p23、p14、p05均为为方程的拟合常系数。
作为本发明的进一步改进，还包括将双参数机舱传递函数进行精确性验证的步骤，包括：
通过拟合出的双参数机舱传递函数，计算修正的功率有效风速数据集M_P_R中每个时间戳对应的自由来流风速；
计算采用所得双参数机舱传递函数z＝fs(x,y)拟合出的自由来流风速与其测量值的之间的偏差；
将基于曲面拟合的双参数机舱传递函数与基于曲线拟合的单参数机舱传递函数进行对比，分析比较二者计算出的自由来流风速与实测值的总体偏差。
若采用基于曲面拟合得到的双参数机舱传递函数总体平均偏差大于阈值，或者大于曲线拟合机舱传递函数的偏差，则返回修正方程的拟合系数，直至精确性达到要求。
一种双参数风力机机舱传递函数构建系统，包括：
获取模块，用于通过选择标杆风电机组，测量自由来流风速，获取同时段机舱风速和风电机组输出功率数据，形成数据集；
拟合模块，用于对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理，进而基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数机舱传递函数关系。
一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述双参数风力机机舱传递函数构建方法的步骤。
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述双参数风力机机舱传递函数构建方法的步骤。
与现有技术方法相比，本发明的有益效果是：
本发明提供的一种适用于平坦地形风电场的双参数风力机机舱传递函数的构建方法，将自由来流风速视为机舱风速和风力机发电功率的函数，通过在风电场区标杆机组处安装短期便携式测风工具，对风轮前方的自由来流风速进行采集，同时结合同期机舱风速数据与机组输出功率，通过特定的数据处理流程，建立风力机自由来流风速与机舱风速和风力机输出功率的函数映射关系，解决了风力机运行控制和性能评价所需的自由来流风速参数缺失或估算准确性偏低的问题，提高了自由来流风速计算的准确性。
进一步地，本发明所构建的双参数风力机机舱传递函数，相比于常用的单参数机舱传递函数，增加了机组输出功率作为其自变量，从模型机理上增加了传递函数的考虑因素。由于调度的原因和风力机本身的原因，如风力机翼型处于动态失速区的时候，实际风力机运行状态不是来流风速的单值函数，机舱风速与自由来流风速的关系受到风力机运行状态的影响，风电机组的输出功率可以表示风力机的运行状态，因此基于机舱风速和机组输出功率两个参数构建机舱传递函数，在理论上可提高自由来流风速计算的准确性。
进一步地，本发明在风力机没有配置激光雷达等遥感测风装置的条件下，可结合风力机监控及数据采集系统(SCADA)中的机舱风速和风力机出力记录，根据所构造的机舱传递函数计算出风力机上游自由来流风速，进而依据自由来流风速优化风力机的运行参数，减少了风力机对遥感测风装置的依赖，从而降低了风力机系统的复杂性和成本，有利于提高风力机运行的经济性。
进一步地，本发明建立的风力机机舱函数，可在地形起伏较为平坦的风电场中，可替代激光雷达遥感测量设备，根据风电运营商的后评估需求，还原任意一台机组前方自由来流的风速，该方法计算精度高、快速、易于编程实现，便于工程应用，为风力机运行控制优化和能效性能评估提供了一种实用的解决方案。
附图说明
图1为本发明所涉及的一种适用于平坦地形风电场的双参数风力机机舱传递函数的构建方法的流程图；
图2为机舱风速、机组输出功率与自由来流风速的无量纲三维空间关系图；
图3为自由来流风速的双参数机舱传递函数拟合结果；
图4为自由来流风速的单参数机舱传递函数拟合结果；
图5为本发明双参数风力机机舱传递函数构建系统结构示意图；
图6为本发明电子设备结构示意图。
具体实施方式
为了使本发明的目的和技术方案更加清晰和便于理解。以下结合附图和实施例，对本发明进行进一步的详细说明。
参照附图1，本发明第一个目的是提供一种适用于平坦地形风电场的双参数风力机机舱传递函数的构建方法，具体可从以下八个步骤进行实施：
步骤一：选取目标风电场中的标杆风电机组。
标杆风电机组应从风电场区域主来流风向下的第一迎风排机组中，选取运行故障率较低、历史可利用小时数较高的风力机。
步骤二：在步骤一所述标杆风电机组处安装测风装置，并进行校准与自由来流测风。
所述测风装置包括测风塔、地面激光测风雷达和机舱激光测风雷达。测风塔和地面激光测风雷达通常放置于标杆风电机组的主风向上游，距离机组2倍至4倍叶轮直径处，确保能获取风力机轮毂高度处的自由来流风速；机舱激光测风雷达通常安装在风力发电机组机舱顶端，实现由机舱顶端对机舱叶轮面前方几十到几百米风矢量场的精确测量，叶轮前方最远处的风速测量值可被认为是自由来流风速。同时，安装便携式测风装置时，应对数据采集时间戳进行校准，使其与标杆风电机组SCADA系统中的时间戳保持一致，并提取自由来流风速的10分钟平均值。
此实施例中，选择在标杆风电机组的机舱部位安装机舱激光测风雷达，对叶轮前方0m、60m、80m、100m、120m、140m、160m、180m、190m、200m点位处的风速、湍流强度进行测量，将200m空间点位处的风速作为该机组的自由来流风速，并提取其10分钟风速平均值。
步骤三：获取步骤二所述测风装置所采集数据对应的同期风力机运行SCADA数据，主要包括机舱风速10分钟平均值、以及机组的输出功率10分钟平均值。
步骤四：建立标杆机组的自由来流风速-机舱风速-功率有效数据集M。
4.1将测风装置测风期内每一个时间戳对应的叶轮前自由来流风速数据、机舱风速、机组输出功率进行汇总，建立同一时间下不同空间位置测风与相应功率输出的数据集。
此实施例中，同一个时间戳下的数据，包含风力机叶轮前上游200m空间点位的激光雷达测风数据和机舱风速、风电机组输出功率，共计2个空间点位风速数据和1个输出功率数据。
4.2对同一时间戳下的每条风速数据进行有效检测，如果记录中有缺失数据或超范围数据(风速<0或>50m/s)，则将该时间戳对应的整个数据序列视为无效数据。对无效数据记录进行剔除后形成的数据集作为标杆机组的自由来流风速-机舱风速-功率数据集M。
步骤五：将步骤四得到的自由来流风速-机舱风速-功率数据集M中，剔除风电机组输出功率为0(非运行状态下)的数据序列，得到功率有效风速数据集M_P。
步骤六：将步骤五得到的功率有效风速数据集M_P进行修正，即保留风电机组输出功率大于3％额定功率值的数据序列，形成修正的功率有效风速数据集M_P_R。
若在功率有效风速数据集中，有个别数据显示的是风速较大条件下风电机组的功率输出极小，这些数据极可能是数据测量出现问题。因此，以步骤五得到的功率有效风速数据集为基础，进一步将功率有效的条件修正为机组额定功率的3％，剔除掉有效风速数据集中小于机组额定功率值3％的输出功率及其相应的风速测量数据序列。
步骤七：将自由来流风速视为机舱风速、机组输出功率两个参数的函数，通过曲面拟合模式，进而反映三个变量形成的三维空间关系，得到双参数机舱传递函数方程。
7.1将步骤六得到修正的功率有效风速数据集M_P_R中数据进行归一化(无量纲)处理。无量纲参数x为机组SCADA系统中机舱风速与最大机舱风速的比值，无量纲参数y为机组实际输出功率与额定功率的比值，无量纲参数z为自由来流风速与最大自由来流风速的比值。
7.2将机舱风速、机组输出功率的无量纲参数(x、y)作为自变量，以自由来流风速的无量纲参数(z)为因变量，采用二元五次多项式函数，对三者形成的三维空间关系(如附图2所示)进行拟合，得到双参数机舱传递函数方程z＝fs(x,y)。
所述二元五次多项式函数的形式为：
z＝p00+p10x+p01y+p20x2+p11xy+p02y2+p30x3+p21x2y+p12xy2+p03y3+p40x4+p31x3y+p22x2y2+p13xy3+p04y4+p50x5+p41x4y+p32x3y2+p23x2y3+p14xy4+p05y5
其中，p00、p10、p01、p20、p11、p02、p30、p21、p12、p03、p40、p31、p22、p13、p04、p50、p41、p32、p23、p14、p05均为为方程的拟合常系数。
作为优选实施例，pij为方程的拟合系数，具体如下：
p00＝0.2091，p10＝-2.663，p01＝3.567，p20＝35.13，p11＝-55.95，p02＝21.89，p30＝-188.8，
p21＝387.3，p12＝-255.7，p03＝43.4，p40＝494.5，p31＝-1232，p22＝1101，p13＝-386.5，p04＝51.31，
p50＝-477.9，p41＝1406，p32＝-1586，p23＝839，p14＝-221.2，p05＝23。
步骤八：将步骤七得到的双参数机舱传递函数方程z＝fs(x,y)进行精确性验证。
8.1通过步骤七拟合出的双参数机舱传递函数方程z＝fs(x,y)，计算修正的功率有效风速数据集M_P_R中每个时间戳对应的自由来流风速。
若修正后的功率有效风速数据集M_P_R中数据长度为N，即修正的功率有效风速数据集包括N个时间序列，则第i个时间点的机舱风速为Vni，测量的自由来流风速为Vmfi，风电机组输出功率为Pi。N个时间序列中记录的最大机舱风速为Vnmax，最大自由来流风速为Vfmax，风电机组额定输出功率为Prated，则第i个时间点对应的无量纲参数x、y分别为
Figure BDA0003113003200000091
由双参数机舱传递函数，可计算得到无量纲自由来流风速zi＝fs(xi,yi)，进而可计算出自由来流风速为Vcfi＝zi×Vfmax，其结果如附图3所示。
8.2计算采用步骤七所得双参数机舱传递函数z＝fs(x,y)拟合出的自由来流风速与其测量值的之间的偏差(简称机舱传递函数偏差)。
修正的功率有效风速数据集M_P_R中，第i个时间点对应双参数机舱传递函数的偏差计算方法为
Figure BDA0003113003200000092
修正的功率有效风速数据集M_P_R中，双参数机舱传递函数的总体平均偏差计算方法为：
Figure BDA0003113003200000101
8.3将基于曲面拟合的双参数机舱传递函数与基于曲线拟合(常用)的单参数机舱传递函数进行对比，分析比较二者计算出的自由来流风速与实测值的总体偏差。
基于曲线拟合的单参数机舱传递函数，将自由来流风速仅视为机舱风速的函数，即
zL＝fL(x)＝p1x5+p2x4+p3x3+p4x2+p5x+p6
其中pi为方程的拟合系数，具体为p1＝-1.299e-05，p2＝0.0007325，p3＝-0.01566，p4＝0.1587，p5＝0.1465，p6＝2.134。
由单参数机舱传递函数，计算得到第i个时间点对应的无量纲自由来流风速为
Figure BDA0003113003200000102
进而可计算出自由来流风速为
Figure BDA0003113003200000103
其结果如附图4所示，其与真实测量值的偏差为
Figure BDA0003113003200000104
则单参数机舱传递函数的总体平均偏差为：
Figure BDA0003113003200000105
采用基于曲面拟合得到的双参数机舱传递函数总体平均偏差为5.040％，采用基于曲线拟合(常用)的单参数机舱传递函数的总体平均偏差为5.435％，因此，和曲线拟合机舱传递函数相比，曲面拟合机舱传递函数的偏差绝对值下降0.395％，下降的相对幅度为7.27％。
8.4若采用基于曲面拟合得到的双参数机舱传递函数总体平均偏差较大(大于10％)，或者大于曲线拟合机舱传递函数的偏差，则返回步骤七修正方程的拟合系数，直至精确性达到要求。
采用本发明构建形成的风力机机舱传递函数，从模型机理上增加了传递函数的考虑因素，解决了平坦地形条件下风力机运行控制和性能评价所需的自由来流风速参数缺失或估算准确性偏低的问题。同时减少了风电机组后评估时对遥感测风装置的依赖，为风力机运行控制优化和能效性能评估提供了一种实用的解决方案。
如图5所示，本发明第二个目的是提供一种双参数风力机机舱传递函数构建系统，包括：
获取模块，用于通过选择标杆风电机组，测量自由来流风速，获取同时段机舱风速和风电机组输出功率数据，形成数据集；
拟合模块，用于对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理，进而基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数机舱传递函数关系。
还包括验证，其用于将双参数机舱传递函数进行精确性验证，包括：
通过拟合出的双参数机舱传递函数，计算修正的功率有效风速数据集中每个时间戳对应的自由来流风速；
计算采用所得双参数机舱传递函数拟合出的自由来流风速与其测量值的之间的偏差；
将基于曲面拟合的双参数机舱传递函数与基于曲线拟合的单参数机舱传递函数进行对比，分析比较二者计算出的自由来流风速与实测值的总体偏差。
若采用基于曲面拟合得到的双参数机舱传递函数总体平均偏差大于阈值，或者大于曲线拟合机舱传递函数的偏差，则返回修正方程的拟合系数，直至精确性达到要求。
如图6所示，本发明第三个目的是提供一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述双参数风力机机舱传递函数构建方法的步骤。
本发明第四个目的是提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述双参数风力机机舱传递函数构建方法的步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (9)

1.一种双参数风力机机舱传递函数构建方法，其特征在于，包括以下步骤：
通过选择标杆风电机组，测量自由来流风速，获取同时段机舱风速和风电机组输出功率数据，形成数据集；
对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理，进而基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数机舱传递函数关系；
基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数函数关系具体包括：
将数据集中数据进行归一化处理，无量纲参数x为机组SCADA系统中机舱风速与最大机舱风速的比值，无量纲参数y为机组实际输出功率与额定功率的比值，无量纲参数z为自由来流风速与最大自由来流风速的比值；
将机舱风速、机组输出功率的无量纲参数(x、y)作为自变量，以自由来流风速的无量纲参数(z)为因变量，采用二元五次多项式函数，对三者形成的三维空间关系进行拟合，得到双参数机舱传递函数方程：
z＝fs(x,y)＝p00+p10x+p01y+p20x2+p11xy+p02y2+p30x3+p21x2y+p12xy2+p03y3+p40x4+p31x3y+p22x2y2+p13xy3+p04y4+p50x5+p41x4y+p32x3y2+p23x2y3+p14xy4+p05y5
其中，p00、p10、p01、p20、p11、p02、p30、p21、p12、p03、p40、p31、p22、p13、p04、p50、p41、p32、p23、p14、p05均为方程的拟合常系数。
2.根据权利要求1所述的方法，其特征在于，
选择标杆风电机组是从风电场区域主来流风向下的第一迎风排机组中，选取风力机。
3.根据权利要求1所述的方法，其特征在于，
测量自由来流风速是在所述标杆风电机组处安装测风装置，并进行校准与自由来流测风；
所述测风装置包括测风塔、地面激光测风雷达和机舱激光测风雷达；安装测风装置时，应对数据采集时间戳进行校准，使其与标杆风电机组SCADA系统中的时间戳保持一致，并提取自由来流风速的平均值。
4.根据权利要求1所述的方法，其特征在于，
形成数据集具体包括：
将测风装置测风期内每一个时间戳对应的叶轮前自由来流风速数据、机舱风速、机组输出功率进行汇总，建立同一时间下不同空间位置测风与相应功率输出的数据集；
对同一时间戳下的每条风速数据进行有效检测，如果记录中有缺失数据或超范围数据，则将该时间戳对应的整个数据序列视为无效数据；对无效数据记录进行剔除后形成的数据集作为标杆机组的自由来流风速-机舱风速-功率数据集M。
5.根据权利要求1所述的方法，其特征在于，
对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理，具体包括：
将得到的自由来流风速-机舱风速-功率数据集中，剔除风电机组输出功率为0的数据序列，得到功率有效风速数据集；
将得到的功率有效风速数据集进行修正，即保留风电机组输出功率大于a％额定功率值的数据序列，a为常数，形成修正的功率有效风速数据集。
6.根据权利要求1所述的方法，其特征在于，
还包括将双参数机舱传递函数进行精确性验证的步骤，包括：
通过拟合出的双参数机舱传递函数，计算修正的功率有效风速数据集中每个时间戳对应的自由来流风速；
计算采用所得双参数机舱传递函数拟合出的自由来流风速与其测量值的之间的偏差；
将基于曲面拟合的双参数机舱传递函数与基于曲线拟合的单参数机舱传递函数进行对比，分析比较二者计算出的自由来流风速与实测值的总体偏差；
若采用基于曲面拟合得到的双参数机舱传递函数总体平均偏差大于阈值，或者大于曲线拟合机舱传递函数的偏差，则返回修正方程的拟合系数，直至精确性达到要求。
7.一种双参数风力机机舱传递函数构建系统，基于权利要求1-6任意一项所述的双参数风力机机舱传递函数构建方法，其特征在于，包括：
获取模块，用于通过选择标杆风电机组，测量自由来流风速，获取同时段机舱风速和风电机组输出功率数据，形成数据集；
拟合模块，用于对数据集进行数据有效性和风电机组运行状态有效性识别和数据处理，进而基于曲面拟合的方法建立自由来流风速与机舱风速、风电机组出力的双参数机舱传递函数关系。
8.一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1-6任一项所述双参数风力机机舱传递函数构建方法的步骤。
9.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1-6任一项所述双参数风力机机舱传递函数构建方法的步骤。
CN202110656398.4A 2021-06-11 2021-06-11 双参数风力机机舱传递函数构建方法、系统、设备及存储介质 Active CN113283035B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110656398.4A CN113283035B (zh) 2021-06-11 2021-06-11 双参数风力机机舱传递函数构建方法、系统、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110656398.4A CN113283035B (zh) 2021-06-11 2021-06-11 双参数风力机机舱传递函数构建方法、系统、设备及存储介质

Publications (2)

Publication Number Publication Date
CN113283035A CN113283035A (zh) 2021-08-20
CN113283035B true CN113283035B (zh) 2023-04-07

Family

ID=77284552

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110656398.4A Active CN113283035B (zh) 2021-06-11 2021-06-11 双参数风力机机舱传递函数构建方法、系统、设备及存储介质

Country Status (1)

Country Link
CN (1) CN113283035B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115358606B (zh) * 2022-08-26 2023-09-05 中国电建集团西北勘测设计研究院有限公司 一种平缓地形下在役风电场能效评估方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091209A (zh) * 2014-06-26 2014-10-08 沈阳工业大学 基于bp神经网络的风电机组功率特性评估方法
CN105279384A (zh) * 2015-11-16 2016-01-27 华北电力科学研究院有限责任公司 一种基于风力机机舱风速的来流风速计算方法及装置
CN108427013A (zh) * 2018-03-28 2018-08-21 中国船舶重工集团海装风电股份有限公司 一种机舱风速仪风速畸变函数的获取方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5022102B2 (ja) * 2007-05-25 2012-09-12 三菱重工業株式会社 風力発電装置、風力発電システムおよび風力発電装置の発電制御方法
CN104794347B (zh) * 2015-04-22 2017-08-29 中国大唐集团新能源股份有限公司 一种风电机组机舱传递函数分区拟合方法
CN105354763A (zh) * 2015-11-16 2016-02-24 华北电力科学研究院有限责任公司 风电机组来流风速测量方法及装置
CN106815456A (zh) * 2015-12-02 2017-06-09 中国电力科学研究院 一种风电机组功率特性评价方法
CN109356800A (zh) * 2018-12-12 2019-02-19 国电联合动力技术有限公司 低风速风电机组机舱风速修正函数的获得方法及其装置
CN110956414A (zh) * 2019-12-18 2020-04-03 中电投电力工程有限公司 一种风电场运行监盘和能效评估系统及方法
CN111311021A (zh) * 2020-03-12 2020-06-19 广西电网有限责任公司 一种风电场理论功率预测方法、装置、设备和存储介质
CN112861429B (zh) * 2021-01-20 2022-08-30 中国长江三峡集团有限公司 一种风力机机舱传递函数计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091209A (zh) * 2014-06-26 2014-10-08 沈阳工业大学 基于bp神经网络的风电机组功率特性评估方法
CN105279384A (zh) * 2015-11-16 2016-01-27 华北电力科学研究院有限责任公司 一种基于风力机机舱风速的来流风速计算方法及装置
CN108427013A (zh) * 2018-03-28 2018-08-21 中国船舶重工集团海装风电股份有限公司 一种机舱风速仪风速畸变函数的获取方法

Also Published As

Publication number Publication date
CN113283035A (zh) 2021-08-20

Similar Documents

Publication Publication Date Title
EP2169218B1 (en) System and method for estimating wind condition for wind turbines
CN105041572B (zh) 用于优化风电场操作的系统和方法
CN111091298B (zh) 一种风电场流场耦合特性评估与智能分群方法及系统
CN104794347B (zh) 一种风电机组机舱传递函数分区拟合方法
US20160084233A1 (en) Systems and methods for validating wind farm performance measurements
US20120211986A1 (en) Wind sensor system using blade signals
US20170335827A1 (en) System and Method for Determining Wind Farm Wake Loss
US20170022974A1 (en) Operating wind turbines
CN109784563B (zh) 一种基于虚拟测风塔技术的超短期功率预测方法
CN108317040B (zh) 偏航对风矫正的方法、装置、介质、设备和风力发电机组
CN113205210B (zh) 复杂地形风电场风速与功率预测方法、系统、设备及存储介质
CN107832899B (zh) 风电场出力的优化方法、装置和实现装置
CN108431404B (zh) 用于控制多个风力涡轮机的方法和系统
CN113283035B (zh) 双参数风力机机舱传递函数构建方法、系统、设备及存储介质
CN108062722B (zh) 基于风速变异系数的山地风电场模型风机机械功率计算方法
CN110991701A (zh) 一种基于数据融合的风电场风机风速预测方法及系统
CN113642884A (zh) 一种电网失电情况下风电场发电量损失统计方法及系统
CN108593967A (zh) 风速值修正方法和装置、计算机可读存储介质
CN115898787A (zh) 一种风电机组静态偏航误差动态识别方法及装置
CN113586336B (zh) 风力发电机组的控制方法及其控制装置及计算机可读存储介质
CN114398842A (zh) 一种在运行风电场发电量评估方法
CN113269430A (zh) 基于运行测风塔实测数据的风电项目发电量后评估方法
CN113705126B (zh) 一种基于二维迹线追踪的风电场发电量校核方法、系统、计算机设备、存储介质
CN112861301A (zh) 一种基于风机实时数据的风电场理论功率智能计算方法
CN117145701A (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
GR01 Patent grant
GR01 Patent grant