放電聲發(fā)射波檢測(cè)中數(shù)據(jù)預(yù)處理的小波分析實(shí)現(xiàn)

放電聲發(fā)射波檢測(cè)中數(shù)據(jù)預(yù)處理的小波分析實(shí)現(xiàn)

摘　要：研究了放電聲發(fā)射波檢測(cè)中用小波變換進(jìn)行數(shù)據(jù)預(yù)處理的實(shí)現(xiàn)問(wèn)題。內(nèi)容包括小波分析預(yù)處理的步驟，選擇小波，選擇小波分解層數(shù)，提取有用頻段的小波系數(shù)，小波消噪過(guò)程，小波消噪中的閾值選取，對(duì)提取的小波系數(shù)進(jìn)行小波消噪等。兩組實(shí)測(cè)信號(hào)的小波分析數(shù)據(jù)預(yù)處理結(jié)果證明了數(shù)據(jù)預(yù)處理的有效性。
　　關(guān)鍵詞：局部放電，聲發(fā)射檢測(cè)，小波變換，數(shù)據(jù)預(yù)處理
1　引言
　　對(duì)于放電聲發(fā)射在線監(jiān)測(cè)裝置來(lái)說(shuō)，如何排除干擾的影響是實(shí)現(xiàn)準(zhǔn)確測(cè)量放電聲發(fā)射信號(hào)的主要問(wèn)題。在裝置的實(shí)際應(yīng)用中，所采集的信號(hào)往往包含許多尖峰或突變部分，并且噪聲也不是平穩(wěn)的白噪聲，對(duì)這種信號(hào)進(jìn)行分析，首先需要作信號(hào)的預(yù)處理，提取有用信號(hào)。
　　小波變換以其優(yōu)秀的時(shí)頻局部特性被廣泛用來(lái)從干擾環(huán)境中提取有用信號(hào)。在局部放電信號(hào)提取中小波變換的效果顯著，本文重點(diǎn)研究在局部放電聲發(fā)射波檢測(cè)中用小波變換進(jìn)行數(shù)據(jù)預(yù)處理的實(shí)現(xiàn)問(wèn)題。
2 小波變換與數(shù)據(jù)預(yù)處理
　　小波分析對(duì)非平穩(wěn)信號(hào)消噪有著傅立葉分析不可比擬的優(yōu)點(diǎn)，運(yùn)用小波分析進(jìn)行一維信號(hào)消噪處理是小波分析的一個(gè)重要應(yīng)用之一。由于小波分析能同時(shí)在時(shí)、頻域內(nèi)對(duì)信號(hào)進(jìn)行分析，它能有效地區(qū)分信號(hào)中的突變部分和噪聲，可以很好地保存有用信號(hào)中的尖峰和突變部分，從而實(shí)現(xiàn)信號(hào)的消噪。
　　另外，我們能夠用小波分析將合成信號(hào)中的單純正弦信號(hào)的頻率提取出來(lái)。因?yàn)樵谛〔ǚ纸庀?，不同的尺度具有不同的時(shí)間和頻率分辨率，因而小波分解能將信號(hào)的不同頻率成份分開(kāi)。
　　在局部放電信號(hào)提取中，用小波分析進(jìn)行數(shù)據(jù)預(yù)處理的步驟如下：
　　1) 選定小波函數(shù)和小波分解層數(shù)；
　　2) 用小波分解函數(shù)對(duì)數(shù)據(jù)進(jìn)行一維離散小波（/包）分解；
　　3) 提取有用頻段的小波系數(shù)；
　　4) 對(duì)提取的小波系數(shù)進(jìn)行小波消噪處理；
　　5) 對(duì)消噪后的信號(hào)進(jìn)行小波重構(gòu)。
3 提取有用頻段的小波系數(shù)
　　在小波分解下，不同的尺度具有不同的時(shí)間和頻率分辨率，因而，小波分解能把信號(hào)的不同頻率成分分開(kāi)，因此，可以利用小波分解提取出特定頻率范圍內(nèi)的信號(hào)。
3.1 選擇小波
　　小波分析在工程應(yīng)用中，一個(gè)十分重要的問(wèn)題是最優(yōu)小波基的選擇。這是因?yàn)樾〔ǖ姆N類較多，不同的小波往往具有不同的時(shí)頻特性，用不同的小波基分析同一個(gè)問(wèn)題會(huì)產(chǎn)生不同的結(jié)果。因此從眾多的小波中選取合適的小波，是能否有效地利用小波變換提取干擾環(huán)境中局部放電聲發(fā)射信息的關(guān)鍵。
　　對(duì)于任意信號(hào)f(x)∈L2(IR),其連續(xù)小波變換定義為：
　　小波變換系數(shù)Wf(a,b)實(shí)際上是函數(shù)f(x)和小波 的相關(guān)系數(shù)，它反映的是相應(yīng)時(shí)段的信號(hào)和選定的小波之間的相似程度，Wf(a,b)的值越大，信號(hào)和選定的小波之間相似度就越高。這也正是對(duì)于同一信號(hào)選擇不同的小波進(jìn)行分解，其結(jié)果差異很大的原因。
由此，我們可以得出選擇小波的依據(jù)：被分析信號(hào)的時(shí)域和頻域特性與所用小波的時(shí)域和頻域特性分別相接近，二者的相似程度越高，分析的效果越佳。
分析比較之后發(fā)現(xiàn)：在眾多的小波簇中，bior3.7小波和局部放電聲發(fā)射波的特征最接近（另文專述）。其分解用小波函數(shù)和尺度函數(shù)分別具有很好的對(duì)稱性、緊支性和正則性，中心頻率在0.9336處。
3.2 選擇分解層數(shù)
　　總體上來(lái)說(shuō)，一個(gè)一維離散的信號(hào)，它的高頻部分影響的是小波分解的高頻第一層，低頻部分影響的是小波分解的最深層及其低頻層。如果對(duì)一個(gè)只是由白噪聲組成的信號(hào)進(jìn)行小波分解，則可以看出：高頻系數(shù)的幅值隨著分解層次的增加而很快地衰減，并且，高頻系數(shù)的方差也很快地衰減。
　　在小波變換中，對(duì)應(yīng)每一尺度a，都有一偽頻率和它對(duì)應(yīng)：
式（2）中，a是尺度；△是采樣周期；Fc是小波的中心頻率，單位為Hz；是尺度a對(duì)應(yīng)的偽頻率，單位為Hz。
　在小波多分辨率分解中，若將信號(hào)中的最高頻率成分看作是1，則各層小波分解便是帶通或低通濾波器，且各層所占的具體頻帶為：
　　從式（3）可見(jiàn)，當(dāng)小波函數(shù)和采樣周期選定之后，信號(hào)多分辨率分解的各層所占的頻帶是一定的。根據(jù)實(shí)測(cè)結(jié)果，變壓器局部放電聲發(fā)射波的頻譜一般在20－200kHz內(nèi)，將各層頻帶與實(shí)際信號(hào)的有效頻帶相比較，就可以確定小波分解層數(shù)。
3.3 提取有用頻段的小波系數(shù)
　　提取有用頻段的小波系數(shù)也可以看作是小波濾波，也就是把有用頻帶之外的小波分解系數(shù)強(qiáng)制置0。再進(jìn)行小波重構(gòu)就可得到信號(hào)的有用部分。
　　武漢高壓研究所局部放電實(shí)驗(yàn)室實(shí)測(cè)的兩組變壓器局部放電聲發(fā)射信號(hào)，數(shù)據(jù)長(zhǎng)度為512點(diǎn)，采樣率為1μs。這里用前述方法對(duì)實(shí)測(cè)信號(hào)進(jìn)行數(shù)據(jù)預(yù)處理。
　　取bior3.7小波，用多分辨率分析對(duì)數(shù)據(jù)進(jìn)行2層小波分解，a2層幾乎包括了信號(hào)的所有信息。實(shí)際上經(jīng)計(jì)算可知，a2層所占的頻段為0－233kHz，它完全覆蓋了放電聲發(fā)射波的全部頻段。也就是說(shuō)，放電聲發(fā)射波的有效數(shù)據(jù)都包含在小波分解的a2層內(nèi)，或者說(shuō)小波分解的低頻a2層包含了全部放電聲發(fā)射信號(hào)的有用信息。所以，a2層就是要提取的有用頻段，直接用a2層的系數(shù)進(jìn)行小波重構(gòu)，就可還原出放電聲發(fā)射信號(hào)。
4 小波消噪
4.1 小波消噪過(guò)程
　　通常按時(shí)域信號(hào)特征將干擾分為連續(xù)的周期型干擾、脈沖型干擾和白噪三大類。這里研究用小波分析從已消除周期型和脈沖型干擾的含噪局部放電信號(hào)中提取真實(shí)信息的方法。一個(gè)含白噪聲的一維信號(hào)的模型可以表示成如下的形式：
其中，f(i)為真實(shí)信號(hào)，e(i)為噪聲，s(i)為含噪聲的信號(hào)。
　　對(duì)式（4）表示的含噪信號(hào)的消噪過(guò)程可分為三個(gè)步驟進(jìn)行：
　　(1)一維信號(hào)的小波分解。選擇一個(gè)小波并確定一個(gè)小波分解的層次N，然后對(duì)信號(hào)f(i)進(jìn)行N層小波分解。
　　(2)小波分解高頻系數(shù)的閾值量化。對(duì)第1層到第N層的每一層高頻系數(shù)，選擇一個(gè)閾值進(jìn)行軟閾值量化處理。
　　(3)一維小波的重構(gòu)。根據(jù)小波分解的第N層的低頻系數(shù)和經(jīng)過(guò)量化處理后的第1層到第N層的高頻系數(shù)，進(jìn)行一維信號(hào)的小波重構(gòu)。
4.2 小波消噪中閾值的選取
　　(1) 默認(rèn)閾值消噪處理。該方法利用信號(hào)分析地基本理論，構(gòu)造一個(gè)函數(shù)，利用它可以產(chǎn)生信號(hào)的默認(rèn)閾值，再利用此閾值進(jìn)行小波消噪處理。
　　(2) 給定閾值消噪處理。在局部放電信號(hào)的實(shí)際消噪處理過(guò)程中，閾值可以通過(guò)經(jīng)驗(yàn)公式獲得，更直接地，可以在檢測(cè)裝置地安裝調(diào)試過(guò)程中，測(cè)出現(xiàn)場(chǎng)的噪聲大小，進(jìn)而求出閾值。這種方法考慮了現(xiàn)場(chǎng)的實(shí)際噪聲，具有較高的可信度。
4.3 放電聲發(fā)射信號(hào)的小波消噪
　　提取有用頻段的小波系數(shù)實(shí)質(zhì)上是將小波分解中高頻部分的系數(shù)強(qiáng)制置零，也就是說(shuō)，對(duì)信號(hào)的高頻區(qū)間進(jìn)行了抑制。但是被抑制后的信號(hào)中，仍然存在噪聲。這部分噪聲可以通過(guò)小波消噪來(lái)抑制或衰減。
　　對(duì)提取到的有用頻段小波系數(shù)再用bior3.7小波進(jìn)行兩層分解，然后，選用固定閾值形式進(jìn)行消噪處理：全部消除掉高頻部分（d1層），保留高頻d2層的有用信號(hào)。消噪效果顯著，消噪后的波形中保留有全部放電聲發(fā)射信息，但噪聲幾乎被全部消除。
5 結(jié)論
　　消噪后重構(gòu)表示的是用提取的有用頻段的小波系數(shù)經(jīng)消噪后再進(jìn)行小波重構(gòu)的波形。
　　從其與原始信號(hào)的對(duì)比中可以發(fā)現(xiàn)，用本文所述方法對(duì)局部放電聲發(fā)射信號(hào)進(jìn)行小波數(shù)據(jù)預(yù)