基于相關(guān)分析法的聲發(fā)射信號降噪技術(shù)

摘 要：闡述了相關(guān)分析法在剔除聲發(fā)射信號中環(huán)境噪聲的基本原理。通過LF3鋁合金的腐蝕實驗，將未處理的聲發(fā)射信號與環(huán)境中的噪聲信號進行相關(guān)分析，并以它們相關(guān)系數(shù)的大小作為判斷該信號是否為環(huán)境噪聲信號的依據(jù)。結(jié)果表明該方法能準確地剔除信號中的環(huán)境噪聲信號。
關(guān)鍵詞：聲發(fā)射；腐蝕；相關(guān)分析；降噪；環(huán)境噪聲
0 引言
聲發(fā)射(Acoustic Emission，簡稱AE)是指材料或結(jié)構(gòu)受內(nèi)力或外力作用產(chǎn)生變形或斷裂，以彈性波的形式釋放出應變能的一種物理現(xiàn)象。聲發(fā)射技術(shù)是通過分析聲發(fā)射信號進而提取聲發(fā)射信號特征信息來推斷聲發(fā)射源特性的無損檢測技術(shù)。當前限制這一技術(shù)進一步發(fā)展的一個關(guān)鍵問題是采集到的信號中含有大量的噪聲信號，以致可能淹沒真實的聲發(fā)射信號。同時由于聲發(fā)射信號在傳播過程中受傳播介質(zhì)和傳感器的影響，使得采集到的信號與真實的聲發(fā)射信號有很大的差別。如何準確地從采集到的信號中剔除噪聲信號，提取出真實的聲發(fā)射信號是聲發(fā)射技術(shù)發(fā)展的一個重要問題。
1 相關(guān)分析法的降噪原理
在聲發(fā)射噪聲信號中大部分是環(huán)境中的噪聲信號，它們隨環(huán)境的變化而變化且其特征信息十分豐富。在傳統(tǒng)的聲發(fā)射信號預處理中常采用參數(shù)分析法，尤其在不確定聲發(fā)射信號特征參數(shù)的情況下，噪聲信號的復雜性使得聲發(fā)射的預處理變得十分困難。由于波形分析法比參數(shù)分析法具有更多的信息，本文利用波形的相關(guān)分析法對信號進行預處理。
基于這樣一種設想，用一傳感器檢測金屬腐蝕聲發(fā)射信號，同時另一傳感器置于同樣的環(huán)境中用來檢測環(huán)境噪聲而不檢測金屬腐蝕時的聲發(fā)射信號。由于金屬腐蝕聲發(fā)射信號通常較弱，在工程中可以采取一定的措施使上述設想得以實現(xiàn)。而后將采集到的金屬腐蝕聲發(fā)射信號與采集到的環(huán)境噪聲信號進行實時相關(guān)分析。通過相關(guān)分析，根據(jù)它們的相關(guān)系數(shù)作為判斷信號是否為環(huán)境噪聲信號的依據(jù)。在實際中常用相關(guān)函數(shù)作為兩個信號相關(guān)性的檢測，相關(guān)函數(shù)的定義如下：
信號x(n)和y(n)的互相關(guān)函數(shù)Rxy(m)的定義為：
（1）
的值隨信號x(n)和y(n)相關(guān)性的增強而增大，當信號x(n)和y(n)為同一信號時， 的值達到最大值。實際工作中，信號x(n)、y(n)總是有限長，并且是數(shù)字信號，對于這樣的信號，互相關(guān)公式為：
（2）
自相關(guān)公式即是將公式（1）和（2）中的y(n)替換為x(n)。在實際運用時，為了讓兩個信號的相關(guān)性更為直觀，更有可比性，又引入了相關(guān)系數(shù)的概念。其定義為：
（3）
由許瓦茲(Schwartz)不等式有 ，由公式 (3)可知，當x(n)=y(n)時， =1，兩信號完全相關(guān)(相等)，這時 取得最大值；當x(n)和y(n)完全無關(guān)時， =0， =0；當x(n)和y(n)有某種程度相似時， ，| |在0和1之間取值。因此 和 可用來描述x(n)和y(n)間的相似程度。
2 實驗內(nèi)容
2.1 實驗材料及儀器
本實驗選用硝酸溶液作為腐蝕介質(zhì)，選用LF3鋁合金作為受腐蝕材料。本實驗所用聲發(fā)射儀是德國VALLEN公司生產(chǎn)的16通道AMSY-5型聲發(fā)射系統(tǒng)。該系統(tǒng)能夠同時采集AE信號特征參數(shù)和波形(Transient Recording-TR)數(shù)據(jù)，可檢測到的頻率范圍從5kHz到3MHz。所用傳感器為VS150-RIC型傳感器，其頻率范圍為80～550kHz， 中心頻率為150kHz。實驗時將AMSY-5型聲發(fā)射儀的各通道的門檻值設為28.8dB。
2.2 實驗方案
圖1為整個實驗系統(tǒng)示意圖。硝酸溶液置于玻璃鏡筒內(nèi)，與試樣表面接觸發(fā)生腐蝕。鏡筒尺寸(mm)為φ40×100。鏡筒與試樣之間用凡士林密封，鏡筒內(nèi)加入硝酸溶液后，上面用玻璃片加凡士林密封，防止硝酸揮發(fā)。傳感器耦合劑為凡士林，用彈性帶固定并對傳感器施加一定壓力，保證耦合良好。
試樣表面在被鏡筒限定的區(qū)域內(nèi)發(fā)生腐蝕，產(chǎn)生應力波被傳感器接收并通過電纜傳送至聲發(fā)射儀進行信號的濾波、特征提取、波形記錄等處理，并將處理后的數(shù)據(jù)通過數(shù)據(jù)線送入PC機存儲，再進行進一步的處理、分析和結(jié)果顯示。
安放傳感器之前將試樣表面安放傳感器的部位打磨光滑，鏡筒外部的試樣表面均勻涂上一薄層凡士林，以防止試樣裸露部分在空氣中發(fā)生腐蝕，并降低彈性帶與試樣的摩擦噪聲。為了避免各試樣的腐蝕聲發(fā)射信號被其他傳感器所接受和降低機械噪聲的影響，將四個腐蝕試樣置于鋪有毛氈的木質(zhì)工作臺上并保持約15cm距離。實驗時在玻璃鏡筒1、2、3內(nèi)分別均加入1.4%的硝酸溶液，試樣4上的傳感器用來采集外界環(huán)境中的噪聲信號。
3 實驗數(shù)據(jù)處理及結(jié)果
表1為這段時間各通道采集到的信號數(shù)。
由表1可以看出通道一、二、三采集到信號數(shù)大致一樣，而用來采集環(huán)境噪聲的通道四所記錄的信號數(shù)明顯小于其他三個傳感器的信號數(shù)。因而可以認為通道四只采集到了環(huán)境中的噪聲信號且各試樣的聲發(fā)射信號均只被該試樣上的傳感器所接受。由于四個傳感器處于同一環(huán)境，因而它們對環(huán)境中的噪聲同等敏感，試樣四上的傳感器采集到的環(huán)境中噪聲信號同樣被其他三個傳感器采集。同一噪聲源的噪聲被四個傳感器采集時，它們的波形應大致一致。圖2中a、b、c、d分別為傳感器一、二、三、四幾乎在同一時間采集到的噪聲信號，其開始采集的時間差最長為0.0522 ，最短為0.0009 。從它們的波形圖可以看出，它們的波形十分相似。因而可以利用通道1、2、3采集到的信號與通道4采集到的噪聲信號進行互相關(guān)分析，根據(jù)它們的相關(guān)系數(shù)作為判別通道1、2、3采集到的信號是否是環(huán)境中的噪聲信號的依據(jù)。圖3中e、f、g、h是從采集到的信號中隨機挑選出的四個信號，用于與圖2中的波形進行比較。
表1 各通道采集的信號數(shù)
通道一 通道二 通道三 通道四 總計
信號數(shù) 155 120 165 24 464

a b c d
圖2 同一噪聲源的波形

e f g h
圖3 不同聲源的波形
從圖2可以看出盡管四個波形十分的相近，但是各個波形采集的時間不一樣。這是因為噪聲在傳播途徑中由于噪聲源到達四個傳感器的路徑不同而引起能量的衰減不同，也就是波形的同一點到達傳感器的幅值不一樣，而該聲發(fā)射儀各通道的門檻值均為28.8dB， 從而引起了波形在時間上的平移。因此在進行信號的互相關(guān)分析時只選取其最相似的波形段進行分析。
表2為圖2和圖3中的7個波形中截取的部分波形與中與圖2中的圖d截取的部分波形相關(guān)系數(shù)的最大值。其截取波形數(shù)據(jù)最長為1448個數(shù)據(jù)，最短為1248個數(shù)據(jù)。波形原長為2048個數(shù)據(jù)。
圖4為圖a、b、c、e、f、g、h的波形與圖d中的波形的相關(guān)系數(shù)隨其截取波形數(shù)據(jù)的長度變化而變化的情況。由圖4可以看出來自于同一噪聲源的信號進行自相關(guān)分析時，在選取合適的數(shù)據(jù)長度時它們的互相關(guān)系數(shù)保持在0.7左右，而來自不同噪聲源的信號的互相關(guān)系數(shù)維持在0.25左右，二者有明顯的分界線。因此利用這種辦法可以較客觀地在采集到的信號中剔除來自外界環(huán)境的噪聲信號。利用同樣的辦法，將傳感器四采集的剩余的23個環(huán)境噪聲信號與其它三個通道采集的信號進行相關(guān)分析。同一噪聲源的信號被四個通道幾乎同時采集，為減少計算量我們只對與通道四采集的24個信號的采集時間相差不超過0.5 的信號進行相關(guān)分析。通過計算，通道一、二、三采集的信號中均有24個信號與通道四采集到的環(huán)境噪聲信號的相關(guān)系數(shù)在0.6以上。其結(jié)果見表3。
表2 圖d與其它圖的波形相關(guān)系數(shù)
相關(guān)系數(shù)
圖a 圖b 圖c 圖e 圖f 圖g 圖h
圖d 0.7682 0.6428 0.6885 0.2432 0.2002 0.3028 0.3154
表3 實驗結(jié)果
信號數(shù)（個） 0.6-0.7 0.7-0.8 0.8-0.9 0.9-0.95
通道一 4 10 8 2
通道二 4 9 10 1
通道三 3 12 8 1
4結(jié)