石國杰，王晉國，侯兆陽，田見校，石 剛

( 1.長安大學(xué) 理學(xué)院,陜西 西安 710064，2.長安大學(xué) 公路學(xué)院,陜西 西安 710064)

地質(zhì)雷達(dá)數(shù)據(jù)處理新方法研究

石國杰1，王晉國1，侯兆陽1，田見校2，石 剛1

( 1.長安大學(xué) 理學(xué)院,陜西 西安 710064，2.長安大學(xué) 公路學(xué)院,陜西 西安 710064)

為了改善探地雷達(dá)的目標(biāo)識別能力和分辨率,對雷達(dá)檢測數(shù)據(jù)采用有向小波降噪、速度提取和偏移處理的方法,模型計算和實際工程測試表明，這些方法能夠提高探地雷達(dá)的目標(biāo)識別和分辨能力,可以應(yīng)用于地質(zhì)雷達(dá)數(shù)據(jù)的后期處理,擴(kuò)展探地雷達(dá)在工程領(lǐng)域的實際應(yīng)用范圍。

探地雷達(dá)；二維有向小波降噪；Hough變換；偏移成像

探地雷達(dá)是一種典型的超寬帶電磁脈沖系統(tǒng)，在近地表探測中，探地雷達(dá)具有兼顧探測深度和空間分辨率的優(yōu)點(diǎn)[1-3]。近十年來 ,隨著信號分析[4]和圖像處理技術(shù)的飛速發(fā)展,一些數(shù)據(jù)處理理論和圖像處理方法被廣泛應(yīng)用于地質(zhì)雷達(dá)目標(biāo)識別領(lǐng)域，使得探地雷達(dá)的目標(biāo)識別和分辨能力得到顯著提升，其已經(jīng)被廣泛應(yīng)用于水壩、公路和隧道等結(jié)構(gòu)的無損檢測和質(zhì)量評定中。本文探討應(yīng)用一種從原始數(shù)據(jù)降噪、速度提取和偏移歸位的圖像目標(biāo)特性識別新技術(shù),對雷達(dá)測量數(shù)據(jù)進(jìn)行分析,從而提高探地雷達(dá)的目標(biāo)識別能力和分辨率。

1 地質(zhì)雷達(dá)數(shù)據(jù)處理方法理論研究

1.1 利用二維Cauchy有向小波降噪

由于探地雷達(dá)接收的信號中存在較多天線直耦波和媒質(zhì)表面直接反射波，它們屬于強(qiáng)雜波信號，較大程度降低了系統(tǒng)對深層弱反射目標(biāo)的探測能力。同時，淺層目標(biāo)的反射波與直達(dá)波的到達(dá)時間比較接近，使得能量相對較弱的目標(biāo)信號淹沒在直達(dá)波之中，不利于目標(biāo)的檢測與識別，二維Cauchy小波則是近年來降噪的常用方法。

1.1.1 二維Cauchy有向小波分析理論

根據(jù)小波分析理論，如果小波函數(shù)的有效支撐在傅立葉變換域展布在頂點(diǎn)為原點(diǎn)的一個凸錐上，則說明它具有方向選擇性[5]。大多數(shù)連續(xù)二維有向小波變換具有方向選擇性，但二維Cauchy小波具有一些獨(dú)特的優(yōu)點(diǎn)，并且數(shù)據(jù)處理過程中參數(shù)選擇比較靈活，適合于對雜波進(jìn)行去除處理。二維Cauchy小波時(空)域函數(shù)的實部和虛部分布如圖1所示。

圖1 二維Cauchy 小波時(空)域函數(shù)圖像

從圖1可以看出，二維Cauchy小波在頻率域是緊支撐的,分布范圍很小,因此具有很好的方向選擇性，而且角分辨能力也較高。同時運(yùn)用二維Cauchy小波作為壓制強(qiáng)雜波信號的基函數(shù)進(jìn)行雜波濾除，也可以大大減少數(shù)據(jù)的運(yùn)算時間。

1.1.2 二維Cauchy小波理論降噪基本方法

應(yīng)用Cauchy小波變換理論降噪的基本方法是:設(shè)地平面為xy平面，在y為常數(shù)的測線上一探地雷達(dá)剖面為y(x,t)，其中x表示測線方向，t表示時間。設(shè)測量的有效信號表示為s(x,t)(不考慮各種噪聲)，雜達(dá)波為d(x,t)。則

y(x,t)=s(x,t)+d(x,t)

(1)

其在傅里葉域的表達(dá)式為：

Y(kx,kt)=S(kx,kt)+D(kx,kt)

(2)

式中:Y(kx,kt)、D(kx,kt)及S(kx,kt)分別表示y(x,t)、d(x,t)及s(x,t)的二維傅里葉變換。在(x,t)域，直達(dá)波d(x,t)疊加在有效信號s(x,t)上，兩者無法分離。當(dāng)測量表面較平緩時，在(kx,kt)域，D(kx,kt)的能量主要集中在kx=0附近。對于點(diǎn)目標(biāo)來說，其回波信號的時距曲線為典型的雙曲線，它的能量在(kx,kt)域分布在很大的一個范圍內(nèi)。

(3)

1.2 有效信號處理方法

1.2.1Hough變換速度提取

Hough變換是近年來國內(nèi)外速度提取常用的方法。Hough變換的基本思想是利用圖像域與數(shù)據(jù)域的點(diǎn)到線的對偶性原理，即在圖像空間中共線的點(diǎn)對應(yīng)于參數(shù)空間里相交的線；相反，在參數(shù)空間中相交于同一點(diǎn)的所有曲線在圖像空間中都有共線的點(diǎn)與之對應(yīng)?；谶@樣的理論，Hough變換把在圖像空間需要解決的問題轉(zhuǎn)換到參數(shù)空間，根據(jù)參數(shù)空間里聚焦點(diǎn)的特點(diǎn)來完成目標(biāo)的檢測任務(wù)。

(4)

式中:Δt為雷達(dá)的采樣間隔時間；Δx為雷達(dá)的采樣間隔距離；(i0,j0)為原空間的坐標(biāo)；(i,j)為參數(shù)空間的坐標(biāo)。(i0,j0)的值是由梯度幅度圖像矩陣位置決定的，而(i,j)的取值范圍與梯度幅度圖像矩陣的行與列相同。通過數(shù)據(jù)處理，對于單目標(biāo)的數(shù)據(jù)圖線可提取出單條雙曲線，對于多目標(biāo)的數(shù)據(jù)可提取出多條雙曲線。

根據(jù)介質(zhì)的特性，給定波速的范圍v(min,max)，并選定一個波速步長Δv，以Δv為步進(jìn)單位對波速范圍內(nèi)的每一個波速值，應(yīng)用計算機(jī)對矩陣進(jìn)行迭代運(yùn)算,可以從Hough變換得到所有矩陣元素中最大值在的矩陣中所對應(yīng)的波速v，即波在介質(zhì)中傳播的速度。

1.2.2 探地雷達(dá)數(shù)據(jù)偏移成像

探地雷達(dá)檢測到的信號是時間信號,根據(jù)電磁波成像理論,點(diǎn)目標(biāo)的特征信號在探地雷達(dá)剖面上的展布范圍遠(yuǎn)大于目標(biāo)的尺寸，因此對目標(biāo)成像必須進(jìn)行歸位處理，該處理技術(shù)稱為偏移。探地雷達(dá)檢測得地表記錄可用波場P(x,z,t)表述,z=0表示地面波場。首先對時間t作傅里葉變換得到P(x,z,ω),采用Hough變換提取速度，再對x作傅里葉變換獲得P(kx,0,ω)，然后乘以全通波指數(shù)exp(-ikz,z)獲得深度z的波場P(kx,z,ω)。計算目標(biāo)成像要在下延拓過程中引入一個終結(jié)條件才能完成，過程的終止是由旅行時鐘讀數(shù)為零來控制。因為探地雷達(dá)信號目前多為單發(fā)單收剖面，相當(dāng)于地震勘探中的自激自收剖面,直接給出零偏移距垂直波數(shù)。對ω依次疊加求和,同時在kx方向作反傅里葉變換在此深度上求出清晰的地表場圖像P(x,z,0)。

2 地質(zhì)雷達(dá)實際數(shù)據(jù)圖像處理

通過實測數(shù)據(jù)，對上面提出的方法進(jìn)行驗證。數(shù)據(jù)采集使用某電波所生產(chǎn)的LTD-2100型探地雷達(dá)，選擇400MHz天線進(jìn)行數(shù)據(jù)采集。首先對采集到的數(shù)據(jù)進(jìn)行二維有向小波降噪，再次通過Hough變換進(jìn)行速度提取，最后利用提取的速度進(jìn)行偏移成像處理。采用該系統(tǒng)方法對西安至寶雞高速路面結(jié)構(gòu)變形進(jìn)行評定，并判斷脫空的大小和位置。圖2(a)表示測試路段某位置的原始剖面圖，圖2(b)表示經(jīng)過本文研究的方法首先進(jìn)行降噪、再次經(jīng)過速度提取和偏移成像處理后獲取的剖面圖。處理后的圖像噪聲得到了壓制，分辨率較高，且開挖的結(jié)果與測試的數(shù)據(jù)分析的結(jié)果相吻合。

圖2 西寶高速測試路段某位置數(shù)據(jù)圖像處理

3 結(jié) 語

通過上面的實際數(shù)據(jù)處理分析可知，應(yīng)用有向小波變換可以去除直達(dá)波和壓制噪聲，通過Hough變換提取速度并對雷達(dá)數(shù)據(jù)進(jìn)行偏移處理，能夠提高探地雷達(dá)的分辨能力。實踐證明這些方法可以推廣應(yīng)用于工程質(zhì)量的無損檢測和評定，是探地雷達(dá)數(shù)據(jù)處理行之有效的方法之一。

[1] 吳梅容,曾惠珍.探地雷達(dá)在瀝青路面中的應(yīng)用研究綜述[J].福建建材,2012(12):30-32.

[2] 高拴會,毋光榮,裴少英.探地雷達(dá)在巖溶探測中的應(yīng)用[J].勘察科學(xué)技術(shù),2015(2):59-61.

[3] 李術(shù)才,劉斌,李樹忱,等.基于激發(fā)極化法的隧道含水地質(zhì)構(gòu)造超前探測研究[J].巖石力學(xué)與工程學(xué)報,2011,30(7):1297-1309.

[4] 陸禮訓(xùn)，鄧世坤，冉彌.探地雷達(dá)在隧道施工超前探測中的應(yīng)用[J].工程地球物理學(xué)報，2010，7(2)：201-206.

[5] 鐘雷文.基于MATLAB的煤田核測井信號小波濾波方法研究[J].科學(xué)技術(shù)與工程,2009,9(2):398-401.

[6] 楊富社,徐春龍, 李純.二維小波變換方向的特性[J].長安大學(xué)學(xué)報(自然科學(xué)版),2007,27(4):107-110.

Study on New Methods of Ground Penetrating Radar Data Analysis

SHI Guojie， WANG Jinguo， HOU Zhaoyang， TIAN Jianxiao， SHI Gang

(Chang′an University, Xi′an 710064, China)

In order to improve the target recognition and resolution of ground penetrating radar (GPR), the detect data have been used to speed extraction and migration process after the directional wavelet denoising. Through engineering test, these methods can enhance the target recognition and resolution of GPR effectively, which can be applied in data post processing and extend the practicability of GPR in engineering field.

ground penetrating radar；two dimension directional wavelet denoising；hough transform；migration imaging

2016-11-05

石國杰(1992-)男，江蘇靖江人，在讀碩士研究生，主要從事固體力學(xué)方面的研究.

陜西省工業(yè)攻關(guān)計劃項目(2013K06-27)

10.3969/i.issn.1674-5403.2017.01.012

TU19

A

1674-5403(2017)01-0042-03