基于小波變換的扭振測(cè)量方法

張春海，宋秋艷

(國(guó)核電力規(guī)劃設(shè)計(jì)研究院，北京 100095)

1 概述

汽輪發(fā)電機(jī)組振動(dòng)故障診斷是根據(jù)相關(guān)的數(shù)據(jù)和信息對(duì)故障定性，進(jìn)而對(duì)其產(chǎn)生的原因或機(jī)理做出判斷，并確定解決措施和實(shí)施處理方案。軸類(lèi)部件有三種振動(dòng)形態(tài)：徑向振動(dòng)，縱向振動(dòng)和扭轉(zhuǎn)振動(dòng)。其中扭轉(zhuǎn)振動(dòng)雖然不像徑向振動(dòng)那樣引起大家普遍最多的關(guān)注，但是在70年代初，美國(guó)內(nèi)華達(dá)州南部Mohave電站的一臺(tái)483/426MW機(jī)組發(fā)生大軸斷裂事故，經(jīng)查實(shí)是由扭振疲勞引起的，斷裂的部位發(fā)生在發(fā)電機(jī)和勵(lì)磁機(jī)之間發(fā)電機(jī)的集流環(huán)處。自此以后扭振對(duì)汽輪發(fā)電機(jī)組安全性的影響引起了汽輪機(jī)行業(yè)的密切注意。

眾所周知的臺(tái)灣核電三廠GE公司在1985年發(fā)生的由葉片脫落引起火災(zāi)的事故，具體是900MW一號(hào)機(jī)組由于軸系的扭振固有頻率沒(méi)有避開(kāi)兩倍工頻，導(dǎo)致發(fā)電機(jī)轉(zhuǎn)子要承受兩倍工頻的扭矩，葉片諧振的應(yīng)力比設(shè)計(jì)值高出5倍，最終導(dǎo)致葉片脫落，引起火災(zāi)，直接損失70億臺(tái)幣，修復(fù)時(shí)間約1.5年。同樣是1985年，姚孟電廠300MW 機(jī)組在一次停機(jī)操作中斷路器A相跳不開(kāi)，引起較大負(fù)序電流，導(dǎo)致轉(zhuǎn)子損壞嚴(yán)重。這是由于斷路器的切換操作和電網(wǎng)故障對(duì)機(jī)組的沖擊會(huì)進(jìn)一步導(dǎo)致機(jī)電諧振，很大程度上取決于軸系固有的扭振特性。

同時(shí)，汽輪機(jī)調(diào)速系統(tǒng)的甩負(fù)荷或快關(guān)等擾動(dòng)，均可能導(dǎo)致軸系扭振。汽輪機(jī)扭振得危險(xiǎn)狀況，一般來(lái)說(shuō)有兩種：次同步共振和超同步共振。

2 扭振監(jiān)測(cè)的方法

目前，進(jìn)行扭振監(jiān)測(cè)的方法有扭應(yīng)力監(jiān)測(cè)，扭角測(cè)量，激光測(cè)量等。其中扭應(yīng)力監(jiān)測(cè)有國(guó)內(nèi)設(shè)備廠家研發(fā)使用，該設(shè)備使用傳感器獲得蒸汽壓力、發(fā)電機(jī)組電壓、電流等信號(hào)，再經(jīng)過(guò)掃描、篩選進(jìn)入程序計(jì)算疲勞壽命損失，統(tǒng)計(jì)每次事故達(dá)到的沖擊水平。這種方法結(jié)構(gòu)復(fù)雜，測(cè)點(diǎn)眾多，計(jì)算繁復(fù)，所以設(shè)備價(jià)格昂貴，調(diào)試復(fù)雜。本文介紹的是一種由電機(jī)的定子電流來(lái)間接測(cè)量軸系扭轉(zhuǎn)振動(dòng)的測(cè)量方法。本方法主要有以下幾個(gè)步驟：

(1)軸系建模

現(xiàn)在主要有3種用于軸系扭振分析的模型，其中適合汽輪發(fā)電機(jī)組的軸系扭振的模型可采用和Prohl法類(lèi)似的方式，將轉(zhuǎn)子的轉(zhuǎn)動(dòng)慣量集總到較多的沿軸線的節(jié)點(diǎn)上，可以得到分別只有轉(zhuǎn)動(dòng)慣量和彈性變形的集中質(zhì)量和彈性軸段組成的模型轉(zhuǎn)子。對(duì)由交流電機(jī)驅(qū)動(dòng)的機(jī)械系統(tǒng)進(jìn)行動(dòng)力學(xué)分析，整個(gè)轉(zhuǎn)軸系統(tǒng)可簡(jiǎn)化為圖1所示的模型，

系統(tǒng)振動(dòng)可用下式進(jìn)行描述：

由于J1較小，可忽略，得：

上兩式中c、k為扭振阻尼和剛度，Te為氣隙扭矩，TC為負(fù)載扭矩

上式可得，TC、θ12、ω12、Te有相同的頻率成分。只要可以確定其中任一參數(shù)的頻率成分，即可得其他參數(shù)的頻率成分。

根據(jù)電機(jī)原理，有：

其中：PP為電機(jī)極對(duì)數(shù)，ist為用于產(chǎn)生扭矩得定子電流分量。

又有

其中：ωS= PP(ω0- ωr)為轉(zhuǎn)差， σTrLS分別為電機(jī)常數(shù)，ism為定子電流產(chǎn)生磁場(chǎng)的分量。

當(dāng)系統(tǒng)發(fā)生扭振時(shí)，

其中：f0AT為扭振的頻率和振幅

由(3)(4)(5)，可得

上述兩分量事實(shí)上定義于旋轉(zhuǎn)坐標(biāo)中，是不可測(cè)得的，可通過(guò)將旋轉(zhuǎn)坐標(biāo)向定子線圈坐標(biāo)轉(zhuǎn)換，獲得某一相定子線圈電流為

其中：f 為定子電流頻率。

由上式可以看出，當(dāng)機(jī)械系統(tǒng)發(fā)生扭振，定子電流頻率和扭振頻率差頻和和頻的頻率分量疊加在定子電流中，而且可以看出，該頻率分量的幅值與扭振振幅成正比，所以可得，我們只要分析定子電流的頻譜，得出差頻和和頻的頻率分量，就可確定機(jī)械系統(tǒng)的扭振。所以，系統(tǒng)定子電流于扭振間的傳遞的數(shù)學(xué)模型可由上七式描述。

(2)電流測(cè)量

電機(jī)電流通過(guò)測(cè)取電機(jī)導(dǎo)線電阻的電壓降的方法獲得。本方法的驗(yàn)證階段使用實(shí)驗(yàn)臺(tái)的實(shí)驗(yàn)數(shù)據(jù)，所用的實(shí)驗(yàn)臺(tái)由電機(jī)、齒輪箱和液壓馬達(dá)組成，這種實(shí)驗(yàn)臺(tái)對(duì)于實(shí)際汽輪發(fā)電機(jī)組的應(yīng)用具有非常重要的指導(dǎo)意義。

(3)扭振信號(hào)處理技術(shù)

目前，使用的信號(hào)處理方法不僅有對(duì)穩(wěn)態(tài)信號(hào)處理的各種譜分析方法，而且使用處理瞬態(tài)信號(hào)的小波變換等新的方法。由于機(jī)械調(diào)制的原應(yīng)，還有使用希爾波特(Hilbert)等處理調(diào)制信號(hào)的方法。小波變換是一種信號(hào)的時(shí)間-尺度分析方法，窗口大小可固定不變，時(shí)間窗和頻率窗都可以改變，是一種時(shí)頻局部化分析方法，很適合監(jiān)測(cè)正常信號(hào)中夾帶的瞬態(tài)反?，F(xiàn)象。故本文推薦使用小波變換進(jìn)行主軸瞬態(tài)扭振的估計(jì)。根據(jù)瞬態(tài)定子電流測(cè)量信號(hào)，運(yùn)用小波變換算法對(duì)主軸瞬態(tài)扭振過(guò)程進(jìn)行估計(jì)。以下為小波算法：

設(shè)x(t)是平方可積函數(shù)[記作x(t)∈L2(R)]，φ(t)是被稱(chēng)為基本小波或母小波(mother wavelet)的函數(shù)。則

稱(chēng)為x(t)的小波變換。式中α>0是尺度因子，τ反映位移，其值可正可負(fù)。符號(hào)代表內(nèi)積，它的含義是(上標(biāo)＊代表取共軛)

φaτ(t)前加因子的目的是使不同a值下φaτ(t)的能量保持相等。即，設(shè)ε = ∫|φ(t)|2dt是基本小波的能量，則φaτ(t)的能量是

對(duì)定子電流信號(hào)進(jìn)行連續(xù)小波變換和分析， 通過(guò)小波濾波器在一系列選擇性相同的頻帶上，使之在具有良好濾波器特性的小波時(shí)頻窗中得以體現(xiàn)，以獲得信號(hào)所包含的特征信息。

3 結(jié)論

使用定子電流間接測(cè)量軸的扭振，估計(jì)的扭振的頻譜具有一定的精度和可靠性，對(duì)回轉(zhuǎn)機(jī)械的非穩(wěn)態(tài)扭振，通過(guò)定子電流的小波變換等分析方法可獲得。實(shí)現(xiàn)機(jī)械故障診斷和狀態(tài)監(jiān)測(cè)是扭振測(cè)量的目的之一，通過(guò)定子電流信號(hào)提取機(jī)械系統(tǒng)的狀態(tài)特征和故障特征，最終實(shí)現(xiàn)機(jī)械系統(tǒng)的故障診斷，是本研究的未來(lái)的發(fā)展方向之一，對(duì)于汽輪發(fā)電機(jī)組扭振的故障檢測(cè)具有指導(dǎo)意義。

[1]顧煜炯，何成兵，楊昆.汽輪發(fā)電機(jī)組軸系扭振監(jiān)測(cè)與分析系統(tǒng)[J].儀器儀表學(xué)報(bào)，2001，22增刊.

[2]張俊紅，高文志，付魯華.大型汽輪發(fā)電機(jī)組軸系振動(dòng)及其控制的研究[J].動(dòng)力工程，2003，23(4).

[3]趙華華，童小忠，陸頌元.汽輪發(fā)電機(jī)組軸系扭振固有頻率現(xiàn)場(chǎng)測(cè)定與分析[J].浙江電力，2002(2).

[4]傅忠廣，楊昆，宋之平，等.汽輪發(fā)電機(jī)組軸系扭振的激光測(cè)試技術(shù)[J].現(xiàn)代電力，2001，18(1).

[5]Peng Yong， Zhou Yichen， Hilbert transformation and its application to demodulation， Journal of Wuhan Transportation University， v20，n5，Oct，1996(5).