主权项 |
一种风力发电机组控制器参数辨识与整定方法,其特征在于,针对风力发电机组控制器参数在设计和优化过程中不易计算与整定的问题,利用Bladed软件模型线性化结合模型降阶算法建立了适用于参数整定的机组线性化模型,然后结合风场风速分布概率进行了控制器PI参数最优整定与综合计算,并基于Bladed软件参数辨识结果计算了转矩控制的最优增益系数及自适应桨距控制的增益因子等参数;所述风力发电机组控制器参数辨识与整定方法具体步骤如下:步骤1:基于Bladed的风电机组模型线性化根据机组设备参数由Bladed建立风电机组模型,使用Bladed软件模型线性化进行模型计算,然后通过Matlab软件对计算结果进行转换处理,得到各风速点用于桨距控制器PI参数整定的桨距‑转速传递函数模型和用于转矩控制器PI参数整定的转矩‑转速传递函数模型;步骤2:线性化模型降阶处理对转矩‑转速模型和桨距‑转速模型采用最小二乘拟合法降阶为一阶惯性时延环节模型;步骤3:转矩控制器PI参数整定1)利用Bladed软件辨识机组转速稳定在临界值ω<sub>s</sub>=(S1+S4)/2时的风速v<sub>s</sub>,其中,S1为机组最小转速,S4为机组最大转速;2)根据风场风速分布概率密度按照下式计算切入风速区间[v<sub>in</sub>,v<sub>s</sub>]的风速概率P<sub>in</sub>与额定风速区间[v<sub>s</sub>,v<sub>e</sub>]的风速概率P<sub>e</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>in</mi></msub><mo>=</mo><msubsup><mo>∫</mo><msub><mi>v</mi><mi>in</mi></msub><msub><mi>v</mi><mi>s</mi></msub></msubsup><mi>f</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mi>dv</mi></mtd></mtr><mtr><mtd><msub><mi>P</mi><mi>e</mi></msub><mo>=</mo><msubsup><mo>∫</mo><msub><mi>v</mi><mi>s</mi></msub><msub><mi>v</mi><mi>e</mi></msub></msubsup><mi>f</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mi>dv</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000599942310000011.GIF" wi="386" he="273" /></maths>其中,v<sub>in</sub>为切入风速,v<sub>e</sub>为额定风速,f(v)为风速分布概率密度函数;3)根据切入风速点模型和额定风速点模型分别进行最优PI参数整定,得到PI参数[Kp<sub>in</sub>,Ti<sub>in</sub>]和[Kp<sub>e</sub>,Ti<sub>e</sub>];4)对[Kp<sub>in</sub>,Ti<sub>in</sub>]和[Kp<sub>e</sub>,Ti<sub>e</sub>]按照下式进行综合计算,得到转矩控制器的最终PI参数[Kp<sub>Q</sub>,Ti<sub>Q</sub>]:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>Kp</mi><mi>Q</mi></msub><mo>=</mo><msub><mi>Kp</mi><mi>in</mi></msub><mo>·</mo><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><msub><mi>Kp</mi><mi>e</mi></msub><mo>·</mo><msub><mi>w</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>Ti</mi><mi>Q</mi></msub><mo>=</mo><msub><mi>Ti</mi><mi>in</mi></msub><mo>·</mo><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><msub><mi>Ti</mi><mi>e</mi></msub><mo>·</mo><msub><mi>w</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>P</mi><mi>in</mi></msub><mrow><msubsup><mo>∫</mo><msub><mi>v</mi><mi>in</mi></msub><msub><mi>v</mi><mi>e</mi></msub></msubsup><mi>f</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mi>dv</mi></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>w</mi><mn>2</mn></msub><mo>=</mo><mfrac><msub><mi>P</mi><mi>e</mi></msub><mrow><msubsup><mo>∫</mo><msub><mi>v</mi><mi>in</mi></msub><msub><mi>v</mi><mi>e</mi></msub></msubsup><mi>f</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mi>dv</mi></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000599942310000021.GIF" wi="589" he="594" /></maths>其中,w<sub>1</sub>、w<sub>2</sub>为计算权值;步骤4:桨距控制器PI参数整定1)将[v<sub>e</sub>,v<sub>out</sub>]的风速范围分成m个风速区间:[v<sub>e</sub>,v<sub>1</sub>],[v<sub>1</sub>,v<sub>2</sub>],[v<sub>2</sub>,v<sub>3</sub>],…[v<sub>m‑1</sub>,v<sub>out</sub>],其中,v<sub>out</sub>为切出风速;2)根据风速分布概率密度曲线,计算每个区间的风速概率P<sub>i</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><msubsup><mo>∫</mo><msub><mi>v</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>v</mi><mi>i</mi></msub></msubsup><mi>f</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mi>dv</mi></mrow>]]></math><img file="FDA0000599942310000022.GIF" wi="353" he="125" /></maths>其中,f(v)为风速分布概率密度函数;3)基于每个风速区间中点<img file="FDA0000599942310000024.GIF" wi="363" he="87" />的桨距‑转速模型,进行最优整定得到一组PI参数[Kp<sub>i</sub>,Ti<sub>i</sub>];4)将整定得到的各PI参数[Kp<sub>i</sub>,Ti<sub>i</sub>]按照下式进行综合计算,即可得到桨距控制器最终PI参数[Kp<sub>β</sub>,Ti<sub>β</sub>]:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>Kp</mi><mi>β</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>Kp</mi><mi>i</mi></msub><mo>·</mo><msub><mi>w</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Ti</mi><mi>β</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>Ti</mi><mi>i</mi></msub><mo>·</mo><msub><mi>w</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>P</mi><mi>i</mi></msub><mrow><msubsup><mo>∫</mo><msub><mi>v</mi><mi>e</mi></msub><msub><mi>v</mi><mi>out</mi></msub></msubsup><mi>f</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mi>dv</mi></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000599942310000023.GIF" wi="422" he="573" /></maths>其中,w<sub>i</sub>为计算权值;步骤5:判定系统内部稳定性对由步骤3和步骤4整定得到的闭环转矩和桨距PI控制系统进行内部稳定性判定;步骤6:最优转速‑转矩控制器参数辨识与计算利用Bladed分别辨识最优桨距角、最大风能利用系数和最佳叶尖速比,然后利用辨识参数计算最优增益系数;步骤7:自适应变桨控制参数设置1)利用Bladed软件辨识出额定风速以上各个风速区间中点的桨距角稳态值β<sub>i</sub>;2)将步骤4整定的PI参数及桨距角辨识结果制成如表3所示的自适应变桨控制参数表格;表3 自适应变桨距控制PI参数表<img file="FDA0000599942310000031.GIF" wi="1332" he="523" />3)计算增益因子补入表3中,根据表中增益因子随桨距角变化实现自适应变桨控制。 |