姚顯春， 郭炳煊， 呂　高， 李　寧， 王小佳

(1.西安理工大學(xué) 巖土工程研究所，陜西 西安 710048；2.中石化勝利石油工程有限公司測(cè)井公司，山東 東營(yíng) 257096)

?

地質(zhì)雷達(dá)數(shù)據(jù)介電常數(shù)遞推反演研究

姚顯春1， 郭炳煊1， 呂高1， 李寧1， 王小佳2

(1.西安理工大學(xué) 巖土工程研究所，陜西 西安 710048；2.中石化勝利石油工程有限公司測(cè)井公司，山東 東營(yíng) 257096)

地質(zhì)雷達(dá)進(jìn)行地質(zhì)超前預(yù)報(bào)時(shí)，準(zhǔn)確讀取復(fù)雜地質(zhì)環(huán)境下不均勻介質(zhì)的介電常數(shù)是最為關(guān)鍵的環(huán)節(jié)。本文應(yīng)用地質(zhì)雷達(dá)天線在不同參數(shù)設(shè)置下進(jìn)行采樣，對(duì)采集到的數(shù)據(jù)進(jìn)行分層研究，剖析對(duì)比各地層中電磁波傳播速度，并對(duì)各地層的介電常數(shù)分布散點(diǎn)進(jìn)行歸納。根據(jù)對(duì)地層開(kāi)挖過(guò)程中的地質(zhì)素描驗(yàn)證和巖性比對(duì)，確定了適合該區(qū)域的不均勻混合介質(zhì)的等效介電常數(shù)。研究結(jié)果表明：應(yīng)用現(xiàn)場(chǎng)試驗(yàn)歸納的區(qū)域等效介電常數(shù)值不同于試驗(yàn)室得到的絕對(duì)介電常數(shù)值，前者分布在一定的范圍內(nèi)。與現(xiàn)場(chǎng)巖性的介電常數(shù)值對(duì)比結(jié)果表明，本文方法可準(zhǔn)確預(yù)測(cè)掌子面前方的巖性。

地質(zhì)雷達(dá)； 超前預(yù)報(bào)； 混合介質(zhì)； 介電常數(shù)

地質(zhì)雷達(dá)(Ground Penetrating Rader,簡(jiǎn)稱(chēng)GPR)探測(cè)技術(shù)具有快速、無(wú)損、探測(cè)精度高等優(yōu)點(diǎn)，近年來(lái)在地質(zhì)勘察和工程質(zhì)量檢測(cè)兩大領(lǐng)域的應(yīng)用逐漸廣泛[1-2]。應(yīng)用地質(zhì)雷達(dá)進(jìn)行隧道工程的超前預(yù)報(bào)能夠提前發(fā)現(xiàn)掌子面前方的地質(zhì)變化以及不良地質(zhì)現(xiàn)象，可為及時(shí)調(diào)整施工方法提供技術(shù)基礎(chǔ)，為減少工程事故提供快速準(zhǔn)確的預(yù)警報(bào)告[3]。然而地質(zhì)雷達(dá)在檢測(cè)過(guò)程中，由于地質(zhì)結(jié)構(gòu)復(fù)雜，容易引起誤判而存在多解性[4]。

隧道中巖土體材料介電常數(shù)的差異是地質(zhì)雷達(dá)進(jìn)行超前預(yù)報(bào)的先決條件，雷達(dá)接收到的反射波又是介電常數(shù)的函數(shù)，雷達(dá)圖像的解釋、判讀和反演都依賴(lài)于介質(zhì)的介電性能[5-6]。因此，巖土體介電特性分析是雷達(dá)發(fā)展應(yīng)用的基礎(chǔ)，也是從雷達(dá)圖像數(shù)據(jù)中獲得有用信息的關(guān)鍵。各種單相物質(zhì)的介電常數(shù)可由試驗(yàn)方便地確定。然而，對(duì)于多相混合物來(lái)說(shuō)，介電常數(shù)的確定比較復(fù)雜。因?yàn)槎喾N物質(zhì)混合在一起后，除了有物理方面的變化外，有時(shí)還會(huì)有化學(xué)方面的變化[7-8]，還存在多種因素必須要考慮。不能簡(jiǎn)單地認(rèn)為多相混合物的介電常數(shù)就是各組成物質(zhì)介電常數(shù)的比例組合。巖土體現(xiàn)場(chǎng)的介電常數(shù)值會(huì)因地層的含水量及孔隙率不同而發(fā)生變化[9-12]。為了有效判斷預(yù)測(cè)地層分布情況，需要確定出等效介電常數(shù)。

目前，國(guó)外對(duì)多相混合物的介電常數(shù)研究方面的成果也很突出。例如，Darold Wobschall通過(guò)半離散模型提出如何確定含水土壤中介電常數(shù)的理論，并成功描述如何在高頻范圍(1MHz～1GHz)確定含水土壤中的介電常數(shù)[13]。James R. Wang利用實(shí)驗(yàn)?zāi)Ｐ脱芯亢逝c介電常數(shù)之間的關(guān)系，并提出含水率對(duì)介電常數(shù)影響的合理范圍[14]。Alain Tabbagh等人通過(guò)數(shù)值模型反演介電常數(shù)與土壤含水率之間的物理意義[15]。Yongping Chen等人研究了在溫度和導(dǎo)電率變化的情況下，Maxwell-Wager(麥克斯-瓦格納)極化效應(yīng)對(duì)土壤中介電常數(shù)的影響[16]。

以上成果中，體現(xiàn)了不同混合介質(zhì)對(duì)介電常數(shù)的重要影響，但不能定量確定多相混合物中的介電常數(shù)。本文應(yīng)用地質(zhì)雷達(dá)天線在不同參數(shù)設(shè)置下進(jìn)行采樣，遞推反算得出電磁波速雙程走時(shí)的設(shè)置下電磁波在不同的深度中傳播的修正速度，并對(duì)各地層的介電常數(shù)分布散點(diǎn)進(jìn)行歸納，根據(jù)對(duì)地層開(kāi)挖過(guò)程中的地質(zhì)素描驗(yàn)證和巖性比對(duì)，確定適合該區(qū)域的不均勻混合介質(zhì)的等效介電常數(shù)。

1　檢測(cè)基本原理

1.1反射波幅法

檢測(cè)原理是利用反射波幅法，推算介電常數(shù)。反射波幅法基本依據(jù)是GPR天線的發(fā)射波在遇到多種(本文設(shè)定為2種)不同介質(zhì)的界面時(shí)，不僅波的傳播路徑會(huì)被改變，而且天線接收的反射能量也會(huì)變小，因?yàn)槠渲杏幸徊糠帜芰空凵涞较聦咏橘|(zhì)繼續(xù)傳播。經(jīng)過(guò)這樣不斷的循環(huán)反射，接收天線收到的電磁波就是在兩種介質(zhì)交界面上電磁波的疊加，因此振幅不同的波峰、波谷就是由不斷循環(huán)反射的信號(hào)疊加而成的，具體過(guò)程如圖1所示[17-19]。

圖1　波的反射示意圖Fig.1　The view of reflecting waves

圖1中，Δt1和Δt2可通過(guò)時(shí)間窗信號(hào)得知，A1、A2和A3表示各個(gè)反射波最大振幅。

一般把反射波能量與入射波能量的比稱(chēng)為反射系數(shù)R，把透射波能量與入射波能量的比叫做透射系數(shù)T，具體計(jì)算公式見(jiàn)式(1)、(2)[20]。

(1)

(2)

可以將式(1)改寫(xiě)成如下形式：

(3)

(4)

以上各式中，εr1表示上層介質(zhì)的介電常數(shù)；εr2是下一層介質(zhì)的介電常數(shù)；A表示反射波的波幅；Am表示全反射波的波幅。

(5)

(6)

式中R1是由上層介質(zhì)進(jìn)入下層介質(zhì)的反射系數(shù)，可得：

(7)

(8)

綜上所述，就可推出各層介質(zhì)的介電常數(shù)[21]。

1.2土的相對(duì)介電常數(shù)與含水率關(guān)系

土是由三相組成的：水、土顆粒和空氣。在施工現(xiàn)場(chǎng)，巖土體中存在很多種介質(zhì)，巖性往往并不單一，在含水率和孔隙率不均勻分布的情況下，巖土體的介電常數(shù)值會(huì)產(chǎn)生明顯的變化。這給精確測(cè)量和分析帶來(lái)了困難[22]。

地質(zhì)雷達(dá)主要利用的是電磁波在介質(zhì)中傳播的特性，由于電磁波在不同介質(zhì)中傳播的速度不同，因此可以對(duì)各種介質(zhì)加以區(qū)分。目前較為常用的經(jīng)驗(yàn)公式是：Topp公式，它是根據(jù)介電常數(shù)確定含水量，因此這也是對(duì)介電常數(shù)與含水率關(guān)系進(jìn)行對(duì)比的理論基礎(chǔ)。

(9)

(10)

式中εr是相對(duì)介電常數(shù)，θv是土體的體積含水率。

2　現(xiàn)場(chǎng)巖土介電參數(shù)標(biāo)定

對(duì)于同一掌子面、相同檢測(cè)深度的雷達(dá)數(shù)據(jù)，考慮到現(xiàn)場(chǎng)施工的擾動(dòng)很大，后期處理中存在一定的系統(tǒng)誤差，因此進(jìn)行了多次檢測(cè)，并對(duì)數(shù)據(jù)進(jìn)行歸納處理，保證了現(xiàn)場(chǎng)地質(zhì)介電常數(shù)的準(zhǔn)確性和真實(shí)性。雙程走時(shí)是指雷達(dá)發(fā)射時(shí)電磁波從發(fā)射到接收所用的時(shí)間。為分析各個(gè)圖像之間的差異，特意調(diào)試地質(zhì)雷達(dá)發(fā)射雙程時(shí)長(zhǎng)為：300 ns、340 ns、380 ns、420 ns、450 ns、480 ns、510 ns、540 ns、570 ns、590 ns、610 ns、630 ns、650 ns、670 ns、680 ns。檢測(cè)之后，分析各組雷達(dá)掃描圖以及波列圖，對(duì)比研究各個(gè)圖像之間的差異，并對(duì)雷達(dá)數(shù)據(jù)進(jìn)行交互式解釋?zhuān)謱友芯扛鱾€(gè)地層內(nèi)電磁波傳播速度以及相對(duì)介電常數(shù)，并對(duì)這些介電常數(shù)歸納分析得出一個(gè)趨勢(shì)范圍[23]。

2.1地質(zhì)雷達(dá)掃描圖對(duì)比

當(dāng)對(duì)巖性主要為粘土的隧道進(jìn)行超前預(yù)報(bào)試驗(yàn)時(shí)，其他儀器參數(shù)相同，雙程走時(shí)不同的條件下，檢測(cè)隧道掌子面，現(xiàn)場(chǎng)采集方法照片如圖2所示。

圖2　采集數(shù)據(jù)現(xiàn)場(chǎng)圖Fig.2　The image of data acquisition field

在后期的圖像處理中，因?yàn)楝F(xiàn)場(chǎng)檢測(cè)時(shí)環(huán)境惡劣，干擾很大，為了得到有效的電磁波，必須合理處理步驟，去除多余的干擾波。因此可利用RADAN5.0軟件，對(duì)檢測(cè)的圖片進(jìn)行同樣的自動(dòng)增益、去除背景、水平疊加處理，最終得出地質(zhì)雷達(dá)掃描圖，由于篇幅限制只給出300 ns、450 ns、680 ns的掃描圖，分別如圖3中的(a)、(b)及(c)所示。

圖3　不同雙程走時(shí)設(shè)置下雷達(dá)掃描圖Fig.3　Radar scan figures of different parameters in round-trip time

通過(guò)RADAN5.0軟件對(duì)圖3中的不同雙程走時(shí)掃描圖圖像進(jìn)行wiggle對(duì)比分析(見(jiàn)圖4(a)、(b)、(c))，可以看出隨著雙程走時(shí)的增加，多波列圖出現(xiàn)以下明顯區(qū)別。①隨著探測(cè)深度的增加，探測(cè)精度大幅度降低。如圖4(a)100～150 ns的波形圖，可以明顯預(yù)知前方圍巖結(jié)構(gòu)層松動(dòng)，但(b)(c)誤差范圍太大，不能較為準(zhǔn)確預(yù)測(cè)。②隨著雙程走時(shí)的增加，雖多列波信號(hào)數(shù)明顯增加，但準(zhǔn)確度提高幅度卻不明顯，增加了后期處理的工作量，得不償失。因此為了在保證準(zhǔn)確度條件下不失高效性，在設(shè)定雙程走時(shí)時(shí)必須反復(fù)調(diào)整找到合適值，最終得出最好的多波列圖。

圖4　不同雙程走時(shí)設(shè)置下雷達(dá)波列圖Fig.4　Radar wave columns of different parameters in round-trip time

2.2分層等效介電常數(shù)的遞推反算

通過(guò)RADAN5.0后處理軟件中的交互式解釋進(jìn)行等效介電常數(shù)的推算及結(jié)構(gòu)的分層，其中交互式解釋(interactive interpretation)是利用檢測(cè)數(shù)據(jù)進(jìn)行研究，然后根據(jù)研究結(jié)果分析出即將開(kāi)挖掌子面的圍巖巖性、含水量的一種推算方法。研究基于以下假設(shè)：①地層結(jié)構(gòu)是由不同介質(zhì)體所組成的，各層結(jié)構(gòu)都有明顯的電性差異；②電磁波在傳播過(guò)程中，理想地認(rèn)為傳播深度隨時(shí)間增加而增加。因此通過(guò)記錄雙程走時(shí)，就可以推算出地層結(jié)構(gòu)的深度。

其具體過(guò)程是：利用RADAN5.0后處理軟件中數(shù)據(jù)分析模塊子模塊—“交互式解釋”，將軟件數(shù)據(jù)導(dǎo)出為*.lay文件，然后進(jìn)一步計(jì)算分析以得出相關(guān)的電磁波波速和地層深度。

利用交互式解釋模塊分析掌子面前方結(jié)構(gòu)層時(shí)，通過(guò)電磁波信號(hào)在地下結(jié)構(gòu)層傳播反射的原理，導(dǎo)出電磁波在每個(gè)地層傳播的速度和雙程走時(shí)數(shù)據(jù)文件；通過(guò)對(duì)地層埋深的鉆孔深度確定，計(jì)算電磁波在各層的速度[24]，對(duì)應(yīng)于300 ns、450 ns及680 ns數(shù)據(jù)的相關(guān)參數(shù)分別如表1～3所示。

表1　雙程走時(shí)設(shè)置為300 ns時(shí)交互式解釋數(shù)據(jù)表

表2　雙程走時(shí)設(shè)置為450 ns時(shí)交互式解釋數(shù)據(jù)表

表3　雙程走時(shí)設(shè)置為680 ns時(shí)交互式解釋數(shù)據(jù)表

由表1～3可知，設(shè)置不同雙程走時(shí)，介電常數(shù)的推導(dǎo)值沒(méi)有發(fā)生明顯變化，說(shuō)明檢測(cè)數(shù)據(jù)受外界干擾接近，這就保證數(shù)據(jù)歸納時(shí)能集中，從而可根據(jù)介電常數(shù)的主要分布來(lái)等效該區(qū)域內(nèi)介電常數(shù)?；谶@一思路，本文將各個(gè)參數(shù)設(shè)置的介電常數(shù)值進(jìn)行合并整理，如圖5所示。

圖5　推導(dǎo)介電常數(shù)散點(diǎn)歸納圖Fig.5　Conclude of diffused electrical permittivity

由圖5可知，通過(guò)設(shè)置不同雙程走時(shí)，檢測(cè)以粘性土為主的掌子面時(shí)，介質(zhì)中等效介電常數(shù)并沒(méi)有發(fā)生很明顯的改變，主要分布范圍為4～7。另外圖5中出現(xiàn)了幾個(gè)較大介電常數(shù)值，說(shuō)明出現(xiàn)了較大偏差，這是由于采集數(shù)據(jù)時(shí)受到了嚴(yán)重干擾。但從整體來(lái)看，依然可以根據(jù)其余數(shù)據(jù)的集中程度判斷，該地質(zhì)區(qū)域內(nèi)的巖性等效介電常數(shù)主要為4～7；從而可以進(jìn)一步判斷出掌子面前方25 m處圍巖的巖性基本穩(wěn)定，沒(méi)有發(fā)生明顯變化，進(jìn)而保證了施工人員的安全。

3　介電常數(shù)的遞推算法應(yīng)用

利用某高速公路的淺埋隧洞作為研究對(duì)象，該隧道開(kāi)挖形式為雙洞分離式，其中左線長(zhǎng)1 243 m，最大埋深129 m；右線長(zhǎng)1 215 m，最大埋深120 m。左右兩線地質(zhì)情況較為相似，以灰?guī)r為主，勘察結(jié)果為Ⅳ級(jí)、Ⅴ級(jí)圍巖，并利用ZX 20+734和ZX 20+759這2個(gè)里程之間距離進(jìn)行實(shí)際應(yīng)用和驗(yàn)證。

3.1地層結(jié)構(gòu)推算分析

經(jīng)過(guò)現(xiàn)場(chǎng)檢測(cè)、后期數(shù)據(jù)處理，選取算例是雙程走時(shí)450 ns時(shí)的掃描圖，地質(zhì)雷達(dá)設(shè)置具體參數(shù)如表4所示。

表4　地質(zhì)雷達(dá)試驗(yàn)參數(shù)

預(yù)報(bào)范圍是ZX20+734～ZX20+759之間25 m，掃描圖及波列圖如圖6所示。圖6中(a)是相應(yīng)里程掌子面自左向右探測(cè)的雷達(dá)掃描圖，(b)為自動(dòng)增益處理、疊加處理后的單點(diǎn)波列圖，(c)為波列圖。

圖6　ZX 20+734～ZX 20+759掃描圖、單點(diǎn)波列圖及波列圖Fig.6　Images of radar scanning and oscillogram from ZX 20+734 to ZX 20+759

通過(guò)觀察圖6(c)中的A、B兩點(diǎn)，表示在掌子面前方7 m及20 m處出現(xiàn)了連續(xù)的反射波界面，反射界面明顯，反射面附近波幅顯著增強(qiáng)。由于反射波的同相軸連續(xù)性較好，可以判斷該反射波同相軸連線為破碎帶的位置。

表5　雷達(dá)波傳播走時(shí)與速度

基于地質(zhì)勘探方給出施工掌子面的地質(zhì)資料，即：主體為紅粘土及碎石土，掌子面頂部有滴水現(xiàn)象，結(jié)合檢測(cè)及后處理得到的圖6及表5數(shù)據(jù)可提前預(yù)測(cè)出以下幾點(diǎn)。

1) 距掌子面前方0～4 m之間的圍巖的巖性與掌子面的巖性基本相同，圍巖級(jí)別變化不大，介電常數(shù)基本為4.46。

2) 距掌子面前方4～8 m之間的圍巖介電常數(shù)為3.94，根據(jù)地質(zhì)勘探資料可知巖性基本為紅粘土，因同相軸連續(xù)性較好，預(yù)測(cè)其較之前0～4 m的圍巖結(jié)構(gòu)更為松散。

3) 距掌子面前方8～12 m之間的圍巖介電常數(shù)為10.63，較之前的介電常數(shù)有明顯的變化。由地質(zhì)資料可知，介質(zhì)不會(huì)有較大變化，因此預(yù)測(cè)是由于含水量增高引起的，在開(kāi)挖該范圍掌子面會(huì)有明顯的滴水現(xiàn)象。

4) 距掌子面前方12～20 m之間的圍巖介電常數(shù)基本保持不變。根據(jù)多次檢測(cè)與數(shù)據(jù)分析可知，可能是由于水對(duì)電磁波的吸收削弱所致，因此預(yù)測(cè)該范圍內(nèi)仍有滴水現(xiàn)象。

5) 距掌子面前方20～25 m內(nèi)介電常數(shù)為5，但是由于距離較遠(yuǎn)，傳播的電磁波信號(hào)較弱，保守預(yù)測(cè)該范圍內(nèi)依然含水，但含水率不高。

3.2地質(zhì)素描驗(yàn)證

工程隧道挖掘方式主要為挖掘機(jī)鉆進(jìn)，速度較慢且擾動(dòng)較小，有利于采集到現(xiàn)場(chǎng)圖片。隨著隧道的開(kāi)挖，從2011年12月22日Z(yǔ)X20+734到2012年1月7日Z(yǔ)X20+756，施工掘進(jìn)22 m，屬于檢測(cè)的有效預(yù)測(cè)樁號(hào)范圍之內(nèi)，通過(guò)現(xiàn)場(chǎng)照片(見(jiàn)圖7～10)，可以了解及判斷不同介電常數(shù)雷達(dá)預(yù)報(bào)的準(zhǔn)確性。

圖7　開(kāi)挖3 m處掌子面Fig.7　Excavation 3 m of the tunnel

圖7是開(kāi)挖3 m處掌子面的地質(zhì)現(xiàn)狀圖，從現(xiàn)場(chǎng)觀測(cè)可知：該掌子面的圍巖主要以紅粘土為主，土體緊實(shí)、堅(jiān)硬，但在表層并沒(méi)有明顯看到碎石土，也無(wú)滴水現(xiàn)象。

圖8　開(kāi)挖7 m處掌子面Fig.8　Excavation 7 m of the tunnel

圖8是開(kāi)挖到7 m處的掌子面的地質(zhì)現(xiàn)狀圖，從現(xiàn)場(chǎng)觀測(cè)可知：該掌子面的圍巖主要以紅粘土、碎石土為主，質(zhì)地較為松散，開(kāi)挖時(shí)較為容易，且紅粘土中夾雜著不少碎石土。

圖9　開(kāi)挖10 m處掌子面Fig.9　Excavation 10 m of the tunnel

圖9是開(kāi)挖到10 m處的掌子面的地質(zhì)現(xiàn)狀圖，從現(xiàn)場(chǎng)觀測(cè)可知：該掌子面的圍巖主要以紅粘土、碎石土為主，質(zhì)地松散，雖然施工條件差，但開(kāi)挖較易，開(kāi)挖過(guò)程中圍巖有明顯出水現(xiàn)象，如圖中所示是開(kāi)挖出含水量較大碎石土。

圖10　開(kāi)挖20 m處掌子面Fig.10　Excavation 20 m of the tunnel

圖10是開(kāi)挖到20 m處的掌子面的地質(zhì)現(xiàn)狀圖，從現(xiàn)場(chǎng)觀測(cè)可知：該掌子面的圍巖主要以紅粘土為主，質(zhì)地較硬，紅粘土中夾雜著不少碎石土，現(xiàn)場(chǎng)檢測(cè)時(shí)發(fā)現(xiàn)圍巖中明顯含水，且拱頂、拱腰處有明顯的出水痕跡。

綜上所述，根據(jù)逐步記錄的開(kāi)挖現(xiàn)場(chǎng)照片，可以看出在ZK20+741，即在距開(kāi)始掌子面后隧道開(kāi)挖為7 m處，圍巖地質(zhì)變化比較明顯，巖性變差，開(kāi)挖較易，參考圖6(c)中的A點(diǎn)處，驗(yàn)證了根據(jù)介電常數(shù)判斷較為準(zhǔn)確；另外在距離掌子面20 m處，詳見(jiàn)圖6(c)中的B點(diǎn)處，波列圖有明顯衰減，符合實(shí)際的圍巖含水這一變化情況，綜合對(duì)比，可見(jiàn)在該區(qū)域內(nèi)選取巖性的等效介電常數(shù)為4～7是比較準(zhǔn)確的。

4　結(jié)　論

通過(guò)遞推反算得出電磁波速在不同雙程走時(shí)、不同地層深度中傳播的修正速度，并對(duì)推算出的圍巖介電常數(shù)值分布散點(diǎn)進(jìn)行了歸納，旨在為利用地質(zhì)雷達(dá)數(shù)據(jù)來(lái)遞推反演介電常數(shù)提出一個(gè)支撐。

1) 應(yīng)用地質(zhì)雷達(dá)進(jìn)行隧道現(xiàn)場(chǎng)地質(zhì)等效介電常數(shù)歸納，可以較為準(zhǔn)確快速得出隧道掌子面巖土體的巖性。

2) 巖土體等效介電常數(shù)穩(wěn)定在一定范圍內(nèi)，與標(biāo)準(zhǔn)的介電常數(shù)值并不相等。

3) 為減小雷達(dá)探測(cè)過(guò)程中的干擾因素，如施工臺(tái)車(chē)、鉆孔設(shè)備及電焊機(jī)等的干擾，應(yīng)進(jìn)行多次布置測(cè)線進(jìn)行測(cè)量，才能較為準(zhǔn)確的同以往的等效介電常數(shù)值進(jìn)行比對(duì)。

4) 通過(guò)雷達(dá)檢測(cè)數(shù)據(jù)遞推得到的介電常數(shù)與現(xiàn)場(chǎng)地質(zhì)特征進(jìn)行綜合對(duì)比分析可以提高預(yù)報(bào)準(zhǔn)確度。但由于某些地質(zhì)差異無(wú)法準(zhǔn)確在雷達(dá)掃描圖中呈現(xiàn)，因此對(duì)地質(zhì)雷達(dá)波形圖解析時(shí)存在誤差，所以當(dāng)施工進(jìn)度增加時(shí)，需把預(yù)測(cè)結(jié)果和現(xiàn)場(chǎng)掌子面地質(zhì)相對(duì)比，做進(jìn)一步的修正。

[1] 李大心.探地雷達(dá)方法與應(yīng)用[M].北京:地質(zhì)出版社,1994.

[2]雷宛, 肖宏躍, 鄧一謙.工程與環(huán)境物探教程[M]. 北京:地質(zhì)出版社,2006.

[3]由廣明, 劉學(xué)增, 汪成兵. 地質(zhì)雷達(dá)在公路隧道超前地質(zhì)預(yù)報(bào)中的應(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.

[4]肖宏躍, 雷宛, 楊威. 地質(zhì)雷達(dá)特征圖像與典型地質(zhì)現(xiàn)象的對(duì)應(yīng)關(guān)系[J]. 煤田地質(zhì)與勘探,2008,36(4):57-61.

XIAO Hongyue, LEI Wan, YANG Wei. Correspondence between geological characteristics of radar images and typical geological phenomenon[J]. Coal Geology & Exploration, 2008,36(4):57-61.

[5]楊峰, 彭蘇萍.地質(zhì)雷達(dá)探測(cè)原理與方法研究[M].北京:科學(xué)出版社,2010.

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

[7]陳云敏, 邊學(xué)成, 陳仁朋,等.電磁波在三相土介質(zhì)中的傳播[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2003,24(6):611-618.

CHEN Yunmin, BIAN Xuecheng, CHEN Renpeng, et al. Propagation of electromagnetic wave in the three phases soil media[J]. Applied Mathematics and Mechanics, 2003,24(6):611-618.

[8]陳云敏, 梁志剛, 陳仁朋.電磁波在飽和無(wú)粘性土中的傳播特性及試驗(yàn)研究[J]. 土壤學(xué)報(bào), 2005,42(5):860-865.

CHEN Yunmin, LIANG Zhigang, CHEN Renpeng. Propagation properties of electromagnetic wave in saturated noncohesive soil and laboratory test[J]. Acta Pedologica Sinica, 2005,42(5):860-865.

[9]WALKER P D, Bell M R. Subsurface permittivity estimation from ground-penetrating radar measurements[C]// the IEEE 2000 International Radar Conference,2000.

[10] SPAGNOLINI U. Permittivity measurements of multilayered media with monostatic pulse radar[J]. IEEE Transactions on Geoscience and Remote Sensing, 1997,35(2):454-463.

[11]SCOTT W R, SMITH G S. Measured electrical constitutive parameters of soil as functions of frequency and moisture content[J]. IEEE Transactions on Geoscience and Remote Sensing, 1992,30(3):621-623.

[12]呂高, 李寧,朱才輝,等. 淺埋紅黏土隧洞圍巖含水率遷移研究及地層結(jié)構(gòu)反演分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 2014,33(增刊2):4041-4048.

Lü Gao, LI Ning, ZHU Caihui, et al. Study of moisture content migration of red clay and stratigraphic structure inversion analysis in shallow buried tunnel[J].Chinese Journal of Rock Mechanics and Engineering, 2014,33(S2):4041-4048.

[13]WOBSCHALL D. A theory of the complex dielectric permittivity of soil containing water: the semi-disperse model[J]. IEEE Transactions on Geoscience Electronics, 1977,15(1):49-58.

[14]WANG J R, SCHMUGGE T J. An empirical model for the complex dielectric permittivity of soils as a function of water conten[J]. IEEE Transactions on Geoscience and Remote Sensing, 1980,GE-18(4):288-295.

[15]TABBAGH A, CAMERLYNCK C, COSENZA P. Numerical modeling for investigating the physical meaning of the relationship between relative dielectric permittivity and water content of soils[J]. Water Resources Research, 2000,36(9):2771-2776.

[16]CHEN Yongping, DANI O. Effects of Maxwell-Wagner polarization on soil complex dielectric permittivity under variable temperature and electrical conductivity[J].Water Resources Research,2006,42(6):W06424.

[17]劉海,謝雄耀.上海盾構(gòu)隧道壁后土體介電常數(shù)測(cè)試實(shí)驗(yàn)[J]. 巖土工程學(xué)報(bào),2009,31(2):294-297.

LIU Hai, XIE Xiongyao. Dielectric constant tests on soils behind shield tunnel segment in Shanghai[J]. Chinese Journal of Geotechnical Engineering, 2009,31(2):294-297.

[18]劉敦文,古德生,徐國(guó)元,等.采空區(qū)充填物探地雷達(dá)識(shí)別技術(shù)研究及應(yīng)用[J]. 北京科技大學(xué)學(xué)報(bào),2005,27(1):13-16.

LIU Dunwen, GU Desheng, XU Guoyuan, et al. Exploration study of gob filling by ground penetrating radar and its application[J]. Journal of University of Science and Technology Beijing, 2005,27(1):13-16.

[19]陳權(quán),曾江源,李震,等.遙感監(jiān)測(cè)介電常數(shù)與土壤含水率的關(guān)系模型[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(12):171-175.

CHEN Quan, ZENG Jiangyuan, LI Zhen, et al. Relationship model of soil moisture and dielectric constant monitored with remote sensing[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(12): 171-175.

[20]程久龍,胡克峰,王玉和,等.探地雷達(dá)探測(cè)地下采空區(qū)的研究[J]. 巖土力學(xué),2004，25(增刊1):79-82.

CHENG Jiulong, HU Kefeng, WANG Yuhe, et al. Research on detecting of underground mined-out areas by using GPR[J]. Rock and Soil Mechanics, 2004，25(S1):79-82.

[21]李術(shù)才,薛翊國(guó),張慶松,等.高風(fēng)險(xiǎn)巖溶地區(qū)隧道施工地質(zhì)災(zāi)害綜合預(yù)報(bào)預(yù)警關(guān)鍵技術(shù)研究[J].巖石力學(xué)與工程學(xué)報(bào),2008,27(7):1297-1307.

LI Shucai, XUE Yiguo, ZHANG Qingsong, et al. Key technology study on comprehensive prediction and early-warning of geological hazards during tunnel construction in high risk karst areas[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(7):1297-1307.

[22]蘇立海,李寧,呂高,等.碎石泥土混合體介電常數(shù)研究及隧洞軟弱帶FDTD正演[J]. 巖石力學(xué)與工程學(xué)報(bào), 2015,34(6):1172-1178.

SU Lihai, LI Ning, Lü Gao, et al. Dielectric constant of gravel-soil mixture and FDTD forward analysis of weak zone in tunnel[J]. Chinese Journal of Rock Mechanics and Engineering, 2015,34(6):1172-1178.

[23]譚紹泉, 劉泰生, 徐錦璽,等. 地質(zhì)雷達(dá)技術(shù)在表層結(jié)構(gòu)調(diào)查中的應(yīng)用與研究[J]. 石油物探, 2003,42(1):59-62.

TAN Shaoquan, LIU Taisheng, XU Jinxi, et al. Application of ground-penetrating radar technology in investigation of surface structure[J]. Geophysical Prospecting for Petroleum, 2003,42(1):59-62.

[24]袁明德.探地雷達(dá)檢測(cè)中如何計(jì)算速度[J]. 物探與化探, 2003,27(3):220-222.

YUAN Mingde. The speed calculation in ground-penetrating radar detection[J]. Geophysical & Geochemical Exploration,2003,27(3):220-222.

(責(zé)任編輯王衛(wèi)勛，王緒迪)

Recursive inversion research on dielectric permittivity with data by ground penetrating rader

YAO Xianchun1， GUO Bingxuan1， Lü Gao1， LI Ning1， WANG Xiaojia2

(1.Institute of Geotechnical Engineering， Xi’an University of Technology， Xi’an 710048，China；2.Well Logging Company Sinopec Victory Petroleum Engineering Co., LTD， Dongying 257096，China)

When the ground penetrating radar is used to carry out the advance geological forecast, to accurately take the permittivity of non-homogeneous dielectric under the complex geological environment is the most important key link. In this paper the GPR antenna is employed to take samples under the different parameter setups, and the sampled data are studied and analyzed in strata to contrast the electromagnetic wave transmission velocity in each stratum, and the scattered distribution points of the dielectric permittivity in each stratum are summarized. And the equivalent dielectric permittivity of non-homogeneous mixed mediums suitable for this region in terms of geological sketch validation and lithological comparison in the process of stratum excavation. The research results indicate that the regional equivalent dielectric permittivity values finalized in the applied in-situ tests are different from those of absolute dielectric permittivity obtained from lab tests, whose values are distributed within a certain range. As contrasted with the in-situ lithological dielectric permittivity value, the results reveal that the method in this paper can accurately predict the lithology in front of tunnel face.

ground penetrating radar； advance geological forecast； mixed medium； permittivity

1006-4710(2016)02-0199-08

10.19322/j.cnki.issn.1006-4710.2016.02.012

2015-03-23

陜西省教育廳基金資助項(xiàng)目(15JK1540)；博士創(chuàng)新基金資助項(xiàng)目(207-002J1407)

姚顯春，男，講師，博士，研究方向?yàn)閹r體力學(xué)與工程。E-mail: yxc2049@163.com

U459.2

A