劉 偉，朱玉婷，付平勇

(華東交通大學(xué)道路與鐵道工程重點(diǎn)實驗室，江西南昌330013)

鐵路邊坡的深部水平側(cè)移曲線時常鋸齒狀波動，一方面反映降雨、地震等自然環(huán)境因素或爆破、鉆孔等人為施工因素作用使邊坡內(nèi)部位移產(chǎn)生突變，另一方面，反映監(jiān)測信號由于電壓不穩(wěn)、電磁干擾等不良因素產(chǎn)生了噪聲。為了提高鐵路邊坡穩(wěn)定性監(jiān)測與評價的準(zhǔn)確性，有必要對采集的監(jiān)測信號進(jìn)行降噪處理[1，2]。小波分析作為一種有力的數(shù)學(xué)工具，具有多分辨分析的特點(diǎn)，在時頻和頻域都具有表征信號局部特征的能力，時間窗和頻率窗都可以根據(jù)信號的具體形態(tài)動態(tài)調(diào)整，基于這些特征，小波分析可以廣泛應(yīng)用于監(jiān)測信號降噪[3、4]。

1 小波降噪基本原理

利用小波分析降噪[5]－[7]，即在不同尺度下作小波變換，其實質(zhì)是用不同中心頻率的帶通濾波器對信號進(jìn)行濾波，以取得較好質(zhì)量的信號。在實際的工程中，有用信號通常表現(xiàn)為比較平穩(wěn)的信號或低頻信號，而噪聲信號通常表現(xiàn)為高頻信號。

設(shè)f(x)是觀測點(diǎn)觀測序列，采用Mallat算法，將信號分解成不同頻率成分:

其中，信號在低頻空間Vj上的投影為:

在高頻空間Wj上的投影為:

Ajf(x)是信號f(x)的頻率不超過2－j的成分，Djf(x)是頻率介于 2－j與 2－j+1之間的成分。

上述小波分解式實際可以寫成如下矩陣形式:

式中:j=1，2，…J;H=(hk－2n)，H 為尺度函數(shù)對應(yīng)的低通濾波器;G=(gk－2n)，G為小波函數(shù)對應(yīng)的帶通濾波器;Cj為在2－j分辨率下的離散逼近;Dj為在2－j分辨率下的離散細(xì)節(jié);{hk}k∈Z與{gk}k∈Z是一對離散正交鏡像濾波器(低通和高通)，(4)式即為Mallat塔式分解算法。

Mallat重構(gòu)算法:

其中，j=J，J－1，…，1，H*和 G*，分別為 H 和 G 的共軛。

小波分解與重構(gòu)降噪法就是根據(jù)Mallat分解式(1)至式(3)，將含噪聲的信號分解在不同的頻帶內(nèi)，然后，在信號重構(gòu)時，將與噪聲相應(yīng)的高頻細(xì)節(jié)信號有關(guān)部分Dj置零，重構(gòu)信號便可達(dá)到降噪的目的。是f(x)經(jīng)過濾波后的平滑信號表達(dá)式，即為提取的

?真實信號，這便是小波降噪的基本原理。因此信號降噪可分三步處理:首先對原始信號進(jìn)行小波分解，則噪聲通常含在高頻系數(shù)中;然后對小波分解的高頻系數(shù)以閾值等形式進(jìn)行量化處理;最后再對信號重構(gòu)即可達(dá)到降噪的目的。小波降噪的方法很多，如Mallat算法、小波變換模極大值法、非線性小波變換閾值法和平移不變量小波降噪法等，而小波變換模極大值法降噪和非線性小波變換閾值法降噪最為常用。

1.1 小波變換模極大值法

由于信號和噪聲的小波變換系數(shù)在不同尺度上具有不同的傳播特性，即隨著尺度的增大，噪聲所對應(yīng)的模極大值迅速衰減。而信號的模極大值有三情形:若為緩變信號，則模極大值逐漸增大;若為階越信號，則模極大值保持不變;若為脈沖信號，正、負(fù)極值組成的脈沖對的幅值將同時變小。于是，連續(xù)進(jìn)行若干次小波分解后，綜合各個尺度上模極大值的位置和幅值，可斷定哪些模極大值是由噪聲引起，哪些是由信號產(chǎn)生。去除噪聲引起的模極大值，再重構(gòu)剩余的模極大值，便可實現(xiàn)降噪。

1.2 非線性小波變換閾值法

非線性小波變換閾值法也稱小波收縮。當(dāng)信號的Lipschitz指數(shù)為正時，其小波變換極大值的幅度將隨著尺度而增加;但噪聲的Lipschitz指數(shù)為負(fù)時，則噪聲所對應(yīng)的小波變換模極大值將隨尺度的增加而減少。利用這個特性，在不同的分解尺度上設(shè)置一定的閾值，將小于該閾值的模極大值點(diǎn)認(rèn)為是噪聲的小波變換引起的，將其置為零;大于該閾值的模極大值點(diǎn)則認(rèn)為是信號小波變換引起的而予以保留。常見的閾值規(guī)則有sqtwolog規(guī)則、rigrsure規(guī)則、heursure規(guī)則和minimaxi規(guī)則。工程中常用的簡單閾值處理法有非線性小波變換軟閾值處理法、非線性小波變換硬閾值處理法等。

1.3 小波降噪的實現(xiàn)與性能評價

(1)為了選擇較優(yōu)的小波降噪方案，將涉及到小波函數(shù)、閾值的選擇。實際工程中可根據(jù)自相似原則選取小波基函數(shù)，若選擇的小波基函數(shù)與待處理信號有相似性，小波變換后，利于信號能量集中和減少計算量。土木工程[7]中關(guān)于測斜監(jiān)測曲線，用Symlet、Coflet、Daubechies小波系均可實現(xiàn)曲線的降噪平滑，但以Daubechies小波系效果最優(yōu)。對于閾值的選擇，其實現(xiàn)十分簡單，一個簡單的方法就是用每個小波分解系數(shù)的幅值減去閾值，若差為負(fù)值則將系數(shù)設(shè)置為0;若差是正值則保持系數(shù)不變。在小波降噪方案實施過程中，同時也可以采用Matlab中的小波工具實現(xiàn)。Matlab中獲取閾值的函數(shù)有:ddencmp、thselect、wbmpen、wdcbm;實現(xiàn)信號閾值去噪的函數(shù)有:wden、wdencmp、wthresh、wthcoef、wpthcoef、wpdencmp。

(2)為了對各種降噪方法進(jìn)行比選，可用小波降噪的性能評價標(biāo)準(zhǔn)來評判，工程中常用的小波降噪的性能評價標(biāo)準(zhǔn)有信噪比(SNR)、信噪比增益(A)和均方根誤差(RSME)，這三種評價標(biāo)準(zhǔn)是等價的[8]。其中信噪比是測量信號中噪聲量度的傳統(tǒng)方法，常被用來作為評價降噪效果的指標(biāo)，信噪比越高則濾波效果越好，則降噪性能越好。信噪比單位是分貝，其定義為:

2 工程實例

以峰福鐵路K294+680～K294+900段邊坡CX－6測斜孔，從2010年12月至2011年3月(測量間隔為1周，共間隔16周)的水平位移監(jiān)測數(shù)據(jù)作為小波降噪處理的對象，監(jiān)測數(shù)據(jù)均通過某高科技實業(yè)有限公司CX－806B型數(shù)字顯示滑動式鉆孔測斜儀采集獲得，并使用小波變換模極大值法、非線性小波變換軟閾值法和非線性小波變換硬閾值法對所采集數(shù)據(jù)進(jìn)行降噪處理，小波函數(shù)選擇Db3，尺度選取3，閾值選取常用的hersure閾值，并用信噪比對降噪效果作對比，提取鐵路邊坡體內(nèi)部水平變形的信息，進(jìn)行鐵路邊坡測斜監(jiān)測分析。

(1)邊坡測斜布置

該測區(qū)屬剝蝕低丘，邊坡為土質(zhì)邊坡，山體自然坡度20°～30°，植被發(fā)育，線路沿水口庫區(qū)北岸行走，整個病害長度約為300 m，包括一段路塹及其兩端各一段路堤。根據(jù)路段地形地貌及地質(zhì)情況，布置2排共6個測斜孔，其中一排布置在左側(cè)路塹坡頂以上的山坡上，且位于已有縱向裂縫的下方，測孔間距50 m，布置有 CX－2、CX－6、CX－13 與 CX－7四個測孔;另一排布置在路塹右側(cè)水溝外側(cè)，布置有CX－1和CX－9兩個測孔，布置見圖1所示。

圖1 測斜孔位布置示意圖

(2)數(shù)據(jù)采集與分析

待測斜管安裝完，使用測斜儀進(jìn)行預(yù)通并采集初讀數(shù)，再按時進(jìn)行其他測斜監(jiān)測。下面列出峰福鐵路K294邊坡CX－6測斜孔采集的部分水平位移監(jiān)測數(shù)據(jù)，并用小波變換模極大值法、非線性小波變換軟閾值法與非線性小波變換硬閾值法三種方法對所采集數(shù)據(jù)進(jìn)行降噪處理，降噪結(jié)果見表1。

表1 CX－6測斜監(jiān)測降噪處理結(jié)果

續(xù)表1

CX－6孔在第3、4、8、10與12周采集的原始數(shù)據(jù)有明顯的波動，其中第4、10與12周表現(xiàn)較劇烈。通過降噪前后的數(shù)據(jù)結(jié)果對比，小波變換模極大值法、非線性小波變換軟閾值法與非線性小波變換硬閾值法三種方法降噪效果都較明顯，降噪后的數(shù)據(jù)穿行在原始數(shù)據(jù)之間。第4與12周時，小波變換模極大值法與非線性小波變換軟閾值法降噪后的曲線，較非線性小波變換硬閾值法更好地保留了突變波動。

(3)邊坡測斜降噪評價

由表1可見，非線性小波變換閾值法降噪的信噪比大于小波變換模極大值法;而在非線性小波變換閾值法中，又以非線性小波變換軟閾值法的信噪比最高。由此可見，由于非線性小波變換軟閾值法和小波變換模極大值法都建立在奇異性分析基礎(chǔ)上，故降噪后的曲線對原始曲線中的尖峰和突變進(jìn)行了有效地保持。但從曲線的平滑程度上看，小波變換模極大值法的平滑性優(yōu)于非線性小波變換閾值法。同理，可對K294處邊坡其他測斜孔的監(jiān)測數(shù)據(jù)進(jìn)行降噪。通過對其他各測孔采集得較多數(shù)據(jù)降噪處理效果對比驗證，基于hersure閾值的非線性小波變換軟閾值法更合理有效，能夠敏感識別監(jiān)測數(shù)據(jù)噪聲和有用水平位移信息，更適合鐵路邊坡測斜監(jiān)測分析。

3 結(jié)論

(1)小波變換模極大值法、非線性小波變換軟閾值法與非線性小波變換硬閾值法三種方法降噪效果都較明顯，利于測斜監(jiān)測數(shù)據(jù)反映鐵路邊坡的位移變化趨勢。

(2)基于hersure閾值的非線性小波變換軟閾值法，降噪后的曲線有更高的信噪比，有較好的光滑性，且有效地保留了測斜水平位移監(jiān)測值的原有真實特性，故可以直觀地獲取產(chǎn)生水平位移突變所在深度，有利于鐵路邊坡變形的判定和穩(wěn)定性的評價，更適合鐵路邊坡測斜監(jiān)測分析。

(3)利用小波降噪理論，結(jié)合Matlab處理工具箱，在工程中處理深層水平位移測斜監(jiān)測可節(jié)約測量人員內(nèi)業(yè)處理分析時間，也使降噪后的數(shù)據(jù)更具科學(xué)性與說服性，值得在相關(guān)領(lǐng)域推廣。

[1]閆振，陸陽，田斌鵬.小波分析在邊坡深部水平位移監(jiān)測數(shù)據(jù)降噪中的應(yīng)用[J].中南公路工程，2007，32(2):41－44

[2]黃飛瀾，肖紅.測斜儀在高填方地基側(cè)向水平位移監(jiān)測中的應(yīng)用[J].公路工程，2010，35(5):112－120

[3]馮毅，王香華.小波變換降噪處理及其Matlab實現(xiàn)[J].數(shù)據(jù)采集與處理，2006，21(12):37－39

[4]張鵬，李獻(xiàn)勇，陳建平.基于小波降噪的隧道圍巖監(jiān)測數(shù)據(jù)分析[J].吉林大學(xué)學(xué)報(地球科學(xué)版)，2008，38(6):1010－1014

[5]張鵬.基于小波理論的白鶴隧道圍巖穩(wěn)定性非線性研究[D].長春:吉林大學(xué)，2009

[6]田勝利，周擁軍，葛修潤，等.基于小波分解的建筑物變形監(jiān)測數(shù)據(jù)處理[J].巖石力學(xué)與工程學(xué)報，2004，23(15):2639－2642

[7]Mallt S，Hwang W L.Singularity Detection and Processing with Wavelets[J].IEEE Trans Information Theory，1992，38(2):617－643

[8]Donoho DL.De－ noising by Soft－ thresholding[J].IEEE Transactions on Information Theroy，1995，41(3):613－627