廖翼強(qiáng)

（廣西路建工程集團(tuán)有限公司，廣西 南寧 530000）

0 引言

巖溶隧道施工過程中，隧道突泥的存在一直是制約隧道施工進(jìn)度、造價以及人員財產(chǎn)安全的重要隱患。本文所依托工程為西南某省高速公路隧道，隧址區(qū)域地質(zhì)環(huán)境復(fù)雜，巖石風(fēng)化程度高且較為破碎；雨季施工時，由于巖溶形成地下水系的活動與地表水系聯(lián)系緊密，掌子面涌水的水質(zhì)變化明顯，往往伴隨出泥現(xiàn)象，在此類地區(qū)進(jìn)行隧道施工建設(shè)時極有可能存在大規(guī)模突泥風(fēng)險［1］。巖溶隧道大規(guī)模突水突泥往往具有突發(fā)性，前期預(yù)兆不明顯，對隧道施工人員安全威脅極大。如2018年貴南客運專線在建的朝陽隧道發(fā)生突泥突水事件，致使2名施工人員死亡；2018年湖北廣水的一個在建隧道發(fā)生突水突泥事故，導(dǎo)致6名施工人員被困。這些均為突發(fā)性突泥突水，前期預(yù)兆不明顯，因此，隧道施工前準(zhǔn)確預(yù)報掌子面前方是否存在充泥型溶洞，是隧道工程建設(shè)亟待解決的關(guān)鍵問題［2］。

探地雷達(dá)是一種利用電磁波對地下淺埋目標(biāo)進(jìn)行無損探測的技術(shù)，近些年來已經(jīng)成為隧道建設(shè)過程中不可或缺的重要地質(zhì)預(yù)報技術(shù)手段。因為電磁波在理論和數(shù)據(jù)處理方面與地震波相似，且地震屬性分析得較為完善，探地雷達(dá)屬性也逐漸由地震屬性引入到探地雷達(dá)數(shù)據(jù)解釋［3］，有些學(xué)者直接借用地震屬性分析軟件來分析探地雷達(dá)屬性，如Senecha等［4］直接用地震處理軟件PROMAX處理雷達(dá)數(shù)據(jù)，并取得了不錯的應(yīng)用效果。目前常用的探地雷達(dá)屬性有三瞬屬性、頻譜屬性、時頻屬性等，探地雷達(dá)數(shù)據(jù)解釋也包含了波形、頻率、振幅、瞬時屬性等方法。如岳全貴等［5］通過分析幾種典型巖溶不良地質(zhì)體的雷達(dá)剖面圖的波形、頻率、振幅等特征，推斷隧道開挖掌子面前方潛在不良地質(zhì)體的類型以及位置、構(gòu)造走向等空間分布特征，開挖結(jié)果驗證了探地雷達(dá)超前地質(zhì)預(yù)報的可行性與準(zhǔn)確性；劉東坤等［6］分析了電磁波的反射、吸收（衰減）特性在充填型溶洞、空腔型溶洞、富水破碎帶三種不良地質(zhì)中的變化；劉成禹等［7］提出一種基于探地雷達(dá)原始記錄數(shù)據(jù)的單道信號，通過將特定區(qū)域單道信號進(jìn)行希爾伯特變換計算獲取信號瞬時振幅、相位及頻率的三瞬屬性，從而推斷巖溶位置和大小的方法，并在實際工程中得到良好應(yīng)用。

雖然探地雷達(dá)屬性分析技術(shù)在地質(zhì)探測領(lǐng)域得到成功應(yīng)用，但目前鮮有文獻(xiàn)針對巖溶隧道充泥型溶洞分析其探地雷達(dá)屬性特征。隧道突泥突水是隧道施工過程中的重要安全隱患，對隧道施工危害極大。因此，分析掌子面前方存在充泥型溶洞時探地雷達(dá)的屬性特征，對地質(zhì)超前預(yù)報判定掌子面前方是否存在突泥突水不良地質(zhì)具有重要實際意義。

1 理論分析

本節(jié)簡要介紹探地雷達(dá)基本理論與現(xiàn)場工作方法。

1.1 探地雷達(dá)基本理論

在隧道建設(shè)過程中采用探地雷達(dá)進(jìn)行超前地質(zhì)預(yù)報，探地雷達(dá)從發(fā)射天線T將高頻脈沖發(fā)射到掌子面前方，然后由接收天線R記錄所發(fā)射的脈沖（見圖1a）。這種方法所記錄的簡單單道波形即稱為單道波（見圖1b）。

脈沖信號向掌子面前方輻射時，一部分能量直接通過空氣和淺表面?zhèn)鞑ゲ⒂山邮仗炀€接收，即直達(dá)波；另一部分能量則穿過介質(zhì)傳播到更深的區(qū)域。當(dāng)透射波遇到空洞溶洞、充填型溶洞和溶蝕破碎帶等不連續(xù)面時，接收天線接收并記錄反射回來的數(shù)據(jù)信息，此部分信息稱為反射波。在隧道中，對于整體巖石，大部分電磁波能量是透射的，而在受損巖石，即質(zhì)量差的巖石中，能量隨著填充材料（空氣、水、泥）的介電特性而變?yōu)榉聪蛏⑸?。通過對雷達(dá)記錄的波形、強(qiáng)度、雙向走時和有效信號強(qiáng)度的分析和處理，可以得到巖溶的具體空間位置、結(jié)構(gòu)和幾何形狀［8］。

圖1 探地雷達(dá)工作原理圖和單道波示意圖

1.2 檢測設(shè)備及工作方法

儀器設(shè)備采用美國進(jìn)口勞雷地質(zhì)雷達(dá)SIR－4000型探地雷達(dá)系統(tǒng)，該系統(tǒng)包括儀器主機(jī)和100MHz天線兩部分。本文檢測采用地質(zhì)雷達(dá)法對圍巖情況進(jìn)行檢測，根據(jù)現(xiàn)場實際情況測線布置如圖2所示，測線方向為掌子面從左至右。

圖2 現(xiàn)場實際測線布置示意圖

2 常用探地雷達(dá)屬性介紹

2.1 三瞬屬性

三瞬屬性是目前探地雷達(dá)屬性應(yīng)用中最為常見的一種屬性分析方法，具體指的是單道信號經(jīng)希爾伯特變換構(gòu)造復(fù)信號計算得到瞬時振幅、瞬時相位與瞬時頻率［9］。其中瞬時振幅屬性反映了電磁波在前方不良地質(zhì)中傳播時，反射信號的能量大小以及能量損失特征，因此可以很好地反映掌子面前方介質(zhì)的性質(zhì)，大致推定異常區(qū)域；瞬時相位屬性最大的優(yōu)點在于相位的大小不隨反射波的振幅變化，可以較好地反映掌子面前方深層位置處的介質(zhì)相位變化；瞬時頻率屬性則是相位的時間變化率，因此可以很好地體現(xiàn)掌子面前方斷層等變化情況。

2.2 頻譜屬性

探地雷達(dá)頻譜屬性是根據(jù)波形圖判定的地質(zhì)異常區(qū)域，通過對截取異常區(qū)域的有效單道信號作傅里葉變換以獲取雷達(dá)數(shù)據(jù)在頻率域中的信息，頻譜屬性可以較為直觀地反映電磁波在穿過異常區(qū)域時，巖溶不良地質(zhì)對電磁波的吸收衰減特性［6］。

2.3 時頻屬性

通過頻譜屬性分析可以獲取信號在頻率域中的整體分布情況，但未能得到各時刻的具體頻率成分信息。因此，有必要對雷達(dá)信號進(jìn)行時頻變換以獲取雷達(dá)信號的時頻局部特征［10］，本文時頻屬性的提取方法為對單道信號數(shù)據(jù)進(jìn)行廣義S變換，其數(shù)學(xué)表達(dá)式為：

其中k和p分別為實際各頻率調(diào)節(jié)參數(shù)，具體大小應(yīng)根據(jù)實際需要選擇。

3 工程實例分析

本文以西南某省高速公路隧道為依托，選取經(jīng)現(xiàn)場開挖驗證后的充泥型溶洞探地雷達(dá)數(shù)據(jù)為分析對象，分別從波形圖、三瞬屬性、頻譜屬性以及時頻屬性等方面分析掌子面前方充泥型溶洞的不同特征。

3.1 工程地質(zhì)概況

隧址屬亞熱帶季風(fēng)氣候，夏長多雨，年平均氣溫18.2℃～21.7℃，年降雨量在1 249～1 673mm之間；隧道地貌屬于巖溶峰叢地貌，地形切割劇烈，溝谷發(fā)育；隧道進(jìn)出口主要由微風(fēng)化灰?guī)r組成，結(jié)構(gòu)面發(fā)育。洞身段主要為微風(fēng)化灰?guī)r，巖質(zhì)較硬，局部節(jié)理發(fā)育，多呈閉合狀態(tài)，巖體較完整，局部較破碎，工程地質(zhì)穩(wěn)定性較好。隧址范圍內(nèi)無大的地表水體，巖溶水主要分布于灰?guī)r中的溶洞、溶蝕裂隙及巖溶破碎帶中。

3.2 隧道掌子面情況

掌子面揭露圍巖為微風(fēng)化灰?guī)r，巖質(zhì)堅硬，巖體完整性較好，節(jié)理裂隙弱發(fā)育，巖體結(jié)構(gòu)面結(jié)合較好，掌子面干燥。

3.3 掌子面探地雷達(dá)記錄

根據(jù)本文1.2節(jié)現(xiàn)場測線布置，獲取掌子面處雷達(dá)原始數(shù)據(jù)，并對原始雷達(dá)數(shù)據(jù)進(jìn)行常規(guī)的數(shù)據(jù)處理，如靜態(tài)校正、增益調(diào)整、去直流漂移、濾波、壓制和剔除干擾波及突出有效波等處理，處理后波形圖如圖3所示。圖3中框線圈定區(qū)域存在明顯的同相軸連續(xù)，上方反射界面約為時間深度50ns處，且后方存在多次強(qiáng)反射。結(jié)合掌子面地質(zhì)情況，按照常規(guī)分析可推斷此區(qū)域為充填型溶洞，具體充填物質(zhì)類別不能精確推定。此外，該異常區(qū)域位于隧道掌子面左前方，具體不良地質(zhì)規(guī)模大小未能得到很好體現(xiàn)。因此，在此區(qū)域進(jìn)行隧道掘進(jìn)時應(yīng)謹(jǐn)慎開挖，短進(jìn)尺、弱爆破、早支護(hù)、勤量測，做好預(yù)防突水突泥措施。實際開挖進(jìn)尺約3m后，掌子面左前方出現(xiàn)大量突泥，開挖結(jié)果與預(yù)期一致，驗證了利用探地雷達(dá)對前方不良地質(zhì)進(jìn)行預(yù)報的必要性。

圖3 處理后雷達(dá)圖像示例圖

3.4 充泥型溶洞探地雷達(dá)屬性分析

本文選取掌子面測線0.8m處（見圖3）單道波數(shù)據(jù)為分析對象，單道波形圖如圖4所示，圖中圈定區(qū)域即為掌子面左前方異常區(qū)域，該區(qū)域波形較為規(guī)律，振幅變化明顯且幅值衰減較慢。此類現(xiàn)象可解釋為電磁波傳播在充泥型溶洞不良地質(zhì)區(qū)域上反射界面時介質(zhì)電性差異較大，因此振幅變化明顯；在充泥型溶洞內(nèi)部傳播時電磁波傳播速度明顯降低，使得上下界面反射波出現(xiàn)明顯分離，即在反射界面后方形成多次反射。

圖4 單道波形圖

3.4.1 瞬時屬性分析

將上述選取的雷達(dá)單道信號作希爾伯特變換并作歸一化處理，可獲得瞬時振幅屬性如圖5所示，圖中顯示在異常區(qū)域振幅衰減較為規(guī)律，主峰明顯且存在多個峰值，表明其振幅衰減較慢并呈規(guī)律性變化。這是因為該異常區(qū)域介質(zhì)較為均勻，電磁波能量衰減較慢。

圖5 瞬時振幅示意圖

3.4.2 幅度譜屬性分析

將上述截取的充泥型溶洞區(qū)段雷達(dá)反射波作傅里葉變換并歸一化處理得到幅度譜，如圖6所示。從圖中可以看出充泥型溶洞反射波主頻明顯，中心頻率分布在80～100MHz之間。

圖6 幅度譜示意圖

3.4.3 時頻屬性分析

圖7為上述選定的單道波進(jìn)行廣義S變換的結(jié)果，參數(shù)k取0.4，p取0.7。經(jīng)時頻屬性分析可獲得單道波各時刻頻率信息，從時頻分布圖中可以看出充泥型溶洞上、下反射界面反射波所對應(yīng)的中心頻率接近雷達(dá)反射波的中心頻率，并隨著反射次數(shù)增加逐漸減小，此現(xiàn)象可解釋為雷達(dá)電磁波在泥漿中傳播時高頻成分不斷衰減所致。

圖7 時頻分布示例圖

4 結(jié)語

本文對經(jīng)開挖驗證后的充泥型溶洞探地雷達(dá)數(shù)據(jù)進(jìn)行多種雷達(dá)屬性分析，結(jié)果表明：

（1）充泥型溶洞雷達(dá)上反射截面振幅變化明顯，后方存在多次反射；反射波主頻明顯，且中心頻率分布在80～100MHz之間。

（2）時頻分布表明此類異常區(qū)域中心頻率隨時間衰減較為明顯，反映電磁波在充泥型溶洞內(nèi)存在頻率衰減現(xiàn)象。

（3）本文采用的多種屬性分析技術(shù)及對充泥型溶洞探地雷達(dá)分析結(jié)果對巖溶隧道探地雷達(dá)數(shù)據(jù)解釋與推斷前方是否存在類似突泥不良地質(zhì)具有一定的參考和借鑒意義。