黃嵐清 洪亮 張亮 肖育旗

摘? 要：針對(duì)隧道襯砌結(jié)構(gòu)鋼筋數(shù)量缺失與布置間距不合理等質(zhì)量問題，基于探地雷達(dá)正演模擬試驗(yàn)，對(duì)隧道襯砌雙層正對(duì)鋼筋進(jìn)行探測(cè)研究，分別采用F-K法和Kirchhoff法對(duì)鋼筋檢測(cè)圖像進(jìn)行偏移歸位處理，探究鋼筋回波信號(hào)偏移歸位效果，比較兩種偏移方法的優(yōu)劣。結(jié)果表明：理想無噪聲條件下，F(xiàn)-K法和Kirchhoff法能獲取相同的偏移速度，并使襯砌雙層鋼筋回波信號(hào)準(zhǔn)確歸位，但由F-K偏移法得到的偏移圖像信息熵更小，偏移效果更精確。含噪聲情況下，F(xiàn)-K偏移法具有更強(qiáng)的抗噪性，能獲取穩(wěn)定有效的偏移速度，偏移處理后鋼筋信號(hào)能量集中且易于識(shí)別。

關(guān)鍵詞：隧道襯砌；探地雷達(dá)；偏移成像；F-K偏移法；Kirchhoff偏移法

中圖分類號(hào)：TP39；TN957.52 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：2096-4706（2023）15-0189-06

Research on Migration Imaging of Double-layer Reinforcement in Tunnel Lining Based on F-K and Kirchhoff Methods

HUANG Lanqing， HONG Liang， ZHANG Liang， XIAO Yuqi

（School of Civil Engineering， Hunan City University， Yiyang? 413000， China）

Abstract： In response to quality issues such as missing quantity and unreasonable spacing of reinforcements in the tunnel lining structure， based on the forward modeling test of Ground Penetrating Radar， the detection of double-layer facing reinforcement in the tunnel lining is studied. The F-K method and Kirchhoff method are used to migrate and return the reinforcement detection images， and the effect of reinforcement echo signal migration and return is explored. The advantages and disadvantages of the two migration methods are compared. The results show that the F-K method and Kirchhoff method can obtain the same migration speed and accurately locate the echo signal of the double-layer reinforcement in the lining under the ideal noiseless conditions. However， the information entropy of the migration image obtained by the F-K migration method is smaller， and the migration effect is more accurate. In the presence of noise， the F-K migration method has stronger noise resistance and can obtain stable and effective migration speed. After migration processing， the reinforcement signal energy is concentrated and easy to identify.

Keywords： tunnel lining; Ground Penetrating Radar; migration imaging; F-K migration method; Kirchhoff migration method

0? 引? 言

在經(jīng)濟(jì)快速發(fā)展、科學(xué)實(shí)力和技術(shù)水平不斷提高的背景下，我國(guó)交通設(shè)施網(wǎng)在不斷完善，隧道工程建設(shè)也在快速發(fā)展，規(guī)模不斷擴(kuò)大，隧道襯砌結(jié)構(gòu)的質(zhì)量檢測(cè)也隨之越來越重要。隧道襯砌檢測(cè)的內(nèi)容主要包括襯砌結(jié)構(gòu)的厚度、襯砌密實(shí)度、襯砌結(jié)構(gòu)內(nèi)部鋼筋分布是否滿足設(shè)計(jì)要求等。探地雷達(dá)由于有操作簡(jiǎn)單、成果直觀、適應(yīng)性強(qiáng)、分辨率高且為無損探測(cè)的特點(diǎn)，而被廣泛應(yīng)用于隧道襯砌結(jié)構(gòu)鋼筋分布、質(zhì)量缺陷檢測(cè)中。采用探地雷達(dá)進(jìn)行隧道襯砌探測(cè)與狀態(tài)分析對(duì)襯砌施工質(zhì)量的把控有至關(guān)重要的作用。襯砌結(jié)構(gòu)內(nèi)鋼筋介電常數(shù)遠(yuǎn)高于混凝土介電常數(shù)，電磁波遇到鋼筋會(huì)發(fā)生強(qiáng)烈反射，多層鋼筋引起的反射和繞射問題復(fù)雜，導(dǎo)致檢測(cè)人員難以對(duì)圖像進(jìn)行準(zhǔn)確的解釋和判斷。需探索不同條件下適用于鋼筋數(shù)量檢測(cè)與間距分析的方法，對(duì)鋼筋進(jìn)行準(zhǔn)確的位置還原與識(shí)別，發(fā)現(xiàn)結(jié)構(gòu)中存在的鋼筋錯(cuò)位等問題，及時(shí)檢測(cè)襯砌結(jié)構(gòu)服役性能。

諸多學(xué)者對(duì)隧道襯砌鋼筋探地雷達(dá)檢測(cè)及其回波信號(hào)分析展開研究。陳偉等[1]采用探地雷達(dá)對(duì)隧道襯砌模型進(jìn)行正演模擬，深入分析了不同介質(zhì)對(duì)電磁波反射的影響，得出檢測(cè)結(jié)果中鋼筋分布等圖像特點(diǎn)；吳慶全等[2]在對(duì)工程實(shí)例的地質(zhì)雷達(dá)鋼筋檢測(cè)研究中發(fā)現(xiàn)回波信號(hào)受鋼筋埋深、直徑以及分布情況等因素影響；舒志樂等[3]對(duì)雷達(dá)數(shù)據(jù)進(jìn)行二維和三維的分析研究，發(fā)現(xiàn)鋼筋間距的不同會(huì)產(chǎn)生不同程度的干擾；尹濤等[4]通過FDTD數(shù)值模擬發(fā)現(xiàn)不同類型的襯砌缺陷的識(shí)別難度不同，雙層鋼筋網(wǎng)的屏蔽作用更強(qiáng)；楊宇等[5]通過對(duì)模擬實(shí)際混凝土結(jié)構(gòu)中不同埋深的鋼筋的檢測(cè)及信號(hào)處理，提高了混凝土內(nèi)部深淺埋鋼筋的可視化程度。探地雷達(dá)常用的偏移算法包括柯?；舴颍↘irchhoff）積分偏移法[6]、F-K偏移法[7，8]等。偏移處理的作用在于提高圖像的分辨率，從而恢復(fù)復(fù)雜結(jié)構(gòu)的真實(shí)狀態(tài)和準(zhǔn)確位置。黃玲[9]等對(duì)單層鋼筋及缺陷正演結(jié)果進(jìn)行F-K偏移處理，很好地把鋼筋網(wǎng)的繞射波聚集歸位，有效地改善了圖像的成像效果；林志強(qiáng)等[10]提出基于圖像熵的探地雷達(dá)Kirchhoff偏移成像算法，較準(zhǔn)確地估計(jì)出地下電磁波的傳播速度.

本文通過引入圖像熵的概念，動(dòng)態(tài)估計(jì)可使F-K法和Kirchhoff法偏移成像效果最好的波速參數(shù)，并采用F-K法和Kirchhoff法對(duì)理想無噪及存在強(qiáng)干擾噪聲的襯砌雙層鋼筋正演圖像進(jìn)行偏移成像處理，對(duì)兩種方法的鋼筋偏移歸位效果進(jìn)行比較和判斷，從而為實(shí)際工程中襯砌鋼筋的偏移歸位處理及準(zhǔn)確識(shí)別提供理論參考和依據(jù)。

1? 偏移方法和理論

1.1? kirchhoff偏移法

Kirchhoff偏移法在20世紀(jì)70年代被提出并逐步成熟，其始于繞射偏移原理并依據(jù)Kirchhoff積分公式來實(shí)現(xiàn)偏移成像，并借此發(fā)展為地震波動(dòng)方程積分法。該偏移方法通過對(duì)探測(cè)波波動(dòng)方程進(jìn)行Kirchhoff積分求解，將地面路徑上接收到的地下物體某一點(diǎn)處探測(cè)波繞射信號(hào)匯集起來形成偏移圖像，在實(shí)際工程的探測(cè)中也得到了一定應(yīng)用。

在介質(zhì)均勻的條件下，探地雷達(dá)發(fā)射電磁波的電場(chǎng)分量E（x，y，z，t）滿足標(biāo)量波動(dòng)方程：

式中，v為波速。

式（1）的Kirchhoff積分解為：

其中，S表示包圍點(diǎn)（x0，y0，z0）的任意光滑閉曲面，r表示點(diǎn)（x0，y0，z0）到曲面S上任意點(diǎn)的距離，k = ω / v稱為波數(shù)，ω表示電磁波的角頻率，? / ?n表示沿S的外法線方向的方向?qū)?shù)。

由于在探地雷達(dá)探測(cè)過程中，波場(chǎng)僅存在于地面以下的半空間中，所以選取閉曲面由無限大的地平面S0和一個(gè)半徑無窮大的半球面構(gòu)成。由電磁波的傳播規(guī)律可知，無窮大半球面上的波函數(shù)為零，因此對(duì)式（2）進(jìn)行推導(dǎo)可得：

其中，E - = E（x，y，z，t - r / v）。

式（3）說明P點(diǎn)t時(shí)刻的波場(chǎng)值是由前一時(shí)刻t - r / v地面上的場(chǎng)源激發(fā)的，這符合電磁波“向前”傳播的規(guī)律，而探地雷達(dá)記錄的數(shù)據(jù)是地下目標(biāo)反射至地面的波場(chǎng)函數(shù)值，探地雷達(dá)偏移成像的目的是利用記錄的數(shù)據(jù)反演出目標(biāo)在地下的真實(shí)位置，這要依靠波“倒退”的規(guī)律。事實(shí)證明，波“倒退”也符合惠更斯—菲涅爾原理，也可以用Kirchhoff偏移積分來描述，只是需要將時(shí)間“逆轉(zhuǎn)”，于是有：

其中，E + = E（x，y，z，t + r / v）。

利用探地雷達(dá)的反射成像原理，偏移速度取真實(shí)速度的1/2，并利用在地平面（z = 0）處記錄的數(shù)據(jù)在t = 0時(shí)進(jìn)行成像，可得：

Kirchhoff偏移法具有精確成像、快速計(jì)算、適應(yīng)性強(qiáng)、處理靈活、效率高等特點(diǎn)，能夠解決橫向速度劇烈變化帶來的偏移成像不精確的技術(shù)難點(diǎn)，極大地提高了回波成像的分辨率和信噪比，并且積分法能夠適應(yīng)變化的觀測(cè)系統(tǒng)和起伏的地表，因此積分法能夠適應(yīng)復(fù)雜結(jié)構(gòu)的偏移成像。

1.2? 頻率波數(shù)（F-K）偏移法

設(shè)二維信號(hào)或圖像f （x，z = 0，t）的傅里葉變換為F（kx，z = 0，w），則：

式中：x表示水平坐標(biāo)，z表示垂直坐標(biāo)，向下為正，t表示時(shí)間.在頻率波數(shù)域內(nèi)，深度z的波場(chǎng)可表示為：

F（kx，z，w）=F（lx，0，w）eik，z（7）

設(shè)二維信號(hào)f （x，z，t）{ f （x，t）在z方向的延拓}為F（kx，z，w）關(guān)于kx，w的二維傅里葉逆變換，則：

然后，根據(jù)頻散關(guān)系式：

將式（9）代入式（8），令t = 0，經(jīng)過推導(dǎo)可得：

根據(jù)式（10），可以對(duì)探地雷達(dá)圖像進(jìn)行頻率波數(shù)偏移成像分析。

2? 圖像信息熵評(píng)價(jià)

信息熵用以描述事件或集合的不確定性或無序性，圖像信息熵則反映了圖像平均信息量的大小。若被分析對(duì)象越是規(guī)則有序，所包含的有用信息量越大，對(duì)應(yīng)的信息熵值就越小；若被分析對(duì)象越雜亂無序，則所具備的有用信息量就越小，對(duì)應(yīng)的信息熵值就越大。信息熵測(cè)度是反映被分析對(duì)象規(guī)則性水平的一個(gè)重要度量。

在地質(zhì)雷達(dá)回波圖像熵值的度量中，當(dāng)偏移成像信息越集中，則圖像信息熵測(cè)度值越小，反之，若接收?qǐng)D像中多次波或繞射波越多，則圖像熵值越大。在偏移歸位處理中，恰當(dāng)?shù)钠扑俣仁瞧瞥上衲芊癯晒Φ年P(guān)鍵。當(dāng)選取的偏移速度過小，則雙曲線收斂不完全，反映在對(duì)應(yīng)的圖像熵值偏大；當(dāng)偏移速度過大，此時(shí)會(huì)出現(xiàn)偏移過聚焦問題，繞射波反向沿更淺處延伸，圖像信息熵值也會(huì)偏大，合適的偏移速度會(huì)使圖像中反射界面繞射波盡可能地收斂，表現(xiàn)在回波圖像的信息熵處于較低值。

地質(zhì)雷達(dá)圖像信息可表示為：

式中，aij表示地質(zhì)雷達(dá)圖像第i行第j道采樣數(shù)據(jù)的樣本值。圖像信息熵值定義為：

上式中，E_t表示所有采樣點(diǎn)的總能量值；Pij表示單個(gè)信號(hào)點(diǎn)能量占總能量的比值；D表示圖像信息熵值。

3? 隧道襯砌鋼筋正演模擬圖像偏移處理分析

3.1? 理想不含噪聲正演計(jì)算圖像

采用GPRMAX軟件程序建立隧道襯砌結(jié)構(gòu)探測(cè)正演計(jì)算模型，襯砌結(jié)構(gòu)存在兩層正對(duì)鋼筋，以模擬實(shí)際工程中隧道襯砌結(jié)構(gòu)內(nèi)鋼筋的探測(cè)與識(shí)別。模型尺寸為長(zhǎng)×寬= 2.00 m × 0.95 m，空氣層厚度為0.05 m，天線步距為0.01 m，天線收發(fā)距為0.1 m，信號(hào)道數(shù)為192，天線中心頻率為900 MHz，混凝土的介電常數(shù)為6.0，鋼筋直徑0.04 m，橫向間距0.14 m，排距為0.1 m。上層鋼筋埋深為0.7 m，上層鋼筋埋深0.8 m，每排鋼筋數(shù)量均為10根，如表1所示。正演建模及計(jì)算結(jié)果如圖1所示。

由圖可知，由于鋼筋與混凝土的電性差異大，雷達(dá)波在傳播過程中遇到鋼筋發(fā)生強(qiáng)烈的電磁反射，鋼筋回波信號(hào)呈開口向下的雙曲線弧形（月牙狀），其頂部為鋼筋所在位置，淺層鋼筋的雙曲線反射較為強(qiáng)烈，但同時(shí)也會(huì)削弱下層鋼筋的反射，導(dǎo)致下層鋼筋的反射信號(hào)變?nèi)酰灰追直?。首先?duì)正演計(jì)算圖像進(jìn)行靜校正處理，以去除直達(dá)波。由于正演模型襯砌混凝土介電常數(shù)為εr = 6，根據(jù)介電常數(shù)法，電磁波在襯砌混凝土中傳播速度估計(jì)值為? cm/ns。基于該偏移速度值，設(shè)定偏移速度區(qū)間為[9.25 cm/ns，15.25 cm/ns]，速度掃描間隔為0.2 cm/ns。分別采用F-K法和Kirchhoff法對(duì)鋼筋正演探測(cè)圖像進(jìn)行偏移校正處理，并計(jì)算每個(gè)速度值對(duì)應(yīng)的鋼筋偏移圖像信息熵值，結(jié)果如圖2和圖3所示。

由圖可知，兩種偏移方法的速度-信息熵曲線都非常規(guī)則，有明顯的變化規(guī)律，并且兩種方法的速度-信息熵曲線中均存在極小信息熵值，信息熵的極小值對(duì)應(yīng)的偏移速度均為Vbest = 12.45 cm/ns，表明相對(duì)于其他速度值，采用該最佳速度值對(duì)鋼筋正演圖像進(jìn)行偏移成像處理，鋼筋回波信號(hào)能得到最大程度歸位和收斂?；谠撟罴哑扑俣戎?，分別采用F-K法和Kirchhoff法對(duì)鋼筋探測(cè)正演圖像進(jìn)行偏移歸位處理，結(jié)果如圖4所示。

經(jīng)偏移處理后，鋼筋兩側(cè)繞射信號(hào)得到有效歸位，偏移圖像鋼筋信號(hào)的形狀與其真實(shí)形狀接近，位置也得到還原，圖像分辨率較高，鋼筋回波信號(hào)清晰可見，能量較為聚焦。經(jīng)計(jì)算，F(xiàn)-K法襯砌鋼筋偏移圖像信息熵值為9.31，而Kirchhoff法對(duì)應(yīng)的鋼筋偏移圖像信息熵值為9.60，由F-K法得到的鋼筋偏移圖像信息熵較小。相對(duì)于Kirchhoff法，經(jīng)F-K法偏移處理后，鋼筋回波信號(hào)分布更為集中，繞射信號(hào)收斂更為完全，表明在理想的無噪聲環(huán)境下，F(xiàn)-K偏移法對(duì)襯砌雙層鋼筋的偏移成像效果優(yōu)于Kirchhoff法。

3.2? 含噪聲正演計(jì)算圖像

同樣建立襯砌結(jié)構(gòu)雙層鋼筋正演探測(cè)幾何模型，模型大小、鋼筋及背景介質(zhì)情況與上述理想不含噪聲情況下的襯砌鋼筋檢測(cè)模型一致，并開展模型探測(cè)正演計(jì)算。在正演圖像中加入強(qiáng)度為45 dBW的高斯白噪聲，幾何模型及加噪后的正演計(jì)算圖像如圖5所示。對(duì)比圖5和圖1可知，在其他參數(shù)一致的條件下，由于圖像受噪聲影響，正演計(jì)算圖像中鋼筋回波信號(hào)更為模糊，雙曲線尾部交叉點(diǎn)難以辨別。

對(duì)含噪的正演計(jì)算圖像進(jìn)行去直達(dá)波處理，由于襯砌混凝土介電常數(shù)不變，因此電磁波在襯砌中的傳播速度估計(jì)值同樣為v2 = 12.25 cm/ns，設(shè)定速度掃描區(qū)間為[9.25 cm/ns，15.25 cm/ns]。基于該速度區(qū)間，分別采用F-K法和Kirchhoff法對(duì)去直達(dá)波后的襯砌鋼筋探測(cè)圖像進(jìn)行成像處理，結(jié)果如圖6、圖7所示。

對(duì)比圖6和圖7可知，F(xiàn)-K偏移法速度-圖像信息熵曲線最小信息熵為11.28，最小信息熵對(duì)應(yīng)的偏移速度為12.45 cm/ns；Kirchhoff偏移法速度-圖像信息熵曲線最小信息熵為11.31，對(duì)應(yīng)的最佳偏移速度為13.05 cm/ns；與理想無噪聲情況下經(jīng)偏移處理后得到的速度-信息熵曲線相比較，F(xiàn)_K偏移法速度-信息熵曲線仍然非常規(guī)則，規(guī)律明顯，辨識(shí)度高，存在極小信息熵值，可采用該準(zhǔn)確的偏移速度值進(jìn)行偏移成像處理；而Kirchhoff偏移法速度-信息熵曲線較為凌亂，圖像信息熵值無明顯變化規(guī)律，難以獲取準(zhǔn)確有效的偏移速度值，從而影響襯砌鋼筋回波信號(hào)偏移成像歸位。

由圖8可知，F(xiàn)-K法偏移處理后，圖像信息熵值略小于由克希霍夫法偏移處理得到的圖像信息熵，表明F-K法對(duì)含噪圖像仍然具有很好的偏移歸位效果，較克希霍夫偏移法具有一定的抗噪性，且F-K法處理后，含噪圖像分辨率較高，鋼筋反射信號(hào)清晰可見，能量較為聚焦，而克?；舴蚍ㄆ坪螅捎谄漭^差的抗噪性，圖像鋼筋反射回波易被噪聲干擾信號(hào)所掩蓋，經(jīng)偏移處理后的圖像模糊，較難識(shí)別和辨認(rèn)，無法準(zhǔn)確還原鋼筋位置。

4? 結(jié)? 論

基于探地雷達(dá)電磁波基本理論和GPRMAX正演數(shù)值方法，分別采用F-K法與Kirchhoff法對(duì)隧道襯砌雙層鋼筋探測(cè)正演圖像進(jìn)行偏移成像處理，通過分析含噪及理想不含噪條件下雙層鋼筋的偏移歸位情況，并運(yùn)用圖像熵理論對(duì)偏移效果進(jìn)行定量計(jì)算和評(píng)價(jià)，比較兩種偏移方法的實(shí)際運(yùn)用效果。相對(duì)于Kirchhoff偏移法，F(xiàn)-K偏移法獲取的偏移速度-信息熵曲線更規(guī)則，計(jì)算得到的偏移速度更為準(zhǔn)確合理，偏移處理后鋼筋信號(hào)更聚焦、能量更為集中，兩側(cè)繞射波收斂則更完全，并且F-K偏移法對(duì)含高強(qiáng)度噪聲的鋼筋檢測(cè)圖像也具有很好的歸位效果，較克?；舴蚱品ň哂幸欢ǖ目乖胄浴-K偏移壓制了多次波和繞射波，有利于鋼筋位置的還原和識(shí)別，提高了圖像解釋精度。因此，在隧道襯砌雙層鋼筋的偏移成像處理與識(shí)別中，F(xiàn)-K法具有優(yōu)越性和較好的適用性。

參考文獻(xiàn)：

[1] 陳偉，常軍，鞏文龍.探地雷達(dá)在隧道襯砌缺陷檢測(cè)中的應(yīng)用 [J].蘇州科技大學(xué)學(xué)報(bào)：工程技術(shù)版，2020，33（3）：54-61+71.

[2] 吳慶全，楊子龍，羅偉斌.基于GprMax正演模擬及工程實(shí)例的地質(zhì)雷達(dá)鋼筋探測(cè)研究 [J].礦產(chǎn)與地質(zhì)，2017，31（3）：652-656.

[3] 舒志樂，胥晏，黃強(qiáng)，等.隧道襯砌空洞模型試驗(yàn)與有限元正演模擬 [J].現(xiàn)代電子技術(shù)，2022，45（11）：18-23.

[4] 尹濤，王瑞燕，高鵬，等.基于探地雷達(dá)技術(shù)的隧道二襯缺陷探測(cè)研究 [J].地下空間與工程學(xué)報(bào)，2018，14（S1）：220-226.

[5] 楊宇，凌同華，廖艷程.混凝土構(gòu)件深淺埋鋼筋模擬檢測(cè)試驗(yàn)與偏移分析 [J].土木建筑與環(huán)境工程，2018，40（6）：139-145.

[6] 鄧世坤.克?；舴蚍e分偏移法在探地雷達(dá)圖象處理中的應(yīng)用 [J].地球科學(xué)，1993（3）：303-309+368.

[7] 李文超，王緒本，鄧薇.探地雷達(dá)FDTD正演模擬及F-K偏移處理 [J].勘察科學(xué)技術(shù)，2007（3）：33-36.

[8] 修志杰，陳潔，方廣有，等.基于F-K偏移及最小熵技術(shù)的探地雷達(dá)成像法 [J].電子與信息學(xué)報(bào)，2007（4）：827-830.

[9] 黃玲，曾昭發(fā)，王者江，等.鋼筋混凝土缺陷的探地雷達(dá)檢測(cè)模擬與成像效果 [J].物探與化探，2007（2）：181-185+180.

[10] 林志強(qiáng)，王磊，樊斌斌.基于圖像熵的探地雷達(dá)Kirchhoff偏移成像算法 [J].火力與指揮控制，2020，45（12）：97-100.