李 君

(湖南聯(lián)智科技股份有限公司，湖南 長沙 410200)

0 引言

隨著我國交通行業(yè)高速發(fā)展，運(yùn)營的公路隧道、鐵路隧道等總長度持續(xù)飛速增長。隧道襯砌大多數(shù)為混凝土結(jié)構(gòu)，具備對地質(zhì)條件適應(yīng)性強(qiáng)、易于按需要成型、整體性好、抗?jié)B性強(qiáng)，適用于多種施工條件等優(yōu)點。由于周邊圍巖結(jié)構(gòu)復(fù)雜、隱蔽工程較多，施工難度較大，在隧道建設(shè)和運(yùn)營期暴露出的質(zhì)量安全問題大多由襯砌病變誘發(fā)[1-3]。襯砌施工鋼筋密度不足、襯砌存在脫空病害未及時處理、混凝土離析等隱形缺陷將給隧道安全帶來隱患。所以，在隧道投入運(yùn)營前及運(yùn)營中，需開展隧道襯砌病害檢測，查明病害類型、規(guī)模和分布情況并及時處理，以保障隧道健康安全運(yùn)營。

目前，在隧道襯砌質(zhì)量檢測環(huán)節(jié)中應(yīng)用最廣的當(dāng)屬探地雷達(dá)技術(shù)，其高效率、高分辨率、連續(xù)檢測、直觀、無損等優(yōu)點是其他方法無法替代的。在隧道襯砌病害檢測方面，已有諸多學(xué)者進(jìn)行了研究。杜良等[4]基于時域有限差分法對隧道檢測中常見的鋼筋、鋼拱架、脫空、空洞等典型目標(biāo)體進(jìn)行了探地雷達(dá)正演模擬，并結(jié)合實測雷達(dá)圖像總結(jié)了典型目標(biāo)體雷達(dá)特征圖譜；許德根等[5]基于時域有限差分法和有限單元法，對隧道襯砌單層鋼筋網(wǎng)和無鋼筋襯砌背后的脫空進(jìn)行正演模擬，可提高實測圖譜解釋精度；徐浩等[6]基于時間域有限差分法對隧道襯砌中的脫空、離析、鋼筋、鋼筋下脫空和離析的波場特征進(jìn)行正演模擬，得出鋼筋下出現(xiàn)脫空和離析現(xiàn)象時，從雷達(dá)剖面圖中只能清楚地識別出鋼筋下的脫空現(xiàn)象；葉良應(yīng)等[7]在分析試驗室正演模擬雷達(dá)反射波波形圖的基礎(chǔ)上，進(jìn)一步進(jìn)行了2D切片分析和3D透視成像處理，查明了空洞的位置和尺寸，提高了圖像解釋的準(zhǔn)確度；鄧方進(jìn)[8]利用時域有限差分方法模擬了不同間距、不同深度、不同對應(yīng)位置的單層和雙層鋼筋混凝土結(jié)構(gòu)，并運(yùn)用F-K偏移技術(shù)分別進(jìn)行了成像處理，得出F-K 偏移可較好地壓制上層鋼筋多次波，有利于鋼筋的位置識別。

綜上所述，已有文獻(xiàn)主要側(cè)重于隧道襯砌存在的鋼筋、脫空、空洞等典型目標(biāo)體進(jìn)行單一正演模擬分析，關(guān)于鋼筋混凝土襯砌中鋼筋具體間距對病害識別影響的相關(guān)指標(biāo)研究較少，對比研究素混凝土與鋼筋混凝土襯砌下的脫空病害識別效果亦涉及較少，且對襯砌正演模型設(shè)置方面的研究與工程實例設(shè)計結(jié)合不夠緊密。因此，本文參照相關(guān)隧道襯砌質(zhì)量無損檢測規(guī)程和實際隧道施工設(shè)計圖，建立具有代表性的素混凝土和雙層鋼筋混凝土襯砌下的脫空病害模型，并進(jìn)行探地雷達(dá)正演模擬研究，總結(jié)和歸納出雷達(dá)圖譜特征規(guī)律，并緊密結(jié)合工程實例進(jìn)行應(yīng)用分析，以提高隧道襯砌病害檢測工作效率和解譯準(zhǔn)確度。

1 原理

1.1 探地雷達(dá)工作原理

探地雷達(dá)是利用高頻電磁波探測目標(biāo)體的一種地球物理勘探方法，其工作原理是發(fā)射天線向地下發(fā)射高頻電磁波，電磁波向前方傳播的過程中，遇到介電參數(shù)差異的目標(biāo)體時，會發(fā)生反射，由接收天線接收并記錄，再根據(jù)接收的雷達(dá)波形、電磁場強(qiáng)度、振幅、頻譜特征和雙程走時等參數(shù)來推斷目標(biāo)體的類型和分布情況。探地雷達(dá)檢測隧道襯砌的原理如圖1所示，圖中：h為脫空目標(biāo)體頂面埋深，x為兩天線之間的距離。

圖1 探地雷達(dá)工作原理

根據(jù)電磁波在介質(zhì)中傳播的雙程走時t和式(1)可求得目標(biāo)體頂面埋深h。

(1)

(2)

式中：V為電磁波在介質(zhì)中的傳播速度，C為電磁波在真空中的傳播速度，3×108m/s，為介質(zhì)的相對介電常數(shù)。介電常數(shù)εr通過參考常見介質(zhì)參數(shù)可得出，詳細(xì)參考值如表1所示。

表1 常見介質(zhì)介電常數(shù)參考值材質(zhì)介電常數(shù)εr導(dǎo)電率σ/(S·m-1)襯砌混凝土70.015脫空/空洞10鋼筋300109圍巖50.01

1.2 時域有限差分正演原理

時域有限差分法[9-13]是一種直接求解偏微分方程邊值問題的數(shù)值方法，其思想是將所研究的空間分成一定的網(wǎng)格，將時域麥克斯韋方程用有限差分方程組近似進(jìn)行時間離散化后加上初始條件和邊界條件，再按照時間步進(jìn)方法求解。該方法具有精度高、系統(tǒng)性強(qiáng)等特點，被廣泛應(yīng)用于公路、橋梁和隧道等工程檢測數(shù)值模擬理論研究。

探地雷達(dá)電磁波在介質(zhì)中傳播滿足麥克斯韋(Maxwell)方程，在二維介質(zhì)情況下麥克斯韋旋度方程組 TM模式方程可表示為：

(1)

(2)

(3)

式中：E為電場強(qiáng)度，V/m；H為磁場強(qiáng)度，A/m；μ為介質(zhì)磁導(dǎo)率，H /m；ε為介電系數(shù)；σ為介質(zhì)電導(dǎo)率，S /m；t為時間，s。

在模型邊界進(jìn)行迭代時，為模擬出無線半空間，采用 PML[14-17]邊界條件，而GprMax軟件正是基于 FDTD 算法和 PML 邊界吸收條件的探地雷達(dá)正演模擬開源代碼。所以本文應(yīng)用 GprMax 進(jìn)行探地雷達(dá)正演模擬研究，主要過程包括：首先新建一個空白的.IN格式文本文件，輸入需要植入的模型參數(shù)；然后運(yùn)行Python 腳本環(huán)境并讀取設(shè)置好的.IN文件，再利用GprMax軟件進(jìn)行數(shù)值模擬計算，對完成的所有單道波(A-scan)進(jìn)行數(shù)據(jù)融合輸出二維剖面(B-scan)；接下來再應(yīng)用MATLAB軟件進(jìn)行編程讀取B-scan數(shù)據(jù)體，并對讀取的二維GPR剖面進(jìn)行濾波、背景去噪、增益處理等；最后輸出預(yù)期的二維剖面圖。

2 二維正演模擬

采用 FDTD實現(xiàn)探地雷達(dá)正演模擬，首先要建立相應(yīng)的背景模型，然后植入信號源，接著對各個目標(biāo)體電磁場進(jìn)行賦值，最后進(jìn)行時間和空間數(shù)值模擬。

隧道襯砌正演模擬背景模型如圖2所示，設(shè)計了2層隧道斷面結(jié)構(gòu)，由上至下分別為襯砌和圍巖，斷面結(jié)構(gòu)模型長為 2.2 m、寬為1.1 m，包含了0.1m的 PML 邊界電磁波吸收寬度。上部襯砌混凝土厚度為 0.5 m，相對介電常數(shù)εr=7，電導(dǎo)率 σ=0.015 S/m；下部為相對介電常數(shù)εr=5，電導(dǎo)率σ=0.01 S/m，厚度為 0.5 m的圍巖。植入的信號源為500 MHz的ricker子波，道間距為0.01 m。

圖2 隧道襯砌正演模擬背景模型

2.1 素混凝土襯砌脫空病害正演模擬

素混凝土襯砌脫空病害正演模擬實驗結(jié)果如圖3所示。在隧道襯砌底部設(shè)計1個長方形脫空病害體，其頂面埋深為0.5 m，長為0.4 m，寬為0.1 m，脫空病害體內(nèi)部填充空氣，相對介電常數(shù)εr=1，如圖3(a)所示。

經(jīng)GprMax軟件模擬和MATLAB軟件處理后得到的探地雷達(dá)剖面如圖3(b)所示，從正演特征圖上可以得出：在脫空異常區(qū)頂面出現(xiàn)了明顯的強(qiáng)反射信號，振幅能量強(qiáng)，表現(xiàn)為雙曲線反射波組特征，多次波較發(fā)育，周圍存在明顯的繞射波。

(a)隧道襯砌脫空病害模型

(b)隧道襯砌脫空病害正演模擬圖圖3 素混凝土襯砌脫空病害正演模擬實驗結(jié)果

2.2 雙層鋼筋混凝土襯砌脫空病害正演模擬

雙層鋼筋混凝土襯砌脫空病害正演模擬實驗如圖4所示，第1層鋼筋中心距襯砌表面0.1 m，第2層鋼筋中心距襯砌表面0.3 m，鋼筋的間距分別為0.1、0.2、0.3 m[見圖4(a)～圖4(c)]，鋼筋相對介電常數(shù)εr=300；在襯砌底部分別設(shè)計1個長方形脫空病害體，其距離襯砌表面的頂面埋深為0.5 m，長為0.4 m，寬為0.1 m，脫空病害體內(nèi)部填充空氣，相對介電常數(shù)εr=1。

經(jīng)GPRMAX軟件正演模擬和MATLAB軟件處理后得到的探地雷達(dá)剖面[見圖4(d)～圖4(f)]，從正演特征圖上可以得出：第1層鋼筋的雙曲線反射能量強(qiáng)，分辨率高，均可清晰分辨；第2層鋼筋受到第1層鋼筋屏蔽干擾，隨著鋼筋間距的加密，分辨率逐漸下降。同理，隨著雙層鋼筋間距減小，信號被屏蔽越厲害，鋼筋間距為0.1 m時，幾乎不能分辨出其下方的脫空病害；鋼筋間距為0.2 m時，能分辨出其下方存在脫空病害，但不能準(zhǔn)確判斷其分布情況；鋼筋間距為0.3 m時，可以較清晰地分辨出脫空病害位置和分布情況，異常區(qū)頂面出現(xiàn)了明顯的強(qiáng)反射信號，振幅能量強(qiáng)，表現(xiàn)為似雙曲線反射波組特征，且多次波較發(fā)育，周圍存在較明顯的繞射波干擾。

圖4 雙層鋼筋混凝土襯砌脫空病害正演模擬實驗結(jié)果

3 工程應(yīng)用實例

3.1 素混凝土襯砌脫空病害典型案例

該隧道為貴州省境內(nèi)某公路隧道，工程地質(zhì)情況為：地貌類型為山區(qū)，地形起伏較大，圍巖以層狀構(gòu)造的石灰?guī)r、泥質(zhì)灰?guī)r等為主，夾層含水及泥質(zhì)物充填，溶洞、溶腔較發(fā)育，巖體大部較完整，未發(fā)現(xiàn)有大的斷裂地質(zhì)構(gòu)造經(jīng)過。隧道初期支護(hù)拱墻采用濕噴工藝噴射混凝土，二次襯砌采用素混凝土澆筑，初支與襯砌之間設(shè)置防水層。其中，V級圍巖襯砌設(shè)計厚度為 55 cm；Ⅳ級圍巖襯砌設(shè)計厚度為 50 cm；Ⅲ級圍巖襯砌設(shè)計厚度為 45 cm。

在隧道襯砌病害檢測過程中發(fā)現(xiàn)隧道襯砌的總體質(zhì)量較好，但存在一些質(zhì)量問題。圖5為隧道Ⅲ級圍巖段襯砌背后的異常反射雷達(dá)剖面，由圖可以得出：襯砌厚度整體約為45 cm，滿足設(shè)計要求，異常1為頂部雷達(dá)圖譜反射信號能量明顯增強(qiáng)，表現(xiàn)為似雙曲線反射波組特征，多次波較發(fā)育，兩端存在明顯的繞射波，根據(jù)二維正演模擬特征圖推斷該處異常區(qū)為脫空病害，產(chǎn)生此現(xiàn)象的初步原因分析可能是預(yù)埋件周邊振搗不密實或者防水板接縫處理不當(dāng)，有空氣沒有完全排出；異常2為頂部雷達(dá)圖譜反射信號能量明顯增強(qiáng)，表現(xiàn)為似雙曲線反射波組特征，且多次波發(fā)育，端部出現(xiàn)較明顯的繞射波，根據(jù)二維正演模擬特征圖推斷該處異常區(qū)為脫空病害，產(chǎn)生此現(xiàn)象的初步原因分析可能是在澆筑混凝土施工過程中振搗不到位形成的。

圖5 探地雷達(dá)實測剖面

3.2 雙層鋼筋混凝土襯砌脫空病害典型案例

該隧道為湖南省境內(nèi)的某公路隧道，工程地質(zhì)情況為：區(qū)域地貌類型為丘陵，地表剝蝕較強(qiáng)烈；圍巖以弱-中風(fēng)化花崗巖巖體為主，節(jié)理及裂隙較發(fā)育，裂隙間由鈣質(zhì)物、含水及泥質(zhì)物充填，巖體大部較完整，局部較破碎且未發(fā)現(xiàn)大的斷裂構(gòu)造經(jīng)過。隧道初期支護(hù)拱墻采用濕噴工藝噴射混凝土，二次襯砌采用雙層鋼筋混凝土澆筑，鋼筋間距為20 cm，初支與襯砌之間設(shè)置防水層。其中，Ⅵ級圍巖襯砌設(shè)計厚度為 65 cm；V級圍巖襯砌設(shè)計厚度為 60 cm；Ⅳ級圍巖襯砌設(shè)計厚度為 55 cm。

圖6為隧道V級圍巖段襯砌背后的異常反射雷達(dá)剖面，圖中可以得出：鋼筋表現(xiàn)出的小月牙形狀明顯，可較清晰地分辨雙層鋼筋，但鋼筋的屏蔽信號較強(qiáng)烈，下方出現(xiàn)了明顯的震蕩信號，給解譯帶來一定影響，初步解譯該襯砌整體厚度約為60cm，滿足設(shè)計要求。如圖6(a)所示，異常1為頂部雷達(dá)圖譜反射信號能量明顯增強(qiáng)，表現(xiàn)為似雙曲線反射波組特征，多次波較發(fā)育，端部存在繞射波，根據(jù)二維正演模擬特征圖推斷該處異常區(qū)為脫空病害，產(chǎn)生此現(xiàn)象的初步原因分析可能是澆筑混凝土施工過程中振搗不到位或者防水板接縫處理不當(dāng)，有少量的空氣未完全排出造成；如圖6(b)所示，異常2為頂部雷達(dá)圖譜反射信號能量明顯增強(qiáng)，表現(xiàn)為似雙曲線反射波組特征，且異常區(qū)反射波同相軸不連續(xù)，多次波發(fā)育，存在明顯的繞射波，根據(jù)二維正演模擬特征圖推斷該處異常區(qū)為脫空病害或已經(jīng)達(dá)到空洞的規(guī)模，分析產(chǎn)生此現(xiàn)象的初步原因可能是澆筑混凝土施工過程中振搗不到位所致。

(a)實測GPR異常1剖面

(b)實測GPR異常2剖面圖6 探地雷達(dá)實測剖面

4 結(jié)語

通過探地雷達(dá)正演模擬隧道襯砌病害檢測和工程實例應(yīng)用相結(jié)合的一系列研究，經(jīng)分析總結(jié)可得出以下結(jié)論：

1)理論正演數(shù)值模擬中，GPR對素混凝土襯砌內(nèi)部存在的脫空病害分辨率高，探測效果好，可較準(zhǔn)確圈定病害位置和分布情況；對雙層鋼筋混凝土襯砌內(nèi)部存在的脫空病害探測效果與鋼筋間距有關(guān)。鋼筋間距為10 cm時，鋼筋的震蕩干擾信號劇烈，幾乎不能識別出其下方的異常；鋼筋間距為20 cm時，可分辨出其下方的脫空異常，但對準(zhǔn)確解譯其位置和大小有一定影響；鋼筋間距為30 cm時，其下方的脫空異常信號清晰，可較準(zhǔn)確地解譯其下方的病害位置和分布范圍。

2)隧道襯砌實際工程檢測中，正演模擬特征圖可準(zhǔn)確指導(dǎo)現(xiàn)場襯砌脫空病害的識別。其中，在素混凝土條件下能夠很好地分辨脫空病害位置和大小，探測效果好；值得注意的是，襯砌內(nèi)鋼筋會給GPR解釋工作帶來干擾和難度，常見的雙層鋼筋混凝土襯砌鋼筋間距一般為20 cm左右，在理論正演模擬和工程應(yīng)用實例上均可以成功識別下方的脫空病害，但需要一定的數(shù)據(jù)處理、解譯能力以及豐富的現(xiàn)場經(jīng)驗。

3)現(xiàn)階段的隧道襯砌檢測，兩車道的隧道一般設(shè)計為5條探地雷達(dá)測線，包括拱頂1條，拱腰2條和墩墻2條，探地雷達(dá)二維線性剖面對整體的代表性較差，病害可能分布在兩條測線之間。建議今后選擇多通道天線進(jìn)行斷面全覆蓋探測，最終以三維剖面形式來解譯整個隧道的襯砌質(zhì)量。