发明名称 一种双馈型风力发电机组最大风能捕获方法
摘要 本发明公开了风力发电领域的一种双馈型风力发电机组最大风能捕获方法,构建以风速值、风向值、风轮转速值为输入量,电磁转矩值为输出量的BP神经网络模型,再根据叶尖速比点公式,计算期望的理论最佳风轮转速值,代入所述BP神经网络模型中确定电磁转矩值初始控制量,然后采用变步长的爬山搜索法,随着控制周期数的升高缩小电磁转矩值变化步长,在小范围内改变电磁转矩值大小以寻找最大功率值点,并在每次搜索前测量实时发电机转速值,在发电机转速值相比整个搜索开始时的初始发电机转速值的变化值超出限值时,重新启动搜索,在相邻两个控制周期的发电机转速值或者发电功率值的差值的绝对值小于限值时结束搜索。其技术效果是可以在爬山搜索中避免搜索周期过长或发生振荡,保证最大功率的捕获。
申请公布号 CN104612904A 申请公布日期 2015.05.13
申请号 CN201410743221.8 申请日期 2014.12.08
申请人 上海电气集团股份有限公司 发明人 李卫;孙佳林;杨文斌
分类号 F03D9/00(2006.01)I;G06F17/50(2006.01)I 主分类号 F03D9/00(2006.01)I
代理机构 上海兆丰知识产权代理事务所(有限合伙) 31241 代理人 屠轶凡
主权项 一种双馈型风力发电机组最大风能捕获方法,包括下列步骤:建立BP神经网络模型步骤:采集风机转子平衡状态下多个时间点的风速值,风向值,风轮转速值和电磁转矩值,建立以风速值,风向值和风轮转速值为输入量,电磁转矩值为输出量的BP神经网络模型;电磁转矩值初始控制量计算及控制步骤:依据期望的理论最佳风轮转速值公式,即:<img file="FDA0000627092730000011.GIF" wi="302" he="138" />计算得到期望的理论最佳风轮转速值<img file="FDA0000627092730000012.GIF" wi="88" he="52" />并将其代入所述的BP神经网络模型中,计算得到电磁转矩值初始控制量T(0),并以此时的发电机转速值作为启动爬山搜索前的初始发电机转速值ω<sub>0</sub>,其中,ν为10秒风速平均值,R为风机的风轮半径,λ<sub>opt</sub>为最佳叶尖速比点;电磁转矩值初次扰动步骤:等待T<sub>L</sub>个采样周期,待风机转子处于平衡状态后,记录当前发电机转速值ω<sub>1</sub>及发电功率值P<sub>1</sub>,并根据公式T(n)=T(n‑1)+T<sub>step</sub>×α(n)计算控制周期数n为1时的电磁转矩值T(1),其中α(n)=exp(‑30×(n/N)<sup>p</sup>),N为最大控制周期数,p为1~30间任意一个整数,T<sub>step</sub>为电磁转矩值最大变化步长;发电机转速值判断步骤:等待T<sub>L</sub>个采样周期后,先将控制周期数n加1,再记录控制周期数为n时的发电机转速值ω<sub>n</sub>,并计算ω<sub>n</sub>相比于初始发电机转速值ω<sub>0</sub>的变化值Δω<sub>0</sub>,若|Δω<sub>0</sub>|&gt;Δω<sub>max</sub>,返回电磁转矩值初始值计算及控制步骤,其中Δω<sub>max</sub>为设定的最大发电机转速偏差值,否则进入下一步骤;爬山搜索启动判断步骤:计算相邻两个控制周期间的发电机转速值的差值Δω,以及发电功率值的差值ΔP,其中Δω=ω<sub>n</sub>‑ω<sub>n‑1</sub>,ΔP=P<sub>n</sub>‑P<sub>n‑1</sub>,若|ΔP|≤ΔP<sub>min</sub>或|Δω|≤Δω<sub>min</sub>,则结束本次搜索;若|ΔP|&gt;ΔP<sub>min</sub>且|Δω|&gt;Δω<sub>min</sub>,则进入下一步骤;ΔP<sub>min</sub>为寻优发电功率容许变化值,Δω<sub>min</sub>为寻优发电机转速容许变化值;控制周期数判断步骤:若控制周期n大于最大控制周期数N,结束本次搜索,若控制周期数n小于等于最大控制周期数N,返回发电机转速值判断步骤;爬山搜索电磁转矩值及扰动步骤:根据公式T(n+1)=T(n)+ΔT<sup>*</sup>,计算并输出控制周期数为n+1时的电磁转矩值T(n+1);其中ΔT<sup>*</sup>为电磁转矩值变化步长,爬山搜索确定电磁转矩值的规则如下:当ΔP&gt;0且Δω≥0时,ΔT<sup>*</sup>=‑T<sub>step</sub>×α(n),即T(n+1)=T(n)‑T<sub>step</sub>×α(n);当ΔP&gt;0且Δω&lt;0时,ΔT<sup>*</sup>=T<sub>step</sub>×α(n),T(n+1)=T(n)+T<sub>step</sub>×α(n);当ΔP&lt;0且Δω≥0时,ΔT<sup>*</sup>=T<sub>step</sub>×α(n),T(n+1)=T(n)+T<sub>step</sub>×α(n);当ΔP&lt;0且Δω&lt;0时,ΔT<sup>*</sup>=‑T<sub>step</sub>×α(n),即T(n+1)=T(n)‑T<sub>step</sub>×α(n)。
地址 200336 上海市长宁区兴义路8号30层