地質(zhì)雷達在活動斷裂探測中的應(yīng)用與進展

張 迪，李家存，吳中海，任麗麗

(1.河南工程學(xué)院土木工程學(xué)院，鄭州 451191；2.首都師范大學(xué)資源環(huán)境與旅游學(xué)院，北京 100048；3.中國地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所，北京 100081)

?

地質(zhì)雷達在活動斷裂探測中的應(yīng)用與進展

張 迪1，李家存2，吳中海3，任麗麗2

(1.河南工程學(xué)院土木工程學(xué)院，鄭州 451191；2.首都師范大學(xué)資源環(huán)境與旅游學(xué)院，北京 100048；3.中國地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所，北京 100081)

結(jié)合國內(nèi)外最新研究成果，對地質(zhì)雷達技術(shù)在活動斷裂探測中的應(yīng)用進行了系統(tǒng)分析和總結(jié)，重點對地質(zhì)雷達在城市活斷層探測、古地震探測和斷層識別方法的研究現(xiàn)狀進行了闡述，并結(jié)合理塘活動斷裂的應(yīng)用實例進行評述。最后，提出了地質(zhì)雷達探測活動斷裂現(xiàn)存的問題，對該技術(shù)探測活動斷裂的發(fā)展方向和應(yīng)用前景進行了展望。

地質(zhì)雷達；活動斷裂；城市活斷層；古地震；圖像解譯

0 引言

斷裂出露于地表或隱伏于地下，由于沉積、風(fēng)化等自然動力及人類生產(chǎn)與生活活動影響，地貌形態(tài)會發(fā)生一定程度的退化，而斷裂地下部分受周圍環(huán)境和人類活動影響較小，最大程度上保留了區(qū)域內(nèi)已發(fā)生地震事件的遺跡。傳統(tǒng)探槽方法雖然可以將斷裂附近淺層地下結(jié)構(gòu)直觀反映出來，但探槽位置的合理選擇及開挖將耗費大量人力財力，且檢查頻度較低并對地表環(huán)境破壞不可恢復(fù)。盡管一些地球物理技術(shù)(地震波勘探技術(shù)、電法儀和磁法儀等)已應(yīng)用于活動斷層空間分布的探測，但這些方法不同程度上存在設(shè)備成本昂貴、采集效率慢、分辨率較低和對地表環(huán)境影響等不利因素，使此類方法對活動斷裂探測適用性有限。

地質(zhì)雷達，又名探地雷達(Ground Penetrating Radar, GPR)，作為利用電磁波時頻和振幅特征來確定中淺表地層構(gòu)造的無損地球物理探測技術(shù)，具有操作簡便、檢測范圍廣、探測深度深和分辨率高的特點，已被廣泛應(yīng)用于文物考古[1-2]、資源勘探[3-4]、工程檢測[5-6]、水文水利[7-8]等諸多領(lǐng)域，并取得了較好的應(yīng)用效果。近年來，國外學(xué)者對地質(zhì)雷達探測活動斷裂進行了大量深入的研究工作，并取得了一系列研究成果[9-12]。而國內(nèi)針對地質(zhì)雷達應(yīng)用于斷裂探測也做了一些研究，但尚處于起步階段，主要集中在驗證性研究和定性描述分析方面。本文通過國內(nèi)外大量已有文獻，主要對地質(zhì)雷達在活動斷裂上的國內(nèi)外研究進展進行綜述，并探討其中存在的問題及發(fā)展方向。

1 地質(zhì)雷達技術(shù)簡介

地質(zhì)雷達是采用高頻電磁波(0.1～3 GHz)通過介質(zhì)之間的電性差異(主要是介電常數(shù)和電導(dǎo)率)來確定介質(zhì)結(jié)構(gòu)和物性特征的近地表地球物理探測技術(shù)。地質(zhì)雷達系統(tǒng)主要由中心控制單元、發(fā)射天線、接收天線和測距輪組成。地質(zhì)雷達系統(tǒng)工作時，中心控制單元觸發(fā)信號控制發(fā)射天線發(fā)射電磁波，電磁波在地下介質(zhì)傳播過程中遇到電性差異的界面時發(fā)生反射，接收天線接收反射信號并將異?，F(xiàn)象顯示出來，通過后期數(shù)據(jù)處理并根據(jù)電磁波波形、反射強度和時頻特征等可以推斷地下異常區(qū)域的空間位置、幾何形態(tài)等參數(shù)，以此實現(xiàn)對地下空間結(jié)構(gòu)的有效探測，其工作原理如圖1所示。

圖1 地質(zhì)雷達工作原理示意圖Fig.1 The principle of ground penetrating radar

根據(jù)數(shù)據(jù)采集方式和數(shù)據(jù)采集過程中收發(fā)天線間距是否發(fā)生變化，現(xiàn)有的地質(zhì)雷達系統(tǒng)可分為收發(fā)共置天線對的反射測量、寬角反射測量和透射測量3類[13-14]。針對不同的研究對象，選擇不同測量方式的雷達天線進行數(shù)據(jù)采集。

1.1 收發(fā)共置天線對的反射測量法

數(shù)據(jù)采集過程中，發(fā)射天線和接收天線以固定的間隔沿測線同步移動的方式稱為收發(fā)共置天線對反射測量(見圖2)。當(dāng)收發(fā)天線在地面上同步移動時，即得到由多道雷達道數(shù)據(jù)組成的二維時間剖面圖。由于收發(fā)天線極化方式和天線對的幾何結(jié)構(gòu)保持不變，獲取的二維時間剖面可以直觀準確地反映地下反射信息隨空間位置的變化情況。目前大多數(shù)商用地質(zhì)雷達系統(tǒng)均采用此種反射測量法，如瑞典RAMAC/ALA地質(zhì)系統(tǒng)、加拿大EKKO地質(zhì)雷達系統(tǒng)和中國電波傳播研究所LTD系列雷達系統(tǒng)等。

圖2 收發(fā)共置天線對反射測量方式及雷達圖像Fig.2 The common-offset reflection survey and GPR profile

1.2 寬角反射測量法

與地震探測中的折射和寬角反射方法類似，寬角反射測量主要是通過改變發(fā)射天線和接收天線之間的間距來記錄雷達波傳播速度與深度的變化關(guān)系，測量過程中收發(fā)天線之間的間距是不斷變化的(見圖3)。寬角發(fā)射測量法的優(yōu)點是可以提高圖像信噪比和獲取全域波速信息[15-16]，但其數(shù)據(jù)采集和分析比較復(fù)雜。一般的地質(zhì)雷達系統(tǒng)很少采用寬角反射測量方式，該方法多用于求地下介質(zhì)層位的電磁波速度。

圖3 寬角測量方式Fig.3 The wide-angle reflection survey

1.3 透射測量法

與反射探測相比，透射測量應(yīng)用領(lǐng)域較少。透射探測的雷達稱為鉆孔雷達，通過在2個鉆孔中分別移動發(fā)射和接收天線來實現(xiàn)測量(見圖4)。由于鉆孔探測一般都是在狹長的鉆孔中移動天線，所以一般不考慮天線取向的問題，但也要考慮工作頻率、空間采樣間隔、時窗、時間采樣間隔和鉆孔間距等參數(shù)。

圖4 透射探測的方式及成像剖面示意圖Fig.4 Illustration of a transillumination multioffer gather

2 地質(zhì)雷達在活動斷裂研究中的應(yīng)用

斷裂帶內(nèi)或斷層面附近的介質(zhì)成分比較復(fù)雜，與周圍介質(zhì)的結(jié)構(gòu)和組分存在較大差異，使斷層破碎帶內(nèi)的介質(zhì)與周圍介質(zhì)的相對介電常數(shù)相差較大。當(dāng)雷達天線經(jīng)過斷裂時，與周圍圖像相比，其雷達波的反射強度、波形和時頻特征都會發(fā)生明顯變化，從理論上論證了地質(zhì)雷達探測斷裂的可行性。20世紀90年代以來，國內(nèi)外學(xué)者利用地質(zhì)雷達對世界上不同地區(qū)的斷裂進行了一些研究，主要集中在城市活斷層探測、古地震探測和斷層識別方法3個方面，由最初單頻率天線二維剖面數(shù)據(jù)采集方式到多頻率天線相結(jié)合三維剖面方式，從簡單的定性研究逐漸發(fā)展到半定量、定量研究。

2.1 城市活斷層探測

城市地震和地質(zhì)體的活動是危及人們安全的主要因素之一，準確查明地表附近活斷層的空間分布，確定深部延伸情況，從而揭示地下介質(zhì)的特性和深部構(gòu)造環(huán)境，評估其地震危險性和危害程度，使重要建筑設(shè)施有效避開，最大程度上減少地震災(zāi)害[17-19]。城市活斷層探測常用的技術(shù)手段是地震勘探法，但在第四系蓋層較淺的地質(zhì)環(huán)境，其分層效果較差，提供的上斷點深度往往存在較大偏差。近幾年來，地質(zhì)雷達技術(shù)以其操作簡便、成本低、分辨率高的優(yōu)勢，在覆蓋層薄的地區(qū)，作為一種重要的輔助手段被廣泛應(yīng)用于城市活斷層探測，尤其對于潛伏斷層的位置及活動性，不僅能確定出斷層上部的形態(tài)特征、上斷點埋深和產(chǎn)狀，還能清晰地顯示出斷裂帶附近巖石的變形情況。

Audru 等[20]利用地質(zhì)雷達對經(jīng)過Wellington的走滑斷裂進行探測，通過3個不同位置的二維時間剖面的對比，確定出Wellington 走滑斷裂主斷裂經(jīng)過的區(qū)域，并清晰顯示出斷裂上部的形態(tài)特征。Slatera等[21]利用地質(zhì)雷達評定橫穿Aqaba City 的斷裂，通過50 MHz、100 MHz、200 MHz和400 MHz等4種不同中心頻率天線地質(zhì)雷達的探測效果對比，選擇100 MHz中心頻率的分離式天線共采集7道二維時間剖面，以共中心點測量方式估算出電磁波的平均速度為0.095 m/ns，通過雷達剖面與探槽剖面的效果對比，總結(jié)出經(jīng)過Aqaba City主斷裂的異常特征為傾斜近似45°的直線反射，并用于附近地區(qū)地質(zhì)雷達圖像的解譯。Rashed 等[22-23]為確定Uemachi斷裂區(qū)附近的淺層地質(zhì)構(gòu)造，通過地質(zhì)雷達在Yamato River 區(qū)域采集的二維時間剖面確定出斷層崖的位置和淺表層斷層分布(見圖5)，并進一步分析判斷出斷層呈南北走向。Liberty等[24]結(jié)合地震波、地質(zhì)雷達和磁法儀對Portland Hills 隱伏斷裂進行探測，確定出斷裂的變形區(qū)域至少為400 m，結(jié)合探槽開挖，最后確定出此斷裂歷史上至少發(fā)生過2次較大震級的地震。Khorsandi等[25]利用中心頻率為100 MHz和200 MHz地質(zhì)雷達確定Tehran南部 North Shahre Ray 斷裂的區(qū)域，綜合3處不同位置雷達剖面，不僅準確確定出斷裂的位置，還判斷出斷層主要為逆斷層，斷層傾角呈北方向60°。Carpentier等[26]在確定Christchurch City 淺地表層斷裂的形態(tài)特征中，通過地質(zhì)雷達將由斷裂形成的河流階地不同時期沉積層在雷達圖像上顯示出來(見圖6)，通過分析圖像上異常區(qū)域確定出斷層的分布位置，并判斷其主要以逆斷層為主。

圖5 Uemachi斷裂二維剖面圖[22]Fig.5 The 2D profiles of the Uemachi fault

圖6 Christchurch city 淺層地表二維雷達圖像[26]Fig.6 The profiles of GPR in Christchurch city

國內(nèi)的一些學(xué)者也利用地質(zhì)雷達對城市活動斷層進行了研究。薛建等[27-28]在長春市活動斷層的探測與活動評價中，利用地質(zhì)雷達技術(shù)清晰地顯示出斷層上部的形態(tài)特征、上斷點埋深和巖土分層，并結(jié)合鉆孔資料對雷達分層結(jié)果進行了確認，初步驗證了地質(zhì)雷達在探測城市活動斷層中的有效性。崔國柱等[29]和李征西[30]對不同地球物理方法探測長春市已知斷層的效果進行了對比研究，表明地質(zhì)雷達對活動斷層的剖面形態(tài)研究具有很好的作用，可以提供斷層的形態(tài)特征、近地表的活動規(guī)律，還能為活動斷層的分段性研究提供依據(jù)。李建軍等[31]在東昆侖斷裂帶潛伏斷層的研究中，利用地質(zhì)雷達技術(shù)顯示出潛伏斷層的形態(tài)特征和巖土分層，并結(jié)合沉積序列，對斷層的活動性進行了分析和評價。

2.2 古地震探測

古地震研究可以很大程度上彌補儀器和歷史地震記錄的短暫性和局限性，對進一步揭示地震的發(fā)生時間、空間和強度以及其趨勢性變化規(guī)律，進而研究該地區(qū)的強震活動具有重要意義。地質(zhì)雷達技術(shù)在古地震研究中的應(yīng)用主要在探槽位置的確定和斷層特征參數(shù)的獲取。

探槽是古地震研究的主要環(huán)節(jié)，選擇合理的探槽位置是進行古地震研究的前提。探槽位置選取一般是根據(jù)野外實地調(diào)查，在地震地表遺跡比較明顯的地方開挖，尤其是地面破裂遺跡保存較完整的地區(qū)。但在地質(zhì)調(diào)查中，由于時間比較久遠，或者沉積、風(fēng)化等外界環(huán)境的變化和人類活動的持續(xù)影響，使地表地震遺跡不明顯或已遭受破壞，從而導(dǎo)致古地震研究時探槽位置的選擇具有一定的難度。地質(zhì)雷達作為新型的無損探測技術(shù)，可被用于在地表破裂不明顯的地區(qū)確定探槽的最佳位置。Salvi等[32]采用中心頻率為50 MHz、100 MHz和200 MHz的地質(zhì)雷達天線根據(jù)地表下巖層分布及變形情況，選擇5處不同的地點，結(jié)合探槽剖面確定此處歷史上共發(fā)生3次較大地震事件。Anderson等[33]以南加利福尼亞的逆沖斷層為例證明了地質(zhì)雷達在沉積環(huán)境下確定古地震探槽位置的有效性。Malik等[34]在研究Hajipur 斷裂的活動性時，為選擇合適的探槽位置，選擇SIR 3000型中心頻率為200 MHz的地質(zhì)雷達采集二維和三維地質(zhì)雷達圖像(見圖7)，二維圖像上斷層處的雷達波呈雙曲線反射，連續(xù)時間同相軸發(fā)生中斷，從三維剖面上可以清楚看到斷層傾斜方向及沿斷裂的變形情況。Cahit等[35]在Buyuk Menderes Graben 區(qū)域正斷層探測中，首先采用低中心頻率(250 MHz)天線確定出電磁波異常區(qū)域，然后用高中心頻率(500 MHz)天線對電磁波異常區(qū)域進行重點探測，最后結(jié)合兩種不同頻率天線探測結(jié)果并開挖探槽進行驗證，充分克服了高頻和低頻天線各自的缺點，建立了利用不同頻率天線確定探槽位置的基本方法。

圖7 Sandhwal Village地區(qū)處理后頻率為200 MHz的2D和3D 地質(zhì)雷達圖像[34]Fig.7 Processed 2D and 3D GPR profiles collected using 200 MHz antenna across the HF2 scarp at Sandhwal Village

古地震研究中，地質(zhì)雷達除了用來選擇探槽合適的位置外，在地質(zhì)條件較好的區(qū)域也可以代替探槽，通過獲取斷裂附近的雷達圖像，確定古地震地質(zhì)標(biāo)志及地下層位分布，探測方法由二維圖像逐漸向三維圖像發(fā)展。Chow等[36]利用地質(zhì)雷達和高精度地震反射法對Chihshang 斷裂進行研究，通過中心頻率為200 MHz的地質(zhì)雷達圖像判斷出斷層上斷點、崩積楔和巖土分層，結(jié)合地震勘探法重建Chihshang 斷層的淺層地下介質(zhì)三維模型，為評估地震活動性提供數(shù)據(jù)。Dentith等[37]以1968年 Meckering 地震形成的現(xiàn)已嚴重風(fēng)化的斷層崖為研究對象，結(jié)合附近已知探槽，驗證了地質(zhì)雷達技術(shù)在復(fù)雜地質(zhì)條件下探測斷層崖的可行性。Ercoli等[38-39]提出利用二維和三維地質(zhì)雷達圖像結(jié)合的方式進行古地震探測，首先利用二維地質(zhì)雷達圖像獲取區(qū)域內(nèi)大范圍的異常區(qū)域，然后對電磁波異常集中的區(qū)域進行等間距多道二維剖面數(shù)據(jù)采集，通過后期數(shù)據(jù)處理生成三維數(shù)據(jù)，相對于二維圖像，三維圖像可以將斷層上部形態(tài)特征和產(chǎn)狀等更加形象、直觀地反映出來。

2.3 斷層圖像解譯方法

在利用地質(zhì)雷達探測斷裂的過程中，由于斷裂附近的地質(zhì)構(gòu)造比較復(fù)雜，斷層又分為正斷層、逆斷層和平移斷層，其形態(tài)及發(fā)育方向不相同且分布不均勻，加之電磁波在介質(zhì)傳播過程中的能量衰減和外界因素的干擾，使獲取的雷達剖面中的電磁波特征比較復(fù)雜，并伴有多次反射波、信號振鈴和電磁波繞射等現(xiàn)象，極大地影響了雷達圖像上斷層的正確判讀。斷層圖像解譯方法最初主要依靠目視解譯，與解譯者的經(jīng)驗有很大關(guān)系。為提高斷層解譯準確性，一些學(xué)者采用結(jié)合已知探槽剖面方法進行圖像解譯[40-42]。

隨著計算機技術(shù)的發(fā)展，數(shù)值分析方法被廣泛應(yīng)用于地質(zhì)雷達圖像解譯[43-44]。地質(zhì)雷達圖像數(shù)值模擬的方法較多，但以時間域有限差分法應(yīng)用最為廣泛，基本實現(xiàn)方法是利用計算機以離散差分形式在時間和空間上實現(xiàn)電磁波在地下介質(zhì)中傳播路徑的模擬。Maurizio等[45]利用已知探槽建立斷層的數(shù)值模型，分別模擬出500 MHz和250 MHz的正演圖像，通過與實際剖面對比，總結(jié)出斷層在雷達圖像上的雷達波響應(yīng)特征為：斷裂或斷裂區(qū)域內(nèi)的電磁波反射特征與周圍介質(zhì)的電磁波反射特征差異較大，斷裂兩側(cè)有時會伴有雙曲線繞射現(xiàn)象，但其強度較弱，且連續(xù)層位反射波信號會發(fā)生中斷或錯斷。

二維雷達剖面存在顯示形式單一、無法以多角度方式對斷層的特征參數(shù)進行定量分析等缺點。隨著地質(zhì)雷達硬件及軟件技術(shù)的進一步發(fā)展，地質(zhì)雷達探測斷裂方法開始由二維向三維和二、三維交互方向發(fā)展，對此國外學(xué)者已經(jīng)進行了某些研究，并取得了一定的成果[46-50]。通過三維或者二維和三維相互交互的顯示方式，不僅能將斷裂附近的淺層變形結(jié)構(gòu)以不同的視角顯示出來，而且可以實現(xiàn)對特征參數(shù)的定量分析，例如層位的錯距、上斷點離地面的距離等。

3 應(yīng)用實例

理塘斷裂(禾尼處)位于理塘縣西北部，距縣城57 km(見圖8)。地理位置為東經(jīng)99.89°，北緯30.23°，位于川藏公路318國道旁，從地貌上判斷主要表現(xiàn)為正斷層。正斷層陡坎的左邊為河流，右邊是公路。靠近河流的地方，存在少量較明顯的地表破裂。由于長時間的沉積作用，斷層陡坎右邊的地表形變不太明顯。地面較平坦，夾雜有直徑較大的石塊，位置和地貌概況如圖8所示。

圖8 研究區(qū)位置和地貌概況Fig.8 Location and geomorphologic environment

根據(jù)研究點的地質(zhì)條件和周圍環(huán)境，為了更全面反映斷裂的淺層地下幾何形態(tài)和走向，沿斷裂走向均勻布置4條測線，如圖8上紅色虛線所示，選擇中心頻率為500 MHz和250 MHz的天線沿各測線分別采集水平距離約為70 m的二維剖面。

地質(zhì)雷達圖像的初步解譯結(jié)果如圖9所示。圖9a為測線1位置處中心頻率為500 MHz和250 MHz經(jīng)過數(shù)據(jù)處理后的地質(zhì)雷達圖像，地形校正后的地質(zhì)雷達圖像基本與采集區(qū)域?qū)嶋H地形變化相一致。通過兩種不同頻率地質(zhì)雷達圖像對比，根據(jù)電磁波反射強度，可判斷在水平距離58 m處地表發(fā)生明顯變化(紅色箭頭)，且此處電磁波強度明顯比周圍強度大(紅色陰影區(qū)域)。從500 MHz地質(zhì)雷達圖像可以明顯看出異常區(qū)呈鍥狀分布，且異常區(qū)域右側(cè)電磁波異常自地表至地下呈傾斜狀分布，初步判斷測線58 m處為斷裂經(jīng)過區(qū)域。

圖9b為測線2位置處中心頻率為500 MHz和250 MHz經(jīng)過數(shù)據(jù)處理后的地質(zhì)雷達圖像。對比兩種不同頻率天線圖像，結(jié)合電磁波的強度變化和地形變化情況發(fā)現(xiàn)，水平距離50 m處為電磁波異常區(qū)域，形狀與測線1上紅色區(qū)域相同，初步判斷為斷裂經(jīng)過區(qū)域。水平距離10～30 m之間，存在電磁波異常區(qū)，根據(jù)電磁波變形情況初步斷定為沉積區(qū)域(藍色區(qū)域)。

測線3位置處的雷達圖像如圖9c所示，從不同頻率的雷達圖像上可以清晰地識別水平距離40～70 m之間的電磁波發(fā)射強度與周圍的反射強度有明顯差別(紅色區(qū)域)，且變形區(qū)呈鍥狀分布，此位置與地面上地表破裂位置對應(yīng)，并與測線1和測線2上變形區(qū)域相對應(yīng)，判斷此處為主斷裂經(jīng)過的區(qū)域，并且可以清楚將斷裂的地下形態(tài)反映出來。水平距離30～40 m之間的異常區(qū)域與測線1和測線2(藍色區(qū)域)相對應(yīng)，測線1和測線2僅能反映處異常區(qū)域的局部特征，通過測線3可以清楚判斷出此區(qū)域為水平距離30～70 m異常區(qū)的邊界，結(jié)合正斷層形成的過程，初步推斷為沉積區(qū)。水平距離0～30 m異常區(qū)域在測線3上明顯反映出來，在測線1和測線2上沒有出現(xiàn)，結(jié)合電磁波異常區(qū)域的分布形狀，初步判斷為沉積區(qū)。

圖9 地質(zhì)雷達斷裂區(qū)的二維解譯圖像Fig.9 The interpreted of GPR profiles (Heni)

相對于測線3，測線4位置處的雷達圖像(見圖9d)中可以清晰地識別水平距離30～50 m之間的電磁波發(fā)射強度與周圍的反射強度有明顯的差別(紅色區(qū)域)，此處為主斷裂經(jīng)過的區(qū)域，并與測線1、2和3位置反映出來的一致，藍色異常區(qū)域與測線2反映出來的一致。

從禾尼處4條不同測線雷達圖像的解譯結(jié)果可以得出斷裂帶的寬度約為20 m，并且通過測線2、3、4可以清楚將地下幾何形態(tài)反映出來；在距離斷裂帶10 m和20 m的位置分別存在2個沉積區(qū)，靠近斷裂沉積區(qū)域較大，距離斷裂帶較遠的沉積區(qū)只在測線3上顯示出來，寬度為30 m，深度約為2 m。

為了將斷裂的地下空間分布更形象展現(xiàn)出來，根據(jù)4條二維剖面的探測結(jié)果，選擇在測線1和測線2之間地面比較平坦且碎石分布比較少的區(qū)域，利用中心頻率為500 MHz的天線，以間隔1 m的距離采集10道二維剖面(見圖10a)，測線長20 m左右，在實現(xiàn)雷達圖像和GPS坐標(biāo)信息精確匹配的基礎(chǔ)之上，利用插值的方法實現(xiàn)此區(qū)域內(nèi)雷達圖像的三維顯示(見圖10b)。相對于二維剖面圖，三維效果圖對斷裂位置及其幾何分布形態(tài)的展示更形象。

4 存在問題與展望

4.1 存在問題

縱觀國內(nèi)外已有資料，雖然地質(zhì)雷達技術(shù)在活動斷裂探測中的應(yīng)用取得了較大進展，但仍存在一些需要改進的地方，主要表現(xiàn)在以下幾個方面：

①地質(zhì)雷達的探測效果與研究區(qū)地質(zhì)條件密切相關(guān)，土壤中的水含量對探測效果有較大的影響。不同研究區(qū)斷裂的探測效果存在一定差異，同一研究區(qū)不同地方的探測效果也存在一定差異。

②地質(zhì)雷達技術(shù)具有一定的局限性，需與其他地球物理方法相結(jié)合，才能最大程度上提高探測的效果和圖像判斷的準確性。

③雖然地質(zhì)雷達的硬件和軟件方面都取得較大的發(fā)展，但斷層在雷達圖像上的識別仍然是阻礙了地質(zhì)雷達技術(shù)廣泛應(yīng)用于斷層探測的最大的技術(shù)難題之一。目前斷層在雷達圖像上的解譯仍以人工目視解譯為主，對斷層在探地雷達圖像方面的自動識別和提取有待進一步研究。

4.2 展望

經(jīng)過近幾年的不斷發(fā)展，地質(zhì)雷達技術(shù)已成為一種重要的無損地球物理探測技術(shù)，具有較大的發(fā)展空間，特別是在后期數(shù)據(jù)處理及數(shù)據(jù)解譯方面。隨著地質(zhì)雷達技術(shù)的不斷發(fā)展，三維雷達技術(shù)、雷達圖像解譯和多傳感器集成將成為未來發(fā)展趨勢，并將應(yīng)用于更廣泛的領(lǐng)域。從總體的發(fā)展趨勢上，地質(zhì)雷達在活動斷裂的探測應(yīng)用中仍處于初步階段，在數(shù)據(jù)采集方式、圖像解譯和與其他方法相互結(jié)合應(yīng)用方面仍具有較大的發(fā)展前景。

[1] Zhao W K, Emanuele F, Sara T L, et al. Improved high-resolution GPR imaging and characterization of prehistoric archaeological features by means of attribute analysis [J]. Journal of Archaeological Science, 2015, 54: 77～85.

[2] Bo¨niger* U, Tronicke J. Improving the interpretability of 3D GPR data using target-specific attributes: Application to tomb detection[J]. Journal of Archaeological Science, 2010, 37: 672～679.

[3] Hiroko O, Youngjoo K, Toru T. Depositional and erosional architectures of gravelly braid bar formed by a flood in the Abe River, central Japan, inferred from a three-dimensional ground-penetrating radar analysis[J]. Sedimentary Geology, 2015, 324: 32～46.

[4] Li X J, Hu Z Q, Li S C, et al. Anomalies of mountainous mining paddy in western China[J]. Soil & Tillage Research, 2015, 145: 10～19.

[5] Mercedes S, Higinio G G, Henrique L, et al. Uncertainty evaluation of the 1 GHz GPR antenna for the estimation of concrete asphalt thickness[J]. Measurement, 2013, 46: 3032～3040.

[6] Arias P, Armesto J, Capua D D, et al. Digital photogrammetry, GPR and computational analysis of structural damages in a mediaeval bridge[J]. Engineering Failure Analysis, 2007, 14: 1444～1457.

[7] Brosten T R, Bradford J H, McNamara J P, et al. Estimating 3D variation in active-layer thickness beneath arctic streams using ground-penetrating radar[J]. Journal of Hydrology, 2009, 73: 479～486.

[8] Andrea B. Water content evaluation in unsaturated soil using GPR signal analysis in the frequency domain[J]. Journal of Applied Geophysics, 2010, 71: 26～35.

[9] Pitambar P, Parkash B, Awasthi A K, et al. Concealed thrusts in the Middle Gangetic plain, India: A ground penetrating radar study proves the truth against the geomorphic features supporting normal faulting [J]. Journal of Asian Earth Sciences, 2011, 40: 315～325.

[10] Malik J N, Shah A A, Sahoo A K, et al. Active fault, fault growth and segment linkage along the Janauri anticline (frontal foreland fold), NW Himalaya, India[J]. Tectonophysics, 2010, 483: 327～343.

[11] Mike D, Adam O, Dan C. Ground penetrating radar as a means of studying palaeofault scarps in a deeply weathered terrain, southwestern Western Australia [J]. Journal of Applied Geophysics, 2010, 72: 92～101.

[12] Liberty L, Hemphill-Haley M A, Madin I P. The Portland Hills Fault: Uncovering a hidden fault in Portland, Oregon using high-resolution geophysical methods [J]. Tectonophysics, 2003, 368: 89～103.

[13] 楊峰，彭蘇萍．地質(zhì)雷達探測原理與方法研究[M]．北京：科學(xué)出版社，2006.

YANG Feng, PENG Su-ping. Principles and Methods of Ground penetrating radar [M]. Beijing: Sciences Press, 2006.

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

ZENG Zhao-fa, LIU Si-xin, FENG Xuan, et al. The principle and application of Ground Penetrating Radar[M]. Beijing: Electronics Industry Press, 2010.

[15] Fisher E, McMechan G A, Annan A P. Acquisition and processing of wide-aperture ground penetrating radar data[J]. Geophysics, 1992, 57: 495.

[16] Grasmueck M. 3-D ground-penetrating radar applied to fracture imaging in gneiss[J]. Geophysics, 1996, 61: 1050～1064.

[17] 李樹德．活動斷層分段研究[J]．北京大學(xué)學(xué)報：自然科學(xué)版，1999，35(6)：768～773.

LI Shu-de. Study on segmentation of active faults[J]. Act a Scicentiarum Naturalum Universitis Pekinesis, 1999, 35(6): 768～773.

[18] 徐錫偉．活動斷層、地震災(zāi)害與減災(zāi)對策問題[J]．震災(zāi)防御技術(shù)，2006，1(1)：7～13.

XU Xi-wei. Active faults, associated earthquake disaster distribution and policy for disaster reduction [J]. Technology for Earthquake Disaster Prevention, 2006, 1(1): 7～13.

[19] 沈建文，蔡長青．地震危險性分析與抗震設(shè)防標(biāo)準的確定[J]．地震工程與工程振動，1997，17(2)：27～36.

SHEN Jian-wen, CAI Chang-qing. Seismic hazard analysis and earthquake resistant level [J]. Earthquake Engineering and Engineering Vibration, 1997, 17(2): 27～36.

[20] Audrua J C, Banob M, Beggc J, et al. GPR investigations on active faults in urban areas: The Georisc-NZ project in Wellington, New Zealand[J]. Earth and Planetary Sciences. 2001, 333: 447～454.

[21] Slater L, Niemi T M. Ground-penetrating radar investigation of active faults along the Dead Sea Transform and implications for seismic hazardswithin the city of Aqaba, Jordan[J]. Tectonophysics, 2003, 368: 33～50.

[22] Rasheda M, Kawamuraa D, Nemotoa H, et al. Ground penetrating radar investigations across the Uemachi fault, Osaka, Japan[J]. Journal of Applied Geophysics, 2003, 53: 63～75.

[23] Rasheda M, Kawamuraa D. High-resolution shallow seismic and ground penetrating radarinvestigations revealing the evolution of the Uemachi Fault system, Osaka, Japan[J]. The Island Arc, 2004, 13: 144～156.

[24] Libertya L M, Hemphill-Haley M A, Madinc I P. The Portland Hills Fault: uncovering a hidden fault in Portland,Oregon using high-resolution geophysical methods[J]. Tectonophysics, 2003, 368: 89～103.

[25] khorsandi A, Abdali M, Miyata T, et al. Application of GPR Method Due to Active Faults Determination in Urban Area, Case Study: North Shahre Ray Fault, South of Tehran, Iran[C]. 2011 International Conference on Environment Science and Engineering, 2011.

[26] Carpentier S F, Green A G, Doetsch J, et al. Recent deformation of Quaternary sediments as inferred from GPR images and shallow P-wave velocity tomograms: Northwest Canterbury Plains, New Zealand [J]. Journal of Applied Geophysics, 2012, 81:2～15.

[27] 薛建，賈建秀，黃航，等．應(yīng)用探地雷達探測活動斷層[J]．吉林大學(xué)學(xué)報(地球科學(xué)版)，2008，38(2)：347～350.

XUE Jian, JIA Jian-xiu, HUANG Hang, et al. Application of GPR in Active Fault Detection[J]. Journal of Jilin University(Earth Science Edition), 2008, 38(2):347～350.

[28] 薛建，黃航，張良懷．探地雷達方法探測與評價長春市活動斷層[J]．物探與化探．2009，33(1)：63～66.

XUE Jian, HUANG Hang, ZHANG Liang-huai. The Application of the GPR Method to Detecting and Estimating Active Faults in Changchun[J]. GEOPHYSICAL & GEOCHEM ICAL EXPLORATION. 2009, 33(1): 63～66.

[29] 崔國柱，李恩澤，曾昭發(fā)．活動斷層與地球物理方法[J]．世界地質(zhì)，2003．22(2)：185～190.

CUI Guo-zhu, LI Ee-ze, ZENG Zhao-fa. Active Fault and Geophysical Methods [J]. Global Geology, 2003, 22(2): 185～190.

[30] 李征西，曾昭發(fā)，李恩澤，等．地球物理方法探測活動斷層效果和方法最佳組合分析[J]．吉林大學(xué)學(xué)報(地球科學(xué)版)，2005，35：110～114.

LI Zheng-xi, ZENG Zhao-fa, LI Ee-ze, et al. The Function of Geophysical Method in Active Fault Detection and Discuss of Combining Methods[J]. Journal of Jilin University (Earth Science Edition). 2005, 35: 109～224.

[31] 李建軍，張軍龍．探地雷達在探測潛伏活動斷層中的應(yīng)用[J]．地震，2015，35(4)：83～89.

LI Jian-jun, ZHANG Jun-long. Application of GPR in Surveying Underlied Active Faults [J]. EARTHQUAKE, 2015, 35(4): 83～89.

[32] Salvi S, Cinti F R, Colini L, et al. Investigation of the active Celano-L’Aquila fault system, Abruzzib(central Apennines, Italy) with combined ground-penetrating radarand palaeoseismic trenching[J]. Geophys. J. Int. 2003, 155: 805～811.

[33] Anderson K B, Spotila J A, Hole J A. Application of geomorphic analysis and ground-penetrating radar to characterization of paleoseismic sites in dynamic alluvial environments: an example from southern California[J]. Tectonophysics, 2003, 368: 25～32.

[34] Malik J N, Kumar A, Satuluri S, et al. Ground-Penetrating Radar Investigations along Hajipur Fault: Himalayan Frontal Thrust—Attempt to Identify Near Subsurface Displacement, NWHimalaya, India[J]. International Journal of Geophysics, 2012.

[35] Cahit C Y, Erhan A, Maksim B, et al. Application of GPR to normal faults in the Buyuk Menderes Graben, western Turkey[J]. Journal of Geodynamics, 2013, 65: 218～227.

[36] Chow J, Angelier J, Hua J, et al. Paleoseismic event and active faulting:from ground penetrating radar and high-resolution seismic reflection profiles across the Chihshang Fault. eastern Taiwan[J]. Tectonophysics, 2001, 33: 241～259.

[37] Dentith M A O, Clark D. Ground penetrating radar as a means of studying palaeofault scarps in a deeply weathered terrain, southwestern Western Australia[J]. Journal of Applied Geophysics, 2010, 72: 92～101.

[38] Ercoli M, Pauselli C, Frigeri A, et al. 2D AND 3D GROUND PENETRATING RADAR (GPR) CAN IMPROVE PALEOSEISMOLOGICAL RESEARCHES: AN EXAMPLE FROM THE MT. VETTORE FAULT (CENTRAL APPENNINES, ITALY)[C], GNGTS 2011.

[39] Ercoli M, Pauselli C, Frigeri A,et al. 2D-3D GPR signature of shallow faulting in the Castelluccio di Norcia basin (Central Italy)[C]. EGU General Assembly, 2012.

[40] Gross R, Green A, Holliger K, et al. Shallow geometry and displacements on the San Andreas Fault near Point Arena based on trenching and 3-D georadar surveying[J]. GEOPHYSICAL RESEARCH LETTERS 2002, 29.

[41] Cristina P, Costanzo F, Alessandro F, et al. Ground penetrating radar investigation to study active faults in the Norcia Basin (central Italy) [J]. Journal of Applied Geophysics, 2010, 72: 39～45.

[42] 張迪，李家存，吳中海，等．探地雷達在探測玉樹走滑斷裂帶活動性中的初步應(yīng)用[J]．地質(zhì)通報，2015，34(1)：204～216.

ZHANG Di, LI Jia-cun, WU Zhong-hai, et al. A preliminary application of ground penetrating radar to the detection of active faults along Yushu strike-slip faulted zone. Geological Bulletin of China, 2015, 34(1): 204～216.

[43] 栗毅，黃春琳，雷文太．探地雷達理論與應(yīng)用[M]．北京：科學(xué)出版社，2011.

LI Yi, HUANG Chun-lin, LEI Wen-tai. Ground Penetrating Radar: Theory and Applications [M]. Beijing: Sciences Press, 2011.

[44] Millard S G, Shaw M R, Giannopoulos A, et al. Modeling of subsurface pulsed radar for nondestructive testing of structures[J]. ASCE J Mater Civil Eng 1998, 10: 96～188.

[45] Maurizio E, Cristina P, Alessandro F, et al. "Geophysical paleoses-mology" through high resolution GPR data: A case of shallow faulting imaging in Central Italy[J]. Journal of Applied Geophysics, 2013, 90: 27～40.

[46] Green A G, Gross R, Holliger K, et al. Results of 3-D georadar surveying and trenching the San Andreas fault near its northern landward limit[J]. Tectonophysics, 2003, 368: 7～23.

[47] Gross R, Green A G, Horstmeyer H, et al. 3-D georadar images of an active fault: efficient data acquisition, processing and interpretation strategies[J]. Subsurface Sensing Technologies and Applications, 2003, 4 (1): 19～40.

[48] Vanneste K, Verbeeck K, Petermans T. Pseudo-3D imaging of a low-slip-rate active normal fault using shallow geophysical methods: the Geleen fault in the Belgian Mass River valley[J]. Geophysics, 2008, 73 (1): B1～B9.

[49] McClymont A F, Green A G, Kaiser A, et al. Shallow fault segmentation of the Alpine fault zone, New Zealand revealed from 2-and 3-D GPR surveying[J]. Journal of Applied Geophysics, 2010, 70 (4): 343～354.

[50] Carpentier S F A, Green A G, Langridge R. et al. Flower structures and Riedel shears at a step over zone along the Alpine Fault (New Zealand) inferred from 2-D and 3-D GPR images[J]. Journal of Geophysical Research, 2012:117.

APPLICATION AND PROGRESS OF GROUND PENETRATING RADAR IN ACTIVE FAULT DETECTION

ZHANG Di1, LI Jia-cun2, WU Zhong-hai3, REN Li-li2

(1.College of Civil Engineering, Henan Institute of Engineering, Zhengzhou 451159, China;2.CollegeofResourceEnvironmentandTourism,CapitalNormalUniversity,Beijing100048,China;3.InstituteofGeomechanics,ChineseAcademyofGeologicalSciences,Beijing100081,China)

Combining with the recent research at home and abroad, the application of Ground Penetrating Radar in active fault detection was summed up and concluded. Three aspects on the application of GPR in active fault detection were put forward and reviewed: the application of GPR in urban active fault, the application of GPR in paleoearthquake and the interpretation of fault in the GPR profile, and an example of the method was illustrated to demonstrate the effectiveness of GPR method in Litang active. Finally, the existing problems of the application of GPR in active fault detection was put forward, and the direction of development and the prospect of this technology were also denoted.

Ground Penetrating Radar；active fault, urban active fault; paleoearthquake; data interpretation

1006-6616(2016)03-0733-14

2016-04-13

中國地質(zhì)調(diào)查局地質(zhì)調(diào)查工作項目“長江經(jīng)濟帶活動構(gòu)造與區(qū)域地殼穩(wěn)定性調(diào)查”(121201104000150010)；煤化工資源綜合利用與污染治理河南省工程實驗室開放基金資助項目(502002-B01)

張迪(1987-)，男，博士，講師，主要從事激光和探地雷達相互結(jié)合應(yīng)用研究。E-mail：zhangdi1987228@163.com

李家存(1974-)，男，博士，副教授，主要從事遙感技術(shù)應(yīng)用與研究。E-mail：lijiacun@163.com

P546；P627

A