变速恒频风力发电机动态最优能量捕获方法技术

技术编号:3391102 阅读:278 留言:0更新日期:2012-04-11 18:40
本发明专利技术涉及一种变速恒频风力发电机动态最优能量捕获方法,包括风轮气动转矩测量计算方法和动态最优能量捕获调节方法两部分,风轮气动转矩测量计算方法是:在低速轴两端安装扭转角度测量仪,实时测量出低速轴的扭转角度,将低速轴的扭转刚度和扭转阻尼系数测量出来,根据低速轴的这两个参数以及实时测量的低速轴的扭转角度,可以实时测量计算出低速轴的扭转力矩,根据低速轴的扭转力矩和风轮转速的变化计算出瞬时风轮的气动转矩,所述动态最优能量捕获调节方法是根据测量计算出的瞬时风轮的气动转矩以及功率的变化梯度计算出最优能量捕获转矩,施加到发电机轴上,该发电机轴上转矩再通过齿轮箱施加到低速轴,从而调节风轮转速,使风轮达到最优能量捕获状态。

【技术实现步骤摘要】

本专利技术涉及一种。主要应 用于现代水平轴升力型变速恒频风力发电机,可以提高风力发电系统在低 于额定风速情况下的从湍流风中捕获的能量,增加风力发电机的年发电量。
技术介绍
现代风力发电机一般为水平轴升力型,也就是说利用发电机叶片的特 殊设计的气动外形,在风速和风轮转速共同作用下在叶片的每一段上产生 升力,在升力的作用下,产生旋转力矩,带动发电机发出电能。相对于阻 力型风力发电机,现代水平轴升力型风力发电机具有效率高的特点。为了 能够在不同的风速下,获得最大的能量捕获,现代大型水平轴升力型风力 发电机, 一般采用变速恒频方式操作风力发电机。因此,会在不同的风速 下采用不同的控制策略。根据参考文献l,在风速低于额定风速下工作时, 风力发电机主要是工作在欠功率状态,即叶片的桨距角被调整到最佳桨距角位置,调节风轮转速,达到较高的能量捕获;在风速高于额定风速情况 下,保持转速不变,调节桨距角,减少风轮能量捕获,进入恒功率工作状 态。图1是IEC标准中的三类风场的风速分布,我国大部分风场介于二类风场和三类风场之间,具有很强的代表性。以1.5兆瓦变速恒频风力发电 机为例,风轮直径可以达到70-77米,额定风速在12米/秒左右。风 力发电机主要年工作小时数,从切入风速(4m/s左右)到额定风速(12m/s 左右)为大约6000小时,而额定风速以上的年工作小时数大约为1381小 时。这样,提高在风力发电机工作在低于额定风速情况下的发电效率,将 对于提高风力发电机年发电量产生巨大的影响。对于兆瓦级变速恒频风力发电机来说,为了能够优化地捕获风中的能 量,主要依靠通过静态仿真得出的优化转矩转速曲线进行开环控制。如图 2所示的黑色曲线(摘自参考文献l)。而一般的控制系统框图如图3所示。 根据静态曲线控制有以下问题.-1、 静态曲线考虑的是正面吹来的风作用。而实际上,由于风向时刻改 变,而发电机的偏航控制反应较慢,导致偏航误差,风会在左右方向和上 下方向与风轮所成的夹角,使风力发电机处于非最优能量捕获状态,降低 了能量捕获。2、 由于静态曲线根据稳态风得出,而实际的风场为三维湍流风场。这 样,就造成风轮获得的升力与静态情况有所不同,导致能量捕获下降。3、 静态曲线一般根据叶片数据,有常见的载荷仿真软件,如GHBladed 和Flex5,根据叶片数据直接计算出来,由于实际叶片和叶片数据模型之 间存在误差,经过计算得出的曲线本身就存在误差,使用这个曲线控制时, 也会降低能量捕获。
技术实现思路
本专利技术的目的是针对上述现有技术的不足,克服上述恒频变速风力发 电机在低于额定风速下的静态非最优能量捕获的问题,提供一种变速恒频 风力发电机动态最优能量捕获方法。本专利技术的目的是这样实现的 一种变速恒频风力发电机动态最优能量捕 获方法,其特征在于所述方法包括风轮气动转矩测量计算方法和动态最优 能量捕获调节方法两部分,所述的风轮气动转矩测量计算方法是在低速轴两端安装扭转角度测 量仪,该测量仪测量低速轴在工作状态下,实时测量出低速轴的扭转角度, 扭转角度以弧度为单位,低速轴在锻造的时候通过做实验的方法将低速轴 的扭转刚度和扭转阻尼系数测量出来,根据低速轴的这两个参数以及实时 测量的低速轴的扭转角度,可以实时测量计算出低速轴的扭转力矩,根据 低速轴的扭转力矩和风轮转速的变化计算出瞬时风轮的气动转矩,所述动态最优能量捕获调节方法是根据测量计算出的瞬时风轮的气动 转矩以及功率的变化梯度计算出最优能量捕获转矩,施加到发电机轴上, 该发电机轴上转矩再通过齿轮箱施加到低速轴,从而调节风轮转速,使风 轮达到最优能量捕获状态,达到动态最优能量捕获的目的,具体计算过程 如下,所述瞬时风轮的气动转矩的计算公式如下-r"=力ty二+r =+尺,e+w , 其中,r。为风轮获得的气动转矩,r为低速轴与齿轮箱链接处的扭矩,单 位为牛顿米,,为风轮的转动惯量,单位是公斤米平方,w;为风轮转速的导数,6为低速轴两端扭转角度,单位是弧度,^为低速轴的扭转角度的导 数,《,是低速轴的扭转刚度,单位是牛顿米/弧度,5、为低速轴的扭转阻 尼系数,单位是牛顿米/(弧度/秒)。所述动态最优能量捕获调节方法的控制公式为7; (/t+1)=厶("—-1)) / r、.+^+s, c—维—1)) / r、 + 7(《w - d i)) / ;其中,7;为发电机的电转矩,单位是牛顿米,(t+i), (w,和(/t-i)分 别表示下一时刻,当前时刻,和前一时刻对应变量的值,令7 = /《,《为 测量的电功率。"作为归一化的调节参数,"的调节范围为0.1-0.3。附图说明图1是以往IEC标准中的三类风场的风速分布图。图2是优化转矩转速曲线图。图3是通常的控制系统框图。图4是本专利技术的方案框图。图5是本专利技术的风力发电机简化模型图。具体实施方式本专利技术所涉及的具体的方案框图4所示,方案由两部分组成风轮气动转矩测量计算方法和动态最 优能量捕获调节方法。所述风轮气动转矩测量计算所述的风轮气动转矩测量计算方法包括测 量和计算两部分。第一部分、风轮气动转矩测量风轮气动转矩测量是通过设置在低速轴两端的扭转角度测量仪,该测 量仪测量低速轴在工作状态下,实时直接测量出低速轴的扭转角度入以 弧度为单位。再根据低速轴的扭转刚度和扭转阻尼系数(低速轴在锻造的 时候通过做实验的方法将低速轴的扭转刚度和扭转阻尼系数测量出来),根 据低速轴的这两个参数以及实时测量的低速轴的扭转角度,可以实时测量 出低速轴的扭转力矩。根据低速轴的扭转力矩和风轮角速度的变化计算出 瞬时风轮的气动转矩。风轮的气动转矩直接跟风的湍流密度,偏航误差, 以及风的垂直入射角度有关。因此,直接测量的风轮捕获的气动转矩考虑 到了风速湍流密度,偏航误差,以及风的垂直入射角度。根据测算出的风 轮捕获的气动转矩,叠加功率的梯度信息,设计了动态优化跟踪调节算法。风轮气动转矩计算风轮捕获的气动转矩直接反映了风力发电机的能量捕获状态,因此如 何获得风轮捕获的气动转矩,是最优能量捕获算法的基础。根据参考文献 ,可知风轮捕获的气动转矩与风速和风轮转速以及叶片特性相关7;=上一3~^2 (i)其中,7;为风轮获得的气动转矩,单位为牛顿米,p为空气密度,单 位公斤/立方米,冗为圆周率,/ 为风轮半径,单位为米, 〃为风轮平面上 的等效风速,单位为米/秒,G为功率捕获系数,入为叶尖速比。由于风力发电机的风轮平面很大,例如1.5MW的风力发电机风轮平面面积至少为 3800平方米,因此测量整个风轮平面的等效风速非常困难,所以这个公式在实际的控制中无法使用。根据风力发电机的简化模型(如图5),可以得出如下关系<formula>formula see original document page 9</formula>其中,r。为风轮获得的气动转矩,r为低速轴与齿轮箱链接处的扭矩,z;为 发电机电转矩,转矩的单位为牛顿米,入和^分别为风轮和发电机的转动 惯量,单位是公斤米平方,A和^分别是风轮和电动机的转速,单位是弧 度/秒,w;和《分别是风轮和电动机的转速的导数,e为低速轴两端扭转角 度,单位是弧度,伊为低速轴的扭转角度的导数,K、是低速轴的扭转刚度, 单位是牛本文档来自技高网
...

【技术保护点】
一种变速恒频风力发电机动态最优能量捕获方法,其特征在于所述方法包括风轮气动转矩测量计算方法和动态最优能量捕获调节方法两部分,所述的风轮气动转矩测量计算方法是:在低速轴两端安装扭转角度测量仪,该测量仪测量低速轴在工作状态下,实时测量出低速轴的扭转角度,扭转角度以弧度为单位,低速轴在锻造的时候通过做实验的方法将低速轴的扭转刚度和扭转阻尼系数测量出来,根据低速轴的这两个参数以及实时测量的低速轴的扭转角度,可以实时测量计算出低速轴的扭转力矩,根据低速轴的扭转力矩和风轮转速的变化计算出瞬时风轮的气动转矩,所述动态最优能量捕获调节方法是根据测量计算出的瞬时风轮的气动转矩以及功率的变化梯度计算出最优能量捕获转矩,施加到发电机轴上,该发电机轴上转矩再通过齿轮箱施加到低速轴,从而调节风轮转速,使风轮达到最优能量捕获状态,具体计算过程如下,所述瞬时风轮的气动转矩的计算公式如下:T↓[a]=J↓[r]ω′+T=J↓[r]ω′↓[r]+K↓[s]θ+B↓[s]θ′,其中,T↓[a]为风轮获得的气动转矩,T为低速轴与齿轮箱链接处的扭矩,单位为牛顿米,J↓[r]为风轮的转动惯量,单位是公斤米平方,ω′↓[r]为风轮转速的导数,θ为低速轴两端扭转角度,单位是弧度,θ′为低速轴的扭转角度的导数,K↓[s]是低速轴的扭转刚度,单位是牛顿米/弧度,B↓[s]为低速轴的扭转阻尼系数,单位是牛顿米/(弧度/秒),所述动态最优能量捕获调节方法的控制公式为:T↓[e](k+1)=J↓[r](ω↓[r](k)-ω↓[r](k-1))/T↓[s]+K↓[s]θ+B↓[s](θ(k)-θ(k-1))/T↓[s]+η(P↓[e](k)-P↓[e](k-1))/T↓[s]其中,T↓[e]为发电机的电转矩,单位是牛顿米,(k+1),(k),和(k-1)分别表示下一时刻,当前时刻,和前一时刻对应变量的值,令η=α/P↓[e],其中P↓[e]为测量的电功率。α作为归一化的调节参数,α的调节范围为0.1-0.3。...

【技术特征摘要】

【专利技术属性】
技术研发人员:张雷
申请(专利权)人:江阴远景能源科技有限公司
类型:发明
国别省市:32[]

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1