发明名称 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法
摘要 本发明公开一种变速变桨风力发电机组基于传动系统传动特性补偿控制的最大功率追踪控制方法。该方法通过在主轴低速端与轮毂连接处安装转矩传感器并将检测的阻力矩反馈实现最优阻力矩闭环控制,有效的使轮毂与主轴低速轴连接处阻力矩跟踪最优阻力矩变化曲线。基于该方法实现的最大功率追踪控制方法相对传统的自寻优控制算法有效的对传动系统的非理想传动特性进行了补偿控制,同时有效的消除了由于风速的随机性和不确定性对最优阻力矩的动态影响,成功的解决了传统的自寻优控制算法进行最大功率追踪时难以准确跟踪最优阻力矩跟踪的情况,有效的提升了风轮的最大功率追踪能力。
申请公布号 CN103362735B 申请公布日期 2015.10.28
申请号 CN201210096910.5 申请日期 2012.04.05
申请人 北京能高自动化技术股份有限公司 发明人 贾利民;刘展;雷涛
分类号 F03D7/00(2006.01)I 主分类号 F03D7/00(2006.01)I
代理机构 北京正理专利代理有限公司 11257 代理人 张雪梅
主权项 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法,其特征是该方法基于对风力发电机组传动系统的传动特性进行补偿控制,通过闭环控制方式实现低于额定风速以下轮毂与主轴低速轴连接处最优阻力矩跟踪,实现风轮风能吸收最大功率追踪控制;传动系统可以等效为低速轴和高速轴两质量块环节,按照以下方法建模:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msub><mi>d&omega;</mi><mi>b</mi></msub><mi>dt</mi></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></mtd></mtr><mtr><mtd><mfrac><msub><mi>d&omega;</mi><mi>g</mi></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><mi>&eta;T</mi><msub><mi>NJ</mi><mi>g</mi></msub></mfrac><mo>-</mo><mfrac><msub><mi>T</mi><mi>e</mi></msub><msub><mi>J</mi><mi>g</mi></msub></mfrac></mtd></mtr><mtr><mtd><mfrac><mi>dT</mi><mi>dt</mi></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><msub><mi>d&omega;</mi><mi>b</mi></msub><mi>dt</mi></mfrac><mo>-</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mfrac><msub><mi>d&omega;</mi><mi>g</mi></msub><mi>dt</mi></mfrac><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000702720790000011.GIF" wi="830" he="493" /></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表示传动系统等效到低速轴阻尼系数;设状态变量为<maths num="0002" id="cmaths0002"><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="FDA0000702720790000012.GIF" wi="217" he="231" /></maths>输入变量为<maths num="0003" id="cmaths0003"><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="FDA0000702720790000013.GIF" wi="208" he="162" /></maths>输出为y=T,则传动系表示为一个三阶模型:<img file="FDA0000702720790000027.GIF" wi="918" he="567" />其中<maths num="0004" id="cmaths0004"><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="FDA0000702720790000022.GIF" wi="233" he="149" /></maths>J<sub>t</sub>表示低速轴等效到高速轴后的转动惯量;则系统传递函数表示为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>y</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><msub><mi>G</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>G</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>]</mo><mi>u</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>G</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><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><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><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>]]></math><img file="FDA0000702720790000023.GIF" wi="1143" he="475" /></maths>
地址 100044 北京市海淀区西直门外北京交通大学西门交大知行大厦六层