地質(zhì)雷達波頻率與圍巖等級相關(guān)關(guān)系統(tǒng)計分析

黃 振，劉成禹

(福州大學環(huán)境與資源學院，福州 350108)

0 引言

探地雷達(Ground Penetrating Radar，簡稱 GPR)，又稱地質(zhì)雷達，是近些年發(fā)展起來的高效的淺層地球物理探測新技術(shù)，在隧道地質(zhì)超前預(yù)報中得到了廣泛的應(yīng)用［1－2］。

由于在巖溶地區(qū)隧道施工過程中，經(jīng)常遇到突水、突泥等無法預(yù)料的地質(zhì)災(zāi)害［3］，所以對超前預(yù)報研究的較多。王夢?。?］在對巖溶地區(qū)隧道施工水文地質(zhì)超前預(yù)報的意見中對儀器及方法進行了討論。方建勤等［5］關(guān)于巖溶地質(zhì)對隧道施工穩(wěn)定性的影響進行了分析評價。劉斌等［6］將復(fù)信號分析技術(shù)引入到地質(zhì)雷達探查水體的研究中。但是對于地質(zhì)雷達波頻率與圍巖性質(zhì)相關(guān)關(guān)系的研究較少。

由于地質(zhì)雷達探測存在多解性現(xiàn)象，積累豐富的工程經(jīng)驗是準確解釋雷達探測成果剖面的基礎(chǔ)［7］。場采集地質(zhì)雷達頻率數(shù)據(jù)，統(tǒng)計整體波形頻率變化特征及該區(qū)域內(nèi)不同圍巖情況下的雷達波頻率特征。結(jié)合現(xiàn)場地質(zhì)描述資料，分析了地質(zhì)雷達波頻率與圍巖性質(zhì)的相關(guān)特性，得到了一些定量指標，以便能更好地應(yīng)用于地質(zhì)雷達輔助隧道施工。

1 基本原理

現(xiàn)對掌子面前方圍巖進行超前探測，根據(jù)反射波的穿越時間、幅度和波形，判斷地下目標體的空間位置、結(jié)構(gòu)及其分布［2］;確定隧道前方斷層破碎帶、裂隙風化破碎帶、富水帶等不良地質(zhì)狀況的位置與范圍，進而詳細劃分或界定圍巖類型以便采取不同開發(fā)方式及超前支護措施，防止隧道災(zāi)害的發(fā)生［8－11］。

2 數(shù)據(jù)分析區(qū)域概況

2.1 區(qū)域地質(zhì)情況

贛龍鐵路3標段位于閩、粵、贛3省交界，以低山、丘陵地貌為主，地層巖性以沉積巖和變質(zhì)巖為主。在標段隧道內(nèi)，出露巖層巖性為花崗巖、石英砂巖、粉砂巖、變質(zhì)粉砂巖、變質(zhì)砂巖和炭質(zhì)頁巖;部分區(qū)域出露沖－殘積土和粉質(zhì)黏土;部分隧道地下水發(fā)育，地下水主要為基巖裂隙水和構(gòu)造裂隙水。

2.2 數(shù)據(jù)情況

所采用的數(shù)據(jù)取自2011年3—7月，贛龍鐵路3標地質(zhì)雷達探測報告。地質(zhì)雷達探測選用瑞典MALA地球科學公司(Sweden MALA Geoscience Inc)生產(chǎn)的RAMAC X3M型探地雷達，天線中心頻率100MHz。通過仔細閱讀探測報告，使用Reflexw軟件對原始數(shù)據(jù)重新處理，得到雷達剖面圖和頻譜圖。按照已劃分的不同圍巖段，統(tǒng)計得到了231個不同圍巖分段的雷達波頻率。

3 數(shù)據(jù)處理方法

為了對所得數(shù)據(jù)的總體分布及分布的數(shù)字特征等做出統(tǒng)計推斷，采用區(qū)間估計的數(shù)理統(tǒng)計方法，以得到可信度較高的數(shù)據(jù)區(qū)間。

假設(shè)統(tǒng)計點頻率服從正態(tài)分布N(μ，σ2)，計算統(tǒng)計點頻率的均值 μ 對應(yīng)于置信度為 0.95，0.98，0.99的置信區(qū)間，則計算公式為:

式中:n為計算參數(shù)樣本個數(shù);X為樣本平均值;S為樣本均方差;t為分位點。當統(tǒng)計點數(shù)≤6時，取均值的±4%作為區(qū)間界限。

正態(tài)測試是用于確定一組數(shù)據(jù)(Xi，i=1～n)是否服從正態(tài)分布的非常有用的工具。在正態(tài)測試中統(tǒng)計數(shù)

為了便于分析討論，按以下幾種情況進行數(shù)理統(tǒng)計:1)圍巖等級;2)巖性與圍巖等級;3)水文狀況與圍巖等級。

4 結(jié)果與分析

通過數(shù)理統(tǒng)計，得到3組不同置信度下的頻率分布區(qū)間結(jié)果。頻率分布區(qū)間不宜過大或過小，通過比較，認為選取區(qū)間范圍接近于10 MHz時作為所要求的結(jié)果較合理。

4.1 地質(zhì)雷達波頻率與圍巖等級

該部分共統(tǒng)計了231組數(shù)據(jù)，其中Ⅲ級圍巖15組，Ⅳ級圍巖63組，Ⅴ級圍巖153組。

對Ⅳ級、Ⅴ級圍巖統(tǒng)計數(shù)據(jù)作直方統(tǒng)計圖，并擬合正態(tài)分布曲線，如圖1所示。Ⅳ級、Ⅴ級圍巖統(tǒng)計數(shù)據(jù)較多，得到結(jié)果可信度更高，統(tǒng)計結(jié)果見表1。

圖1 圍巖等級直方統(tǒng)計圖Fig.1 Histograms showing relationship between GPR wave frequency and grade of rock mass

表1 圍巖等級計算結(jié)果Table 1 Calculation results based on grade of rock mass

在對地質(zhì)雷達原始數(shù)據(jù)的處理過程中得到不同圍巖等級的典型波形頻率圖，如圖2和圖3所示。

由圖1—3可知:

1)不同等級圍巖的雷達波頻率范圍有明顯差異，頻率分布范圍交叉較小，隨圍巖工程性質(zhì)變差，頻率明顯下降。

2)對于某一等級圍巖的波形頻率分布區(qū)間范圍較大，但有明顯集中分布的區(qū)間，通過擬合的正態(tài)曲線也能較好地反映其分布規(guī)律。

3)當統(tǒng)計數(shù)據(jù)足夠多時，對統(tǒng)計數(shù)據(jù)進行正態(tài)測試結(jié)果符合正態(tài)分布的假設(shè)，可見不同等級圍巖波形頻率變化為正態(tài)分布。

4.2 不同巖性地質(zhì)雷達波頻率與圍巖等級

該部分共統(tǒng)計了198組數(shù)據(jù)，按花崗巖39組、變質(zhì)粉砂巖75組、變質(zhì)砂巖33組、粉砂巖39組、石英砂巖12組統(tǒng)計，結(jié)果見表2。

對不同巖性統(tǒng)計結(jié)果作直方統(tǒng)計圖，對統(tǒng)計數(shù)據(jù)較多的圖形進行正態(tài)曲線擬合，以便更好地分析變化規(guī)律，如圖4所示。

由圖4可知:

1)粉砂巖等沉積巖的圍巖雷達波頻率較變質(zhì)粉砂巖等變質(zhì)巖的圍巖頻率高。花崗巖波形頻率變化較大，出現(xiàn)了較低的頻率，可能是由于整體性差、地下水發(fā)育較強引起的。

2)對于相同巖性的不同等級圍巖的變化規(guī)律與只考慮圍巖等級的情況大致相同。

表2 巖性與圍巖等級統(tǒng)計結(jié)果Table 2 Statistics on basis of lithology and grade of rock mass

4.3 不同水文狀況地質(zhì)雷達波頻率與圍巖等級

該部分共統(tǒng)計了186組數(shù)據(jù)，按地下水較發(fā)育78組、地下水不發(fā)育108組分別進行統(tǒng)計，結(jié)果見表3。

表3 水文狀況與圍巖等級統(tǒng)計結(jié)果Table 3 Statistics on basis of hydrology and grade of rock mass

對不同水文狀況統(tǒng)計結(jié)果作直方統(tǒng)計圖，并對圖形進行正態(tài)曲線擬合，以便更好地分析變化規(guī)律，如圖5所示。

圖5 Ⅴ級圍巖波形頻率直方統(tǒng)計圖Fig.5 Histograms showing relationship between groundwater and GPR wave frequency in GradeⅤrock mass

雷達波速和頻譜與巖體破碎程度相關(guān)，對巖體中的高電導物體，如含水情況、含泥情況影響也很大。

由圖5可以看出:對比不同水文狀況相同圍巖等級，波形頻率分布區(qū)間變化沒有明顯的規(guī)律性，擬合的正態(tài)分布曲線也基本相同，說明單一考慮地下水的發(fā)育情況，不能對圍巖等級做出很好的判斷。

5 結(jié)論與討論

雖然地質(zhì)描述及地質(zhì)預(yù)報存在很大的不確定性，但是這些現(xiàn)場資料的統(tǒng)計分析對區(qū)域地質(zhì)預(yù)報具有一定的輔助作用。

通過對大量地質(zhì)雷達隧道超前預(yù)報報告和實測掌子面產(chǎn)狀的統(tǒng)計，按照圍巖等級、巖性、水文情況3個方面分別進行分析，得出地質(zhì)雷達波頻率與圍巖等級相關(guān)關(guān)系的統(tǒng)計結(jié)果:1)不同等級圍巖的雷達波頻率范圍存在明顯差異，圍巖工程性質(zhì)變差，頻率變低;2)對于不同巖性圍巖，圍巖等級從Ⅲ級到Ⅴ級，雷達波頻率變低;3)不同水文狀況下，波形頻率分布區(qū)間變化不明顯。

由于區(qū)域性的地質(zhì)特點以及索取數(shù)據(jù)數(shù)量有限，所得規(guī)律性的統(tǒng)計特征具有局限性。雷達波速和頻譜對含泥、含水情況，破碎情況影響都較大，因此，需要更多的統(tǒng)計分析并結(jié)合其他方法(如地震波法、電法等手段)，才能將地質(zhì)預(yù)報的輔助作用發(fā)揮得更好。

［1］ 何發(fā)亮，李蒼松，陳成宗.隧道地質(zhì)超前預(yù)報［M］.成都:西南交通大學出版社，2006.

［2］ 薄會申.地質(zhì)雷達技術(shù)實用手冊［M］.北京:地質(zhì)出版社，2006.

［3］ 李術(shù)才，李樹忱，張慶松，等.巖溶裂隙水與不良地質(zhì)情況超前預(yù)報研究［J］.巖石力學與工程學報，2007，26(2):217 －225.(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 －225.(in Chinese))

［4］ 王夢恕.對巖溶地區(qū)隧道施工水文地質(zhì)超前預(yù)報的意見［J］.鐵路勘查，2004(1):7 －10.(WANG Mengshu.Hydrologic and geological forecast of tunnel construction in the karst district［J］.Railway Investigation and Surveying，2004(1):7－10.(in Chinese))

［5］ 方建勤，胡居義.巖溶地質(zhì)對隧道施工穩(wěn)定性的影響評價［J］.礦業(yè)研究與開發(fā)，2010(1):40－41.(FANG Jianqin，HU Juyi.Evaluation on the influence of karst geology on the excavation stability of tunnel［J］.Mining Research and Development，2010(1):40 －41.(in Chinese))

［6］ 劉斌，李術(shù)才，李樹忱，等.復(fù)信號分析技術(shù)在地質(zhì)雷達預(yù)報巖溶裂隙水中的應(yīng)用研究［J］.巖土力學，2009，30(7):2191 －2196.(LIU Bin，LI Shucai，LI Shuchen，et al.Study of application of complex signal analysis to predicting karst-fractured ground water with GPR［J］.Rock and Soil Mechanics，2009，30(7):2191 －2196.(in Chinese))

［7］ 楊天春，周勇，李好.超前探測中探地雷達應(yīng)用與結(jié)果的處理分析［J］.工程地質(zhì)學報，2010，18(6):971 －975.(YANG Tianchun，ZHOU Yong，LI Hao.Use and analysis of ground penetration radar for forcasting geological conditions in highway tunneling［J］.Journal of Engineering Geology，2010，18(6):971 －975.(in Chinese))

［8］ 白冰，周健.探地雷達測試技術(shù)發(fā)展概況及其應(yīng)用現(xiàn)狀［J］.巖石力學與工程學報，2001，20(4):527－531.(BAI Bing，ZHOU Jian.Advances and applications of ground penetrating radar measure technology［J］.Chinese Journal of Rock Mechanics and Engineering，2001，20(4):527 －531.(in Chinese))

［9］ 由廣明，劉學增，汪成兵.地質(zhì)雷達在公路隧道超前地質(zhì)預(yù)報中的應(yīng)用［J］.公路交通科技，2007，24(8):92 －95.(YOU Guangming，LIU Xuezeng，WANG Chengbing.Application of ground penetrating radar to geological forecast for expressway tunnel construction［J］.Journal of Highway and Transportation Research and Development，2007，24(8):92 －95.(in Chinese))

［10］ 崔芳，高永濤，吳順川.基于地質(zhì)雷達的隧道掌子面前方地質(zhì)預(yù)報研究［J］.中國礦業(yè)，2011，20(3):115－117.(CUI Fang，GAO Yongtao，WU Shunchuan.Research on the geological forecast of work face in tunnel with ground penetrating radar［J］.China Mining Magazine，2011，20(3):115 －117.(in Chinese))

［11］ 薛建，曾昭發(fā)，王者江，等.隧道掘進中掌子面前方巖石結(jié)構(gòu)的超前預(yù)報［J］.長春科技大學學報，2000，30(1):87－89.(XUE Jian，ZENG Zhaofa，WANG Zhejiang，et al.Prediction of rock structure in front of face during tunneling［J］.Journal of Changchun University of Science and Technology，2000，30(1):87 －89.(in Chinese))