李 巖 羅 茜 劉睿浩 云 澤 張冀翔

（中國石油大學(xué)（北京）a.機械與儲運工程學(xué)院；b.過程流體過濾與分離技術(shù)北京市重點實驗室）

氟化工泛指所有含氟元素的產(chǎn)品及其衍生品，可分為有機氟化工和無機氟化工兩類，其品類繁多，廣泛應(yīng)用于國防軍工、建筑、農(nóng)業(yè)、醫(yī)療及新能源等行業(yè)，被稱為“黃金產(chǎn)業(yè)”。 氫氟酸作為氟化工行業(yè)最重要的中間體，高溫會促使HF電離并提高介質(zhì)與金屬間的反應(yīng)速率，進(jìn)而加劇其對金屬管道的腐蝕。 管道發(fā)生腐蝕以后，會表現(xiàn)出管壁變薄、出現(xiàn)蝕損斑等情況，存在一定的安全隱患。 因此，通過科學(xué)有效的方式對管道壁厚進(jìn)行監(jiān)測，及時發(fā)現(xiàn)腐蝕部位以便采取相應(yīng)措施預(yù)防安全事故的發(fā)生極為重要。

目前， 檢測管道內(nèi)腐蝕的方法主要有漏磁檢測、渦流檢測、超聲波和射線檢測技術(shù)［1～3］。漏磁檢測方法適用于中小型管道的快速檢測，成本較低［2］，但是檢測靈敏度低；渦流檢測方法屬于無損檢測技術(shù)，通過檢測交變磁場中的渦流變化進(jìn)而完成腐蝕部位的定位分析［4］；超聲波檢測方法適用于壁厚較厚、口徑大的管道的檢測，能夠高效完成缺陷定位，但無法確定缺陷性質(zhì)［5］；射線檢測方法能夠?qū)艿乐械娜毕莶课贿M(jìn)行有效檢測且靈敏度高，但是其檢測成本較高。

在一些惡劣條件下， 如檢測高溫管道壁厚時，傳統(tǒng)的超聲檢測技術(shù)面臨很大的挑戰(zhàn)，主要表現(xiàn)在直接接觸會導(dǎo)致傳感器失效，其原因在于傳感器的耦合特性和耐受性［6］。 基于超聲導(dǎo)波的檢測技術(shù)是一種比較新的超聲檢測方法，可以應(yīng)用于極端工況下的材料評估和表征［7］。因此，基于超聲導(dǎo)波檢測技術(shù)開發(fā)出一套非侵入式的壁厚檢測裝備，并利用COMSOL Multiphysics多物理場仿真軟件對波導(dǎo)桿和試件內(nèi)超聲波場進(jìn)行了模擬分析，為氟化工裝置中高溫管道壁厚在線檢測提供了理論指導(dǎo)。

1 波導(dǎo)桿結(jié)構(gòu)設(shè)計

波導(dǎo)桿是用來引導(dǎo)波傳播的結(jié)構(gòu)， 通常由銅、鋁及不銹鋼等材質(zhì)制成。 它也能作為溫度緩沖結(jié)構(gòu)，率先提出采用波導(dǎo)桿作為高溫試件的溫度緩沖結(jié)構(gòu)并申請了專利的是美國學(xué)者CROSS N O［8］，他采用不完全焊接的方式將一根圓柱形波導(dǎo)桿焊接在工件表面。 隨后，很多學(xué)者也采用了與被測管道相同材質(zhì)的圓柱形波導(dǎo)桿開展了相關(guān)試驗研究［3,9］。 但是，使用這種長圓柱形波導(dǎo)桿的一個關(guān)鍵問題是，在超聲波傳播過程中發(fā)生在桿側(cè)表面的反射波和模態(tài)轉(zhuǎn)換波會引起拖尾信號的產(chǎn)生。 FOUDZI F M和IHARA I［10］曾基于有限差分法的三維數(shù)值模擬軟件Wave 3000， 探究具有不同多邊形截面的波導(dǎo)桿的超聲波脈沖回波行為。 IZMAILOVA E V等開發(fā)了一種通過聲學(xué)方法檢測管道腐蝕失效的技術(shù)，根據(jù)回波信號判斷管壁減薄等缺陷及其幾何尺寸［11］。

采用波導(dǎo)桿作為溫度緩沖結(jié)構(gòu)使被檢工作表面的高溫降至超聲探頭正常工作的溫度范圍，不僅與波導(dǎo)桿的長度有關(guān)， 還與其材料屬性、幾何形狀等因素有關(guān)，因此合理設(shè)計波導(dǎo)桿結(jié)構(gòu)是非侵入式壁厚檢測技術(shù)的關(guān)鍵之一。

筆者采用COMSOL Multiphysics軟件分別對銅、 鋁和不銹鋼3種材質(zhì)的不同截面形狀和尺寸進(jìn)行優(yōu)化， 最后給出適用于500 ℃高溫工況下的波導(dǎo)桿結(jié)構(gòu)設(shè)計方案。 該溫度緩沖過程包含熱傳導(dǎo)、熱對流和熱輻射3種傳熱方式。 在COMSOL中選擇輻射和固體傳熱模塊進(jìn)行穩(wěn)態(tài)分析，仿真模型和網(wǎng)格劃分如圖1所示。 圓柱形波導(dǎo)桿和矩形波導(dǎo)桿的截面面積和長度均相同。 在被檢工件下表面施加500 ℃高溫?zé)嵩矗?以模擬管道內(nèi)的高溫介質(zhì)，環(huán)境溫度設(shè)為25 ℃，添加對流換熱及表面對表面輻射邊界條件。

圓柱形和矩形波導(dǎo)桿的溫度緩沖曲線如圖2a所示， 可見矩形波導(dǎo)桿的溫度緩沖效果更好；316L不銹鋼、鋁和銅質(zhì)矩形波導(dǎo)桿的溫度緩沖曲線如圖2b所示， 相對而言，316L不銹鋼是制作波導(dǎo)桿的最佳材質(zhì)， 且300 mm即可滿足降溫需求；采取控制變量法研究了矩形波導(dǎo)桿截面不同寬度w和厚度d對溫度緩沖效果的影響（圖2c、d），可以發(fā)現(xiàn)，矩形波導(dǎo)桿的厚度是影響其降溫效果的重要尺寸，寬度影響不大。 對于波導(dǎo)桿的結(jié)構(gòu)設(shè)計，應(yīng)在滿足降溫長度的前提下，盡可能選用大寬厚比的316L不銹鋼矩形形狀波導(dǎo)桿。 考慮到加工難易程度， 設(shè)計一個長寬厚分別為300 mm×18 mm×1 mm的矩形波導(dǎo)桿可滿足500 ℃高溫測量需求。

圖2 影響波導(dǎo)桿降溫效果的主要因素

2 非侵入式管道壁厚檢測系統(tǒng)設(shè)計

2.1 原理

圖3為本研究所用基于波導(dǎo)結(jié)構(gòu)的管道壁厚檢測示意圖，利用超聲脈沖反射原理對被測管道進(jìn)行壁厚測量。 超聲探頭發(fā)射的脈沖信號經(jīng)波導(dǎo)桿傳播到工件表面，發(fā)生反射和透射。 透射波進(jìn)入被檢工件并在下表面發(fā)生反射，反射信號沿波導(dǎo)返回并由超聲探頭接收。 在聲速已知的情況下，通過測量相鄰兩底面反射信號間的時間間隔便可計算出壁厚。

圖3 管道壁厚檢測示意圖

本研究采用的超聲導(dǎo)波為水平剪切波（SH波），其振動方向與波傳播的方向垂直，其頻率大小需要滿足頻寬積大于準(zhǔn)SH0模態(tài)下的非頻散臨界頻寬積（不銹鋼的頻寬積為15 MHz·mm），頻厚積小于SH1模態(tài)下的截止頻厚積（不銹鋼的頻厚積為1.6 MHz·mm）［12］，如下：

式中 d——波導(dǎo)桿的厚度；

f——激勵信號頻率；

w——波導(dǎo)桿的寬度。

模擬過程中采用頻率為1 MHz的高斯脈沖信號作為激勵信號，其函數(shù)表達(dá)式為：

2.2 超聲導(dǎo)波傳播特性

衡量超聲波作用效果的一個重要參數(shù)是聲場中的聲壓分布， 因此采用COMSOL軟件對超聲導(dǎo)波在波導(dǎo)桿和待測工件中的傳播過程進(jìn)行瞬態(tài)求解，根據(jù)聲壓分布進(jìn)而研究超聲導(dǎo)波的傳播特性與機理。 圖4為不同時刻下的聲壓分布云圖，根據(jù)聲壓分布情況可以看出超聲導(dǎo)波在待測工件上表面發(fā)生了反射和透射現(xiàn)象。 反射信號沿波導(dǎo)桿返回并由超聲探頭接收，透射進(jìn)入待測工件的超聲導(dǎo)波在工件內(nèi)部的有限邊界發(fā)生反射，并伴有模態(tài)轉(zhuǎn)換，反射信號同樣經(jīng)波導(dǎo)桿返回并由超聲探頭接收。

圖4 不同時刻的聲壓分布

圖5為超聲換能器表面的平均聲壓曲線，除初始信號，波峰較高的信號是發(fā)生在待測工件上表面的反射信號；波峰較低的是發(fā)生在下表面的反射信號，分別記為一次波、二次波、三次波。 觀察聲壓曲線可以發(fā)現(xiàn)，接收到的脈沖回波信號彼此分離純度較高，波形穩(wěn)定，并且清晰可辨。 由于矩形波導(dǎo)桿的厚度遠(yuǎn)小于水平剪切波的波長，當(dāng)這種水平剪切波入射時理論上可以將其視為一種反平面剪切線源，其在待測工件內(nèi)部只會激勵出水平剪切波，并沿各方向均勻擴散。 因此，水平剪切波可以作為一種用于壁厚檢測的理想的激勵源模式。

圖5 換能器表面平均聲壓曲線

2.3 試驗結(jié)果分析

圖6所示為搭建的非侵入式超聲波壁厚檢測裝備，該裝備由信號發(fā)生器、超聲探頭、耦合劑、固定夾具、波導(dǎo)桿、待測工件、示波器和計算機組成。 與Tektronix-AFG信號發(fā)生器連接的奧林巴斯V154-RB橫波探頭通過夾具固定在波導(dǎo)桿的一端作為激勵端，另一個橫波探頭連接波導(dǎo)桿作為接收端， 接收到的脈沖回波信號通過Tekteonix-TBS1102B-EDU示波器顯示。 試驗時，由信號發(fā)生器產(chǎn)生一個五周期正弦脈沖信號，該脈沖信號接到奧林巴斯橫波探頭作用于波導(dǎo)桿，并沿波導(dǎo)桿傳播到待測工件表面，在待測工件中激勵出水平剪切波，并在厚度方向上發(fā)生多次反射，反射回波由示波器進(jìn)行采集和顯示。

圖6 非侵入式超聲波壁厚檢測裝備

為準(zhǔn)確計算待測工件厚度值，首先應(yīng)確定回波信號中波峰的具體位置。 將示波器采集到的回波信號導(dǎo)入至MATLAB進(jìn)行Hilbert變換，圖7是經(jīng)過Hilbert變換后得到的10 mm待測工件的回波信號包絡(luò)曲線，波峰波谷位置清晰可辨。

圖7 回波信號包絡(luò)曲線

對兩個相鄰波峰鄰域內(nèi)的包絡(luò)曲線求極值點，分別記為tA和tB，根據(jù)下式即可計算得到待測試件厚度δ：

式中 c——超聲波波速。

試驗所用波導(dǎo)桿和待測工件材質(zhì)均為316L不銹鋼，已知常溫狀態(tài)下（25 ℃），超聲波在該材質(zhì)中的剪切波速為3 154.4 m/s， 因而得到試件測量厚度值為10.09 mm，相對誤差為0.9%。 利用該設(shè)備對不同溫度下的試件進(jìn)行了測量，測量結(jié)果如圖8所示。 從圖中可以看出，隨著溫度上升，測量厚度幾乎線性增長， 導(dǎo)致測量誤差越來越大，試件溫度高于250 ℃時， 誤差已經(jīng)達(dá)到5%以上。這是由于溫度升高，超聲波速會發(fā)生變化，而當(dāng)前是按照固定波速計算，因此會導(dǎo)致測量誤差增大。 因此針對更高溫度工件檢測時，必須要考慮溫度的影響，后續(xù)研究需要對測量數(shù)據(jù)計算進(jìn)行溫度補償修正，從而得到更準(zhǔn)確的結(jié)果。

圖8 不同溫度下的試件測量厚度

3 結(jié)束語

基于超聲導(dǎo)波無損檢測理念，設(shè)計了一種非侵入式的管道壁厚檢測裝備， 其結(jié)構(gòu)簡單輕便，可通過夾具進(jìn)行壓力耦合固定，實現(xiàn)壁厚定點檢測。 采用有限元仿真軟件對波導(dǎo)桿進(jìn)行熱仿真分析， 設(shè)計了一種大寬厚比的矩形波導(dǎo)桿結(jié)構(gòu)， 該結(jié)構(gòu)既能實現(xiàn)溫度緩沖， 又能夠有效地傳遞非頻散水平剪切波。 聲場模擬結(jié)果也說明， 使用該矩形波導(dǎo)對待測工件進(jìn)行厚度測量， 超聲換能器接收到的各次反射回波信號清晰可辨， 利用相鄰回波信號之間的時間差可計算得到厚度值。 試驗測量了室溫下10 mm厚的316L不銹鋼板， 接收到的脈沖回波信號波峰波谷清晰可見， 信號彼此分離， 并且回波信號經(jīng)處理后較為準(zhǔn)確地計算出了工件厚度。 利用該設(shè)備對不同溫度下試件厚度進(jìn)行測量， 發(fā)現(xiàn)測量誤差隨溫度幾乎線性增加， 因此高溫工件檢測時需要進(jìn)行溫度修正， 從而得到更精準(zhǔn)的結(jié)果。 筆者基于波導(dǎo)結(jié)構(gòu)并利用SH超聲導(dǎo)波的厚度測量方法， 將常規(guī)超聲換能器與高溫管道分離開來， 實現(xiàn)了對高溫管道壁厚的間接測量，若結(jié)合無線傳輸技術(shù)將能夠?qū)崿F(xiàn)對高溫管道腐蝕狀況的長期在線監(jiān)測。