周云耀 王海斌 向 涯,3,4 呂永清,3,4 齊軍偉,3,4 吳 歡,3,4 何 彬

1 中國地震局地震研究所，武漢市洪山側(cè)路40號，430071 2 引力與固體潮國家野外觀測研究站，武漢市洪山側(cè)路40號，430071 3 中國地震局地震大地測量重點(diǎn)實(shí)驗(yàn)室，武漢市洪山側(cè)路40號，430071 4 湖北省地震局，武漢市洪山側(cè)路40號，430071

隨機(jī)噪聲具有強(qiáng)烈的隨機(jī)性，在地震資料中普遍存在[1]，因此無論是天然地震還是人工地震，隨機(jī)噪聲的處理都是提高地震資料信噪比的關(guān)鍵之一。常用的地震信號隨機(jī)噪聲壓制方法有拉東變換、經(jīng)驗(yàn)?zāi)B(tài)分解、曲波變換與小波變換等[2]，其中小波變換具有多分辨率分析的特點(diǎn)，能夠同時表征信號在時域和頻域的局部特性[3]。經(jīng)小波分解得到小波系數(shù)后，根據(jù)有效信號與隨機(jī)噪聲小波系數(shù)的不同而設(shè)定閾值，對隨機(jī)噪聲的小波系數(shù)進(jìn)行處理后再對小波系數(shù)進(jìn)行重構(gòu)，可達(dá)到信號去噪的目的。

提升小波算法被稱為“第二代小波”，是一種基于時域運(yùn)算的信號分析方法，可克服經(jīng)典小波降噪方法局部信息丟失的缺點(diǎn)[4]。對于小波閾值去噪，閾值函數(shù)與閾值的選取對地震資料去噪的結(jié)果具有很大影響。經(jīng)典的閾值函數(shù)包含硬閾值和軟閾值，其中硬閾值函數(shù)是當(dāng)小波系數(shù)的絕對值大于設(shè)定的閾值時系數(shù)保持不變、小于閾值時系數(shù)為零，這會導(dǎo)致其在閾值處不連續(xù)，從而使重構(gòu)信號產(chǎn)生一定的震蕩；而軟閾值函數(shù)可克服硬閾值函數(shù)在閾值處不連續(xù)的缺點(diǎn)，當(dāng)小波系數(shù)的絕對值大于閾值時減去閾值，小于閾值時則進(jìn)行置零處理，但經(jīng)過軟閾值處理的小波系數(shù)與實(shí)際小波系數(shù)存在恒定偏差，從而影響重構(gòu)信號的精度。

本文在經(jīng)典軟、硬閾值函數(shù)的基礎(chǔ)上提出一種改進(jìn)的閾值函數(shù)，從理論上證明該閾值函數(shù)可克服經(jīng)典閾值函數(shù)的不連續(xù)與存在恒定偏差的缺點(diǎn)，并通過信號仿真證實(shí)經(jīng)過該閾值函數(shù)處理后可很好地壓制高斯白噪聲，提高含噪信號的信噪比。最后將改進(jìn)閾值函數(shù)的提升小波算法運(yùn)用于氣槍震源，分析其在人工震源去噪方面的可行性。

1 改進(jìn)閾值函數(shù)的提升小波變換

1.1 提升小波

Mallat算法利用多分辨率分析分解信號，從而去除信號間的相關(guān)性，而提升小波算法則是將信號分解為偶數(shù)采樣點(diǎn)和奇數(shù)采樣點(diǎn)[5]。由于相鄰數(shù)據(jù)的相關(guān)性較強(qiáng)，可通過奇數(shù)采樣點(diǎn)估計偶數(shù)采樣點(diǎn)，再去除包含在偶數(shù)采樣點(diǎn)中關(guān)于奇數(shù)采樣點(diǎn)的信息，從而消除信號間的相關(guān)性。提升小波算法的實(shí)現(xiàn)包含分裂(spilit)、預(yù)測(predict)和更新(update)。其中，分裂是將原始信號Si按照序列的奇偶性分裂為2個互不相關(guān)的子集，即偶數(shù)序列Se與奇數(shù)序列So。預(yù)測是利用原始信號相鄰數(shù)據(jù)間的相關(guān)性，用偶數(shù)序列Se與1個預(yù)測算子P來預(yù)測奇數(shù)序列So，得到P(So)，并將奇數(shù)序列So與奇數(shù)序列預(yù)測值P(So)的差值di-1作為小波系數(shù)，其中di-1=So-P(So)。更新則是為使分裂后得到的偶數(shù)序列Se保持原始信號的某些特性，構(gòu)造一個更新算子U(di-1)對其進(jìn)行更新，得到尺度系數(shù)ai-1[6]：

ai-1=Se+U(di-1)

(1)

1.2 閾值去噪

對含噪地震信號進(jìn)行一維離散小波分解，得到一系列小波系數(shù)。由于有效信號在時間域或空間域具有一定的連續(xù)性，而噪聲信號則表現(xiàn)為隨機(jī)性，從而導(dǎo)致小波域中有效信號對應(yīng)的小波系數(shù)值較大，而噪聲信號對應(yīng)的小波系數(shù)值較小。根據(jù)該特點(diǎn)設(shè)置一個閾值，將小于該閾值的小波系數(shù)去除，即將噪聲主導(dǎo)對應(yīng)的小波系數(shù)去除，再對剩余的小波系數(shù)進(jìn)行重構(gòu)，達(dá)到去除噪聲的目的，因此閾值與閾值函數(shù)的選擇具有重要意義：

(2)

(3)

研究表明[3,6-7]，隨著分解層數(shù)的增加，各層小波細(xì)節(jié)系數(shù)中噪聲含量會逐漸減小，因此閾值應(yīng)當(dāng)有所變化。經(jīng)過一系列仿真實(shí)驗(yàn)，本文使用一種改進(jìn)的閾值，具體公式為：

(4)

式(2)～式(4)中，λ為去噪閾值，σ為噪聲強(qiáng)度，N為信號長度，cdn為小波分解后第n層高頻小波系數(shù)，j為分解層數(shù)。

在得到閾值后，可使用硬閾值函數(shù)和軟閾值函數(shù)對小波系數(shù)進(jìn)行處理。硬閾值函數(shù)表達(dá)式為：

(5)

圖1為硬閾值函數(shù)圖像，從圖中可以看出，硬閾值函數(shù)在閾值±λ處不連續(xù)。

圖1 硬閾值函數(shù)圖像Fig.1 Function image of hard threshold

軟閾值函數(shù)表達(dá)式為：

(6)

圖2 軟閾值函數(shù)圖像Fig.2 Function image of soft threshold

1.3 改進(jìn)的閾值函數(shù)

為彌補(bǔ)傳統(tǒng)閾值函數(shù)存在的缺點(diǎn)，本文提出一種改進(jìn)的閾值函數(shù)，具體表達(dá)式為：

(7)

本文將從數(shù)學(xué)方面來分析改進(jìn)的閾值函數(shù)在±λ處的連續(xù)性、改進(jìn)閾值函數(shù)的漸近線及克服軟閾值函數(shù)的恒定偏差。

從圖3可以直觀看出，改進(jìn)的閾值函數(shù)可彌補(bǔ)硬閾值函數(shù)在閾值±λ處不連續(xù)的缺點(diǎn)，并且在閾值之后能快速趨近硬閾值函數(shù)，從而彌補(bǔ)了軟閾值函數(shù)存在恒定偏差的缺點(diǎn)。

圖3 改進(jìn)的閾值函數(shù)與軟硬閾值函數(shù)對比Fig.3 Comparison of improved threshold function and soft and hard threshold function

2 信號仿真

本文選用MATLAB中4種標(biāo)準(zhǔn)測試信號(Blocks，Bumps，Doppler，Heavy sine)分別加入標(biāo)準(zhǔn)偏差為8的高斯白噪聲，再使用本文改進(jìn)的閾值函數(shù)及傳統(tǒng)的軟、硬閾值函數(shù)進(jìn)行去噪。

從圖4～7可以看出，軟閾值函數(shù)與硬閾值函數(shù)去噪均會使原始信號嚴(yán)重失真，甚至變形。相比于軟閾值函數(shù)，硬閾值函數(shù)去噪結(jié)果存在較多“毛刺”凸起，這是由硬閾值函數(shù)在λ處不連續(xù)造成的。相比于硬閾值函數(shù)，雖然軟閾值函數(shù)可彌補(bǔ)其不連續(xù)的缺點(diǎn)，使得去噪后的信號較為光滑，但軟閾值函數(shù)在小波系數(shù)大于λ時減去λ，存在恒定偏差，從而造成去噪后的結(jié)果過于平滑，影響去噪效果。另外，改進(jìn)的閾值函數(shù)能在很好地去除高斯白噪聲的同時，改進(jìn)硬閾值函數(shù)的“毛刺”現(xiàn)象及軟閾值函數(shù)的過度光滑，使信號不失真。

圖4 Doppler信號去噪仿真Fig.4 Doppler signal denoising simulation

圖5 Blocks信號去噪仿真Fig.5 Blocks signal denoising simulation

本文將仿真信號的信噪比(SNR)、均方根誤差(RMSE)及能量比例(Per)作為去噪后的評價指標(biāo)：

(8)

(9)

(10)

圖6 Heavy sine信號去噪仿真Fig.6 Heavy sine signal denoising simulation

圖7 Bumps信號去噪仿真Fig.7 Bumps signal denoising simulation

在4種信號中加入標(biāo)準(zhǔn)偏差為8的高斯白噪聲，利用控制變量法，選擇提升小波基為bior6.8，分解與重建尺度為5，分別使用改進(jìn)的閾值函數(shù)與傳統(tǒng)的閾值函數(shù)進(jìn)行去噪處理，評價指標(biāo)結(jié)果見表1。通過SNR與RMSE可以看出，相比于傳統(tǒng)的閾值函數(shù)，改進(jìn)的閾值函數(shù)可大大提高去噪后信號的信噪比，有效去除高斯白噪聲，同時保證仿真信號不失真，去噪效果明顯優(yōu)于傳統(tǒng)方案。從去噪前后的Per可以看出，改進(jìn)的閾值函數(shù)的提升小波在4種仿真信號中均保持較高比例，說明相對于軟閾值函數(shù)與硬閾值函數(shù)，改進(jìn)的閾值函數(shù)可保留較多的原始信號能量成分，而軟閾值函數(shù)與硬閾值函數(shù)在去噪過程中不僅會去除噪聲，同時也會去除較多的有效信號。

表1 仿真信號去噪評價指標(biāo)

3 實(shí)際氣槍震源去噪

3.1 改進(jìn)閾值函數(shù)與傳統(tǒng)閾值函數(shù)實(shí)際去噪對比

地震波是研究地球結(jié)構(gòu)的有效工具之一[8]，天然地震學(xué)與勘探地震學(xué)都面臨著提高地震資料信噪比的難題[9]。氣槍震源是應(yīng)用于地球內(nèi)部結(jié)構(gòu)探測的一種新型人工震源[10]，相比于其他人工震源具有可控性好、低頻成分豐富、信號重復(fù)性高等優(yōu)點(diǎn)。但氣槍震源信號較弱，受外界噪聲干擾較大，且傳播距離有限，衰減與噪聲干擾會使遠(yuǎn)處臺站接收到的氣槍信號信噪比極低，難以分辨。

提高氣槍信號信噪比的常用方法是利用氣槍震源的重復(fù)性進(jìn)行信號疊加[11]，并在此基礎(chǔ)上進(jìn)行帶通濾波處理。由于氣槍信號的有效信號集中在2～5 Hz之間，通常使用2～6 Hz的帶通濾波處理氣槍信號[9]，但該方法存在一定問題，即不能有效濾除與有效信號具有相同頻率的噪聲。本文使用云南賓川氣槍信號發(fā)射臺2016年接收到的實(shí)驗(yàn)激發(fā)數(shù)據(jù)，信號接收臺站是震中距為151.13 km的53252臺站和震中距為115.89 km的53253臺站。在疊加氣槍信號和進(jìn)行2～6 Hz帶通濾波的基礎(chǔ)上對改進(jìn)閾值函數(shù)的提升小波進(jìn)行去噪處理(圖8～9)，結(jié)果表明，改進(jìn)的閾值函數(shù)不僅可克服傳統(tǒng)軟、硬閾值函數(shù)不連續(xù)和存在恒定偏差的缺點(diǎn)，還可極大地提高氣槍震源地震資料的信噪比。

圖8 53252臺站接收到的信號去噪效果Fig.8 Denoising effect of signal received by station 53252

圖9 53253臺站接收到的信號去噪效果Fig.9 Denoising effect of signal received by station 53253

3.2 實(shí)際氣槍震源信號評價指標(biāo)分析

為直觀對比改進(jìn)的閾值函數(shù)的去噪信號與原始信號，將2個信號進(jìn)行疊加，結(jié)果見圖10～11。結(jié)合圖8～11可以明顯看出，波軟閾值函數(shù)與硬閾值函數(shù)的提升小波在去噪過程中會濾除部分氣槍信號的相位信息，導(dǎo)致有效信號嚴(yán)重失真；而改進(jìn)閾值函數(shù)的提升小波閾值去噪則在壓制噪聲的同時幾乎未造成有效信息的丟失。上述分析表明，改進(jìn)的閾值函數(shù)在實(shí)際的人工地震信號去噪中比傳統(tǒng)的閾值函數(shù)表現(xiàn)更優(yōu)。

圖10 53252臺站改進(jìn)閾值函數(shù)去噪信號與原始信號對比Fig.10 Comparison of modified threshold function denoising signal and original signal of station 53252

圖11 53253臺站改進(jìn)閾值函數(shù)去噪信號與原始信號對比Fig.11 Comparison of modified threshold function denoising signal and original signal of station 53253

從表2可以看出，相比于軟、硬閾值函數(shù)，改進(jìn)閾值函數(shù)的提升小波在實(shí)際的氣槍信號去噪中信噪比有大幅提高，均方根誤差減小且能量比例接近1。傳統(tǒng)的閾值函數(shù)去噪能量比例較小，表明在去噪過程中傳統(tǒng)閾值函數(shù)的提升小波會去除許多有效信號，這與圖8～9的結(jié)果相符。分析結(jié)果表明，在實(shí)際的氣槍震源信號處理中，改進(jìn)閾值函數(shù)的提升小波算法可克服軟閾值函數(shù)存在恒定偏差與硬閾值函數(shù)不連續(xù)的缺點(diǎn)，提升遠(yuǎn)處臺站氣槍震源信號的信噪比是一種可用于氣槍震源信號數(shù)據(jù)處理的有效方法。

表2 遠(yuǎn)處臺站信號去噪評價指標(biāo)

3.3 改進(jìn)閾值函數(shù)去噪與帶通濾波對比

圖8～9均為在進(jìn)行900次氣槍信號疊加的基礎(chǔ)上通過2～6 Hz帶通濾波進(jìn)行提升小波的閾值去噪處理，為比較提升小波閾值去噪與帶通濾波去噪的效果，對53253臺站接收到的信號進(jìn)行500次疊加，并在此基礎(chǔ)上分別進(jìn)行2～6 Hz與2～8 Hz的帶通濾波。

從圖12可以看出，氣槍信號在進(jìn)行疊加后，改進(jìn)的閾值函數(shù)的提升小波可較好地壓制氣槍信號的隨機(jī)噪聲，且未造成氣槍信號失真；而帶通濾波雖然可以濾除隨機(jī)噪聲，但同時也會濾除部分有效信號，如55～60 s信號，從而造成一定程度上的信號失真。

圖12 改進(jìn)閾值函數(shù)去噪與帶通濾波對比Fig.12 Comparison of improved threshold function denoising and bandpass filtering

4 結(jié) 語

本文在傳統(tǒng)軟、硬閾值函數(shù)基礎(chǔ)上提出一種新的閾值函數(shù)，并從理論上證明其可在一定程度上彌補(bǔ)軟閾值函數(shù)存在恒定偏差及硬閾值函數(shù)不連續(xù)的缺點(diǎn)，同時將其運(yùn)用到Blocks、Doppler、Heavy sine和Bumps四種加入高斯白噪聲的仿真信號中進(jìn)行去噪實(shí)驗(yàn)。從評價指標(biāo)SNR和RMSE可以看出，改進(jìn)的閾值函數(shù)去噪效果優(yōu)于傳統(tǒng)閾值函數(shù)，并且仿真信號不失真。將改進(jìn)的閾值函數(shù)運(yùn)用到氣槍人工震源中發(fā)現(xiàn)，相比于傳統(tǒng)閾值函數(shù)，改進(jìn)的閾值函數(shù)能在壓制隨機(jī)噪聲并提高遠(yuǎn)處臺站信噪比的同時，保證有效氣槍信號的相位和幅值不失真。將氣槍信號進(jìn)行疊加后發(fā)現(xiàn)，改進(jìn)的閾值函數(shù)在保留有效信號方面優(yōu)于帶通濾波。

致謝：云南省地震局“主動源創(chuàng)新團(tuán)隊”提供氣槍震源數(shù)據(jù)，在此表示感謝。