发明名称 一种降低风电机组载荷的非线性桨距角控制方法
摘要 本发明涉及一种降低风电机组载荷的非线性桨距角控制方法,属于风力发电技术领域。该方法将传动轴看成一个以气动转矩为控制输入的伪线性系统,采用比例积分和前馈算法计算气动转矩的指令值,然后采用逆方法获得用于转速调节的桨距角指令值;将塔架看成一个以气动推力为控制输入的伪线性系统,气动推力指令值为稳态指令值和小信号指令值之和,然后采用逆方法得到考虑塔架振动阻尼后的桨距角指令值;所得到的桨距角指令值经过串联超前校正环节后,得到最终的桨距角控制指令值,用于风电机组的桨距角控制。本发明在控制风电机组转速恒定的同时,有效降低风电机组的载荷,有助于降低风电机组运维成本和延长风电机组的使用寿命,且适于工程应用。
申请公布号 CN103629047B 申请公布日期 2016.02.10
申请号 CN201310540834.7 申请日期 2013.11.05
申请人 清华大学 发明人 肖帅;耿华;杨耕
分类号 F03D7/00(2006.01)I 主分类号 F03D7/00(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 一种降低风电机组载荷的非线性桨距角控制方法,其特征在于,包括以下各步骤:1)获取风电机组的气动转矩特性函数M<sub>a</sub>(Ω,β,v<sub>0</sub>)、气动推力特性函数F<sub>a</sub>(Ω,β,v<sub>0</sub>)、齿轮箱的变比i、发电机额定电磁转矩<img file="FDA0000793199980000011.GIF" wi="106" he="86" />风轮额定转速Ω<sup>*</sup>和塔架前后振动模态的自然频率f<sub>T</sub>,其中:Ω为风轮转速测量值,β为桨距角测量值,v<sub>0</sub>为风轮有效风速测量值;2)采用编码器测量风轮转速,并经过截止频率为3Ω<sup>*</sup>/2π的低通滤波器后得到风轮转速测量值Ω,其中Ω<sup>*</sup>为步骤1)中的风轮额定转速;3)采用速度传感器测量塔顶速度,并经过以频率f<sub>T</sub>为中心频率的带通滤波器后得到塔顶速度测量值v<sub>T</sub>,其中f<sub>T</sub>为步骤1)中的塔架前后振动模态的自然频率;4)采用风速测量装置获得风轮有效风速测量值v<sub>0</sub>;5)将传动轴看作一个以气动转矩M<sub>a</sub>为控制输入的伪线性系统,采用线性比例积分和前馈算法,根据步骤1)中的齿轮箱的变比i、发电机额定电磁转矩<img file="FDA0000793199980000012.GIF" wi="80" he="84" />和风轮额定转速Ω<sup>*</sup>,以及步骤2)中的风轮转速测量值Ω,计算气动转矩M<sub>a</sub>的控制指令<img file="FDA0000793199980000013.GIF" wi="109" he="82" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>M</mi><mi>a</mi><mo>*</mo></msubsup><mo>=</mo><msub><mi>K</mi><mrow><mi>&Omega;</mi><mo>,</mo><mi>p</mi></mrow></msub><mrow><mo>(</mo><msup><mi>&Omega;</mi><mo>*</mo></msup><mo>-</mo><mi>&Omega;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mrow><mi>&Omega;</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&Integral;</mo><mrow><mo>(</mo><msup><mi>&Omega;</mi><mo>*</mo></msup><mo>-</mo><mi>&Omega;</mi><mo>)</mo></mrow><mi>d</mi><mi>t</mi><mo>+</mo><msubsup><mi>M</mi><mi>g</mi><mo>*</mo></msubsup><mo>&CenterDot;</mo><mi>i</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793199980000014.GIF" wi="1462" he="100" /></maths>其中:K<sub>Ω,p</sub>表示比例系数,K<sub>Ω,i</sub>表示积分系数;6)根据步骤2)中的风轮转速测量值Ω,步骤4)中的风轮有效风速测量值v<sub>0</sub>,以及步骤5)中的气动转矩的控制指令<img file="FDA0000793199980000015.GIF" wi="107" he="79" />计算用于转速调节的桨距角控制指令<img file="FDA0000793199980000016.GIF" wi="96" he="80" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&beta;</mi><mi>&Omega;</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>M</mi><mi>a</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mi>&Omega;</mi><mo>,</mo><msubsup><mi>M</mi><mi>a</mi><mo>*</mo></msubsup><mo>,</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793199980000017.GIF" wi="1313" he="85" /></maths>其中<img file="FDA0000793199980000018.GIF" wi="320" he="85" />为步骤1)中的气动转矩特性函数M<sub>a</sub>(Ω,β,v<sub>0</sub>)的逆,选择满足稳定解条件dM<sub>a</sub>/dΩ<0的稳定解作为控制指令<img file="FDA0000793199980000019.GIF" wi="100" he="82" />7)将塔架看作一个以气动推力F<sub>a</sub>为控制输入的伪线性模型,根据步骤1)中的气动推力特性函数F<sub>a</sub>(Ω,β,v<sub>0</sub>),步骤2)中的风轮转速测量值Ω,步骤4)中的风轮有效风速测量值v<sub>0</sub>,以及步骤6)中的桨距角控制指令<img file="FDA0000793199980000021.GIF" wi="99" he="80" />计算气动推力指令值<img file="FDA0000793199980000022.GIF" wi="65" he="84" />的稳态分量<img file="FDA0000793199980000023.GIF" wi="103" he="84" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>F</mi><mrow><mi>a</mi><mo>,</mo><mn>0</mn></mrow><mo>*</mo></msubsup><mo>=</mo><msub><mi>F</mi><mi>a</mi></msub><mrow><mo>(</mo><mi>&Omega;</mi><mo>,</mo><msubsup><mi>&beta;</mi><mi>&Omega;</mi><mo>*</mo></msubsup><mo>,</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793199980000024.GIF" wi="1110" he="84" /></maths>8)根据步骤3)中塔顶速度测量值v<sub>T</sub>,计算气动推力指令值<img file="FDA0000793199980000025.GIF" wi="65" he="76" />的小信号分量<img file="FDA0000793199980000026.GIF" wi="108" he="82" /><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>F</mi><mrow><mi>a</mi><mo>,</mo><mi>&Delta;</mi></mrow><mo>*</mo></msubsup><mo>=</mo><mo>-</mo><msub><mi>C</mi><mrow><mi>&Delta;</mi><mi>T</mi></mrow></msub><msub><mi>v</mi><mi>T</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793199980000027.GIF" wi="1072" he="87" /></maths>其中C<sub>ΔT</sub>表示期望的塔架阻尼增量;9)根据步骤7)中的气动推力指令值<img file="FDA0000793199980000028.GIF" wi="64" he="88" />的稳态分量<img file="FDA0000793199980000029.GIF" wi="88" he="86" />和步骤8)中的气动推力指令值<img file="FDA00007931999800000210.GIF" wi="62" he="82" />的小信号分量<img file="FDA00007931999800000211.GIF" wi="110" he="83" />计算气动推力指令值<img file="FDA00007931999800000212.GIF" wi="89" he="78" /><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>F</mi><mi>a</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>F</mi><mrow><mi>a</mi><mo>,</mo><mn>0</mn></mrow><mo>*</mo></msubsup><mo>+</mo><msubsup><mi>F</mi><mrow><mi>a</mi><mo>,</mo><mi>&Delta;</mi></mrow><mo>*</mo></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007931999800000213.GIF" wi="1074" he="85" /></maths>10)根据步骤2)中的风轮转速测量值Ω,步骤4)中的风轮有效风速测量值v<sub>0</sub>,以及步骤9)中的气动推力指令值<img file="FDA00007931999800000214.GIF" wi="85" he="82" />计算考虑塔架振动阻尼后的桨距角指令值<img file="FDA00007931999800000215.GIF" wi="91" he="82" /><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>&beta;</mi><mi>T</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>F</mi><mi>a</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mi>&Omega;</mi><mo>,</mo><msubsup><mi>F</mi><mi>a</mi><mo>*</mo></msubsup><mo>,</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007931999800000216.GIF" wi="1096" he="84" /></maths>其中<img file="FDA00007931999800000217.GIF" wi="290" he="85" />为步骤1)中的气动推力函数F<sub>a</sub>(Ω,β,v<sub>0</sub>)的逆;11)根据上述考虑塔架振动阻尼后的桨距角指令值<img file="FDA00007931999800000218.GIF" wi="87" he="86" />经过传递函数为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>G</mi><mrow><mi>l</mi><mi>e</mi><mi>a</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&alpha;</mi><mi>T</mi><mi>s</mi><mo>+</mo><mn>1</mn></mrow><mrow><mi>T</mi><mi>s</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>,</mo><mi>&alpha;</mi><mo>&gt;</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007931999800000219.GIF" wi="1112" he="133" /></maths>的超前补偿环节后,得到最终的桨距角控制指令值β<sup>*</sup>,并输入到风电机组的变桨距执行器,实现对风电机组桨距角的控制。
地址 100084 北京市海淀区清华园1号