吳寶杰，高岳洲，董 躍

(1.浙江省建筑科學(xué)設(shè)計(jì)研究院有限公司，浙江 杭州 310012；2.浙江省大成建設(shè)集團(tuán)有限公司，浙江 杭州 310012；3.浙江省建設(shè)工程質(zhì)量檢驗(yàn)站有限公司，浙江 杭州 310012)

1 引 言

地下廢棄化糞池年代較遠(yuǎn)、無(wú)資料可查，難以知道其準(zhǔn)確位置，給工程界帶來(lái)了一定的困擾。為了準(zhǔn)確定位廢棄化糞池的位置，通過(guò)各種探測(cè)方法的嘗試，由于地質(zhì)雷達(dá)具有探測(cè)效率高、分辨率高、能夠?qū)崟r(shí)顯示地下探測(cè)結(jié)果且屬于無(wú)損探測(cè)等特性，故優(yōu)先選用該方法。地質(zhì)雷達(dá)憑借其優(yōu)點(diǎn)，在工程中得到越來(lái)越廣泛的應(yīng)用。薛飛等[1]利用探地雷達(dá)探測(cè)地下障礙物。岳崇旺等[2]利用探地雷達(dá)探測(cè)地下空洞；趙明堂[3]利用地質(zhì)雷達(dá)調(diào)查城市地面塌陷隱患；許澤善等[4]利用探地雷達(dá)檢測(cè)道路脫空空洞病害。

地質(zhì)雷達(dá)信號(hào)的瞬時(shí)分析能夠更真實(shí)地反應(yīng)目標(biāo)物的信息，有助于提高數(shù)據(jù)解釋的準(zhǔn)確性[5-8]。Hilbert變換利用實(shí)部與虛部的關(guān)系，定義出任意時(shí)刻的瞬時(shí)幅度、瞬時(shí)相位和瞬時(shí)頻率[9]。這三個(gè)瞬時(shí)參數(shù)可以表征信號(hào)的動(dòng)力學(xué)特征，獲取信號(hào)中更多有效的信息。

2 地質(zhì)雷達(dá)探測(cè)技術(shù)

地質(zhì)雷達(dá)是利用超高頻電磁波探測(cè)地下介質(zhì)分布的一種探測(cè)技術(shù)。地質(zhì)雷達(dá)可對(duì)隧道襯砌質(zhì)量[10]、機(jī)場(chǎng)跑道質(zhì)量、各類(lèi)地下管線[11]、地基、路基、道橋質(zhì)量、工程地質(zhì)、地下埋設(shè)物、放射性廢棄物等進(jìn)行非破壞性探測(cè)。地質(zhì)雷達(dá)發(fā)射的電磁波頻率越高，空間分辨率就越高，電磁波在地下介質(zhì)中衰減就越強(qiáng)，探測(cè)距離就越小。

化糞池是處理糞便并加以過(guò)濾沉淀的設(shè)備，是埋在地下的容器。地質(zhì)雷達(dá)探測(cè)地下廢棄化糞池，其工作原理是發(fā)射天線向地表以下發(fā)射電磁波信號(hào)，當(dāng)遇到電磁阻抗發(fā)生變化的目標(biāo)體(如化糞池)時(shí)，電磁波便發(fā)生反射，由接收天線接收反射波并成圖，探測(cè)原理如圖1所示。

圖1 地質(zhì)雷達(dá)探測(cè)原理示意圖Fig.1 Schematic diagram of ground penetrating radar detection principle

3 Hilbert變換原理

假定一連續(xù)的時(shí)間信號(hào)f(t)，其Hilbert變換的定義如下[9]

假設(shè)窄帶信號(hào)μ(t)=a(t)cos(ω0t+θ0)，其共軛信號(hào)υ(t)為：

將它們組成一個(gè)復(fù)信號(hào)：

z(t)=μ(t)+jυ(t)

瞬時(shí)幅度:

瞬時(shí)相位：

瞬時(shí)頻率：

4 數(shù)據(jù)采集

工程位于杭州地鐵5號(hào)線江暉路站附近興業(yè)嘉園西側(cè)，據(jù)當(dāng)?shù)鼐用窕貞洠@個(gè)地方以前有個(gè)化糞池，但具體在哪個(gè)位置資料缺失無(wú)從考證，為了確保地鐵車(chē)站施工過(guò)程的安全性，需要找到該化糞池的具體位置。本次探測(cè)使用加拿大Sensors & Software 公司生產(chǎn)的pulse EKKO PRO型地質(zhì)雷達(dá)，根據(jù)場(chǎng)地環(huán)境和探測(cè)目的選用100 MHz頻率天線探測(cè)。采集參數(shù)如下：天線間距1.0 m，采樣間隔0.8 ns，時(shí)窗180 ns，疊加次數(shù)32次。本次地質(zhì)雷達(dá)探測(cè)采用剖面法，即地質(zhì)雷達(dá)發(fā)射天線和接收天線沿測(cè)線同步移動(dòng)的測(cè)量方法。這

種方式能準(zhǔn)確描述測(cè)線下方各反射界面的形態(tài)。根據(jù)場(chǎng)地條件，本次探測(cè)興業(yè)嘉園西側(cè)化糞池共布置6條測(cè)線，編號(hào)為1#～6#測(cè)線，見(jiàn)圖2。

圖2 測(cè)線布置及孔洞異常范圍Fig.2 Survey line layout and abnormal range of holes

5 數(shù)據(jù)處理與解釋

地質(zhì)雷達(dá)采集的數(shù)據(jù)包含各種噪聲和干擾因素，因此需要對(duì)原始數(shù)據(jù)做進(jìn)一步的處理。數(shù)據(jù)處理的目的就是提高信噪比，突出目標(biāo)物的有效信號(hào)。數(shù)據(jù)處理主要有濾波、增益。濾波可以把信號(hào)中的低頻或高頻等干擾成分濾掉。由于信號(hào)傳播過(guò)程的衰減性，增益可以對(duì)信號(hào)的衰減進(jìn)行一定的補(bǔ)償，從而增強(qiáng)信號(hào)能量。通過(guò)對(duì)原始數(shù)據(jù)進(jìn)行Dewow濾波、時(shí)間域?yàn)V波、空間域?yàn)V波、信號(hào)增益等一系列處理，其中1#、2#測(cè)線雷達(dá)剖面圖見(jiàn)圖3、圖4。廢棄化糞池的地質(zhì)雷達(dá)信號(hào)呈現(xiàn)出頻率較低、同相軸缺失錯(cuò)斷、繞射等特征，如圖3、圖4紅線所劃異常范圍。圖3中1#測(cè)線異常位于測(cè)線0.6～6.6 m，埋深1.0～3.0 m；圖4中2#測(cè)線異常位于測(cè)線0.8～6.8 m，埋深1.0～3.0 m，結(jié)合3#～6#測(cè)線探測(cè)結(jié)果得到圖2藍(lán)線所劃孔洞異常平面分布范圍。

圖3 1#測(cè)線地質(zhì)雷達(dá)剖面Fig.3 Ground penetrating radar profile of line 1

圖4 2#測(cè)線地質(zhì)雷達(dá)剖面Fig.4 Ground penetrating radar profile of line 2

通過(guò)對(duì)地質(zhì)雷達(dá)信號(hào)進(jìn)行Hilbert變換，得到各測(cè)線瞬時(shí)振幅圖，瞬時(shí)振幅反映出給定時(shí)刻反射信號(hào)能量的大小，其中1#、2#測(cè)線瞬時(shí)振幅圖見(jiàn)圖5、圖6，從圖中可以較為明顯地反應(yīng)廢棄化糞池的圖像特征，如圖5、圖6黑線所劃異常范圍。圖5中1#測(cè)線異常位于測(cè)線0.6～6.8 m，埋深1.0～3.1 m；圖6中2#測(cè)線異常位于測(cè)線0.8～6.8 m，埋深1.0～3.0 m。

為了避免廢棄化糞池影響工程施工，根據(jù)地質(zhì)雷達(dá)所探測(cè)廢棄化糞池位置和埋深指導(dǎo)開(kāi)挖，現(xiàn)場(chǎng)開(kāi)挖照片見(jiàn)圖7。根據(jù)開(kāi)挖情況可知，廢棄化糞池的位置和埋深與探測(cè)結(jié)果基本一致。

圖5 1#測(cè)線瞬時(shí)振幅Fig.5 Instantaneous amplitude of line 1

圖6 2#測(cè)線瞬時(shí)振幅Fig.6 Instantaneous amplitude of line 2

圖7 開(kāi)挖驗(yàn)證Fig.7 Excavation verification

6 結(jié) 語(yǔ)

隨著城市的高速建設(shè)，對(duì)地下空間的利用率越來(lái)越高，由于地下障礙物的隱蔽性，將影響建設(shè)工程施工，而且還可能帶來(lái)巨大的經(jīng)濟(jì)損失。地質(zhì)雷達(dá)憑借其探測(cè)效率高、準(zhǔn)確率高、結(jié)果直觀，在地下障礙物探測(cè)領(lǐng)域的應(yīng)用越來(lái)越多。

地質(zhì)雷達(dá)電磁波的傳播特性反映了地下媒質(zhì)電磁性參數(shù)的差異。廢棄化糞池與周?chē)寥来嬖诿黠@的電磁性差異，因此地質(zhì)雷達(dá)能夠有效地探測(cè)地下廢棄化糞池的位置。通過(guò)對(duì)地質(zhì)雷達(dá)信號(hào)進(jìn)行Hilbert變換，提取信號(hào)的瞬時(shí)振幅，能夠更直觀地顯示廢棄化糞池的圖像特征，從而指導(dǎo)工程施工，避免經(jīng)濟(jì)損失。