束 裕，程博文，于曉賀 ，李 沖，羅 蓉

(1.武漢理工大學(xué)，湖北 武漢 430063；2.湖北省公路工程技術(shù)研究中心，湖北 武漢 430063)

0 引言

自20世紀(jì)80年代起步以來，伴隨我國(guó)經(jīng)濟(jì)的增長(zhǎng)，公路干線網(wǎng)迅速發(fā)展升級(jí)。隨著公路總里程的增加，道路在表面和內(nèi)部都將出現(xiàn)不同形式的損傷，如表面裂縫、車轍、基層反射裂縫和空洞等病害。相比較來說，路表病害比路面內(nèi)部病害更容易被檢測(cè)出來，但路面內(nèi)部病害對(duì)路面造成損傷更大。路面含水空洞病害是路面早期水損害形式一種，其出現(xiàn)頻率非常高，具有普遍性、難以發(fā)現(xiàn)等特點(diǎn)，常用檢測(cè)儀器和評(píng)價(jià)指標(biāo)體系無法判斷其發(fā)展、影響程度。因此，需要一種有效檢測(cè)手段來判定其損傷程度。

探地雷達(dá)(Ground Penetrating Radar，簡(jiǎn)稱GPR)檢測(cè)是一項(xiàng)用于淺層地質(zhì)探測(cè)的新技術(shù)，又稱為地質(zhì)雷達(dá)，具有高效、簡(jiǎn)便、分辨率高、無損的特點(diǎn)。GPR圖像正演技術(shù)是一種通過計(jì)算機(jī)模擬電磁波在道路介質(zhì)中傳播過程的方法[1]。通過模擬電磁波在道路瀝青混合料介質(zhì)的傳播過程，可以提高對(duì)探地雷達(dá)探測(cè)剖面的認(rèn)識(shí)，提高實(shí)際雷達(dá)檢測(cè)數(shù)據(jù)解譯、分析結(jié)果的準(zhǔn)確性。

陶連金[2]等采用對(duì)離散元軟件對(duì)城市道路地下空洞病害進(jìn)行機(jī)理研究，從細(xì)觀的角度揭示了地下空洞發(fā)展至路面塌陷的破壞機(jī)制；張青山[3]通過現(xiàn)場(chǎng)雷達(dá)檢測(cè)與室內(nèi)試驗(yàn)結(jié)合分析瀝青路面水損害的狀況，并提出分級(jí)養(yǎng)護(hù)措施；劉俊[4]采用二維時(shí)域有限差分法對(duì)路面雷達(dá)電磁波進(jìn)行模擬，建立了路面脫空病害的二維正演模型；姜守東[5]等通過探地雷達(dá)對(duì)瀝青路面內(nèi)部的檢測(cè)效果進(jìn)行研究，發(fā)現(xiàn)探地雷達(dá)對(duì)材料的孔隙率和含水率反應(yīng)強(qiáng)烈，對(duì)層間積水有良好的檢測(cè)效果。目前研究均明確：路面中含水空洞病害對(duì)路面的影響較大，但由于路面含水空洞病害為路面內(nèi)部病害，要實(shí)現(xiàn)準(zhǔn)確測(cè)量其含水率比較困難。

經(jīng)研究發(fā)現(xiàn)，含水空洞中的含水率與雷達(dá)波形圖之間存在著定量關(guān)系。因此，本文通過分析實(shí)際雷達(dá)數(shù)據(jù)，建立不同含水率空洞模型，建立含水率與反射波形的數(shù)值關(guān)系，進(jìn)而通過峰值反向計(jì)算空洞含水率，為分析典型病害發(fā)展趨勢(shì)和高速公路預(yù)防性養(yǎng)護(hù)提供數(shù)據(jù)支撐，這對(duì)于精細(xì)化分析路面內(nèi)部含水空洞病害等瀝青路面早期水損害問題具有較大的實(shí)際工程意義。

1 理論推導(dǎo)

1.1 探地雷達(dá)在道路結(jié)構(gòu)中的傳播規(guī)律

在確定路面空洞含水率的研究中，構(gòu)建典型病害識(shí)別圖譜是第1步，而其基礎(chǔ)便是通過程序準(zhǔn)確模擬電磁波在三維介質(zhì)中的傳播過程。模擬電磁波在瀝青混合料面層的傳播過程則是通過有限元法求解基于電磁場(chǎng)參數(shù)的Maxwell方程組。

根據(jù)電磁波理論，高頻電磁波在介質(zhì)中的傳播規(guī)律服從Maxwell方程組[6]，該方程組用于描述磁場(chǎng)和電場(chǎng)相互激發(fā)、相互聯(lián)系形成電磁場(chǎng)中的關(guān)系，并結(jié)合介質(zhì)的本構(gòu)方程可得到式(1)。式(1)滿足電磁波在介質(zhì)傳播規(guī)律，被稱為電磁場(chǎng)的Helmholtz方程[7]。基于式(1)可知：①發(fā)射天線在瀝青混合料中形成的電磁場(chǎng)以波的形式運(yùn)動(dòng)；②電磁波是可以脫離波源獨(dú)立傳播的。

(1)

式中：E為電場(chǎng)強(qiáng)度；H為磁感應(yīng)強(qiáng)度；μ為磁導(dǎo)率，由于瀝青混合料是非磁性路面材料，因此瀝青混合料的磁導(dǎo)率μ為0；δ為電導(dǎo)率，其物理意義為瀝青混合料中的自由電荷在外加電場(chǎng)作用下產(chǎn)生瞬變電流的能力；ε為介電常數(shù)，表示為瀝青混合料在外部電場(chǎng)的作用下電能存儲(chǔ)、損耗的能力。

1.2 激勵(lì)源的設(shè)置

有限元模型的后續(xù)迭代計(jì)算需要選取合適的入射波與其數(shù)值變化規(guī)律。本文將選取1 GHz的雷克子波作為本次正演過程的激勵(lì)源。雷克子波被廣泛應(yīng)用于地表淺層探測(cè)工作中。在實(shí)際檢測(cè)過程中，初始場(chǎng)值始終為零，因此在雷達(dá)波數(shù)值模擬過程中，激勵(lì)源的初始值為零。為使雷克子波波形的初始值為零，通常在時(shí)域中保持其頻譜不變條件下，通過移動(dòng)雷克子波來實(shí)現(xiàn)。延時(shí)雷克子波的表達(dá)式如下式所示[8]：

(2)

式中：t0代表延遲時(shí)間。一般認(rèn)為t0=1/fm。

圖1為探地雷達(dá)天線的中心頻率為1GHz的雷克子波時(shí)域波形圖。

圖1 1 GHz雷克子波的時(shí)域波形圖

1.3 三維GPR波動(dòng)方程有限元求解

本文選取基于單行波算子的吸收邊界條件(ABC)作為本文正演的邊界條件。ABC吸收邊界是通過在網(wǎng)格剖分外部添加一層具有與邊界散射電磁波頻率入射角無關(guān)的波阻抗的媒質(zhì)。與傳統(tǒng)的分裂式PML吸收邊界條件相比，ABC邊界具有計(jì)算量小，算法穩(wěn)定和吸收效果好等優(yōu)點(diǎn)，可以提高三維波動(dòng)方程數(shù)值模擬的計(jì)算效率[9]。

本文將采用伽遼金法[10]推導(dǎo)電磁波在瀝青混合料中探地雷達(dá)有限元方程GPR波動(dòng)方程把場(chǎng)矢量隨空間、時(shí)間的變化聯(lián)系在一起，說明電磁場(chǎng)是以波動(dòng)方式在空間的傳播的。式(3)表示為GPR的波動(dòng)方程的余量。

(3)

(4)

式中：R為加權(quán)余量;v和v1為權(quán)函數(shù)；n為電場(chǎng)的切線方向向量，公式中其他含義與上式相同。將公式(4)中的右邊第1項(xiàng)展開為笛卡爾坐標(biāo)系的表示方法，并利用函數(shù)乘積的分部積分公式進(jìn)行化簡(jiǎn)得到式(5)：

(5)

將式(5)代回式(4)，選擇試函數(shù)v=-v1，得到式(6)，再通過正六面體將三維區(qū)域劃分成多個(gè)單元，得到單元矩陣方程。

(6)

假設(shè)網(wǎng)格區(qū)域有N個(gè)節(jié)點(diǎn)，采用六面體單元將計(jì)算區(qū)域進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分，設(shè)置插值函數(shù)后，進(jìn)行化簡(jiǎn)，得到求解方程，再將方程進(jìn)行擴(kuò)展至整體區(qū)域、所有節(jié)點(diǎn)，得到式(7)：

(7)

將式(7)采用中心差分法進(jìn)行求解，得到可以依次遞推電磁任意時(shí)間的電場(chǎng)值的式(8)：

(8)

2 數(shù)據(jù)采集

路面結(jié)構(gòu)層的厚度參數(shù)和混合料的相對(duì)介電常數(shù)是進(jìn)行計(jì)算的重要參數(shù)。為了使構(gòu)建的模型符合路面結(jié)構(gòu)真實(shí)情況，保證驗(yàn)證計(jì)算模型的真實(shí)性，因此采用真實(shí)路面結(jié)構(gòu)層的厚度參數(shù)構(gòu)建模型，采用瑞利(Rayleigh)模型來計(jì)算相對(duì)介電常數(shù)。

秦明月正要猶豫著要不要告訴她所敬愛的毛德君已經(jīng)死掉的事實(shí)時(shí)，正在檢查辦公室的趙大剛指著房間角落內(nèi)一個(gè)紙箱說：“這是什么？”

本文的數(shù)據(jù)采集地點(diǎn)為湖北省漢十高速公路孝襄段，主要進(jìn)行數(shù)據(jù)采集的3個(gè)部分：首先，通過探地雷達(dá)檢測(cè)設(shè)備對(duì)漢十高速公路孝襄段進(jìn)行雷達(dá)檢測(cè)，以獲取道路縱向雷達(dá)檢測(cè)二維剖面圖像和所有檢測(cè)點(diǎn)的單道反射波形數(shù)據(jù)數(shù)據(jù)；其次，本文對(duì)漢十高速公路孝襄段開展鉆芯取樣工作，測(cè)量面層芯樣結(jié)構(gòu)層厚度，以建立更加符合實(shí)際路況的道路三維模型；最后，由于相對(duì)介電常數(shù)測(cè)量的不穩(wěn)定性，本文需要通過計(jì)算模型進(jìn)行計(jì)算混合料的相對(duì)介電常數(shù)。由于介電常數(shù)數(shù)值較大，并且難以準(zhǔn)確得到，因此通常采用介電常數(shù)ε與真空介電常數(shù)ε0的比值表示，即相對(duì)介電常數(shù)ε1=ε/ε0。該參數(shù)是描述材料的極化特性的物理參數(shù)，即材料在外加電場(chǎng)的作用下，材料在電場(chǎng)中存儲(chǔ)電荷的能力。

2.1 路面雷達(dá)檢測(cè)

本次檢測(cè)使用的是由武漢理工大學(xué)和美國(guó)德克薩斯農(nóng)機(jī)大學(xué)聯(lián)合開發(fā)的WB1-21型路用探地雷達(dá)檢測(cè)設(shè)備，如圖2所示。本文設(shè)計(jì)的探地雷達(dá)檢測(cè)設(shè)備的參數(shù)如下：①工作頻率：1 GHz；②探測(cè)深度：0～60 cm；③豎向分辨率0.1 mm；④發(fā)射點(diǎn)間距：0.5 m。由于探地雷達(dá)檢測(cè)深度為60 cm，此次檢測(cè)目標(biāo)為道路面層和基層。并通過設(shè)備等距發(fā)射電磁波，接受道路反射的回波，以得到道路雷達(dá)檢測(cè)數(shù)據(jù)。

圖2 WB1-21型探地雷達(dá)檢測(cè)裝置

2.2 現(xiàn)場(chǎng)鉆芯取樣

路面取芯作為傳統(tǒng)的檢測(cè)手段，鉆取的芯樣不僅可以進(jìn)行病害判斷，還可以在路面使用中準(zhǔn)確的得到路面各結(jié)構(gòu)層的厚度，為后續(xù)建立道路三維模型奠定基礎(chǔ)[11]。由于本文鉆芯的目的是得到各結(jié)構(gòu)層的厚度，因此本次所取芯樣厚度為整個(gè)瀝青面層。為了不影響高速公路的正常運(yùn)營(yíng)，本次取芯選擇慢車道進(jìn)行取芯，取芯過程如圖3所示，一共取了8個(gè)芯樣(見圖4)，芯樣的厚度數(shù)據(jù)如表1所示。2.3相對(duì)介電常數(shù)的計(jì)算

圖3 漢十高速公路現(xiàn)場(chǎng)取芯

圖4 現(xiàn)場(chǎng)芯樣圖

采用儀器所測(cè)量的相對(duì)介電常數(shù)具有不穩(wěn)定性，因此采用計(jì)算的方法得到相對(duì)介電常數(shù)。本文采用的相對(duì)介電常數(shù)的計(jì)算模型為瑞利(Rayleigh)模型，該模型相對(duì)于線形模型、均方根模型等復(fù)合材料介電模型來說誤差較小，比較適用于計(jì)算路面各結(jié)構(gòu)層的相對(duì)介電常數(shù)[12]。

表1 芯樣厚度數(shù)據(jù)Table 1 Core sample thickness datacm上面層中面層下面層總厚度芯樣一3.64.77.916.2芯樣二3.44.38.215.9芯樣三3.84.68.016.4芯樣四3.65.17.916.6芯樣五3.65.58.217.3芯樣六3.16.17.516.7芯樣七3.666.616.2芯樣八3.86.28.218.2平均3.65.37.816.7

瀝青混合料是由空氣、瀝青結(jié)合料和骨料組成的三相復(fù)合介質(zhì)，其瑞利(Rayleigh)模型可表示為：

(9)

式中：θa、θas、θs分別為空氣、瀝青、集料的體積率；εac、εa、εas、εs分別為混合料試件、空氣、瀝青、集料的相對(duì)介電常數(shù)。

表2為常用路面使用材料的相對(duì)介電常數(shù)值，表3為漢十高速公路孝襄段所取的芯樣各結(jié)構(gòu)層的材料體積率的平均值，將其代入到式(9)中，可以得到各結(jié)構(gòu)層瀝青混合料的相對(duì)介電常數(shù)。

表2 常用路用材料介電常數(shù)范圍Table 2 Dielectric constant range of common circuit materials路面材料介電常數(shù)范圍空氣1水81瀝青5.1石灰?guī)r8.4

表3 各結(jié)構(gòu)層的平均材料體積率和相對(duì)介電常數(shù)Table 3 Mean material volume ratio and relative dielectric constant of each structure laye結(jié)構(gòu)層厚度/cm組成材料體積率/%空氣瀝青集料相對(duì)介電常數(shù)上面層(SUP-12.5)3.604.711.284.16.90中面層(AC-20C)5.304.49.286.47.03下面層(AC-25C)7.805.59.586.07.00

由分析可知：上面層的厚度為3.60 cm，相對(duì)介電常數(shù)平均值為6.90；下面層的厚度為5.30 cm，相對(duì)介電常數(shù)平均值為7.03，下面層的厚度為7.80 cm，相對(duì)介電常數(shù)平均值為7.00。

3 基于GPR正演技術(shù)的道路空洞含水率計(jì)算

3.1 網(wǎng)格劃分

根據(jù)前文各結(jié)構(gòu)層厚度、相對(duì)介電常數(shù)的結(jié)果，構(gòu)建真實(shí)道路三維模型，并將模型導(dǎo)入程序中計(jì)算含水率與雷達(dá)反射波形的數(shù)學(xué)關(guān)系。由于瀝青混合料試件由瀝青和集料拌和而成，因此將模型中的每一層的瀝青混合料視為定值。所建立的道路三維模型是一個(gè)10 m×0.22 m×1 m的長(zhǎng)方體，如圖5所示。上部為空氣層，厚度為3.30 cm，相對(duì)介電常數(shù)為1.00；第1層介質(zhì)為上面層，厚度為3.60 cm，相對(duì)介電常數(shù)為6.90；第2介質(zhì)為中面層，厚度為5.30 cm，相對(duì)介電常數(shù)7.03；第3層介質(zhì)為下面層，厚度為7.80 cm，相對(duì)介電常數(shù)為7.00。整個(gè)區(qū)域用單元1 cm×0.1 cm×1 cm的矩形單元進(jìn)行單元?jiǎng)澐?，共劃分? 000×220×100的網(wǎng)格空間。

圖5 真實(shí)道路三維模型

3.2 構(gòu)建模型

將上述真實(shí)路面三維模型導(dǎo)入MATLAB計(jì)算程序中，將數(shù)值結(jié)果圖像化處理，可得到單道反射波形圖，如圖6所示。

圖6 單道反射波形圖

圖中A1為空氣層與上面層產(chǎn)生的反射波幅；A2為上面層與中面層交界處產(chǎn)生的反射波幅；A3為中面層與下面層交界處產(chǎn)生的反射波幅。

由于雨水會(huì)從路表滲入和道路基層的水氣會(huì)向上擴(kuò)散，會(huì)導(dǎo)致在空洞處存在積水現(xiàn)象，且不易消除。由于水與空氣的相對(duì)介電常數(shù)有這較大的差異，導(dǎo)致雷達(dá)在發(fā)射高頻電磁波時(shí)會(huì)在兩者接觸面發(fā)生反射與折射現(xiàn)象[13]。因此，本文建立了不同含水率的三維空洞模型，來確定典型含水空洞雷達(dá)檢測(cè)特征圖像和單道反射波形圖，從而得到含水率與反射波幅間的模型。

3.3 不同含水率空洞病害正演模擬

通過查閱歷年工程養(yǎng)護(hù)情況和相關(guān)文獻(xiàn)，可以得到漢十高速公路孝襄段的空洞病害一般出現(xiàn)在層間交界處，其中下面層的空洞病害最嚴(yán)重。因此本文為了提高空洞識(shí)別能力，假設(shè)空洞形狀為5 cm×5 cm×2 cm的長(zhǎng)方體空洞，且位于道路三維模型的中面層與下面層的交界處，其中心點(diǎn)空間坐標(biāo)為(5 m，0.5 m，0.078 m)。同時(shí)選取了10%、20%、35%、50%、75%、100%的6種不同的含水率來建立三維空洞模型，如圖7所示。

(a)10%含水空洞

將模型導(dǎo)入到MATLAB中，并將模擬結(jié)果進(jìn)行圖像處理，可得到電磁波在空洞正上方發(fā)射、接收所得的中心單聲道反射波形圖，由于該模型反射波峰較多，但在400 ps至450 ps內(nèi)僅存在1個(gè)反射波峰，因此提取400 ps至450 ps間反射波谷進(jìn)行繪圖，如圖8所示。將反射波谷的峰值與含水率進(jìn)行線性擬合，如圖9所示。

圖8 反射波谷變化圖

圖9 空洞含水率與反射波谷擬合

根據(jù)含水率和該處的峰值可建立如下關(guān)系，即X表示含水率，Y表示該處峰值，單位為V，擬合公式為式(7)，相應(yīng)R2=0.992 2。該式可以反映反射波谷峰值隨空洞含水率變化的數(shù)學(xué)關(guān)系，可通過峰值反算空洞含水率的大小。

y=176.83x-121.6

(10)

3.4 實(shí)際應(yīng)用

根據(jù)現(xiàn)有的研究表明，存在含水空洞的雷達(dá)檢測(cè)二維剖面圖容易在中下面層反射界面出現(xiàn)“紅+藍(lán)+紅”特征圖像，藍(lán)色部分代表空洞中水的位置，紅色部分代表空洞中上部的空氣，通過“紅+藍(lán)”的位置準(zhǔn)確定位空洞位置，與此同時(shí)隨著這空洞含水率的增加，下部分紅色線條將不斷遠(yuǎn)離中下面層反射界面[14]。含水空洞病害的典型雷達(dá)檢測(cè)二維剖面圖見圖10。

圖10 含水空洞病害的典型雷達(dá)檢測(cè)二維剖面圖

基于這個(gè)典型特征，通過對(duì)照漢十高速公路實(shí)際雷達(dá)檢測(cè)二維剖面圖，發(fā)現(xiàn)漢十高速公路路面存在較多的含水空洞典型病害。圖11、圖12分別為漢十高速公路孝襄段實(shí)際含水空洞病害的雷達(dá)檢測(cè)二維剖面圖和反射波形圖。

圖11 實(shí)際雷達(dá)檢測(cè)二維剖面圖

圖12 雷達(dá)檢測(cè)反射波形圖

通過反射波形圖可以得到反射波谷峰值為-50.15 mV，代入之前所得到的反射波谷峰值與含水率的計(jì)算模型中，可以得到該處的含水空洞的含水率為：40.41%，其余含水空洞病害的圖像所在位置和含水率的大小如表4所示。

從表4可知，漢十高速公路孝襄段含水空洞病害共有6處，病害程度中等，特別是在K1156～K1158路段，如圖13所示。該路段可能出現(xiàn)連續(xù)病害，養(yǎng)護(hù)部門應(yīng)重點(diǎn)關(guān)注，及時(shí)養(yǎng)護(hù)，減少對(duì)路面造成的損傷。

表4 孝襄段隱藏含水空洞病害含水率Table 4 Moisture content of hidden water-bearing cavity disease in Xiaoxiang section樁號(hào)反射波谷峰值/mV含水率/%K1068+466-50.1540.41K1143+442-62.3733.50K1156+135-41.2345.45K1157+227-49.6240.71K1158+932-38.9446.75K1162+664-51.4739.66

(a) 雷達(dá)二維剖面圖

4 結(jié)語

本文通過對(duì)典型路面含水空洞模型的雷達(dá)正演模擬，得到典型含水空洞病害的雷達(dá)檢測(cè)反射波形圖，探究路面空洞病害含水率與反射波谷峰值之間的數(shù)學(xué)關(guān)聯(lián)。主要結(jié)論如下：

a.本文通過推導(dǎo)求解電磁場(chǎng)的Maxwell方程組，并設(shè)置激勵(lì)源和邊界條件，得到GPR三維波動(dòng)方程。基于該方程進(jìn)行編寫MATLAB代碼，并將路面結(jié)構(gòu)層厚度和混合料的相對(duì)介電常數(shù)導(dǎo)入，得到了真實(shí)道路三維模型。本正演模型是基于漢十高速公路實(shí)際路面數(shù)據(jù)建立的，但本文所提供的模型構(gòu)建方法仍適用于其他高速公路，推動(dòng)了探地雷達(dá)無損檢測(cè)技術(shù)在道路結(jié)構(gòu)病害的精細(xì)化檢測(cè)的應(yīng)用，為道路無損檢測(cè)技術(shù)的未來發(fā)展奠定了基礎(chǔ)。

b.本文發(fā)現(xiàn)在道路中下面層層間黏結(jié)處易含水空洞，因此在此處插入含水空洞模型，通過分析不同含水率的漢十空洞模型對(duì)電磁波在路面結(jié)構(gòu)傳播的影響，發(fā)現(xiàn)路面含水空洞含水率與反射波谷峰值呈線性關(guān)系，具體關(guān)系如式(10)所示，此數(shù)學(xué)關(guān)系式可以應(yīng)用于漢十高速公路路面空洞病害含水率的計(jì)算。這對(duì)其他高速公路準(zhǔn)確測(cè)量路面空洞含水率提供了參考依據(jù)，為路面雷達(dá)回波的正演模擬提供了一個(gè)有效的方法，也為進(jìn)一步開展路面雷達(dá)電磁波反演打下了理論基礎(chǔ)，具有極大的工程實(shí)際價(jià)值。