发明名称 变速变桨风力发电机组基于前馈解耦控制的最大功率追踪控制方法
摘要 本发明公开一种变速变桨风力发电机组基于前馈解耦控制的最大功率追踪控制方法。本发明通过估算风轮气动力矩,根据估算气动力矩结果基于风力发电机组电磁力矩响应模型及传动系统模型设计气动力矩前馈解耦控制器,有效的消除了由于外部气动力矩的随机性和不确定性对最优阻力矩的动态跟踪影响,使轮毂与主轴低速轴连接处的阻力矩完全由电磁力矩闭环控制,高效的实现了最优阻力矩精确动态跟踪,有效的提升了风轮的最大功率追踪能力。
申请公布号 CN103362738B 申请公布日期 2016.03.30
申请号 CN201210103854.3 申请日期 2012.04.11
申请人 北京能高自动化技术股份有限公司 发明人 贾利民;刘展;雷涛
分类号 G05B13/00(2006.01)I 主分类号 G05B13/00(2006.01)I
代理机构 北京正理专利代理有限公司 11257 代理人 张雪梅
主权项 变速变桨风力发电机组基于前馈解耦控制的最大功率追踪控制方法,其特征是该方法考虑风力发电机组传动系统为柔性环节且存在摩擦的情况下,对风力发电机组传动系统进行建模,通过对风电机组气动力距的前馈解耦控制方式消除了风轮气动力矩的波动性和不确定性对阻力矩控制的影响,高效的实现额定风速以下轮毂与主轴低速轴连接处最优阻力矩跟踪,有效的实现风轮风能吸收最大功率追踪;其中风力发电机组传动系统等效为低速轴和高速轴两质量块环节,按照以下方法建模:<math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><msub><mi>d&omega;</mi><mi>b</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><msub><mi>T</mi><mi>b</mi></msub><msub><mi>J</mi><mi>b</mi></msub></mfrac><mo>-</mo><mfrac><mi>T</mi><msub><mi>J</mi><mi>b</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>d&omega;</mi><mi>g</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>&eta;</mi><mi>T</mi></mrow><mrow><msub><mi>NJ</mi><mi>g</mi></msub></mrow></mfrac><mo>-</mo><mfrac><msub><mi>T</mi><mi>e</mi></msub><msub><mi>J</mi><mi>g</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><mi>d</mi><mi>T</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mi>K</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>b</mi></msub><mo>-</mo><mfrac><msub><mi>&omega;</mi><mi>g</mi></msub><mi>N</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>d&omega;</mi><mi>b</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mfrac><mrow><msub><mi>d&omega;</mi><mi>g</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000829656940000011.GIF" wi="828" he="486" /></maths>式中ω<sub>b</sub>表示低速轴转速;T<sub>b</sub>表示风轮气动力矩;J<sub>b</sub>表示低速轴转动惯量;T表示轮毂与主轴低速轴连接处与旋转方向相反的阻力矩;ω<sub>g</sub>表示高速轴转速;J<sub>g</sub>表示高速轴转动惯量;T<sub>e</sub>表示发电机电磁力矩;N表示传动系统增速比;η表示传动系统传动效率;K表示传动系统等效到低速轴刚性系数;B表示传动系统等效到低速轴阻尼系数;设状态变量为<math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>&omega;</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&omega;</mi><mi>g</mi></msub></mtd></mtr><mtr><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000829656940000012.GIF" wi="214" he="230" /></maths>输入变量为<math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>T</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>e</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000829656940000013.GIF" wi="203" he="151" /></maths>输出为y=T,则传动系统表示为一个三阶模型:<math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><msub><mi>J</mi><mi>b</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mfrac><mi>N</mi><msub><mi>J</mi><mi>t</mi></msub></mfrac></mtd></mtr><mtr><mtd><mi>K</mi></mtd><mtd><mrow><mo>-</mo><mfrac><mi>K</mi><mi>N</mi></mfrac></mrow></mtd><mtd><mrow><mo>-</mo><mfrac><mi>B</mi><msub><mi>J</mi><mi>b</mi></msub></mfrac><mo>-</mo><mfrac><mi>B</mi><msub><mi>J</mi><mi>t</mi></msub></mfrac></mrow></mtd></mtr></mtable></mfenced><mi>x</mi><mo>+</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mfrac><mn>1</mn><msub><mi>J</mi><mi>b</mi></msub></mfrac></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><msub><mi>J</mi><mi>g</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mfrac><mi>B</mi><msub><mi>J</mi><mi>b</mi></msub></mfrac></mtd><mtd><mfrac><mi>B</mi><mrow><msub><mi>NJ</mi><mi>g</mi></msub></mrow></mfrac></mtd></mtr></mtable></mfenced><mi>u</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mi>x</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000829656940000021.GIF" wi="909" he="563" /></maths>其中<math><![CDATA[<mrow><msub><mi>J</mi><mi>t</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>N</mi><mn>2</mn></msup><msub><mi>J</mi><mi>g</mi></msub></mrow><mi>&eta;</mi></mfrac></mrow>]]></math><img file="FDA0000829656940000025.GIF" wi="228" he="150" /></maths>J<sub>t</sub>表示低速轴等效到高速轴后的转动惯量;则系统传递函数表示为:<math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mi>y</mi><mo>(</mo><mi>s</mi><mo>)</mo><mo>=</mo><mo>&lsqb;</mo><msub><mi>G</mi><mn>1</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>,</mo><msub><mi>G</mi><mn>2</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>&rsqb;</mo><mi>u</mi><mo>(</mo><mi>s</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>G</mi><mn>1</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>=</mo><mfrac><mrow><msub><mi>BJ</mi><mi>t</mi></msub><mi>s</mi><mo>+</mo><msub><mi>KJ</mi><mi>t</mi></msub></mrow><mrow><msub><mi>J</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mi>B</mi><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow><mi>s</mi><mo>+</mo><mi>K</mi><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>G</mi><mn>2</mn></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>=</mo><mfrac><mrow><msub><mi>BJ</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub><mi>s</mi><mo>+</mo><msub><mi>KJ</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub></mrow><mrow><msub><mi>NJ</mi><mi>g</mi></msub><msub><mi>J</mi><mi>b</mi></msub><msub><mi>J</mi><mi>t</mi></msub><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><msub><mi>BNJ</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow><mi>s</mi><mo>+</mo><msub><mi>KNJ</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>J</mi><mi>b</mi></msub><mo>+</mo><msub><mi>J</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>.</mo></mrow>]]></math><img file="FDA0000829656940000022.GIF" wi="1167" he="471" /></maths>
地址 100044 北京市海淀区西直门外北京交通大学西门交大知行大厦六层
您可能感兴趣的专利