CN106979126B - 基于svr的风力发电机组高风速段有效风速估计方法 - Google Patents

Abstract

本发明公开了一种基于SVR的风力发电机组高风速段有效风速估计方法。该方法包括SVR模型训练和模型在线使用两步。在SVR模型训练的过程中，使用传感器获取训练特征集和目标集，对特征集进行归一化，得到SVR的训练集，使用GA算法选择惩罚参数和核函数参数，得到训练好的SVR模型；在模型在线使用过程中，实时获得机组的输出数据，归一化后输入到训练好的SVR模型中，经过低通滤波器之后，得到最终的有效风速估计值。该方法合理利用了机组的输出数据，能够针对高风速段的风电机组进行有效风速估计，设计过程简单，易于实施，可代替LIDAR测风装置，所得有效风速估计值可用于为减小机组机械载荷提供前馈控制信息和风电场风资源评估，从而提高风电场的经济效益。

Description

基于SVR的风力发电机组高风速段有效风速估计方法

技术领域

本发明涉及风力发电机组控制技术领域，特别涉及风力发电机组的高风速段有效风速估计。

背景技术

风能是一种清洁、成本较低、商业潜力巨大的可再生能源，风力发电技术在近几年得到了突飞猛进的发展。世界风能协会在2015年的世界风能报告中指出，到2020年，全球风电装机容量将达到792.1GW。然而，风电技术的发展仍面临大型机组运维成本高、机械载荷偏大导致机组寿命减少、大规模风电并网难度大等挑战。因此，开发一种风电机组有效风速估计方法，从而为降低风电机组的运行载荷提供前馈控制信息，延长机组使用寿命具有重要的实践意义。

风电机组的有效风速定义为桨叶扫掠面积对应的风速矢量场的空间平均值。有效风速的获取是风力发电系统的关键技术，对于实现最大风能捕获、减少机组各部件的机械载荷及风电场风资源评估具有重要意义。风速具有很强的随机性和间歇性，每一个瞬时的风速大小都不相同，且向大型化趋势发展的风电机组的桨叶扫掠面积日益增大，因此，风电机组的有效风速估计是一个极具挑战性的课题。

目前，工业界有效风速的获取方法通常有两种。一种是在机舱尾部安装风速计，然而该方法只能获得桨叶下风向空间中某一点的风速，且测量误差比较大；另一种是在机舱顶部安装LIDAR(LIght Detection and Ranging)测风装置，该方法虽然能比较精确的获得某一范围内的平均风速，但是LIDAR设备价格十分昂贵，若为风电场的每台机组都安装该设备，则会极大增加风电场的建设和运维成本。

为了解决上述问题，学者们提出了许多风电机组的有效风速估计方法，这些方法大致可以分为两类。一类是基于卡尔曼滤波的方法，该类方法的基本思路是：将气动转矩看成系统状态，在假设风电系统的模型参数精确已知且系统的过程及测量噪声符合高斯分布的前提下，建立系统过程方程和测量方程，使用卡尔曼滤波算法获得气动转矩状态的值，再根据气动转矩与有效风速风速之间的数值关系，使用牛顿迭代法获得有效风速的值。然而，实际中风电机组的模型参数很难准确获得，且系统的噪声也不一定满足高斯分布。另一类方法是基于机器学习的方法，这类方法不需要使用系统的数学模型，而是将机组本身看成测量装置，在离线训练阶段，使用预处理后的历史数据训练选定的机器学习模型，比如神经网络(NN)、支持向量机(SVM)、极限学习机(ELM)等，建立机组输出与有效风速之间的非线性关系，进一步则运用训练好的模型以机组实时输出为模型输入，实时获得机组的有效风速。但是，目前已有的基于机器学习的风速估计方法，在模型输入中通常包含风轮转速、发电机转速和发电功率。需要注意的是，在现代大型风电机组的控制系统中，高风速段的控制目标是维持机组的发电机转速和发电功率为定值，因此，高风速段机组的风轮转速、发电机转速和发电功率并不能反应风速信息，因此将风轮转速、发电机转速和发电功率作为机组高风速段风速估计模型的输入显然不合理，现有的基于机器学习的风速估计方法无法应用于风电机组的高风速段的风速估计。

发明内容

为了合理利用风电机组的输出数据，解决现有风电机组风速估计方法估计误差较大且无法应用于机组高风速运行阶段的问题，本发明提供一种针对高风速段、不需要使用系统数学模型、简单易行的风电机组有效风速估计方法，能够比较准确地建立机组输出数据与有效风速之间的非线性关系，获取的有效风速估计值能够为减小机组的机械载荷提供前馈控制信息，同时可应用于风电场的风资源评估。

本发明解决其技术问题所采用的技术方案是：一种基于SVR的风电机组高风速段有效风速估计方法，包括以下步骤：

(1)使用LIDAR测风装置获得一段时间内的有效风速信息，使用SCADA系统和载荷传感器获得相应时间段内风电机组的相关输出数据，机组的相关输出数据用X'表示，(X'＝[x'(i,j)],i＝1,...,l,j＝1,...,6)。用x'(i,:)表示机组的一次采样输出，x'(i,:)表达式为：

x'(i,:)＝[β,β_r,B₁,a_f-d,a_f-v,a_f-a]

其中，β是桨距角，β_r是桨距角加速度，B₁是叶片1的前后偏移量，a_f-d是塔架前后偏移量，a_f-v是塔架前后速度，a_f-a是塔架前后加速度；

(2)将步骤1获得的机组输出数据进行归一化处理，作为SVR模型的训练特征集X(X＝[x(i,j)],i＝1,...,l,j＝1,...,6)，步骤1获得有效风速信息作为SVR模型的训练目标值，将训练特征集和训练目标值作为SVR的训练集；

(3)使用步骤2获得的训练集求解SVR的原始优化问题，为求解该优化问题，引入拉格朗日函数，然后得到对偶优化问题；

(4)使用GA算法选择惩罚参数和核函数参数，求解步骤3中的对偶优化问题，得到训练好的SVR模型；

(5)在线使用时，将某一控制周期内的机组输出数据进行归一化处理，然后输入到步骤4得到的训练好的SVR模型中，得到每一个采样周期的初步风速估计值。

(6)将步骤5得到的初步风速估计值输入到低通滤波器中，得到最终的风速估计值。

进一步地，所述步骤2中，归一化处理指的是：

其中，用x'(:,j)表示X'中的列分量，max(x'(:,j))和min(x'(:,j))分别是x'(:,j)的最大值和最小值，x(:,j)是X中的列分量。

进一步地，所述步骤2中，SVR模型指的是

y＝<w,φ(x)>+b

其中，是模型输出，是模型输入，是将x从n维映射到N维的函数，是偏置项。

进一步地，所述步骤3中，SVR的原始优化问题是

s.t.y_i-<w,φ(x(i,:))>-b≤ε+ξ_i,i＝1,2,...,l

ξ_i≥0,i＝1,2,...,l

其中，C是惩罚参数，l是SVR训练集中的样本个数，ξ_i和是松弛变量，ε是ε-不敏感函数的参数。

进一步地，所述步骤3中，拉格朗日函数的形式为：

其中，是拉格朗日乘子。

进一步地，所述步骤3中，对偶优化问题的形式是：

其中，K(x(i,:),x(j,:))是核函数，本发明中采用高斯核函数，即

其中σ²是核函数参数。

进一步地，所述步骤4中，GA算法的适应度函数选取为训练SVR模型时产生的均方误差，该算法包括个体编码、产生初始群体、适应度计算、选择运算、交叉运算和变异运算六个步骤。

进一步地，所述步骤4中，训练好的SVR模型，其形式为

其中，和是对偶最优问题的解，x_new是机组的实时输出,其所包含的物理量与x(i,:)相同。

进一步地，所述步骤6中，低通滤波器的形式为：

其中，τ是滤波器参数。

本发明的有益效果是：合理利用机组的输出数据，针对现代风电机组高风速段采用桨距角控制策略的现状，选择桨距角及其变化率，叶片前后偏移量，塔架前后偏移量、前后速度以及前后加速度作为样本特征，设计了针对风电机组高风速段的有效风速估计方法，能够比较准确地建立机组输出与有效风速之间的非线性关系；该有效风速估计方法设计过程简单，使用GA算法选择全局最优参数，所得有效风速估计值能够为减小机组机械载荷提供前馈控制信息，同时，该有效风速估计值可用于风电场的风资源评估。实际中，该有效风速估计方法可代替LIDAR测风设备，极大减小风电场的建设和运维成本，提高风电场的经济效益。

附图说明

图1为基于SVR的风力发电机组高风速段风速估计方法框架；

图2为18m/s湍流风示意图；

图3为基于SVR的风力发电机组高风速段风速估计方法设计流程图；

图4为有效风速真实值及其估计值对比图；

图5为测试阶段1000s-2000s有效风速估计误差。

具体实施方式

下面结合附图和具体实施例对本发明作进一步详细说明。

本发明提供的一种基于SVR的风电机组高风速段有效风速估计方法，包括下述步骤：

步骤1，使用LIDAR测风装置获得一段时间内的有效风速信息，使用SCADA系统和载荷传感器获得相应时间段内风电机组的相关输出数据，机组的相关输出数据用X'表示，(X'＝[x'(i,j)],i＝1,...,l,j＝1,...,6)。用x'(i,:)表示机组的一次采样输出，x'(i,:)表达式为：

x'(i,:)＝[β,β_r,B₁,a_f-d,a_f-v,a_f-a]

其中，β是桨距角，β_r是桨距角加速度，B₁是叶片1的前后偏移量，a_f-d是塔架前后偏移量，a_f-v是塔架前后速度，a_f-a是塔架前后加速度。现代风电机组在高风速段(风速大于额定风速小于切出风速)，普遍采用的控制策略是维持发电机转矩为额定值，然后通过控制桨距角从而使得发电功率和发电机转速维持在额定转速附近，工业上普遍采用如下的变增益PI控制器：

其中，调节误差定义为e＝ω_r-ω_d，ω_d为额定风轮转速，K_i比例控制增益常数，K_p是基于当前桨距角值，查表得到的变化的积分控制增益参数。风电机组在高风速段运行时，桨距角是机组的控制信号，其变化是控制策略对有效风速变化的响应，因此，要将β纳入到SVR的训练特征集中。

步骤2，将步骤1获得的机组输出数据进行归一化处理，具体方式为

其中，用x'(:,j)表示X'中的列分量，max(x'(:,j))和min(x'(:,j))分别是x'(:,j)的最大值和最小值，x(:,j)是X中的列分量。X将作为SVR模型的训练特征集。此处的SVR模型，其具体形式为

y＝<w,φ(x)>+b

其中，是模型输出的有效风速信息，是模型输入，是将x从n维映射到N维的函数，n＝6，N是一个非常大的数，是偏置项，<w,φ(x)>表示向量w和φ(x)之间的内积。值得注意的是，步骤1获得的有效风速信息并不需要进行归一化处理，而是直接作为SVR模型的训练目标值，因为训练好的SVR模型在在线使用时无法进行反归一化操作。训练特征集和训练目标值构成了SVR的训练集。

步骤3，使用步骤2获得的训练集求解如下的SVR原始优化问题：

s.t.y_i-<w,φ(x(i,:))>-b≤ε+ξ_i,i＝1,2,...,l

ξ_i≥0,i＝1,2,...,l

其中，C是惩罚参数，l是SVR训练集中的样本个数，ξ_i和是松弛变量，ε是ε-不敏感函数的参数。可见，原始优化问题优化变量较多，求解过程复杂，且约束条件中的φ(·)未知。为简化SVR的训练过程，且能比较自然地引入核函数，通过引入拉格朗日函数，得到原始优化问题的对偶优化问题。引入的拉格朗日函数为

其中，是拉格朗日乘子。根据拉格朗日函数L对原始优化变量的偏导为零，得到如下对偶优化问题：

其中，K(x(i,:),x(j,:))是核函数，本发明中采用高斯核函数，即

其中σ²是核函数参数。可见，在对偶优化问题中，只需要求解α_i和减小了计算量，且引入了核函数技巧，实现了将训练特征集从低维空间映射到高维空间。在对偶优化问题中，还有两个参数需要选择，一个是惩罚参数C，另一个是核函数参数σ²。

步骤4，使用GA算法寻找步骤3中参数C和σ²的最优值，选取训练SVR模型时产生的均方误差作为GA算法的适应度函数，该算法包括个体编码、产生初始群体、适应度计算、选择运算、交叉运算和变异运算六个步骤。选好参数C和σ²后，求解对偶优化问题，得到解α_i和α_i ^*，得到训练好的SVR模型，其形式为：

其中x_new是归一化后的机组实时输出，其包含的物理量与x(i,:)相同。参数b可以根据KKT条件求得。

步骤5，在线使用步骤4获得的训练好的SVR模型，将某一控制周期内的机组输出数据x'_new x'_new包含的物理量与x'(i,:)相同，进行归一化处理，得到x_new，将x_new输入训练好的SVR模型中，得到每一个采样周期的初步风速估计值

步骤6，设计带宽合适的低通滤波器对初步风速估计值进行处理，滤除其高频噪声，得到最终的有效风速估计值

其中，τ是低通滤波器参数。最终的有效风速估计值与有效风速真实值之间的误差较小。

实施例

本实施例使用风电技术开发软件GH Bladed和Matlab仿真平台，对本发明方法的有效性进行验证。

图1所示为基于SVR的风力发电机组低风速段风速估计方法框架。实施例中使用1.5MW三叶片水平轴变速风力发电机组模型，其主要参数如下表所示：

风力发电机组基本参数	数值范围
		额定功率	1500KW
功率因数	-0.95～+0.95
		切入风速	3m/s
额定风速	11m/s
		切出风速	25m/s
风轮直径	77m
		扫掠面积	4654㎡
叶片数	3
		齿轮箱传动比	104.494
高速轴惯量	12Kg·m
		发电机惯量	123Kg·m
风轮额定转速	1.803rad/s

控制器采用浙江某风电研究院开发的已在工业上推广的变增益PI控制器，采样周期是0.04s，机组运行时间设置为2000s，用前1000s的数据作为训练数据，后1000s的数据作为测试集。GA算法选择的最优参数分别为：σ²＝4.3758，C＝0.5367，低通滤波器的参数取值为τ＝3.96。图2是实施例中使用的18m/s湍流风，该湍流风由GH Bladed产生，其纵向、横向和垂直方向的湍流密度分别为：10％、8％和5％。

图3是基于SVR的风力发电机组低风速段风速估计方法设计流程图。流程图中的细线箭头表示模型训练过程，粗线箭头表示模型在线使用过程。在模型训练的过程中，首先使用传感器获取SVR模型的训练特征集和目标集，对特征集进行归一化，得到SVR的训练集，在模型训练的过程中，使用GA算法选择惩罚参数和核函数参数，进而得到训练好的有效风速估计模型；在模型在线使用过程中，实时获得机组的输出数据，归一化后输入到训练好的SVR模型中，经过低通滤波器之后，得到最终的有效风速估计值。

图4是测试阶段1000s-2000s有效风速真实值与有效风速估计值之间的对比图。测试阶段的MSE＝1.4094，MAPE＝5.2757％。

图5是测试阶段1000s-2000s有效风速估计误差。

Claims (9)

1.一种基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，该方法包括以下步骤：

(1)使用LIDAR测风装置获得一段时间内的有效风速信息，使用SCADA系统和载荷传感器获得相应时间段内风电机组的相关输出数据X'，X'＝[x'(i,j)],i＝1,...,l,j＝1,...,6；用x'(i,:)表示机组的一次采样输出，x'(i,:)表达式为：

x'(i,:)＝[β,β_r,B₁,a_f-d,a_f-v,a_f-a]

其中，l是SVR训练集中的样本个数，β是桨距角，β_r是桨距角加速度，B₁是叶片1的前后偏移量，a_f-d是塔架前后偏移量，a_f-v是塔架前后速度，a_f-a是塔架前后加速度；

(2)将步骤(1)获得的机组输出数据进行归一化处理，作为SVR模型的训练特征集X，X＝[x(i,j)],i＝1,...,l,j＝1,...,6；步骤(1)获得有效风速信息作为SVR模型的训练目标值，训练目标值不需要归一化；将训练特征集和训练目标值作为SVR的训练集；

(3)使用步骤(2)获得的训练集求解SVR的原始优化问题，为求解该优化问题，引入拉格朗日函数，然后得到对偶优化问题；

(4)使用GA算法选择惩罚参数和核函数参数，求解步骤(3)中的对偶优化问题，得到训练好的SVR模型；

(5)在线使用时，将某一控制周期内的机组输出数据进行归一化处理，然后输入到步骤(4)得到的训练好的SVR模型中，得到每一个采样周期的初步风速估计值；

(6)将步骤(5)得到的初步风速估计值输入到低通滤波器中，滤除高频噪声，得到最终的风速估计值。

2.根据权利要求1所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，所述步骤(2)中，归一化处理指的是：

其中，用x'(:,j)表示X'中的列分量，max(x'(:,j))和min(x'(:,j))分别是x'(:,j)的最大值和最小值，x(:,j)是X中的列分量。

3.根据权利要求1所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，所述步骤(2)中，SVR模型指的是：

y＝<w,φ(x)>+b

其中，是模型输出，是模型输入，φ(·)为将x从n维映射到N维的函数，是偏置项。

4.根据权利要求3所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，所述步骤(3)中，SVR的原始优化问题是

s.t.y_i-<w,φ(x(i,:))>-b≤ε+ξ_i,i＝1,2,...,l

ξ_i≥0,i＝1,2,...,l

其中，C是惩罚参数，ξ_i和是松弛变量，ε是ε-不敏感函数的参数。

5.根据权利要求4所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，所述步骤(3)中，拉格朗日函数的形式为：

其中，η_i,α_i,是拉格朗日乘子。

6.根据权利要求5所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，所述步骤(3)中，对偶优化问题的形式是：

其中，K(x(i,:),x(j,:))是高斯核函数，即

其中σ²是核函数参数。

7.根据权利要求1所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，所述步骤(4)中，GA算法的适应度函数选取为训练SVR模型时产生的均方误差，该算法包括个体编码、产生初始群体、适应度计算、选择运算、交叉运算和变异运算六个步骤。

8.根据权利要求6所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于所述步骤(4)中，训练好的SVR模型，其形式为：

其中，和是对偶最优问题的解，x_new是机组的实时采样输出，其所包含的物理量与x(i,:)相同。

9.根据权利要求1所述的基于SVR的风力发电机组高风速段有效风速估计方法，其特征在于，所述步骤(6)中，低通滤波器的形式为：

其中，τ是滤波器参数。