姚森，楊柳青

(1.東方電氣風(fēng)電有限公司，四川 德陽(yáng)，618000； 2.東方電氣自動(dòng)控制工程有限公司，四川 德陽(yáng)，618000)

0 引言

塔架振動(dòng)問(wèn)題一直是風(fēng)電機(jī)組頻發(fā)的振動(dòng)故障之一。塔架振動(dòng)過(guò)大嚴(yán)重危害機(jī)組的壽命和運(yùn)行安全，進(jìn)而影響機(jī)組的有效發(fā)電量和可利用率。所以，減少風(fēng)機(jī)塔架振動(dòng)具有重要意義。塔架頂部的橫向振動(dòng)一般由傳動(dòng)鏈扭矩的反作用或者外界激勵(lì)引起[1]，屬于風(fēng)輪與塔架耦合振動(dòng)的一類(lèi)振動(dòng)問(wèn)題。風(fēng)輪轉(zhuǎn)速引起塔架振動(dòng)主要有2個(gè)頻率：（1）風(fēng)輪旋轉(zhuǎn)頻率，即1 P，此激勵(lì)主要是由于葉片不平衡導(dǎo)致；（2）風(fēng)輪旋轉(zhuǎn)頻率乘以葉片數(shù)目倍數(shù)，如3 P，6 P等，主要是氣動(dòng)力不平衡所致[2]。風(fēng)輪與塔架耦合振動(dòng)分析可以采用坎貝爾圖進(jìn)行輔助分析，通過(guò)坎貝爾圖觀察各模態(tài)頻率在風(fēng)輪變速運(yùn)行范圍內(nèi)，是否和1 P、3 P、6 P等包絡(luò)線相交來(lái)進(jìn)行判斷。本文從分析帶阻濾波器頻率特性入手，并結(jié)合風(fēng)機(jī)共振坎貝爾圖，提出了在風(fēng)機(jī)扭矩-轉(zhuǎn)速控制和變槳-轉(zhuǎn)速控制中，合理地使用帶阻濾波器對(duì)風(fēng)機(jī)轉(zhuǎn)速頻譜特定頻率點(diǎn)處進(jìn)行陷波處理，可以實(shí)現(xiàn)抑制風(fēng)機(jī)塔架振動(dòng)的方法。

1 帶阻濾波器的頻率特性

帶阻濾波器是能通過(guò)大多數(shù)頻率分量，但某些范圍的頻率分量被衰減到極低水平的濾波器。其中陷波濾波器（notch filter）是一種特殊的帶阻濾波器，它的阻帶范圍極小，有著很高的Q值（quality factor）。

一個(gè)通用的二階濾波器連續(xù)傳遞函數(shù)形式如式（1）所示：

當(dāng) ω1＝ω2＝ω， 且 ζ1＜ζ2時(shí)， 為二階帶阻濾波器，可以濾除特定頻率ω。阻尼系數(shù)ζ表明了這個(gè)帶阻濾波器的寬度或強(qiáng)度。它應(yīng)該調(diào)整直至在特定頻率上出現(xiàn)明顯的濾波效果為止，并且不應(yīng)在較低頻率上對(duì)系統(tǒng)產(chǎn)生過(guò)大的影響[3]。

下面對(duì)該二階濾波器作頻率特性分析，取ω1=ω2=2.915 rad/s， ζ1=0.01，ζ2=0.1.該陷波濾波器對(duì)數(shù)頻率特性（伯德圖）如圖1所示。

圖1 陷波濾波器伯德圖

由經(jīng)典自動(dòng)控制理論知識(shí)可知，二階振蕩環(huán)節(jié)在轉(zhuǎn)折頻率ωT處的對(duì)數(shù)幅頻特性的幅值為[4]：

因此，傳遞函數(shù)具有式（1）形式的陷波濾波器在轉(zhuǎn)折頻率ωT處的對(duì)數(shù)幅頻特性的幅值為：

由式（3）可知，式（1）所示陷波濾波器的深度由阻尼系數(shù) ζ1與 ζ2的比值決定， 將本例 ζ1＝0.01， ζ2＝0.1 的值代入式（3）得， L（ω＝ωT=2.91）＝-20 dB。對(duì)比圖1所示ω＝2.915 rad/s處對(duì)數(shù)幅頻特性，二者幅值相等，因而驗(yàn)證了陷波濾波器的幅值濾波性能。實(shí)際使用過(guò)程中，可通過(guò)對(duì)阻尼比ζ1與ζ2進(jìn)行調(diào)節(jié)，來(lái)決定期望頻率處的幅值陷波深度。

2 風(fēng)電機(jī)組坎貝爾圖分析

Bladed是一個(gè)用于風(fēng)力發(fā)電機(jī)組性能和載荷計(jì)算的集成化的軟件包，是Garrad Hassan公司編制的機(jī)組性能和載荷計(jì)算、評(píng)估的軟件，可以用于機(jī)組的測(cè)試、認(rèn)證等環(huán)節(jié)。其所采用的模型和理論方法已得到多家公司的機(jī)組設(shè)計(jì)數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)的比對(duì)驗(yàn)證，是風(fēng)電機(jī)組設(shè)計(jì)和載荷驗(yàn)證的可靠工具。

在使用Bladed軟件對(duì)機(jī)組模型完成仿真搭建后，通過(guò)對(duì)風(fēng)力機(jī)的模態(tài)分析參數(shù)進(jìn)行設(shè)定，可得風(fēng)機(jī)運(yùn)行風(fēng)況下各耦合模態(tài)頻率。某風(fēng)力機(jī)穩(wěn)態(tài)風(fēng)速10.5 m/s風(fēng)況下各個(gè)耦合模態(tài)的自然頻率值和阻尼系數(shù)如表1所示。其中，自然頻率值的計(jì)算以P為單位，此處將風(fēng)輪旋轉(zhuǎn)一周的頻率定義為1 P。

表1 風(fēng)力機(jī)耦合頻率示例（穩(wěn)定風(fēng)速10.5 m/s）

續(xù)表1

在Bladed軟件中可運(yùn)行計(jì)算生成風(fēng)電機(jī)組風(fēng)輪轉(zhuǎn)速坎貝爾圖，用以表征風(fēng)機(jī)各個(gè)耦合模態(tài)頻率。某1.5 MW 77型風(fēng)力發(fā)電機(jī)組風(fēng)輪坎貝爾圖如圖2所示。

圖2 某1.5 MW 77型風(fēng)電機(jī)組坎貝爾圖

如圖2所示，機(jī)組塔架一階左右模態(tài)，頻率為0.460 9 Hz，在風(fēng)輪轉(zhuǎn)速13.8 r/min時(shí)，與2 P斜線有交點(diǎn)。塔架二階左右模態(tài)在17 r/min時(shí)與15 P模態(tài)有交點(diǎn)。傳動(dòng)鏈頻率為1.817 Hz，在12 r/min時(shí)與9 P頻率有交點(diǎn)，風(fēng)輪面內(nèi)三階模態(tài)，頻率為3.757 Hz，在15 r/min時(shí)與15 P有交點(diǎn)。

為防止風(fēng)輪與塔架耦合所產(chǎn)生的振動(dòng)問(wèn)題的影響，在機(jī)組部件設(shè)計(jì)過(guò)程中就要避開(kāi)3 P，6 P，9 P，12 P等3 P的整數(shù)倍的振動(dòng)頻率。如果在風(fēng)力發(fā)電機(jī)組運(yùn)行全風(fēng)速段范圍內(nèi)，出現(xiàn)了與3 P整數(shù)倍頻率相交的耦合頻率，就要看其相應(yīng)的阻尼比，阻尼比越小就越容易引起共振[5]。這是設(shè)計(jì)中需要考慮的問(wèn)題，如果不可避免，就要求控制系統(tǒng)進(jìn)行相應(yīng)調(diào)整，使機(jī)組運(yùn)行轉(zhuǎn)速避開(kāi)共振區(qū)。

3 帶阻濾波器在風(fēng)機(jī)塔架振動(dòng)控制中的應(yīng)用

變速風(fēng)力發(fā)電機(jī)組扭矩-轉(zhuǎn)速控制器和變槳-轉(zhuǎn)速控制器都使用發(fā)電機(jī)測(cè)量轉(zhuǎn)速作為主要反饋輸入，通過(guò)穩(wěn)定轉(zhuǎn)速來(lái)實(shí)現(xiàn)風(fēng)能的吸收與電能的輸出控制。當(dāng)發(fā)電機(jī)轉(zhuǎn)速反饋進(jìn)入控制器前使之經(jīng)過(guò)濾波器濾波，不僅能夠增強(qiáng)控制系統(tǒng)在特定頻率域的功能，并且能夠防止轉(zhuǎn)矩和槳距角在不期望的頻率處動(dòng)作，如葉片穿越頻率。如圖3所示扭矩-轉(zhuǎn)速控制環(huán)反饋通道加入濾波器結(jié)構(gòu)圖。

圖3 扭矩-轉(zhuǎn)速控制環(huán)反饋通道加入濾波器結(jié)構(gòu)圖

在風(fēng)機(jī)塔架橫向振動(dòng)問(wèn)題處理中，因?yàn)樗軝M向振動(dòng)主要與主軸轉(zhuǎn)速及傳動(dòng)鏈轉(zhuǎn)矩有關(guān)，如在扭矩-轉(zhuǎn)速控制環(huán)和變槳-轉(zhuǎn)速控制環(huán)對(duì)發(fā)電機(jī)轉(zhuǎn)速信號(hào)加以主要影響塔架橫向振動(dòng)的塔架一階模態(tài)頻率，3 P及6 P頻率處陷波濾波器，可有效抑制塔架的橫向振動(dòng)問(wèn)題。下面以某風(fēng)場(chǎng)實(shí)際運(yùn)行風(fēng)機(jī)振動(dòng)故障解決情況為例進(jìn)一步分析。

某在運(yùn)行風(fēng)場(chǎng)，1.5 MW 77型風(fēng)機(jī)經(jīng)常在10～12 m/s風(fēng)速運(yùn)行區(qū)間報(bào)出塔架橫向振動(dòng)超限故障，通過(guò)對(duì)故障數(shù)據(jù)中轉(zhuǎn)速及槳距角頻譜分析來(lái)看，塔架一階模態(tài)頻率峰值很高，如圖4、圖5所示。

圖4 某77型風(fēng)機(jī)塔架橫向振動(dòng)故障時(shí)轉(zhuǎn)速頻譜圖

圖5 某77型風(fēng)機(jī)塔架橫向振動(dòng)故障時(shí)槳距角頻譜圖

對(duì)扭矩-轉(zhuǎn)速控制環(huán)發(fā)電機(jī)轉(zhuǎn)速反饋回路串聯(lián)加入塔架一階模態(tài)頻率陷波器，3 P頻率陷波器及6 P頻率陷波器。3個(gè)帶阻濾波器參數(shù)信息如表2所示。

表2 加入帶阻濾波器參數(shù)信息

在加入3個(gè)帶阻濾波器后，扭矩-轉(zhuǎn)速開(kāi)環(huán)伯德圖如圖6所示，在塔架一階模態(tài)，3 P及6 P頻率處，對(duì)數(shù)頻率特性幅值均有了明顯的陷波作用，同時(shí)保證了系統(tǒng)穩(wěn)定域量仍然滿足設(shè)計(jì)要求。

圖6 加入帶阻濾波器后扭矩-轉(zhuǎn)速伯德圖

在相同運(yùn)行風(fēng)速工況下，控制算法中扭矩-轉(zhuǎn)速控制環(huán)加入帶阻濾波器后，對(duì)風(fēng)機(jī)實(shí)際運(yùn)行轉(zhuǎn)速及槳距角頻譜進(jìn)行對(duì)比分析，發(fā)現(xiàn)塔架一階振動(dòng)模態(tài)頻率處幅值有了明顯衰減抑制效果，且在風(fēng)機(jī)實(shí)際運(yùn)行時(shí)發(fā)現(xiàn)大幅度減小了塔筒的橫向振動(dòng)加速度值，如圖7、圖8所示。

圖7 扭矩-轉(zhuǎn)速控制環(huán)加入帶阻濾波器后運(yùn)行轉(zhuǎn)速頻譜圖

圖8 扭矩-轉(zhuǎn)速控制環(huán)加入帶阻濾波器后運(yùn)行槳距角頻譜圖

4 結(jié)語(yǔ)

本文主要研究探討了通過(guò)在風(fēng)機(jī)控制算法轉(zhuǎn)速反饋控制中串聯(lián)使用帶阻濾波器，可以有效解決風(fēng)力發(fā)電機(jī)組經(jīng)常報(bào)出的一類(lèi)塔架振動(dòng)問(wèn)題的方法。先是從通用型帶阻濾波器伯德圖入手，系統(tǒng)分析了陷波濾波器的頻率特性，驗(yàn)證了通過(guò)調(diào)節(jié)濾波器中阻尼比系數(shù)的值來(lái)實(shí)現(xiàn)對(duì)特定頻率處幅值陷波深度的控制。然后結(jié)合Bladed軟件生成風(fēng)力機(jī)的坎布爾圖分析介紹了風(fēng)機(jī)葉輪旋轉(zhuǎn)頻率及其倍頻與機(jī)組固有頻率相互耦合的影響。最后以某風(fēng)場(chǎng)在運(yùn)行1.5 MW 77型風(fēng)電機(jī)組塔架振動(dòng)故障處理情況為例，通過(guò)對(duì)比算法扭矩-轉(zhuǎn)速控制環(huán)反饋通道中加入塔架一階模態(tài)頻率，3 P及6 P陷波濾波器前后，對(duì)風(fēng)機(jī)運(yùn)行轉(zhuǎn)速及槳距角頻譜進(jìn)行分析，詳細(xì)闡述了串聯(lián)帶阻濾波器后對(duì)塔架一階模態(tài)頻率處的陷波作用，并通過(guò)風(fēng)場(chǎng)實(shí)際運(yùn)行情況來(lái)看，取得了較好的振動(dòng)抑制效果，對(duì)于風(fēng)機(jī)塔架振動(dòng)問(wèn)題處理具有實(shí)際參考意義。