劉財，商耀達，鹿琪，等. GPR信號去噪的變分模態(tài)分解. 吉林大學(xué)學(xué)報（地球科學(xué)版），2024，54（3）：10421053. doi：10.13278/j.cnki. jjuese.20230154.

Liu Cai， Shang Yaoda， Lu Qi， et al. Variational Mode Decomposition for GPR Signal Denoise. Journal of Jilin University （Earth Science Edition）， 2024， 54 （3）： 10421053. doi：10.13278/j.cnki.jjuese.20230154.

摘要：

為了進一步提高探地雷達（ground penetrating radar， GPR）數(shù)據(jù)的信噪比，壓制由隨機擾動引起的隨機繞射能量，將二維變分模態(tài)分解（two-dimensional variational mode decomposition， 2DVMD）引入二維GPR數(shù)據(jù)的噪聲壓制處理中。首先，對GPR數(shù)據(jù)進行2DVMD處理，并分析各階本征模態(tài)函數(shù)（intrinsic mode function， IMF）分量及其對應(yīng)的頻率波數(shù)域譜來確定雷達剖面中的各回波類型。然后，計算IMF分量與原始數(shù)據(jù)的互相關(guān)系數(shù)來確定信號模態(tài)和噪聲模態(tài)，并對信號模態(tài)進行重構(gòu)得到降噪后的數(shù)據(jù)。理論數(shù)據(jù)和實測數(shù)據(jù)測試表明，相比于傳統(tǒng)的1DVMD法，2DVMD濾波后的含噪正演記錄峰值信噪比由6.44 dB增加到7.72 dB；經(jīng)2DVMD降噪處理后的雷達剖面在保留有效信號的基礎(chǔ)上，可以有效壓制隨機擾動帶來的噪聲，并且得到的雷達剖面同相軸連續(xù)性更好。

關(guān)鍵詞：

探地雷達；二維變分模態(tài)分解；頻率波數(shù)譜；互相關(guān)系數(shù)；去噪

doi：10.13278/j.cnki.jjuese.20230154

中圖分類號：P631.4

文獻標(biāo)志碼：A

收稿日期：20230626

作者簡介：劉財（1963-），男，教授，博士生導(dǎo)師，主要從事地震波場正反演理論、綜合地球物理等方面的研究，E-mail： liucai@jlu.edu.cn

基金項目：國家自然科學(xué)基金項目（41874125）；吉林省科技發(fā)展計劃項目（20200201045JC）

Supported by the National Natural Science Foundation of China （41874125） and the Project" of Jilin Province Science and Technology Development Plan （20200201045JC）

Variational Mode Decomposition for GPR Signal Denoise

Liu Cai1，2，3， Shang Yaoda1， Lu Qi1，2，3， Xu Yangyang1

1. College of GeoExploration Science and Technology， Jilin University， Changchun 130026， China

2. Key Laboratory of Applied Geophysics， Ministry of Natural Resources， Changchun 130026， China

3. National Engineering Research Center of Offshore Oil and Gas Exploration，" Changchun 130026， China

Abstract：

In order to improve the signal-to-noise ratio （SNR） of ground penetrating radar （GPR） data and reduce the random diffraction energy caused by random perturbation， two-dimensional variational mode decomposition （2DVMD） is introduced into the noise reduction processing of 2D GPR data. First， the GPR data is processed by 2DVMD， and the intrinsic mode function （IMF） components and their corresponding frequency and wave number spectra are analyzed to determine the type of each echo showed in the radar profile. Then， the cross-correlation coefficients between the IMF components and the original data are calculated to determine the signal mode and noise mode， and the signal mode is reconstructed to obtain data after noise reduction. Synthetic and practical data tests demonstrate that compared with the traditional 1DVMD method， the peak SNR of the forward recording with noise after 2DVMD filtering increases from 6.44 dB to 7.72 dB. The newly developed method can get significantly improvement on the SNR of GPR data， and obtain radar profiles with better event continuity.

Key words：

ground penetrating radar （GPR）; two-dimensional variational modal decomposition （2DVMD）; frequency and wave number spectrum; cross-correlation coefficient; noise reduction

0" 引言

探地雷達（ground penetrating radar， GPR）是一種快速、高效、無損的地球物理探測技術(shù)，因其具有超淺層快速無損勘探的特點，在地表以及以下幾十米深度范圍可對地下結(jié)構(gòu)進行高分辨率成像，已廣泛應(yīng)用于國防，公安，鐵路、橋梁、隧道、機場建設(shè)、礦產(chǎn)勘查，地質(zhì)、考古等方面［1］。實際探測過程中，因工作對象及環(huán)境的特殊性，GPR回波信號往往會受到強背景雜波及各種隨機噪聲的嚴(yán)重干擾，導(dǎo)致數(shù)據(jù)信噪比較低，嚴(yán)重影響了GPR的目標(biāo)檢測能力。因此，有效壓制GPR回波信號中的噪聲，突出有效信號，具有十分重要的意義。

GPR回波信號中的噪聲信號具有隨機、非平穩(wěn)的特點，常規(guī)的時域分析或頻域分析往往無法達到理想的去噪效果。因此，國內(nèi)外學(xué)者提出多種去噪技術(shù)來提高GPR信號的信噪比［2］。小波因具有良好的時頻局部化特性，常用于GPR信號去噪處理。但該方法降噪效果往往受小波基函數(shù)選取的影響，缺乏自適應(yīng)性；此外小波變換缺乏平移不變性，輸入信號微小的變化都會對不同尺度下小波系數(shù)的能量分布產(chǎn)生較大的影響［3］。針對這一問題，Huang等［4］于1999年提出一種適用于非線性非平穩(wěn)信號的時頻分析新方法——希爾伯特黃變換（HHT），其核心技術(shù)為經(jīng)驗?zāi)B(tài)分解（empirical mode decomposition， EMD），即基于信號自身時間尺度特點，將其按頻率分解為不同的本征模態(tài)函數(shù)（intrinsic mode function， IMF）分量，并通過去掉高IMF分量來壓制噪聲，該方法具有自適應(yīng)時頻處理信號的優(yōu)勢。楊建軍等［5］將EMD應(yīng)用于GPR信號的降噪處理，并驗證了EMD方法的降噪效果優(yōu)于小波閾值去噪。馮德山等［6］采用EMD實現(xiàn)了低信噪比GPR信號的降噪處理，有效突出了雷達剖面異常體的特征，提高了數(shù)據(jù)解釋精度。然而，EMD在進行信號分解時往往會出現(xiàn)模態(tài)混疊或端點效應(yīng)等問題。為此，后人提出一系列類EMD方法對EMD存在的不足進行改進［711］。為了將EMD方法應(yīng)用于二維信號的時頻分析中，Nunes等［1213］提出了二維經(jīng)驗?zāi)B(tài)分解（bi-dimensional empirical mode decomposition， BEMD），并將該方法用于紋理分析和自然圖像處理。王紅振等［14］采用基于集合經(jīng)驗?zāi)B(tài)分解（ensemble empirical mode decomposition， EEMD）的BEMD方法嘗試處理斷裂帶地震數(shù)據(jù)，獲得了較好的處理效果。任青青［15］針對3D路面裂縫的數(shù)據(jù)特點，提出了一種基于BEMD的三維路面裂縫檢測算法。但二維EMD與一維EMD類似，BEMD中也普遍存在模態(tài)混疊現(xiàn)象，并且基于EMD的模態(tài)分解方法均缺乏精確的數(shù)學(xué)模型［16］。

針對EMD方法的不足，Dragomiretskiy等［17］于2013年提出了變分模態(tài)分解（variational mode decomposition， VMD），將信號分量的獲取過程轉(zhuǎn)移到變分框架內(nèi)，通過構(gòu)造并求解約束變分問題實現(xiàn)對原始信號的分解；該方法具有堅實的數(shù)學(xué)理論支撐，并且能夠有效解決EMD方法中存在模態(tài)混疊、端點效應(yīng)等問題。而后，Dragomiretskiy等［18］在一維VMD理論基礎(chǔ)上，設(shè)計并提出一種全新的完全非遞歸的二維變分模態(tài)分解（2Dvariational mode decomposition， 2DVMD），并將該分解方法應(yīng)用于圖像的去噪處理中。VMD法在故障診斷、特征提取、地震信號和GPR數(shù)據(jù)處理中都取得了成功的應(yīng)用［1922］。Xue等［23］將VMD方法應(yīng)用于地震信號的降噪處理中，相較于EMD等方法，VMD具有更好的局部分解能力且對噪聲魯棒性更好。張杏莉等［24］將能量熵引入到VMD方法中，提出了一種自適應(yīng)微震噪聲去除方法。Zhang等［25］將VMD方法應(yīng)用于雷達數(shù)據(jù)噪聲壓制，通過模擬數(shù)據(jù)與實測數(shù)據(jù)驗證，經(jīng)VMD去噪后的雷達剖面能夠更加清晰、準(zhǔn)確識別地下異常體。許軍才等［26］將樣本熵引入到VMD法中，用于篩選分解后的IMF分量并進行數(shù)據(jù)重構(gòu)，成功去除了GPR數(shù)據(jù)中的高斯白噪聲。塔拉［27］利用核密度估計來計算VMD各階IMF分量的概率密度函數(shù)，將所有模態(tài)分為有效信號模態(tài)和噪聲模態(tài)，并利用連續(xù)小波變換和Savitzky-Golay濾波等方法對高頻噪聲模態(tài)做進一步去噪處理，有效提高了2D地震數(shù)據(jù)的信噪比。戴前偉等［28］通過引入能量損失比來實現(xiàn)VMD最優(yōu)模態(tài)數(shù)的自適應(yīng)選擇，并結(jié)合奇異譜分析法對VMD的低頻模態(tài)進行二次濾波以解決中低頻模態(tài)的振蕩現(xiàn)象。2DVMD多應(yīng)用于圖像處理、生物化學(xué)和納米科學(xué)領(lǐng)域［2930］。

GPR數(shù)據(jù)采集通常以單道或多道二維方式掃描，二維數(shù)據(jù)在包含信號時間維度信息的同時，也包含水平方向距離信息和相鄰道之間的共深度點信息。根據(jù)GPR數(shù)據(jù)采集特點，相鄰道之間共深度點的有效信號具有較強的相關(guān)性，而噪聲往往不具備相關(guān)性。以隨機噪聲和介質(zhì)隨機性帶來的噪聲為例，噪聲往往同時存在于時間維度和空間維度上，特別是地下隨機介質(zhì)帶來的隨機繞射波往往具有較強能量且在空間方向存在嚴(yán)重的拖尾現(xiàn)象，使得有效信號受到嚴(yán)重干擾，進而影響目標(biāo)體的檢測判別。常規(guī)二維頻率波數(shù)域濾波方法對含噪數(shù)據(jù)直接進行二維傅里葉變換，并根據(jù)有效信號域噪聲在頻率波數(shù)域上的差異設(shè)計相應(yīng)濾波器來壓制噪聲，但該方法本質(zhì)上仍屬于信號的頻域分析方法。而基于EMD和VMD等時頻域去噪方法對GPR數(shù)據(jù)進行降噪處理時，前人的成果多從時間維度上考慮，僅針對一維時間信號進行噪聲壓制，忽略了存在于空間維度上的噪聲，往往難以獲得理想的噪聲壓制效果。

針對上述問題，本文將2DVMD引入二維GPR數(shù)據(jù)的噪聲壓制處理中。首先，對GPR數(shù)據(jù)直接進行2DVMD處理，得到二維GPR數(shù)據(jù)的IMF分量以及各分量對應(yīng)的頻率波數(shù)譜。然后，對時間空間域的IMF分量及其頻率波數(shù)譜進行分析，確定其所對應(yīng)的回波類型。最后，通過計算各階IMF分量與原始數(shù)據(jù)的互相系數(shù)來定量確定信號模態(tài)和噪聲模態(tài)，并對信號模態(tài)進行重構(gòu)得到降噪后的數(shù)據(jù)。

1" 2DVMD基本原理

2DVMD為1DVMD的自然擴展，可將二維輸入信號自適應(yīng)且完全非遞歸地分解為多個子模態(tài)。

1.1" 二維解析信號

定義頻率域的二維解析信號為

u^as，k（ω）=（1+sgn（〈ω，ωk〉））u^k（ω）=

2u^k（ω），〈（ω），ωk〉＞0;

u^k（ω），〈（ω），ωk〉=0;

0，""" 〈（ω），ωk〉＜0。

（1）

式中：k=1，2，…，K（K為分解模態(tài)數(shù)）；u^k（ω）為第k階IMF的頻率域表示；ωk為第k階IMF分量的中心角頻率；〈·〉表示對多維分量的平均，以提取特定的模態(tài)分量。為了獲得解析信號的單方向性，將頻率域中的一個半平面設(shè)置為0。u^as，k（ω）在時間空間域表示為

uas，k（x）=uk（x）δ（〈x，ωk〉）+jπ〈x，ωk〉δ（〈x，ωk，⊥〉）。（2）

式中：x=（x1，x2）表示二維信號；uk（x）為uas，k（x）的第k階IMF分量；δ（〈x，ωk〉）為狄拉克脈沖響應(yīng)函數(shù)；ωk，⊥為與ωk正交的單位向量；uas，k（x）為由uk（x）與其希爾伯特變換構(gòu)成的復(fù)解析信號；表示卷積。uas，k（x）沿著ωk的方向計算，且是線性獨立的。

1.2" 約束變分問題

建立2DVMD約束方程：

minuk，ωk∑Kk=1αkSymbolQC@［uas，k（x）e－j〈ωk，x〉］22s.t." x：∑Kk=1uk（x）=f（x）。（3）

式中：αk為Tikhonov正則化參數(shù)；f（x）為待分解的二維GPR數(shù)據(jù)。復(fù)解析信號通過與指數(shù)項相乘混合調(diào)制到估算的中心頻率上，并通過計算解調(diào)的L2范數(shù)平方來估計IMF分量的帶寬，使每個分量的估計帶寬之和最小。

引入二次懲罰項及算子求解約束變分問題，采用交替乘子法（alternating direction method of multipliers， ADMM）進行優(yōu)化計算。

1.3" 2DVMD的ADMM優(yōu)化

將約束變分問題轉(zhuǎn)化為非約束問題，其增廣Lagrange表達式為

L（{uk}，{ωk}，λ）=∑Kk=1αkSymbolQC@［uas，k（x）e－j〈ωk，x〉］22+

f（x）－∑Kk=1uk（x）22+〈λ（x），f（x）－∑Kk=1uk（x）〉。（4）

式中，λ為Lagrange算子。非約束問題求解變成

minuk，ωk" maxλ L（{uk}，{ωk}，λ）。（5）

將λ納入二次懲罰項中，得到

∑Kk=1αkSymbolQC@［uas，k（x）e－j〈ωk，x〉］22+

f（x）－∑Kk=1uk（x）

+λ（x）222-

λ（x）422

。（6）

采用ADMM對IMF分量進行優(yōu)化計算，第k階IMF分量的第n+1次頻率域ADMM迭代格式為：

u^n+1k（ω）=f^（ω）－∑Ki=1，i≠ku^n+1i（ω）+λ^（ω）21+2αkω－ωk2，ω∈Ωk，Ωk=ω〈ω，ωk〉≥0；（7）

ωn+1k=∫Ωkωu^n+1k（ω）2dω∫Ωku^n+1k（ω）2dω。（8）

式中：u^n+1k（ω）、f^（ω）、λ^（ω）分別為un+1k（x）、 f（x）、λ（x）的二維傅里葉變換；f^（ω）－∑i≠ku^i（ω）+λ^（ω）2表示信號的第k個殘差。

2DVMD具體更新過程如下。

1）初始化{u^0k}、{ω^0k}、{λ^0k}。

2）在頻率域內(nèi)交替更新u^n+1k、ω^n+1k。

3）更新λ：

λ^n+1（ω）=λ^n（ω）+τf^（ω）－∑Kk=1u^n+1k（ω）。（9）

式中，τ為噪聲容限參數(shù)。

4）設(shè)置收斂誤差ε，重復(fù)上述步驟，直至滿足以下迭代停止條件：

∑Kk=1u^n+1k－u^nk22u^nk22lt;ε。 （10）

2" GPR數(shù)據(jù)的2DVMD去噪

利用2DVMD可將二維GPR數(shù)據(jù)分解為K個IMF分量。不同于圖像的2DVMD處理，對二維GPR數(shù)據(jù)進行2DVMD處理之后，得到的IMF分量包含了從低頻低波數(shù)到高頻高波數(shù)的分量信息。因此，通過IMF分量及其頻率波數(shù)域譜，可定性分析IMF分量所對應(yīng)的回波類型，并進一步計算各階IMF分量與原始數(shù)據(jù)的二維互相關(guān)系數(shù)確定噪聲模態(tài)，最后重構(gòu)有效信號模態(tài)，完成降噪處理。

2.1" GPR數(shù)據(jù)的2DVMD分解

對二維GPR數(shù)據(jù)實施2DVMD，可將原始數(shù)據(jù)分解為若干個IMF和一個殘差：

f（t，x）=∑Kk=1fIMF，k（t，x）+r（t，x）。（11）

式中：f（t，x）為原始二維時間空間域GPR信號；fIMF，k（t，x）為f（t，x）分解后的第k個IMF分量；r（t，x）為殘差。

2.2" 分解模態(tài)數(shù)選取

對二維GPR數(shù)據(jù)進行2DVMD處理時，需要預(yù)設(shè)K，而K的選取對數(shù)據(jù)的自適應(yīng)分解效果起決定性作用。K選取結(jié)果的好壞容易造成欠分解或過分解現(xiàn)象，直接影響后續(xù)的信號重構(gòu)。常規(guī)1DVMD中，K的選取方法主要包括粒子群算法、中心頻率觀察法、樣本熵等。熵是隨機性的統(tǒng)計度量，可用來表征輸入信號（或圖像）的復(fù)雜性。本文對GPR數(shù)據(jù)進行2DVMD分解后，計算各階IMF分量的熵，將熵最小的IMF分量作為VMD的趨勢項。當(dāng)取適當(dāng)?shù)腒值時，趨勢項的熵變小。隨著K的增大，熵逐漸趨于穩(wěn)定。將熵趨于穩(wěn)定的轉(zhuǎn)折點所對應(yīng)的K作為2DVMD的分解次數(shù)，以避免過度分解。

二維圖像熵的計算公式為

H=－∑iPilog2Pi。（12）

式中：H為二維圖像熵；Pi為i個直方圖的概率。

2.3" 信號重構(gòu)

對經(jīng)2DVMD分解后的IMF分量，計算其與原始信號的二維互相關(guān)系數(shù)，確定噪聲分量。皮爾遜相關(guān)系數(shù)R是衡量兩個變量之間線性相關(guān)程度的統(tǒng)計指標(biāo)，其值的范圍為［－1，1］。當(dāng)R接近0時，表示幾乎沒有線性相關(guān)。通過分析R值，可以對兩個變量如何相互關(guān)聯(lián)有一個直觀的理解。其計算公式為

R=∑x∑t（f－fmean）（f^－f^mean）∑x∑t（f－fmean）2∑x∑t（f^－f^mean）2。 （13）

式中：f為原始含噪數(shù)據(jù)；f^為噪聲壓制后的數(shù)據(jù)；f^mean為f^的均值。

3" 數(shù)值實驗

為了考察2DVMD對壓制GPR數(shù)據(jù)中隨機擾動產(chǎn)生噪聲的有效性，設(shè)計不含隨機擾動的參考水平層狀模型和相應(yīng)的隨機等效介質(zhì)模型。采用時域有限差分（FDTD）生成正演數(shù)據(jù)。發(fā)射天線和接收天線均處于地表，采取自激自收的收發(fā)方式。天線主頻為500 MHz，時間采樣間隔為0.02 ns，時窗為50 ns。模型大小設(shè)定為3 m×3 m，網(wǎng)格為301×301，道間距為0.01 m，共采集151道。

不含隨機擾動的參考水平層狀模型（圖1a）正演模擬記錄如圖1b所示，對其進行自動增益處理后，計算對應(yīng)的二維頻率波數(shù)域譜，結(jié)果如圖1c所示。對應(yīng)隨機等效介質(zhì)模型（圖2a）的正演模擬記錄（信噪比為6.44 dB）和二維頻率波數(shù)域譜分別如圖2b、c所示。從正演記錄可以看出，當(dāng)?shù)叵陆橘|(zhì)含隨機擾動時，會產(chǎn)生能量較強的隨機繞射波且分布無規(guī)律性；隨機繞射波在空間上互相干涉疊加后，繞射能量的拖尾會造成反射波對應(yīng)回波信號同相軸間斷，從而影響地下目標(biāo)體的識別判斷（圖2b）。從對應(yīng)的二維頻率波數(shù)譜圖（圖2c）上可以看出，隨機擾動產(chǎn)生的繞射波能量在波數(shù)域?qū)?yīng)高波數(shù)，而反射波能量在波數(shù)域?qū)?yīng)低波數(shù)，水平層的反射回波信號主要集中在0波數(shù)附近。因此，由隨機擾動引起的隨機繞射波主要分布在空間方向上，利用針對時間維度的1DVMD很難對其進行有效壓制。對水平層狀等效介質(zhì)模型的正演模擬記錄增益后進行2DVMD處理，取K=4，分解后得到IMF分量及其對應(yīng)的二維譜（圖3）。從圖3中可以看出，經(jīng)2DVMD處理后，得到的IMF分量包含了從低頻低波數(shù)到高頻高波數(shù)的分量信息。其中：IMF1對應(yīng)反射波信號，IMF3對應(yīng)來自水平層的反射回波

信號，它們的二維譜能量主要集中在低波數(shù)和0波數(shù)附近且二者頻率較低；IMF2和IMF4對應(yīng)于來自不

同傳播方向的繞射波，且頻率較高。

計算各階IMF分量與分解前數(shù)據(jù)的互相關(guān)系數(shù)（表1），確定隨機繞射波所對應(yīng)的噪聲分量。去掉相關(guān)系數(shù)小于0.3的分量IMF2、IMF4，重構(gòu)得到壓制噪聲后的GPR數(shù)據(jù)，結(jié)果如圖4所示。經(jīng)1DVMD分解去噪后的數(shù)據(jù)信噪比（SNR）為

6.28 dB，且對由隨機擾動引起的隨機繞射噪聲并沒有明顯的壓制效果（圖4a）；而經(jīng)2DVMD去噪后

0.091 3的剖面反射波同相軸連續(xù)性更好、沒有間斷假象，信噪比為7.72 dB，有效提高了含噪數(shù)據(jù)的信噪比（圖4b）。

4 "實測數(shù)據(jù)處理

實測數(shù)據(jù)（圖5a）采用某地區(qū)地下管線GPR數(shù)據(jù)中的一條測線。天線主頻為500 MHz，每道采樣點數(shù)為1 178，時窗長度為40.8 ns，道間距為0.019 m，共采集916道，測線長度為17.41 m。對GPR數(shù)據(jù)進行預(yù)處理，包括去直流漂移、自動增益、帶通濾波等，結(jié)果如圖5b所示。雖然預(yù)處理后的數(shù)據(jù)質(zhì)量已經(jīng)得到明顯改善，壓制了部分高頻噪聲和低頻雜波，但仍可以看出剖面中不僅包含目標(biāo)管線對應(yīng)的能量較強且有規(guī)律分布的繞射波，還包含地層界面粗糙和其他隨機擾動帶來的繞射能量，且繞射波干涉疊加造成空間方向同相軸連續(xù)性較差。對預(yù)處理后的數(shù)據(jù)進行2DVMD處理，取K=8，得到的IMF分量及其對應(yīng)的二維頻率波數(shù)域的譜如圖6、圖7所示。同樣，IMF分量包含低頻低波數(shù)到高

頻高波數(shù)的

分量信息。觀察IMF分量及其二維譜，可以直觀地分析各IMF分量所對應(yīng)的回波類型。其中：IMF5和IMF7對應(yīng)于背景反射層的反射波；IMF6和IMF8對應(yīng)于來自不同方向的主要繞射波，可能來自于目標(biāo)體；而IMF2、IMF3、IMF4則對應(yīng)于更高頻的高波數(shù)繞射波，推測包含隨機擾動產(chǎn)生的隨機繞射噪聲。各階IMF分量與分解前含噪數(shù)據(jù)的互相關(guān)系數(shù)見表2，可以看出，IMF2、IMF3、IMF4分量的互相關(guān)系數(shù)均小于0.3，有效驗證了綜合IMF分量及其二維譜的定性分析結(jié)果。

去掉互相關(guān)系數(shù)小于0.3的IMF分量，并進行數(shù)據(jù)重構(gòu)。如圖8所示，經(jīng)2DVMD去噪后雷達剖面空間方向同相軸連續(xù)性得到明顯改善。

對去噪前后的數(shù)據(jù)隨機抽取某一道和某一時刻的空間采樣信號進行詳細對比，結(jié)果如圖9所示?？梢钥闯?，由于預(yù)處理時已進行過帶通濾波，時間維度上的噪聲已得到較好的壓制，而進一步經(jīng)2DVMD去噪處理之后，空間采樣信號的連續(xù)性得到明顯改善，由隨機擾動所引起的信號震蕩已被壓制。

5" 結(jié)論

1）本文利用2DVMD方法對二維GPR數(shù)據(jù)進行處理，通過計算各階IMF分量與原始信號的二維互相關(guān)系數(shù)及分析各階IMF分量及對應(yīng)的二維圖譜，完成模態(tài)篩選過程。

2）與其他二維譜分析方法對數(shù)據(jù)整體變換后設(shè)計相應(yīng)濾波器進行噪聲壓制造成部分有效信息丟失的情況不同，2DVMD可以分別對不同的IMF分量做圖譜分析，保留更多有效信息的同時結(jié)果更加直觀。

3）與圖像處理不同，經(jīng)2DVMD信號分解后各模態(tài)對應(yīng)的二維圖譜具備尺度分離性、方向選擇性、頻帶分離性等物理性質(zhì)。

4）2DVMD可在頻率波數(shù)域和時間空間域共同分析，直接在時間空間域篩選IMF分量，完成信號重構(gòu)，更加準(zhǔn)確、便捷地壓制信號空間域的隨機噪聲。

參考文獻（References）：

［1］" 曾昭發(fā)，劉四新，馮晅，等.探地雷達原理與應(yīng)用［M］.北京：電子工業(yè)出版社，2010.

Zeng Zhaofa， Liu Sixin， Feng Xuan， et al. Principles and Application of Ground Penetrating Radar［M］. Beijing： Publishing House of Electronics Industry， 2010.

［2］" 李磊.探地雷達信號去噪方法的研究與應(yīng)用［D］.北京：中國地質(zhì)大學(xué)（北京），2013.

Li Lei. Research and Application of Denoising Method for Ground Penetrating Radar Signal［D］. Beijing： China University of Geosciences （Beijing）， 2013.

［3］" 石宏理，胡波.雙樹復(fù)小波變換及其應(yīng)用綜述［J］.信息與電子工程，2007，5（3）：229234.

Shi Hongli， Hu Bo.Survey of Dual-Tree Complex Wavelet Transform and Its Applications［J］. Information and Electronic Engineering， 2007， 5（3）： 229234.

［4］" Huang N E， Shen Z， Long S R. A New View of Nonlinear Water Waves：The Hilbert Spectrum［J］. Annual Review of Fluid Mechanics， 1999， 31（1）： 417457.

［5］" 楊建軍，劉鴻福.基于經(jīng)驗?zāi)B(tài)分解的探地雷達信號去噪處理［J］.工程地質(zhì)計算機應(yīng)用，2008（3）：3742.

Yang Jianjun， Liu Hongfu. Ground Penetrating Radar Signal Denoising Based on Empirical Mode Decomposition［J］. Engineering Geology Computer Application， 2008（3）： 3742.

［6］" 馮德山，戴前偉，余凱.基于經(jīng)驗?zāi)B(tài)分解的低信噪比探地雷達數(shù)據(jù)處理［J］.中南大學(xué)學(xué)報（自然科學(xué)版），2012，43（2）：596604.

Feng Deshan， Dai Qianwei， Yu Kai.GPR Signal Processing Under Low SNR Based on Empirical Mode Decomposition［J］. Journal of Central South University （Science and Technology）， 2012， 43（2）： 596604.

［7］" Abd-El-Malek M B， Hanna S S. Using Filter Bank Property to Simplify the Calculations of Empirical Mode Decomposition［J］. Communications in Nonlinear Science and Numerical Simulation， 2018， 62： 429444.

［8］" Wu Z， Huang N E. Ensemble Empirical Mode Decomposition： A Noise Assisted Data Analysis Method［J］. Advances in Adaptive Data Analysis， 2009， 1（1）： 141.

［9］" Yeh J R， Shieh J S， Huang N E. Complementary Ensemble Empirical Mode Decomposition： A Novel Noise Enhanced Data Analysis Method［J］. Advances in Adaptive Data Analysis， 2010， 2（2）： 135156.

［10］" Torres M E， Colominas M A， Schlotthauer G， et al. A Complete Ensemble Empirical Mode Decomposition with Adaptive Noise［C］//2011 IEEE International Conference on Acoustics， Speech and Signal Processing.［S.1.］： IEEE， 2011： 41444147.

［11］" 李廣，丁迪，石福升，等.基于支持向量機的可控源電磁數(shù)據(jù)智能識別方法［J］.吉林大學(xué)學(xué)報（地球科學(xué)版），2022， 52（3）：725736.

Li Guang， Ding Di， Shi Fusheng， et al. Intelligent Recognition of Controlled-Source Electromagnetic Data by Support Vector Machine［J］. Journal of Jilin University （Earth Science Edition）， 2022， 52（3）： 725736.

［12］" Nunes J C， Niang O， Bouaoune Y， et al. Texture Analysis Based on the Bidimensional Empirical Mode Decomposition with Gray-Level Co-Occurrence Models［C］//Seventh International Symposium on Signal Processing and Its Applications.［S.1.］： IEEE， 2003： 633635.

［13］" Nunes J C， Guyot S， Deléchelle E. Texture Analysis Based on Local Analysis of the Bidimensional Empirical Mode Decomposition［J］. Machine Vision and Applications， 2005， 16： 177188.

［14］" 王紅振，李振春，王姣.二維經(jīng)驗?zāi)B(tài)分解在地震斷裂識別中的應(yīng)用［C］//2016中國地球科學(xué)聯(lián)合學(xué)術(shù)年會論文集.青島：中國石油大學(xué)（華東），2016：12591260.

Wang Hongzhen， Li Zhenchun， Wang Jiao. Application of Two-Dimensional Empirical Mode Decomposition in Seismic Fault Identification［C］//Proceedings of the 2016 China Earth Science Joint Academic Annual Conference. Qingdao： China University of Petroleum （East China）， 2016： 12591260.

［15］" 任青青.基于二維EMD的三維路面裂縫檢測算法研究［D］.西安：長安大學(xué)，2017.

Ren Qingqing. The Study of Three-Dimensional Pavement Crack Detection Algorithm Based on Two-Dimensional EMD［D］. Xi’an： Chang’an University， 2017.

［16］" 周義.快速二維經(jīng)驗?zāi)B(tài)分解和相位追蹤方法及其在導(dǎo)波無損檢測中的應(yīng)用［D］.上海：上海交通大學(xué)，2014.

Zhou Yi. Fast Bi-Dimensional Empirical Mode Decomposition and Phase Tracking Method and Their Applications in the Guided-Wave-Based Nondestructive Testing［D］. Shanghai： Shanghai Jiao Tong University， 2014.

［17］" Dragomiretskiy K， Zosso D. Variational Mode Decomposition［J］. IEEE Transactions on Signal Processing， 2013， 62（3）： 531544.

［18］" Dragomiretskiy

K D， Zosso D.Two-Dimensional Variational Mode Decomposition［C］//Proceedings of the International Workshop on Energy Minimization Methods in Computer Vision and Pattern Recognition. Hong Kong： ［s.n.］， 2015： 1316.

［19］" 喬云，李瓊，錢浩東，等.基于VMD與改進小波閾值的地震信號去噪方法研究［J］.物探化探計算技術(shù)， 2021， 43（6）：690696.

Qiao Yun， Li Qiong， Qian Haodong，et al. Seismic Signal Denoising Method Based on VMD and Improved Wavelet Threshold［J］. Computing Techniques for Geophysical and Geochemical Exploration， 2021， 43（6）： 690696.

［20］" 劉沖，馬立修，潘金鳳，等.聯(lián)合VMD與改進小波閾值的局放信號去噪［J］.現(xiàn)代電子技術(shù)， 2021， 44 （21）：4550.

Liu Chong， Ma Lixiu， Pan Jinfeng， et al. PD Signal Denoising Based on VMD and Improved Wavelet Threshold［J］. Modern Electronics Technique， 2021， 44（21）： 4550.

［21］" 劉建昌，權(quán)賀，于霞，等.基于參數(shù)優(yōu)化VMD和樣本熵的滾動軸承故障診斷［J］.自動化學(xué)報，2022， 48 （3）：808819.

Liu Jianchang， Quan He， Yu Xia， et al. Rolling Bearing Fault Diagnosis Based on Parameter Optimization VMD and Sample Entropy［J］. Acta Automatica Sinica， 2022， 48（3）： 808819.

［22］" 冉利民，李健偉，杜娟，等.基于變分模態(tài)分解算法的探地雷達信號去噪研究［J］.世界地質(zhì)，2022，41（1）：198207.

Ran Limin， Li Jianwei， Du Juan， et al. Ground Penetrating Radar Signal Denoising Based on Variational Mode Decomposition［J］. World Geology， 2022， 41（1）： 198207.

［23］" Xue Y J， Cao J X， Wang D X， et al. Application of the Variational-Mode Decomposition for Seismic Time-Frequency Analysis［J］. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing， 2016， 9（8）： 38213831.

［24］" 張杏莉，盧新明，賈瑞生，等.基于變分模態(tài)分解及能量熵的微震信號降噪方法［J］.煤炭學(xué)報，2018，43（2）：356363.

Zhang Xingli， Lu Xinming， Jia Ruisheng， et al. Micro-Seismic Signal Denoising Method Based on Variational Mode Decomposition and Energy Entropy［J］. Journal of China Coal Society， 2018， 43（2）： 356363.

［25］" Zhang X， Nilot E， Feng X， et al. IMF-Slices for GPR Data Processing Using Variational Mode Decomposition Method［J］. Remote Sensing， 2018， 10（3）：476.

［26］" 許軍才，任青文，黃臨平.基于變分模態(tài)分解的探地雷達信號分析方法［J］. 河海大學(xué)學(xué)報（自然科學(xué)版），2018，46（6）：545550.

Xu Juncai， Ren Qingwen， Huang Linping. GPR Signal Analysis Method Based on Variational Mode Decomposition［J］. Journal of Hohai University （Natural Sciences）， 2018， 46（6）： 545550.

［27］" 塔拉.地震數(shù)據(jù)去噪中的幾種變分模態(tài)分解方法［D］.哈爾濱：哈爾濱工業(yè)大學(xué)，2020.

Tara. Several Mode Decomposition Methods in Seismic Data Denoising［D］. Harbin： Harbin Institute of Technology， 2020.

［28］" 戴前偉，丁浩，張華，等.基于變分模態(tài)分解和奇異譜分析的GPR信號去噪［J］.吉林大學(xué)學(xué)報（地球科學(xué)版），2022，52（3）：701712.

Dai Qianwei， Ding Hao， Zhang Hua， et al. Noise Reduction Method of GPR Signal Based on VMD-SSA［J］. Journal of Jilin University （Earth Science Edition）， 2022， 52（3）： 701712.

［29］" 高宏宇，馬立園，董宏麗，等.基于2DVMD與CC 結(jié)合的模態(tài)重構(gòu)算法［J］.吉林大學(xué)學(xué)報（信息科學(xué)版），2019，37（6）：603609.

Gao Hongyu， Ma Liyuan， Dong Hongli， et al. Reconstruction Algorithm Based on 2DVMD and CC［J］. Journal of Jilin University （Information Science Edition）， 2019， 37（6）： 603609.

［30］" 樊博，金旭榮，田瑞，等. 一種2DVMD與NLM結(jié)合的芯片圖像去噪算法研究［J］.計算機測量與控制， 2021， 29（6）：199204.

Fan Bo， Jin Xurong， Tian Rui， et al. Chip Image Denoising Algorithm Based on Bi-Dimensional Variational Mode and Non-Local Mean Filtering［J］. Computer Measurement Control， 2021， 29（6）： 199204.

［31］" Das M， Choudhary T， Bhuyan M K， et al. Non-Contact Heart Rate Measurement from Facial Video Data Using a 2DVMD Scheme［J］. IEEE Sensors Journal， 2022， 22（11）： 1115311161.

［32］" Messagier M， Meguellati S， Mahgoun H. Fringe Pattern Denoising Using Two-Dimensional Variational Mode Decomposition （2DVMD） Method for Inspection of Flatness of Reduced Surfaces［J］. Experimental Techniques，2022， 46（1）： 2741.

［33］" 汪冬華.多元統(tǒng)計分析與SPSS應(yīng)用［M］. 2版. 上海：華東理工大學(xué)出版社，2018.

Wang Donghua. Multivariate Statistical Analysis and SPSS Application ［M］. 2nd ed. Shanghai： East China University of Science and Technology Press， 2018.