基于彈性波映像識(shí)別法的鋼管隧洞缺陷檢測(cè)

武 鵬，袁剛烈，車愛蘭

(上海交通大學(xué) 船舶海洋與建筑工程學(xué)院，上海 200040)

0 引 言

鋼管隧洞是一種由鋼管和混凝土構(gòu)成的組合構(gòu)件，鋼與混凝土的協(xié)調(diào)作用使得鋼管隧洞具有承載能力高、抗災(zāi)性能好、施工方便等特點(diǎn)，已廣泛應(yīng)用于輸水管道和橋梁等土木工程結(jié)構(gòu)中[1-3]。由于周圍土、水等復(fù)雜荷載的長(zhǎng)期作用，鋼管隧洞的鋼管出現(xiàn)變形、連接處開裂、連接層脫空等病害，大大降低了鋼管隧洞的抗剪能力[4]，嚴(yán)重影響工程安全[5]。

由于彈性波測(cè)試方法適用于復(fù)雜環(huán)境、可實(shí)現(xiàn)定量評(píng)價(jià)等優(yōu)勢(shì)，近些年來(lái)，在土木工程檢測(cè)、監(jiān)測(cè)等領(lǐng)域得到了廣泛的應(yīng)用。彈性波的傳播速度與介質(zhì)的密度、強(qiáng)度等參數(shù)有著密切的關(guān)系[6]。張獻(xiàn)民等[7]通過(guò)建立土石混合料的剪切波速模型，研究了土石混合料剪切波速與瑞雷面波波速的關(guān)系，完成了土石混填地基壓實(shí)質(zhì)量的評(píng)價(jià)；朱仁杰等[8]探究了彈性波沖擊響應(yīng)能量與結(jié)構(gòu)缺陷類型、大小、位置之間的相關(guān)性，對(duì)高位井結(jié)構(gòu)下部缺陷開展了現(xiàn)場(chǎng)試驗(yàn)，評(píng)價(jià)了缺陷分布和可能存在的滲透路徑；B.ZIMA等[9]基于導(dǎo)波技術(shù)和橢圓二值損傷成像算法，實(shí)現(xiàn)了通過(guò)彈性波對(duì)鋼絞線、鋼筋、錨固螺栓等張拉構(gòu)件的無(wú)損監(jiān)測(cè)；王平等[10]建立了高速道岔有限元模型，對(duì)在健康尖軌和不同斷軌位置傳播的彈性波信號(hào)進(jìn)行了時(shí)頻分析，建立了斷軌位置和彈性波信號(hào)間的關(guān)聯(lián)。

隨著計(jì)算機(jī)學(xué)、信息學(xué)的發(fā)展，根據(jù)樣本特征將樣本歸類的模式識(shí)別法得到了發(fā)展。在土木工程中，模式識(shí)別法常用于結(jié)構(gòu)的缺陷識(shí)別，如將模式匹配神經(jīng)網(wǎng)絡(luò)用于鋼梁的損傷定性和定量分析[11]。k-means算法是模式識(shí)別中的重要算法，是一種經(jīng)典的無(wú)監(jiān)督聚類和學(xué)習(xí)的算法，其因簡(jiǎn)潔、運(yùn)算速度快，并能有效地處理大數(shù)據(jù)集而被廣泛應(yīng)用于模式化問(wèn)題的解析中[12]。

筆者針對(duì)鋼管隧洞缺陷檢測(cè)問(wèn)題，基于彈性波傳播基本原理，考慮鋼管隧洞的結(jié)構(gòu)特點(diǎn)以及缺陷的模式化特征，結(jié)合聚類分析方法，提出了一種鋼管隧洞缺陷快速檢測(cè)評(píng)價(jià)方法。首先，詳細(xì)闡述了彈性波映像識(shí)別方法的工作原理、采集方法及算法思路；然后，對(duì)某鋼管隧洞開展了彈性波測(cè)試并對(duì)測(cè)試結(jié)果進(jìn)行了分析；最后，將試驗(yàn)結(jié)果與3D激光掃描結(jié)果及施工數(shù)據(jù)進(jìn)行對(duì)比分析，驗(yàn)證了彈性波映像識(shí)別法對(duì)缺陷檢測(cè)的有效性。

1 彈性波映像識(shí)別方法

1.1 鋼管隧洞中彈性波傳播特性

鋼管隧洞由空心鋼管外包混凝土構(gòu)成，混凝土填充鋼殼體外部以提高管節(jié)結(jié)構(gòu)的整體剛度和強(qiáng)度。將鋼管隧洞視為多層層狀介質(zhì)，當(dāng)介質(zhì)中某一層的彈性參數(shù)發(fā)生變化時(shí)，彈性波場(chǎng)的傳播特性會(huì)隨著彈性參數(shù)的改變而改變。如圖1，在介質(zhì)表面激勵(lì)i次(i=1，2，3，…，n)，當(dāng)激勵(lì)位置位于Si時(shí)，振動(dòng)沿著介質(zhì)深度方向傳播，形成波動(dòng)，彈性波信號(hào)Ri由固定在激勵(lì)點(diǎn)旁的檢波器接收。根據(jù)接收到的彈性波波形特征，如振幅大小、衰減快慢以及持續(xù)時(shí)間長(zhǎng)短等，分析彈性波響應(yīng)變化，判斷鋼管隧洞洞壁受腐蝕程度。

圖1 彈性波探測(cè)原理Fig. 1 Detecting principle of elastic wave

彈性波場(chǎng)以體波及面波形式在介質(zhì)內(nèi)部傳播，介質(zhì)表面的彈性波場(chǎng)分布與介質(zhì)內(nèi)部的構(gòu)造以及介質(zhì)的物理性質(zhì)密切相關(guān)。彈性波場(chǎng)在兩層介質(zhì)中傳播時(shí)，彈性波的反射系數(shù)按式(1)計(jì)算：

(1)

式中：RPP、RSS分別為縱波、橫波的反射系數(shù)；ρ1(ρ2)、VP1(VP2)、VS1(VS2)分別為介質(zhì)1(2)的密度、縱波速度、橫波速度。

可見，反射波的反射系數(shù)由2種介質(zhì)的波速V與密度ρ的乘積即波阻抗ρV決定。2種介質(zhì)的波阻抗相差越大，反射系數(shù)也越大，產(chǎn)生的反射波響應(yīng)能量也就越強(qiáng)。當(dāng)介質(zhì)內(nèi)部結(jié)構(gòu)發(fā)生變化時(shí)，如內(nèi)部產(chǎn)生裂縫或空洞，介質(zhì)密度減小將導(dǎo)致波速降低[13]，形成反射界面，接收到的波形響應(yīng)能量就會(huì)變強(qiáng)。裂縫或空洞大小不同，檢波器接收到的彈性波響應(yīng)Ri特征就不同。

根據(jù)彈性波響應(yīng)能量E，結(jié)合圖像處理技術(shù)及拓?fù)鋵W(xué)方法，可以評(píng)價(jià)響應(yīng)能量分布，實(shí)現(xiàn)缺陷分布的定量評(píng)價(jià)。

彈性波響應(yīng)能量E按式(2)計(jì)算：

(2)

式中：f(t)為所接收彈性波隨時(shí)間t變化的函數(shù)；t1、t2分別為響應(yīng)能量計(jì)算起始、終止時(shí)間。

1.2 基于k-means聚類的彈性波映像識(shí)別方法

彈性波信號(hào)在鋼管隧洞的空隙、空洞、裂縫等缺陷處的傳播特征不同于在完好結(jié)構(gòu)處的。在缺陷處接收到的彈性波信號(hào)及響應(yīng)能量分別用R′i、E′表示。將彈性波響應(yīng)能量進(jìn)行聚類分類，確定缺陷響應(yīng)能量閾值，從而可判斷檢測(cè)區(qū)域內(nèi)屬于正常或缺陷模式。

在對(duì)彈性波響應(yīng)能量進(jìn)行k-means聚類分析時(shí)，任意一個(gè)彈性波響應(yīng)的位置屬于正常位置處的Ri或缺陷位置處的R′i。對(duì)所有的彈性波響應(yīng)數(shù)據(jù)以空間的2個(gè)點(diǎn)進(jìn)行聚類，對(duì)最靠近這2個(gè)點(diǎn)的其余響應(yīng)能量值進(jìn)行歸類，通過(guò)迭代法，循環(huán)更新2聚類中心的值，直至得到目標(biāo)聚類結(jié)果。

假設(shè)有n個(gè)待分析彈性波響應(yīng)能量值，記為E={E1，E2，E3，…，En}，初始化2個(gè)聚類中心{C1，C2}，則每個(gè)響應(yīng)能量值Ei(i=1，2，3，…，n)到聚類中心Cj(j=1，2)的距離按式(3)計(jì)算：

(3)

式中：Eit為第x個(gè)對(duì)象的第t個(gè)屬性；Cjt為第j個(gè)聚類中心的第t個(gè)屬性。

依次比較每一個(gè)響應(yīng)能量值Ei到每一個(gè)聚類中心Cj距離，將對(duì)象分配到距離最近的聚類中心類簇中，得到2個(gè)類簇{S1，S2}。

待所有點(diǎn)分配完畢之后，對(duì)類簇內(nèi)的所有點(diǎn)取平均值，作為該類簇新的中心點(diǎn)，再重新將所有點(diǎn)進(jìn)行分配并更新類簇中心點(diǎn)，直至類簇中心點(diǎn)的變化很小，形成與2類彈性波響應(yīng)能量相對(duì)應(yīng)的劃分閾值，即正常位置處的Ri和缺陷位置處的R′i，兩者的分界閾值即是對(duì)彈性波響應(yīng)能量進(jìn)行狀態(tài)分類的標(biāo)準(zhǔn)。

2 鋼管隧洞缺陷現(xiàn)場(chǎng)檢測(cè)

圖2為現(xiàn)場(chǎng)檢測(cè)的輸水鋼管隧洞截面。鋼管內(nèi)徑4.8 m，鋼管壁厚0.02 m；外部回填厚0.6 m C20混凝土；每節(jié)管段長(zhǎng)100 m。目測(cè)在內(nèi)部水壓力及圍壓作用下，隧洞已出現(xiàn)了明顯的表面變形及鋼管后部積水等病害。

圖2 鋼管隧洞截面Fig. 2 Section of steel tube tunnel

檢測(cè)工序?yàn)椋菏占摴芩矶吹氖┕?shù)據(jù)→對(duì)鋼管隧洞進(jìn)行3D掃描→分析鋼管的變形特性→在鋼管變形區(qū)域開展彈性波測(cè)試→對(duì)比分析彈性波測(cè)試結(jié)果與3D掃描分析結(jié)果以及施工數(shù)據(jù)。

2.1 切割及鉆孔

對(duì)鋼管進(jìn)行局部切割及鉆孔，可見鋼管與后方混凝土形成嚴(yán)重脫空現(xiàn)象，鋼管表面出現(xiàn)鼓突，變形區(qū)域最大直徑約0.46 m，最大脫空深度約0.30 m，已有滲漏水現(xiàn)象。

2.2 3D激光掃描檢測(cè)

采用測(cè)量距離精度為2 mm，角度精度為12″的瑞士徠卡三維脈沖式激光掃描儀Scan Station C10，在隧洞內(nèi)開展3D激光掃描檢測(cè)，以獲取鋼管隧洞內(nèi)部表面變形數(shù)據(jù)。

將鋼管隧洞的圓柱內(nèi)表面沿母線展開呈現(xiàn)長(zhǎng)22 m、寬12 m的矩形區(qū)域，即為檢測(cè)區(qū)域。對(duì)檢測(cè)區(qū)域進(jìn)行掃描得到的等值線如圖3?？梢?，檢測(cè)區(qū)域存在明顯的變形，最大相對(duì)位移達(dá)0.28 m。

圖3 鋼管隧洞3D掃描等值線Fig. 3 3D scanning contour of steel tube tunnel

2.3 彈性波檢測(cè)

2.3.1 檢測(cè)系統(tǒng)

檢測(cè)系統(tǒng)包括：Geometrics公司生產(chǎn)的24通道高精度寬頻帶地震儀；固有頻率為100 Hz的速度檢波器；筆者團(tuán)隊(duì)開發(fā)的手持式泡沫耦合裝置及激勵(lì)錘；筆記本電腦、電纜、激發(fā)裝置、電源等。

2.3.2 檢測(cè)步驟

1)將速度檢波器置于耦合裝置內(nèi)，確保速度檢波器與鋼管隧洞內(nèi)壁良好接觸，實(shí)現(xiàn)彈性波數(shù)據(jù)的快速采集。

2)根據(jù)激勵(lì)錘的重量、錘頭形狀等確定震源頻率。由于本次檢測(cè)深度淺，要求精度高，因此選用250 g圓頭激勵(lì)錘激發(fā)800～1 300 Hz的彈性波信號(hào)。圖4為檢測(cè)現(xiàn)場(chǎng)。

圖4 現(xiàn)場(chǎng)檢測(cè)Fig. 4 On-site inspection

3)選取變形顯著的區(qū)域(圖5)作為研究區(qū)域；在隧洞側(cè)壁、頂、底部共布設(shè)16條測(cè)線，每條測(cè)線長(zhǎng)22 m，測(cè)線間距0.8 m，沿測(cè)線方向每隔0.25 m布設(shè)1個(gè)檢測(cè)點(diǎn)，形成0.25 m×0.8 m測(cè)試網(wǎng)格，共計(jì)1 584個(gè)測(cè)點(diǎn)。

圖5 測(cè)線布置(單位:m)Fig. 5 Layout of survey lines

4)數(shù)據(jù)采集偏移距為0.2 m。采樣間隔20.833 μs，記錄時(shí)長(zhǎng)0.1 s，采樣延遲0.004 s。圖6為測(cè)線3采集的波形。

圖6 測(cè)線3采集波形Fig. 6 Acquisition waveform of survey line 3

3 鋼管隧洞缺陷分布分析

3.1 數(shù)據(jù)分析及映像識(shí)別

3.1.1 彈性波波形可視化

采集到的彈性波以時(shí)間序列數(shù)據(jù)的形式儲(chǔ)存，通過(guò)彈性波波形可視化計(jì)算，可實(shí)現(xiàn)鋼管隧洞缺陷分布的定性分析。將波形網(wǎng)格化，網(wǎng)格中彈性波數(shù)值的大小代表相應(yīng)網(wǎng)格的能量，網(wǎng)格的能量與網(wǎng)格到波形中軸的距離成正比。將所有網(wǎng)格賦值后，可以按照各網(wǎng)格的數(shù)值給波形圖填充顏色。

以測(cè)線3第14～22 m為例，采集的彈性波信號(hào)經(jīng)可視化處理后得到的波形如圖7。

圖7 測(cè)線3彈性波波形可視化效果Fig. 7 Elastic waveform visualization results of survey line 3

從圖7可以看到，不同位置彈性波響應(yīng)差別較大。與圖3的3D掃描結(jié)果相比，在變形較大區(qū)域，彈性波場(chǎng)的響應(yīng)放大更明顯。

3.1.2 彈性波聚類分析

由于鋼管隧洞的管壁很薄，彈性波在鋼材中傳播速度很快，所以，初始時(shí)刻接收到的彈性波波形為直達(dá)波，隨后接收到的是反射波、透射波等復(fù)雜場(chǎng)的疊加，震源的能量主要集中于直達(dá)波。因此，在計(jì)算彈性波響應(yīng)能量E時(shí)，為了忽略震源對(duì)計(jì)算結(jié)果的影響，確定時(shí)間t1= 0.002 5 s，t2=0.010 0 s，利用k-means聚類方法對(duì)1 167個(gè)檢測(cè)點(diǎn)的彈性波響應(yīng)能量值進(jìn)行聚類劃分，如圖8。

圖8 響應(yīng)能量值聚類分析結(jié)果Fig. 8 Cluster analysis results of response energy values

研究區(qū)內(nèi)閾值取值為9，對(duì)于0

3.2 鋼管隧洞缺陷評(píng)價(jià)及修補(bǔ)

3.2.1 鋼管隧洞缺陷評(píng)價(jià)

將研究區(qū)域鋼管隧洞內(nèi)側(cè)壁沿軸線方向展開，繪制出鋼管隧洞缺陷分布，如圖9。

圖9 鋼管隧洞缺陷分布Fig. 9 Defect distribution of steel tube tunnel

由圖9可見：

1)鋼管隧洞研究區(qū)域存在大面積連續(xù)脫空，尤其是右側(cè)壁及頂部，脫空較嚴(yán)重。對(duì)比發(fā)現(xiàn)：大面積連續(xù)脫空區(qū)域與圖3的3D掃描結(jié)果具有很好的一致性。

2)其它區(qū)域表面雖未見變形，但亦存在脫空缺陷，可認(rèn)為這些部位鋼管與其后方混凝土間存在空洞，由于空洞面積及深度均較小，未導(dǎo)致鋼管變形發(fā)生，但存在極大的安全隱患。

3.2.2 鋼管隧洞補(bǔ)強(qiáng)

由于鋼管變形嚴(yán)重，鋼管隧洞已無(wú)法正常使用，需對(duì)變形嚴(yán)重的右側(cè)壁進(jìn)行補(bǔ)強(qiáng)處治。先確定補(bǔ)強(qiáng)區(qū)域長(zhǎng)度范圍L=[16 m，22 m]，如圖9，再整體切割鋼管，對(duì)鋼管后部混凝土壁進(jìn)行了重新澆筑，并對(duì)切割開的鋼管進(jìn)行了重新焊接。圖10為補(bǔ)強(qiáng)前后鋼管隧洞混凝土壁照片。

圖10 補(bǔ)強(qiáng)前后鋼管隧洞混凝土壁照片F(xiàn)ig. 10 Pictures of concrete wall of steel tube tunnel before and afterreinforcement

3.2.3 鋼管隧洞補(bǔ)強(qiáng)效果檢驗(yàn)

對(duì)修復(fù)后的區(qū)域隨機(jī)選取10個(gè)測(cè)點(diǎn)進(jìn)行彈性波測(cè)試，結(jié)果如圖11?？梢姡瑥椥圆憫?yīng)能量E均低于評(píng)價(jià)閾值9 cm/s，可判斷鋼管隧洞補(bǔ)強(qiáng)效果良好。

圖11 補(bǔ)強(qiáng)后測(cè)點(diǎn)彈性波響應(yīng)能量值Fig. 11 Response energy value of elastic wave at measuringpoints after reinforcement

4 結(jié) 論

1)針對(duì)鋼管隧洞及其缺陷的模式化特點(diǎn)，基于彈性波在層狀介質(zhì)中的傳播特性，結(jié)合聚類分析方法，提出了一種彈性波映像識(shí)別方法。

2)針對(duì)映像方法中缺陷分布及評(píng)價(jià)問(wèn)題，開展了現(xiàn)場(chǎng)檢測(cè)。通過(guò)波形映像定性評(píng)價(jià)了缺陷分布；通過(guò)對(duì)彈性波響應(yīng)能量值進(jìn)行聚類分析，得到了劃分彈性波狀態(tài)的閾值。

3)將彈性波映像法評(píng)價(jià)結(jié)果與3D掃描結(jié)果進(jìn)行比較，兩者的結(jié)果具有較好的一致性。表明彈性波映像識(shí)別法可以快速評(píng)價(jià)鋼管隧洞的缺陷狀態(tài)，具有廣闊的應(yīng)用前景。