一种基于数据驱动的风机切入风速计算方法
技术领域
本发明涉及风机切入风速计算技术领域,特别是涉及一种基于数据驱动的风机切入风速计算方法。
背景技术
在当前的能源环境背景下,新能源发电特别是风力发电的发展受到了广泛的关注。发展风电成为了减少国民经济对化石能源的依赖、解决能源生产与消费之间的矛盾和减少温室气体排放保持生态平衡的重要途径。随着资源与环境双重压力的持续增大,发展风力发电已成为我国甚至是国际未来能源利用的发展方向。
随着我国近年风电事业的迅猛发展,早年投入使用的老旧风机已经接近其设计寿命甚至超期服役,零部件老化现象也逐渐显现,原先制定的风机发电切入风速设定值已难以适应实际运行环境与工况的改变,因此,对老旧风机发电的切入风速设定值开展修订计算方法研究,对合理利用老旧风机,提升风场经济效益具有重要的理论指导与使用价值。
目前,对于风机切入风速计算的研究,主要通过现场试验的方法,但这种方法需要长时间进行试验,不仅耗费大量时间,也耗费大量的人力,并且由于风速的不确定性,很难通过现场的试验得出有效的风速结果。
发明内容
本发明提供了一种基于数据驱动的风机切入风速计算方法。
本发明提供了如下方案:
一种基于数据驱动的风机切入风速计算方法,包括:
采集待计算风机的秒级运行数据组,并通过风向角的余弦计算,校正获得对应功率下的理论风速数据组;
通过厂家提供的所述待计算风机型号的风功率曲线表数据,使用最小二乘方法建立风机功率曲线的数学模型,去除所述理论风速数据组中风速与对应功率偏差较大的数据获得待筛选风速数据组;
根据所述待计算风机的运行状态值,建立风机切入状态模型,对所述待筛选风速数据组包含的风速数据进行进一步筛选,筛选出符合所述待计算风机切入状态的切入风速数据组;
根据切入风速下的功率对所述切入风速数据组中包含的风速数据进行分类,保留主动切入的风速数据获得主动切入风速数据组;
将所述主动切入风速数据组中包含的风速数据进行排序,通过统计分析学中的分位数选择方法,选出风机能够启动的风速作为所述待计算风机的切入风速。
优选的:所述余弦计算所用公式为:
S'=S·cosα
上式中S'为校正后风速,S为风机运行秒级数据组中的风速,α为风速的风向角。
优选的:所述最小二乘方法使用的拟合函数为:
y=∑θixi
上式中θi为标准功率曲线的多项式参数,x为风机运行秒级数据组中的风速,i为多项式拟合的阶数,y为标准功率曲线下计算的标准功率。
优选的:所述去除所述理论风速数据组中风速与对应功率偏差较大的数据获得待筛选风速数据组;包括:
去除风速大于5米每秒的风速数据获得待筛选风速数据组。
优选的:去除风速大于5米每秒的风速数据时的规则为:遍历每一条运行数据,若第i条数据风速大于5米每秒,则去掉第i条数据;若风速小于5米每秒,则保留第i条数据。
优选的:在去掉大于5米的风速数据之后,根据标准功率曲线对数据进行进一步筛选获得筛选风速数据组,对于功率大于0的情况,筛选公式为:
上式中error为实际功率与标准功率曲线的偏差,P为实际功率,y为标准功率曲线下计算的标准功率,θi为标准功率曲线的多项式参数,x、S均为风机运行秒级数据组中的风速,i为多项式拟合的阶数;筛选规则为去除error大于0.5的数据,即实际功率与标准功率曲线功率偏差超过50%的数据,保留error小于0.5的数据。
优选的:所述运行状态包括待机状态、正常发电状态、维护状态、故障状态;将与待机状态以及正常发电状态对应的待筛选风速数据组中的风速数据保留,获得所述切入风速数据组。
优选的:所述将所述主动切入风速数据组中包含的风速数据进行排序,通过统计分析学中的分位数选择方法,选出风机能够启动的风速作为所述待计算风机的切入风速;包括:
排序的规则为从小到大排序;再通过四分位数的规则选取第3个四分位数作为所述待计算风机的切入风速。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
通过本发明,可以实现一种基于数据驱动的风机切入风速计算方法,在一种实现方式下,该方法可以包括采集待计算风机的秒级运行数据组,并通过风向角的余弦计算,校正获得对应功率下的理论风速数据组;通过厂家提供的所述待计算风机型号的风功率曲线表数据,使用最小二乘方法建立风机功率曲线的数学模型,去除所述理论风速数据组中风速与对应功率偏差较大的数据获得待筛选风速数据组;根据所述待计算风机的运行状态值,建立风机切入状态模型,对所述待筛选风速数据组包含的风速数据进行进一步筛选,筛选出符合所述待计算风机切入状态的切入风速数据组;根据切入风速下的功率对所述切入风速数据组中包含的风速数据进行分类,保留主动切入的风速数据获得主动切入风速数据组;将所述主动切入风速数据组中包含的风速数据进行排序,通过统计分析学中的分位数选择方法,选出风机能够启动的风速作为所述待计算风机的切入风速。该方法结合了风机标准功率曲线数据,大大提高了切入风速筛选和选择的效率。基于数据驱动的优势,充分利用了风机运行的历史数据,而不必在现场长时间对风机运行进行观测和记录,节省了切入风速计算的时间和成本。本方法考虑了风向角对切入风速的影响,使计算结果更为准确。
当然,实施本发明的任一产品并不一定需要同时达到以上所述的所有优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种基于数据驱动的风机切入风速计算方法的流程图;
图2为导出WT02093风机的秒级运行数据;
图3为cos校正后的风机秒级运行数据;
图4为5阶多项式拟合标准功率曲线;
图5为筛选出小于5米的风速数据;
图6为根据功率曲线筛选后数据;
图7为切入逻辑筛选后结果;
图8为主动切入风速列表。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员所获得的所有其他实施例,都属于本发明保护的范围。
实施例
参见图1,为本发明实施例提供的一种基于数据驱动的风机切入风速计算方法,如图1所示,该方法包括采集待计算风机的秒级运行数据组,并通过风向角的余弦计算,校正获得对应功率下的理论风速数据组;
通过厂家提供的所述待计算风机型号的风功率曲线表数据,使用最小二乘方法建立风机功率曲线的数学模型,去除所述理论风速数据组中风速与对应功率偏差较大的数据获得待筛选风速数据组;
根据所述待计算风机的运行状态值,建立风机切入状态模型,对所述待筛选风速数据组包含的风速数据进行进一步筛选,筛选出符合所述待计算风机切入状态的切入风速数据组;
根据切入风速下的功率对所述切入风速数据组中包含的风速数据进行分类,保留主动切入的风速数据获得主动切入风速数据组;
将所述主动切入风速数据组中包含的风速数据进行排序,通过统计分析学中的分位数选择方法,选出风机能够启动的风速作为所述待计算风机的切入风速。
该方法在实际应用时,可以包含以下步骤:
步骤1:风机运行数据获取及校正,具体包括如下步骤:
步骤1.1:获取风机运行秒级数据
在对风机切入风速进行计算时,需要获取风机运行时的数据,由于风速具有不确定性,因此通过高频率的数据采样能够尽可能保留真实的风速信息。风机的秒级数据是风场现场采样频率最高的数据,对于风电运营企业来说,可以从公司服务器的数据库中获取某风场某台风机的秒级数据。
步骤1.2:余弦计算校正风速数据
风场现场采集的数据为一次数据,即直接从传感器获取的数据。风速传感器只能采集到风速的大小,风向的采集由风向仪完成。本方法为了消除计算的偏差,将实时风速校正为投影到风向角为0的方向的风速。
通过余弦计算来校正风速,使用公式为:
S'=S·cosα
上式中S'为校正后风速,S为风机运行秒级数据中的风速,α为风速的风向角。
步骤2:风功率曲线的拟合与风机运行秒级数据筛选,具体包括如下步骤:
步骤2.1:最小二乘法拟合风功率曲线
由于风机生产商标准功率曲线以列表的形式提供,精度为0.5米,为了更好筛选数据,通过最小二乘拟合标准功率曲线。本方法使用最小二乘法多项式拟合,对于给定的离散风功率曲线数据,恰当选取通用模型中的函数类别与具体形式,以保证拟合效果。功率上升部分符合多项式函数的特征,通过多项式拟合函数,实现风功率曲线的拟合。
多项式拟合可以选择不同阶数的拟合函数:
y=Σθixi
其中,θi为标准功率曲线的多项式参数,x为风机运行秒级数据组中的风速,i为多项式拟合的阶数,y为标准功率曲线下计算的标准功率。
本方法比较了2至6阶的拟合函数,选择拟合效果最好的高阶多项式作为标准功率曲线的数学模型。
步骤2.2:数据初步筛选
风机运行秒级数据的采样周期为1至5秒,即每隔1至5秒采集一组风机运行数据,而短时间内风机的启动和切入次数十分有限,为了计算符合统计规律的切入风速,需要通过3至5个月的风机运行数据进行统计。风机出厂设置的切入风速为3米每秒,在风机的使用过程中切入风速会发生变化,但不会超过5米每秒。因此可以去除风速大于5米每秒的数据来减少数据量,对数据进行初步筛选。筛选规则为:
遍历每一条运行数据,若第i条数据风速大于5米每秒,则去掉第i条数据;若风速小于5米每秒,则保留第i条数据。
步骤2.3:通过标准功率曲线筛选数据
在去掉大于5米的风速数据之后,根据标准功率曲线对数据进行进一步筛选。对于功率大于0的情况,筛选公式为:
其中error为实际功率与标准功率曲线的偏差,P为实际功率,y为标准功率曲线下计算的标准功率,θi为标准功率曲线的多项式参数,x为风速S。筛选规则为去除error大于0.5的数据,即实际功率与标准功率曲线功率偏差超过50%的数据,保留error小于0.5的数据。
步骤3:建立基于风机运行状态的风机切入模型,具体包括如下步骤:
步骤3.1:建立风机切入判断逻辑
风机运行秒级数据中包含了风机状态数据,风机状态包括待机、正常发电、维护、故障等。根据不同的风机协议,风机的状态对应着不同的状态码。通过遍历风机运行数据,假设第i组数据的运行状态为正常发电状态,第i-1组数据的运行状态为待机状态,则这两组数据符合风机启动切入的逻辑,保留这两组数据;运行状态不改变或者状态改变逻辑不符合切入判断逻辑的数据则去除。
步骤3.2:切入风速排序及切出风速选择
通过以上步骤,可以获得满足风机切入情况的风速数据,通常根据4个月的风机数据进行计算可用获得120组切入风速数据,而在统计学中普遍使用四分位数对样本进行定性分析,因此本方法选择第3个分位数,即75%以上情况能主动启动的风速作为风机的切入风速。排序及选择的步骤为:首先将筛选得到的风机切入风速进行排序,排序的规则为从小到大排序;然后通过四分位数的规则选取第3个四分位数作为风机的切入风速。
风力发电在我国的发展已有长达二十多年的发展,距离第一批风电机组的投产运行也已过了20余年,设备老化等引发的故障停运事故层出不穷,因此,老旧风机的性能与出厂时的标准有较大的偏差。若仍按照出厂制定的切入风速启动,则可能在此风速下只能从电网取电来启动风机,这不仅增加了风电场的场用电,而且也大大降低了风机的经济效益。目前风机切入风速的校正还没有太好的研究方法,主要通过风场现场试验校正,不仅费时而且费力。本发明提出的基于数据驱动的风机切入风速计算方法,通过对风机运行秒级数据进行筛选与分析,最后根据统计分位点风速确定风机切入风速的方法,具有快速高效、并且能针对不同风机单独计算其切入风速的特点,有着重要的现实意义。为了验证前面提出的方法的有效性和合理行,采取某风电场某风机的数据进行分析验证。具体方法步骤如下所示:
步骤1:风机运行数据获取及校正,具体包括如下步骤:
步骤1.1:获取风机运行秒级数据
在对风机切入风速进行计算时,需要获取风机运行时的数据。以某风场的风机编号为WT02093的风机为例,从龙源电力集团数据库中导出WT02093风机2017年1月1日0:0到2017年5月15日14:38分的秒级数据,见附图2,数据总量为982895条数据。
步骤1.2:余弦计算校正风速数据
风场现场采集的数据为一次数据,即直接从传感器获取的数据。从附图2可以看到,风机秒级运行数据中包含了风速信息和风向信息,为了消除计算的偏差,将实时风速校正为投影到风向角为0的方向的风速。
通过余弦计算来校正风速,使用公式为:
S'=S·cosα
上式中S'为校正后风速,S为风机运行秒级数据中的风速,α为风速的风向角。校正风速后的风机数据如附表3所示。
步骤2:风功率曲线的拟合与风机运行秒级数据筛选,具体包括如下步骤:
步骤2.1:最小二乘法拟合风功率曲线
由于风机生产商标准功率曲线以列表的形式提供,精度为0.5米,风机标准功率曲线见附表1。
表1 WT02093风机标准功率曲线表
从表中可以看出,出厂制定的切入风速为3米每秒,在风速到达12米每秒时进入满发状态,直到风速大于25米每秒就切出。从3米每秒到12米每秒的风速中,对应的功率并不是线性增加,在使用最小二乘法对这部分风功率曲线进行多项式拟合。多项式拟合可以选择不同阶数的拟合函数:
y=Σθixi
其中,θi为标准功率曲线的多项式参数,x为风机运行秒级数据组中的风速,i为多项式拟合的阶数,y为标准功率曲线下计算的标准功率。本方法比较了2至6阶的拟合函数,其中5阶多项式拟合的效果最好,见附图4。拟合的表达式为:
y=330.8-356.4x+137.5x2-23.8x3+2.142x4-0.07091x5
步骤2.2:数据初步筛选
风机出厂设置的切入风速为3米每秒,在风机的使用过程中切入风速会发生变化,但不会超过5米每秒。因此可以去除风速大于5米每秒的数据来减少数据量,对数据进行初步筛选。筛选规则为:
遍历每一条运行数据,若第i条数据风速大于5米每秒,则去掉第i条数据;若风速小于5米每秒,则保留第i条数据。筛选后的结果如附图5所示,此时数据量减少为230494条。
步骤2.3:通过标准功率曲线筛选数据
在去掉大于5米的风速数据之后,根据标准功率曲线对数据进行进一步筛选。对于功率大于0的情况,筛选公式为:
其中error为实际功率与标准功率曲线的偏差,P为实际功率,y为标准功率曲线下计算的标准功率,θi为标准功率曲线的多项式参数,x为风速S’。筛选规则为去除error大于0.5的数据,即实际功率与标准功率曲线功率偏差超过50%的数据,保留error小于0.5的数据。这一步骤通过编程实现,筛选结果见附图6,筛选后剩余202767条数据。
步骤3:建立基于风机运行状态的风机切入模型,具体包括如下步骤:
步骤3.1:建立风机切入判断逻辑
风机运行秒级数据中包含了风机状态数据,风机状态包括待机、正常发电、维护、故障等。WT02093风机采用的协议为11205,具体协议见附表2。
表2 WT02093风机协议
风机切入过程状态变化为状态代码从其他变化为4或者5,通过程序进行筛选,筛选结果见附图7,共有283组切入风速数据。
步骤3.2:切入风速排序及切出风速选择
通过以上步骤,可以获得满足风机切入情况的风速数据,再根据功率数据将切入情况分成主动切入和被动切入两种。其中主动切入时功率大于0,被动切入时功率小于0。分类后获得的主动切入数据如附图8所示,共有92组主动切入风速数据。然后将风速进行提取,并排序,排序的结果如附表3所示。
表3 主动启动风速排序
在统计学中普遍使用四分位数对样本进行定性分析,因此本方法选择第3个四分位数,即75%以上情况能主动启动的风速作为风机的切入风速。92组数据的第三个四分位数为第69个,因此选择第69组数据的4.7683米每秒作为WT02093风机的校正切入风速。
本申请提供的方法,本方法只要获得某风机的运行秒级数据,以及风机制造厂商提供的风功率曲线表,就能够计算出该风机的主动切入风速,并且在获取的秒级数据量越大的情况下计算越准确。本发明能够利用统计数据计算得出风机的切入风速,方法切实可行,结论合理,有助于老旧风场风机切入风速的调整、主动启动时机的选择等相关应用。结合了风机标准功率曲线数据,大大提高了切入风速筛选和选择的效率。基于数据驱动的优势,充分利用了风机运行的历史数据,而不必在现场长时间对风机运行进行观测和记录,节省了切入风速计算的时间和成本。本方法考虑了风向角对切入风速的影响,使计算结果更为准确。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。