趙寧雨，潘金秋，喻海軍，梁 波

(1.重慶交通大學 土木建筑學院，重慶 400074；2.重慶航天職業(yè)技術學院，重慶 400021)

?

巖溶裂隙水的地質雷達信號分析與應用

趙寧雨1，潘金秋1，喻海軍2，梁 波1

(1.重慶交通大學 土木建筑學院，重慶 400074；2.重慶航天職業(yè)技術學院，重慶 400021)

從巖溶地質發(fā)展規(guī)律的角度分析了富水巖溶裂隙水的發(fā)展規(guī)律；研究了地質雷達電磁信號對巖溶裂隙水的反射特征和基于雷達信號的波速-介電常數(shù)關系的含水量估算方法；具體工程實例應用分析表明：地質雷達對隱埋在低含水量值的巖體中的裂隙水及裂隙巖體環(huán)境的分辨力良好，巖溶裂隙水含水量的估算方法可行。驗證了基于地質規(guī)律認識的信號分析方法對巖溶裂隙水不良地質圈定和含水量估算方法的有效性。

隧道工程；巖溶裂隙水；地質雷達；地質規(guī)律；地質預報；含水量

0 引 言

巖溶地區(qū)隧道建設中遇到的地質災害多為不同程度發(fā)育的巖溶問題引發(fā)。限于地質條件的復雜性，在勘察設計階段，尤其是在隧道深埋的情況下，往往難以較準確查清不良地質體的性質。在巖溶不良地質中，巖溶裂隙水是施工中的重要風險來源。目前三車道大斷面隧道越來越多，因巖溶水引發(fā)的突水、突泥和塌方等工程事故更易發(fā)生，給施工安全和掘進速度造成嚴重影響[1-3]。近年，國內外在隧道施工中運用TSP、地質雷達(GPR)，瞬變電磁法(TEM)、陸地聲納法、紅外探水等技術進行超前地質預報工作，積累了豐富的經驗。但出現(xiàn)預報失敗情況的幾率還是相當?shù)母?，尤其是巖溶裂隙含水量的預報[4]。隧道的巖溶突水預報是國內外工程建設中還沒有很好解決的難題，其整體的技術水平還不能滿足指導隧道工程施工控制安全的要求。相對于TSP、TRT、陸地聲納法等超前地質預報物探方法，探地雷達是目前隧道地質預報短距離探水的主要手段。由于巖溶裂隙水通常隱埋在低含水量的基巖中，而電磁波具有在基巖中的傳播速度快、地質雷達分辨力高和裂隙巖體含水層的水文特性異質性高等特點，使得地質雷達在巖溶裂隙水巖體環(huán)境中的分辨力比其他地球物理探測方法要高。

目前，在采用地質雷達技術進行針對巖溶裂隙水的預報及相關研究中，結合巖溶地質規(guī)律分析的雷達信號解譯和含水量估算等方面還缺乏準確有效的系統(tǒng)預報方法。因此，開展巖溶富水區(qū)域巖溶裂隙水及其含水量估算的高精度超前預報技術具有重要的理論意義和實用價值。

1 巖溶裂隙水的地質發(fā)育特征

巖溶裂隙水是指賦存于可溶性巖石裂隙中的地下水，可包括在廣義的基巖裂隙水范疇內[5]。在巖溶地區(qū)，巖石類型、空隙類型和地下水類型之間有一定的對應關系。可溶性巖石(或含有可溶質成分的碎屑巖)中往往存在裂隙或空隙，在巖溶作用下，發(fā)育成一個良好的儲水構造和水流通道，在一定的補給和排泄條件下形成不同的發(fā)育形態(tài)[6]。

1)褶皺和斷裂構造對巖溶裂隙水的發(fā)育起主要控制作用[7]。通常褶皺和斷裂構造產生的裂隙非常發(fā)育，形成構造溶蝕裂隙密集帶；在斷裂構造節(jié)理裂隙的交匯處容易形成小型溶蝕通道，成為地下水的主要儲存空間和運移通道；巖溶含水層的富水程度與巖溶裂隙發(fā)育程度密切相關，巖溶含水層的水量往往比較豐富。

2)巖體成分對巖溶和巖溶裂隙水的發(fā)育有重要影響。塑性較強的高泥質含量巖石的節(jié)理裂隙發(fā)育程度低，不利于巖溶的發(fā)育，且透水性差而常成為良好的隔水層；成分復雜的巖層一般沉積間斷比較多，節(jié)理裂隙較發(fā)育，有利于形成溶蝕比較均勻的巖溶含水層，巨厚的純灰?guī)r的構造裂隙呈寬而長的稀疏狀，有利于形成大型巖溶洞穴。

3)近代巖溶后期的巖溶裂隙水分布不均的活躍性，對隧道施工安全影響較大。近代巖溶后期溶蝕發(fā)育，多呈溶隙狀出現(xiàn)，當彼此貫通時呈串珠狀，進而可形成較為寬大的溶蝕裂隙。巖溶裂隙發(fā)育的寬度大小和聯(lián)通程度的差異，使巖溶裂隙水流速和流量狀態(tài)在空間上分布不均，進而強化了裂隙空間的方向性和透水的各向異性，造成巖溶裂隙水的分布不均。

2 巖溶裂隙水電磁信號的特征

2.1 地質雷達信號對水體的反射特征

在假設為均勻、線性、各向同性的低耗損介質中，電磁波以電場表達的麥克斯韋方程可寫為[8]：

(1)

式(1)中第2項為能量耗散，第3項為能量存儲。在低損耗介質中，能量耗散往往較能量儲存較小，地質雷達能進行有效的探測。

求解式(1)，無界均勻有耗介質中的電磁波的電場表達式為：

E(r)=E0e-jKr

(2)

式中：E0為場源處的電場強度；r為空間坐標；K=β-jα(α為衰減系數(shù)，β為相位常數(shù))。

在地質雷達的應用中，可將介質在空間中的介電常數(shù)(ε)和電導率(σ)從矢量形式簡化為與場無關的標量，其衰減系數(shù)見式(3)：

(3)

對電導率很大的地下水體，可假設σ/(ωε)>>1，化簡式(3)為[5]：

(4)

在式(4)中磁導率(μ)在介質中變化一般可取1，角頻率(ω)與雷達天線有關。由式(4)可知電磁波在水體中傳播時其高頻成分的衰減系數(shù)較大。同時，經V.I.Arkhipov[9]的研究，水的介電常數(shù)實部可寫為：

(5)

其中，水在室溫附近的弛豫時間t=8.28×10-12s，靜態(tài)相對介電常數(shù)εs=81，高頻相對介電常數(shù)ε∞=5.6，可知電磁波在水中的顯著損耗產生于高頻段。從式(4)和式(5)可以看出，在電導率較大的水體探測中，宜使用低頻天線(中心頻率)。

地質雷達反射信號的能力衰減問題可由波阻抗的表達式描述：

(6)

由于磁導率(μ)在介質中變化較小，波阻抗主要由相對介電常數(shù)決定。同時，功率反射系數(shù)R表達式為：

(7)

其中，入射波電場(Ei)與反射波電場(Er)之間有如式(8)的關系：

(8)

對高電阻介質，在從介電常數(shù)為ε1的介質向介電常數(shù)為ε2的介質垂直入射情況下，有：

(9)

常溫下，地層中可溶性巖石(石灰?guī)r、白云巖等)的介質相對介電常數(shù)在4～15的范圍內，而水的介質相對介電常數(shù)約為81，因此，雷達反射信號在濕潤灰?guī)r與水的界面上的反射強烈，且雷達波在含水層界面的反射波相位與入射波相差180°。

2.2 基于地質雷達信號的地層含水量估算原理

電磁波在介質中的傳播速度主要由地層的介電常數(shù)決定，介電常數(shù)和波速的計算表達式為：

(10)

其中損耗因子為：

(11)

式中：c為電磁波在空氣中的傳播速度，c=3×108m/s； f為電磁波的頻率；ε′(f)和ε″(f)分別為介電常數(shù)ε(f)=ε′(f)-jε″(f)的實部和虛部；u為相對磁導率，u=1；σdc為直流電導率；ε0為空氣中的介電常數(shù)，ε0=8.854×10-12F/m。

由于介電常數(shù)嚴重依賴于電磁波所傳播的介質頻率特性，對純凈的極化介質(水)，通過極化消除介質頻率特性影響，消除影響后的Debye極化介電常數(shù)表達式為[10]：

(12)

圖1 由Debye確定的介電常數(shù)的實部和虛部

由圖1可見，在頻率較低或較高的時候，實部和虛部分別趨于常數(shù)和0。在隧道的超前預報工作中，天線的主頻一般較小(20～200 MHz)，相對介電常數(shù)的虛部趨于0，實部趨于20。因此，式(10)可以簡化為：

(13)

由G.C.Topp，等[11]的模型實驗和現(xiàn)場測試工作可知，介質含水量與介電常數(shù)間有如式(14)經驗關系：

θ=-0.053+0.029 2ε-5.5×10-4ε2+4.3×10-6ε3

(14)

另外，根據(jù)電磁波在介質中的傳播和反射特點，其在介質中的波速可表示為：

(15)

式中：x為異常體距離天線的水平距離；d為異常體的埋置深度；t為電磁波在介質中的傳播時間。

當天線位于異常體正上方時，式(15)可簡化為：

(16)

對于目前應用的大多數(shù)地質雷達信號分析軟件[12]，通過測試信號的分析可以求得υ，將其帶入式(13)，可求得異常體的相對介電常數(shù)，再利用式(14)可推斷異常體的含水量，進而估算前方含水體的水量大小。

3 工程實例應用

3.1 工程概況

成渝高速公路復線云霧山隧道全長3 355m，橫穿揚子準臺地重慶陷褶束華鎣山穹褶束瀝鼻峽背斜，有3條斷層通過隧道洞身。地貌單元屬于構造剝蝕、溶蝕顯著的中低山區(qū)，表現(xiàn)為新生代巖溶地貌特征。區(qū)域年平均降水量為1 000～1 200mm，大氣降水通過斷裂構造和裂隙滲入補給地下水。掌子面在接近雙圣寺斷層(F5)時，巖體節(jié)理裂隙逐漸密集，掌子面滲水較嚴重，溶蝕發(fā)育逐漸增強，在ZK20+350位置，采用中心頻率為40MHz的天線，進行了地質雷達探測。

3.2 不良地質體的圈定

圖2是雷達信號經過基本信號分析處理后(解震蕩→時間零點校正→濾波→時間增益→二維成像)的剖面。在圖2(a)中，圈定位置在圖像上為深色區(qū)，表明該處反射波能量強，其余淺色位置為正常信號反射情況；圖2(b)為異常位置處的單道反射波?？梢娫趫D2(a)中圈定的異常位置處雷達反射信號在不同介質的界面上分別出現(xiàn)兩處強烈反射。同時，在該兩處的反射波相位翻轉了180o。

圖2 雷達剖面及其異常位置的單道反射波

圖3為圖2(b)處某道反射信號的頻譜圖，從圖3中可以看出在異常位置處主頻約為20 MHz，在其他巖層均勻處于天線主頻接近，約為43 MHz，可見異常位置主頻明顯降低，電磁波的高頻成分出現(xiàn)顯著損耗。

圖3 異常位置雷達信號頻譜

在隧道開挖過程中積累的掌子面地質發(fā)展規(guī)律和勘察設計地質資料分析的基礎上，可以判斷隧道巖溶整體處于新生代巖溶發(fā)育期，該里程掌子面處于雙圣寺斷層(F5)的影響區(qū)內，且為巖溶較發(fā)育的區(qū)域。結合雷達信號2維剖面和不良體位置的單道反射波的振幅、相位、主頻變化特征，分析推斷在掌子面前方8～13 m范圍為一主要巖溶水含水層構造，在圈定位置有一含水量豐富的巖溶裂隙(溶縫)水通道。

3.3 含水量的估算

為估算該里程掌子面前方地層含水量的大小，根據(jù)40 MHz雷達天線信號在探測范圍類的傳播時間和剖面圖上圈定的不良地質體的厚度進行了兩個剖面點的地層波速的計算(如圖4)，并利用式(14)估算了掌子面前方地層的含水量大小變化情況。計算過程為：①將探測的30 m范圍按2 m長度進行等分，取測試系統(tǒng)讀波速的平均值；②將每段的平均波速代入式(13)，其中c=3×105km/s，可計算出分段區(qū)間內的相對介電常數(shù)的實部ε′；③因天線主頻較小，可取ε=ε′，將ε′代入式(14)求得分段內的含水量平均值。對1線的估算結果表明，在通過雷達波反射信號特征和地質資料綜合分析基礎上推斷的掌子面前方不良地質體處(圈定范圍)地層介質含水量高達約65%，2線的含水量平均約為4%～10%，與推斷的巖溶裂隙水情況較為吻合。

圖4 基于介電常數(shù)的地層含水量估算

通過開挖驗證，在剖面成像圖中圈定的巖溶裂隙水發(fā)育區(qū)域內巖體節(jié)理裂隙也較發(fā)育。在地質雷達測試掌子面里程前方約8 m，圖4中1線剖面處(掌子面右側，距中心線約3 m)存在一主要的溶縫為巖溶水通道〔圖5(a)〕，且?guī)r溶發(fā)育，初期涌水量約為139 L/min；2線剖面附近的巖層有少量裂隙水滲出，如圖5(b)。驗證了綜合分析推斷的結果，對隧道施工起到了指導作用。

圖5 開挖后揭露的掌子面溶蝕發(fā)育及裂隙水溶縫

4 結 論

1)巖溶地區(qū)隧道地質的復雜性，加深了地質雷達信號的多解性難題。在充分認識探測目標的地質發(fā)展規(guī)律，對揭露地層與探測體間的關系進行正確地質認識的基礎上，有機結合地質雷達信號對不同不良地質體的反射特征，對提高隧道巖溶裂隙水區(qū)域預報的準確性有重要作用。

2)采用介電常數(shù)與地層含水量的經驗關系的巖溶裂隙水量估算方法在超前地質預報中是可行的，且能滿足隧道施工中對施工風險的精度要求。地層含水體的空間大小在剖面圖上圈定的不準確是估算含水體水量的主要誤差來源。

3)新生代巖溶地區(qū)的地下水一般較發(fā)育，地層含水量的大小和圍巖裂隙的發(fā)育程度有著緊密的聯(lián)系。充滿水的裂隙通常隱埋在含水量較小的巖體中。巖溶裂隙水的探測結果可作為進一步評估圍巖完整性的依據(jù)。

[1] 王夢恕.對巖溶地區(qū)隧道施工水文地質超前預報的意見[J].鐵道勘查,2004,30(1):7-9. Wang Mengshu.Hydrologic and geological forecast of tunnel construction in the karst district[J].Railway Investigation and Surveying,2004,30(1):7-9.

[2] 鐘世航,孫宏志,王榮,等.隧道掌子面前方地質預報的進展[J].隧道建設,2007,27(4):7-11. Zhong Shihang,Sun Hongzhi,Wang Rong,et al.Progress of geology forecasting ahead of tunnel face by means of geophysical technology[J].Tunnel Construction,2007,27(1):7-9.

[3] 羅利銳,劉志剛.巖溶地區(qū)超前地質預報方法對比分析[J].巖土工程學報,2011,33(增刊1):351-355. Luo Lirui,Liu Zhigang.Comparative analysis of geological prediction methods in karst areas[J].Geotechnical Engineering,2011,33(Sup1):351-355.

[4] 李術才,李樹忱 ,張慶松,等.巖溶裂隙水與不良地質情況超前預報研究[J].巖石力學與工程學報,2007,26(2):217-215. Li Shucai,Li Shuchen,Zhang Qingsong,et al.Forecast of karst fractured ground water and defective geological conditions [J].Chinese Journal of Rock Mechanics and Engineering,2007,26(2):217-215.

[5] 鐘世航,孫宏志,李術才,等.隧道及地下工程施工中巖溶裂隙水及斷層、溶洞等隱患的探查、預報[J].巖石力學與工程學報,2012,31(增刊1):3298-3327. Zhong Shihang,Sun Hongzhi,Li Shucai,et al.Detection and forecasting for hidden danger of karst fissure water and other geological disasters during construction of tunnels and underground projects [J].Chinese Journal of Rock Mechanics and Engineering,2012,31(Sup1):3298-3327.

[6] 陳宇,溫忠輝,束龍倉.基巖裂隙水研究現(xiàn)狀與展望[J].水電能源科學,2010,28(4):62-65. Chen Yu,Wen Zhonghui,Shu Longcang.Status and prospect of research on bedrock fissure[J].Water Resources and Power,2010,28 (4):62-65.

[7] 徐開禮,朱志澄.構造地質學[M].北京:地質出版社,1987. Xu Kaili,Zhu Zhicheng.Structural Geology[M].Beijing:Geological Publishing House,1987.

[8] Harry M J.探地雷達理論與應用[M].北京:電子工業(yè)出版社,2011. Harry M J.Ground Penetrating Radar:Theory and Applications [M].Beijing:Publishing House of Electronics Industry,2011.

[9] Arkhipov V I.Hierarchy of dielectric-relaxation times in water [J].Journal of Non-Crystalline Solids,2002,305(1/2/3):127-135.

[10] Annan A P.Ground Penetrating Radar:Workshop Notes[M].Ontario,Canada:Sensors & Software Inc.,1999.

[11] Topp G C,Yanuka M,Zebchuk W D,et al.Determination of electrical conductivity using time domain reflectometry:soil and water experiments in coaxial lines[J].Water Resources Research,1988,24(7):945-952.

[12] 左昌群,呂喬森,黃彬.隧道施工地質預報方法模糊優(yōu)選與應用[J].重慶交通大學學報:自然科學版,2012,31(5):948-952. Zuo Changqun,Lv Qiaosen,Huang Bin.Fuzzy optimization and its application of geologic prediction methods in tunnel construction[J].Journal of Chongqing Jiaotong University:Natural Science,2012,31(5):948-952.

GPR for Detecting Karst-Fissure Aquifer in Tunnel

Zhao Ningyu1, Pan Jinqiu1, Yu Haijun2, Liang Bo1

(1. School of Civil Engineering & Architecture, Chongqing Jiaotong University, Chongqing 400074, China; 2. Chongqing Aerospace Polytechnic, Chongqing 400021, China)

The geological development principles of karst-fissure aquifer were analyzed. Then the characteristics of GPR electromagnetic reflection signal on karst-fissure aquifer and water content estimation method based on wave velocity-dielectric constant relationship were studied. Through specific engineering application, the geological radar has good resolution to fracture water buried in low water content rock and in fractured rock environments. And the estimation method of Karst fissure water moisture is feasible. The validity of the signal analysis method based on the geological understanding of the law to delineate adverse geological of karst fissure water and to estimate water content was verified.

tunnel engineering; karst-fissure aquifer; GPR; geological principle; geological detection; water content

10.3969/j.issn.1674-0696.2015.03.07

2014-08-12；

2014-12-30

國家自然科學基金項目(51178490)；重慶市自然科學基金項目(cstc2013jcyjA30019)；重慶市教育科學技術研究項目(KJ1400303)

趙寧雨(1981—)，男，四川南充人，講師，博士研究生，主要從事巖土工程試驗方面的研究。E-mail：zny2008@163.com。

U452.1+1；TU94+3.1

A

1674-0696(2015)03-032-04