基于风能利用系数与多元多项式回归的风力机功率预测方法与流程

1.本发明属于新能源领域，更具体地，涉及一种基于风能利用系数与多元多项式回归的风力机功率预测方法。


背景技术：

2.化石燃料储量的减少与能源需求的稳步上升给世界带来了新的难题。风能是可再生能源，也是清洁能源，风力发电是风能的主要使用方式。功率曲线是风电机组对输入风能的利用效率和机组整体发电性能优劣的直接反映，体现了在轮毂高度测量的不同风速下风力机的输出功率水平。对风电机组功率曲线进行建模和监测能够及时发现风电机组的运行异常和早期故障，提高机组的可利用率，降低机组维护费用。通常，风电功率预测方法是使用制造商提供的风机功率曲线模型。但是，这条曲线没有考虑到风力机的具体安装位置以及运行后的磨损情况，理论风电数据与实际数据之间的差异会导致附加误差。功率曲线建模使用了许多不同的技术：iec61400
‑
12标准中给出的“bin”方法、参数函数、逻辑回归模型等。但是传统的风速
‑
功率曲线只显示功率与风速的关系，没有考虑其他影响因素，例如风向，转子速度和桨距角等，具有一定的局限性。
3.目前，预测风机输出功率的多变量回归模型主要有人工神经网络模型、高斯过程模型等。上述多变量模型的输入特征主要是通过专家对风机运行特性的分析得出的，但是影响风机功率的因素在不同的风电场环境中是不同的。此外，人工神经网络模型的训练过程通常要花费很多时间，在模型的透明度和解释性上有一定的局限性。


技术实现要素：

4.本发明的目的是基于风机scada数据，采用风能利用系数与多元多项式回归来预测风机输出功率。
5.实现本发明目的的技术解决方案为：
6.一种基于风能利用系数与多元多项式回归的风力机功率预测方法，包括以下步骤：
7.步骤1，对scada采样数据集d0以风速、叶片转速、桨距角以及功率作为特征提取出来，组成数据集d1。
8.步骤2，根据风机特性曲线将风力发电机的运行过程分为恒功率、恒转速、最大功率点跟踪三个阶段。以风速作为分割线，将采样数据集d1分割成三部分，为d11，即风速大于等于额定风速，d12，即风速大于等于额定转速风速且小于额定风速，和d13，即风速小于额定转速风速，表示风机所处的三种状态。对数据集d11按一定比例分为两组数据，记为d11
‑
1和d11
‑
2；对数据集d12按一定比例分为两组数据，记为d12
‑
1和d12
‑
2；对数据集d13按一定比例分为两组数据，记为d13
‑
1和d13
‑
2。d11
‑
1，d12
‑
1和d13
‑
1作为训练集，即建模样本；d11
‑
2，d12
‑
2和d13
‑
2作为验证集，即预测样本。
9.步骤3，根据风机理论基础，将风能利用系数c
p
作为功率建模的中间变量。风机的
桨叶具有空气动力学形状，使风轮绕其轴转动，可将风能转化为机械能。气动方程计算如下：
[0010][0011]
式中，p
w
是可利用的空气动能；ρ是空气密度；r是叶片半径；v是风速。理想情况下，风力发电机从风能中获取的功率可表示为：
[0012][0013][0014]
式中，c
p
(λ,θ)是风能利用系数；λ是叶尖速比，是叶尖线速度与风速的比值；θ是桨距角；ω是叶片角速度。c
p
(λ,θ)是关于叶尖速比λ和桨距角θ的非线性函数，表示为：
[0015][0016]
式中，c1‑
c9为功率参数，由风力发电机的类型决定。利用实际输出功率p和风速v可以计算出实际风能利用系数c
p
：
[0017][0018]
步骤4，分别对恒功率段数据集d11
‑
1,恒转速段数据集d12
‑
1和最大功率点跟踪段数据集d13
‑
1建立多元多项式模型：
[0019]
数据集d11
‑
1：
[0020]
采用多元多项式回归模型，将风能利用系数c
p
(λ,θ)表示为叶尖速比λ和桨距角θ的m1阶多项式，如下式：
[0021]
y＝∑a
ij
x
1i
x
2j
；
[0022]
式中，y表示风能利用系数c
p
，x1表示叶尖速比λ，x2表示桨距角，a
ij
表示系数，其中i,j＝0,1,2,
…
,m1；i j≤m1。采用最小二乘法找到使得代价函数j最小的阶次m1和参数a
ij
，使得多项式模型的精度达到最高。代价函数为：
[0023][0024]
式中，n是样本数量，和y
i
分别是风能利用系数c
p
的预测值和实际值。
[0025]
合适的m1值是使模型预测的均方根误差rmse最小对应的阶数，公式如下：
[0026][0027]
式中，n为验证样本个数；y
i
为验证样本风能利用系数c
p
实际值；是多项式模型c
p
预测值。
[0028]
数据集d12
‑
1：
[0029]
由于风机处于恒转速阶段，认为转速恒为额定转速，用λ
′
表征此状态下风机的叶尖速比：
[0030][0031]
其中转速n
n
取额定转速。
[0032]
同样采用多元多项式回归模型，将风能利用系数c
p
(λ,θ)表示为叶尖速比λ
′
和桨距角θ的m2阶多项式，也用最小二乘法找到最佳阶次m2和参数a
ij
，合适的m2值是使模型预测的rmse最小对应的阶数。
[0033]
数据集d13
‑
1：
[0034]
考虑到低风速时由于变桨距控制技术风机桨距角变化极小，对该数据集建立风能利用系数c
p
关于叶尖速比λ的一元多项式模型，经最小二乘拟合找到多项式的最佳阶次m3和参数a
ij
，合适的m3值是使模型预测的rmse最小对应的阶数。
[0035]
步骤5，利用步骤4中数据集d11
‑
1，d12
‑
1和d13
‑
1建立好的风能利用系数c
p
多项式模型，在验证数据集d11
‑
2，d12
‑
2和d13
‑
2上进行风机输出功率预测，功率表达式为：
[0036][0037]
式中，为风能利用系数c
p
的预测值。
附图说明
[0038]
图1为本发明风机运行状态风速
‑
功率曲线示意图。
[0039]
图2为本发明风机运行状态风速
‑
转速曲线示意图。
[0040]
图3为本发明风机运行状态转速
‑
功率曲线示意图。
[0041]
图4为本发明恒功率段均方根误差随多项式阶次变化关系示意图。
[0042]
图5为本发明恒转速段均方根误差随多项式阶次变化关系示意图。
[0043]
图6为本发明最大功率点跟踪段均方根误差随多项式阶次变化关系示意图。
[0044]
图7为本发明多元多项式回归建模流程示意图。
[0045]
图8为本发明scada验证数据上预测时间
‑
功率曲线示意图。
[0046]
图9为本发明scada验证数据上预测风速
‑
功率散点示意图。
具体实施方式
[0047]
基于多元多项式回归和最大信息系数的风机输出功率预测，包括以下步骤：
[0048]
s1、对scada采样数据集d0将风速、叶片转速、桨距角以及功率作为特征提取出来，组成数据集，记为d1。某台风力发电机为例，额定功率为2mw，scada采样时间为60s，共5915条记录d0。该风机额定风速为14m/s，额定转速风速为8.5m/s，额定转速为15.1rpm，叶片直径105m，空气密度取1.25kg/m3。对5915条记录d0以风速、叶片转速、桨距角以及功率作为特征变量提取，组成数据集d1。
[0049]
s2、根据风机特性曲线将风力发电机的运行状态分为恒功率ab、恒转速bc、最大功率点跟踪cd三个阶段，如图1
‑
图3所示。图1是风速
‑
功率曲线。图2是风速
‑
转速曲线。图3是
转速
‑
功率曲线。以风速作为分割线，将数据集d1分为三部分。处于恒功率阶段的数据集d11共2675条，即风速大于等于额定风速；处于恒转速阶段的数据集d12共2481条，即风速大于等于额定转速风速且小于额定风速；处于最大功率点跟踪阶段的数据集d13共759条，即风速小于额定转速风速，表示风机所处的三种状态。将数据集d11分为d11
‑
1共2275条，和d11
‑
2共400条；d12分为d12
‑
1共2131条，和d12
‑
2共350条；d13分为d13
‑
1共609条和d13
‑
2共150条。数据集d11
‑
1，d12
‑
1和d13
‑
1作为训练集；d11
‑
2，d12
‑
2和d13
‑
2作为验证集。
[0050]
s3、分别对恒功率段数据集d11
‑
1,恒转速段数据集d12
‑
1和最大功率点跟踪段数据集d13
‑
1建立多元多项式模型：
[0051]
恒功率段：
[0052]
在数据集d11
‑
1上训练多元多项式回归模型，将预测量表示为叶尖速比λ和桨距角θ的m1阶多项式，如下式：
[0053]
y＝∑a
ij
x
1i
x
2j
；
[0054]
式中，y表示风能利用系数c
p
，x1表示叶尖速比λ，x2表示桨距角，a
ij
表示系数，其中i,j＝0,1,2,
…
,m1；i j≤m1。采用最小二乘法找到使得代价函数j最小的阶次m1和参数a
ij
，使得多项式模型的精度达到最高。代价函数为：
[0055][0056]
式中，n是样本数量，和y
i
分别是风能利用系数c
p
的预测值和实际值。
[0057]
进一步地，均方根误差rmse公式如下：
[0058][0059]
式中，n为验证样本个数；y
i
为验证样本风能利用系数c
p
实际值；是多项式模型c
p
预测值。训练过程rmse随阶数m1变化如图4所示。当多项式阶次m1从1变化到2时，模型的rmse迅速减小后几乎不发生变化。当m1＝2时，训练样本d11
‑
1的rmse为0.0037，验证样本d11
‑
2的rmse为0.0022，此时模型是最佳的，确定适当的m1值为2。建立好的风能利用系数c
p
模型为：
[0060][0061]
式中，y表示风能利用系数c
p
的预测值，x1表示叶尖速比λ，x2表示桨距角θ。
[0062]
恒转速段：
[0063]
由于风机处于恒转速阶段，认为转速恒为额定转速，该状态下叶尖速比λ
′
为：
[0064][0065]
其中转速n
n
取额定转速。
[0066]
将风能利用系数c
p
(λ,θ)表示为叶尖速比λ
′
和桨距角θ的m2阶多项式，也用最小二乘法找到最佳阶次m2和参数a
ij
。合适的m2值是使模型预测的rmse最小对应的阶数。训练过程中rmse随阶数m2变化如图5所示。训练样本d12
‑
1和测试样本d12
‑
2的rmse随多项式阶次
增大先增大后减小，当多项式阶次大于5时，模型出现过拟合的现象。因此当阶次m2＝5时，模型效果最好，训练样本的rmse为0.0215，验证样本的rmse为0.0344。建立好的风能利用系数c
p
模型为：
[0067][0068]
式中，y表示风能利用系数c
p
预测值，x1表示叶尖速比λ
′
，x2表示桨距角θ。
[0069]
最大功率点跟踪段：
[0070]
考虑到低风速时由于变桨距控制技术，风机的桨距角变化极小，在数据集d13
‑
1上建立风能利用系数c
p
关于叶尖速比λ的一元多项式模型，经最小二乘拟合找到多项式的最佳阶次m3和参数a
ij
，合适的m3值是使模型预测的rmse最小对应的阶数。训练过程中rmse随阶数m3变化如图6所示。当多项式阶次m3＝3时，验证样本d13
‑
2的rmse最低，为0.0107。建立好的风能利用系数c
p
模型为：
[0071]
y＝2.9798
‑
0.9899x 0.1091x2‑
0.37
·
10
‑2x3；
[0072]
式中，y表示风能利用系数c
p
预测值，x表示叶尖速比λ。
[0073]
s4、利用s3建立好的风能利用系数c
p
多项式模型，在验证数据集d11
‑
2，d12
‑
2和d13
‑
2上进行风机输出功率预测,功率表达式为：
[0074][0075]
式中，为预测功率值。上述具体建立多元多项式回归模型的步骤如图7所示。结果如图8和图9所示。图8是测试集上预测时间
‑
功率图，实线表示实际的scada数据，虚线表示预测值。图9是测试集上预测风速
‑
功率散点图，圆圈表示真实数据值，加号表示预测值。预测结果显示风速和功率呈“带状”关系。从两个图中可以看出，该模型可以预测大多数功率值。
[0076]
该验证数据预测值的平均绝对百分比误差mape为6.29％，公式如下：
[0077][0078]
式中n是验证集总数量，y
i
是验证样本的实际输出功率，是模型功率预测值。
[0079]
本发明优点如下：
[0080]
1)本发明采用风能利用系数来选择影响风机输出功率的变量，减少了专家分析影响因素的困难。
[0081]
2)本发明在风机输出功率预测方面具有较高的准确率和较低的模型复杂度。