一种风力发电机组优化启动控制方法技术

本发明专利技术公开了一种风力发电机组优化启动控制方法，该方法利用时间序列分析法建立风速序列的AR(p)模型，利用卡尔曼滤波算法得到未来10min风速的预测值，再将前10min的平均风速与后10min预测风速相结合来作为执行风机启动动作的依据，充分保证了风机启动后的发电能力，减少了小风速季节风机的不合理启动，延长了设备寿命，提高了设备可用率，减少了自耗电，提高了发电量。

【技术实现步骤摘要】

：本专利技术属于电力系统风力发电
，具体涉及一种风力发电机组优化启动控 制方法。
技术介绍
： 随着我国风电技术的发展和对开发利用可再生能源的重视，风力发电近几年在我 国得到了飞速发展，其中2012年全国累计安装风电机组53764台，装机容量75324.2Mff，2013 年全国累计装机容量91412.89MW，截止2014年12月，全国累计安装风电机组76241台，累计 装机容量114609MW。 风电装机容量的增大势必对风电机组的利用效率和可靠性提出更高的要求。由于 风速、风向总是频繁变化的且具有随机性，风电机组长期工作在交变载荷的作用下，工作环 境恶劣，从而使机组在运行过程中故障率增高，性能下降，发电量减少。为提高机组发电量、 降低故障率，各种有效地风电机组控制策略、算法被提出，但多集中于风力机控制、发电机 控制、并网控制和偏航控制，风力发电机的启动控制却未能取得有效发展。 现行的风力发电机组启动并网的控制策略为:首先，风电机组在自检正常的情况 下，叶轮处于自由运动状态，当连续lOmin风速在风力发电机组运行风速的范围内（3m/ S〈v〈 25m/s)且叶轮正对风向时，桨距角控制器调节叶片桨距角，使风力发电机组启动，发电机开 始加速，当发电机空载转速保持在切入转速时，变流器对励磁电流幅值、相位和频率进行控 制，使得发电机定子空载电压和电网电压同频率、同相位、同幅值，此时闭合并网接触器实 现并网。 上述的控制机组启动的依据是前lOmin的平均风速，并未考虑未来的风速。在小风 速季节，有可能前lOmin的平均风速已经达到切入风速，在机组启动后又降到切入风速以 下，从而出现机组没有并网就停机，或频繁启停、频繁并网脱网的情况。这类启动没有带来 多少发电量，反而增加了机组的自耗电和设备疲劳损耗，使故障率升高。因此需合理规划小 风期风机启动的时机，减少不合理启动，取得在发电量、设备寿命和自耗电三者之间的平 衡。
技术实现思路
： 本专利技术的目的在于针对上述风机启动控制策略在小风期的不足，提出了一种风力 发电机组优化启动控制方法，该方法利用时间序列分析法建立风速序列的AR(p)模型，利用 卡尔曼滤波算法得到未来lOmin风速的预测值，再将前lOmin的平均风速与后lOmin预测风 速相结合来作为执行风机启动动作的依据，充分保证了风机启动后的发电能力，减少了小 风速季节风机的不合理启动，延长了设备寿命，提高了设备可用率，减少了自耗电，提高了 发电量。 为达到上述目的，本专利技术采用如下技术方案来实现的： -种风力发电机组优化启动控制方法，用于小风速季节风电机组的优化启动控 制，包括以下步骤： 1)采集风速历史数据库，每隔lOmin求一次风速平均值，将此平均值设为样本点， 共取得N个样本点{x(k)}，k=l，2,…，N. 2)基于时间序列分析法建立风速数据的AR(p)模型，AR(p)模型的表达式如下： (p{q ')x(k) - w(k) Cl) 其中： p(q^1) = 1 - f>if-···- fpqrP (2:) w(k)为零均值且方差为〇2的白噪声；只为所建模型的待估系数;cT'q1分别为单位 滞后算子和P阶滞后算子;利用AIC准则确定模型阶数p，利用矩估计法确定参数終，i = 1，2， 3,…，p; 3)基于AR(p)模型，建立卡尔曼滤波算法的状态方程和测量方程，如下： 其中，X(k+1)为k+Ι时刻的状态向量;F为系统状态转移矩阵；Γ为激励转移矩阵;Z (k+Ι)为k+Ι时刻的观测向量;Η为测量系统的观测状态矩阵;w(k+l)和v(k+l)分别为过程白 噪声向量和测量白噪声向量，协方差矩阵分别为Q和R; 根据式（1)和式(2)可得：由于X2(k+l)=xi(k)，X3(k+l)=X2(k)，···，xP(k+l)=xP-i(k)，因此可得：根据上述状态方程，建立如下观测方程：其中，v(k+l)为观测噪声;式(7)和式(8)即为卡尔曼滤波的状态方程和观测方程； 4)利用卡尔曼滤波算法对风速序列进行在线滤波，其中卡尔曼滤波的递推方程如 下： X(k+1 |k)=F · X(k|k) (9) P(k+1 I k) =F · P(k | k) · FT+ Γ Q Γ T (10) K(k+l)=P(k+l|k)HT-1 (11) X(k+1 |k+l)=X(k+l |k)+K(k+l) (12) P(k+1 |k+l) = · P(k+1 |k) (13) 其中，X(k+1 |k)表示在k时刻对k+1时刻的状态预测，P(k+1 |k)表示预测值X(k+1 k) 的预测误差协方差矩阵;K(k+1)为卡尔曼增益矩阵;X(k+l|k+l)为k+1时刻的测量值Z(k+ l) 对状态预测值X (k+11 k)的更新，为系统状态的最优估计;P (k+11 k+1)为X (k+11 k+1)的估 计误差协方差矩阵；I为单位矩阵； 5)卡尔曼滤波中X(k+1 |k)中的分量xi(k+l |k)即为下一个lOmin的平均风速预测 值，将其与k时刻前lOmin的实际平均风速求平均值vave3(k)，若v ave3(k)>3m/s，风机在k时刻开 桨启动;若(k) <3m/s，风机在k时刻不执行启动动作。 本专利技术进一步的改进在于:在步骤4)中，为了使卡尔曼滤波尽快收敛，取初值X为滤波开始前P个采样点构成的向量，P = l〇I，Q=I，R=I。 本专利技术进一步的改进在于:在步骤5)中，首先判断卡尔曼滤波是否收敛，判断方法 为:若k>ks时，| |X(k|k)_Z(k) | | < εηι/s，则认为卡尔曼滤波已经收敛，ε表示状态估计值X(k k)与实际测量值Z(k)的允许偏差范围，取值为(0，1]。与现有技术比较，本专利技术根据风速序列的AR(p)模型建立卡尔曼滤波方程的状态 方程和测量方程后，根据步骤4)和步骤5)判断风机在时刻k是否应该启动。经过测量和计算 得到下一个lOmin的平均风速后，依据式(8)建立Z(k+1)，通过式（11)~式（13)对X(k+l|k) 和P(k+11 k)进行更新，得到第k+1个样本点的最优状态估计X(k+11 k+1)和协方差矩阵P(k+1 k+1)，再根据式(9)和式(10)计算X(k+2 | k+1)和协方差矩阵P(k+2 | k+1)，将X(k+2 | k+1)中 的分量X当前第1页1 2 本文档来自技高网...

【技术保护点】
一种风力发电机组优化启动控制方法，其特征在于，用于小风速季节风电机组的优化启动控制，包括以下步骤：1)采集风速历史数据库，每隔10min求一次风速平均值，将此平均值设为样本点，共取得N个样本点{x(k)}，k＝1,2,…,N.2)基于时间序列分析法建立风速数据的AR(p)模型，AR(p)模型的表达式如下：其中：w(k)为零均值且方差为σ2的白噪声；为所建模型的待估系数；q‑1、q‑p分别为单位滞后算子和p阶滞后算子；利用AIC准则确定模型阶数p，利用矩估计法确定参数i＝1,2,3,…,p；3)基于AR(p)模型，建立卡尔曼滤波算法的状态方程和测量方程，如下：X(k+1)=F·X(k)+Γw(k+1)Z(k+1)=H·X(k+1)+v(k+1)---(3)]]>其中，X(k+1)为k+1时刻的状态向量；F为系统状态转移矩阵；Γ为激励转移矩阵；Z(k+1)为k+1时刻的观测向量；H为测量系统的观测状态矩阵；w(k+1)和v(k+1)分别为过程白噪声向量和测量白噪声向量，协方差矩阵分别为Q和R；根据式(1)和式(2)可得：将式(4)变形为：令x1(k)＝x(k)，x2(k)＝x(k‑1)，…，xp(k)＝x(k‑p+1)，则有由于x2(k+1)＝x1(k)，x3(k+1)＝x2(k)，…，xp(k+1)＝xp‑1(k)，因此可得：根据上述状态方程，建立如下观测方程：Z(k+1)=[10...00]·x1(k+1)x2(k+1)...xp-1(k+1)xp(k+1)+v(k+1)---(8)]]>其中，v(k+1)为观测噪声；式(7)和式(8)即为卡尔曼滤波的状态方程和观测方程；4)利用卡尔曼滤波算法对风速序列进行在线滤波，其中卡尔曼滤波的递推方程如下：X(k+1|k)＝F·X(k|k)  (9)P(k+1|k)＝F·P(k|k)·FT+ΓQΓT  (10)K(k+1)＝P(k+1|k)HT[H·P(k+1|k)·HT+R]‑1  (11)X(k+1|k+1)＝X(k+1|k)+K(k+1)[Z(k+1)‑H·X(k+1|k)]  (12)P(k+1|k+1)＝[I‑K(k+1)·H]·P(k+1|k)  (13)其中，X(k+1|k)表示在k时刻对k+1时刻的状态预测，P(k+1|k)表示预测值X(k+1|k)的预测误差协方差矩阵；K(k+1)为卡尔曼增益矩阵；X(k+1|k+1)为k+1时刻的测量值Z(k+1)对状态预测值X(k+1|k)的更新，为系统状态的最优估计；P(k+1|k+1)为X(k+1|k+1)的估计误差协方差矩阵；I为单位矩阵；5)卡尔曼滤波中X(k+1|k)中的分量x1(k+1|k)即为下一个10min的平均风速预测值，将其与k时刻前10min的实际平均风速求平均值vave(k)，若vave(k)>3m/s，风机在k时刻开桨启动；若vave(k)<3m/s，风机在k时刻不执行启动动作。...

【技术特征摘要】

【专利技术属性】
技术研发人员：姚玲玲，王靖程，马剑民，牛瑞杰，郭辰，冯笑丹，
申请(专利权)人：西安热工研究院有限公司，华能集团技术创新中心，
类型：发明
国别省市：陕西;61