林 敏,黃 劼,甘芳吉,廖俊必,耿鵬武2,師根泰
(1.四川大學 制造科學與工程學院,成都 610065;2.空軍駐雅安地區(qū)軍事代表室,雅安 625000)
油氣管道的腐蝕情況嚴重威脅著石化企業(yè)的安全生產(chǎn),有效地監(jiān)測腐蝕引起的管道壁厚減薄程度,可以及時地對管道進行維護,延長其使用壽命,同時大幅減少因管道腐蝕而發(fā)生的安全事故,降低不必要的經(jīng)濟損失[1-3]。超聲波厚度檢測具有測量速度快,檢測精度高,對管道無影響的特點,因而在管道腐蝕監(jiān)測中應用較多[4-6]。
在超聲波長期在線監(jiān)測過程中,其硬件系統(tǒng)可能出現(xiàn)如下3方面的問題:
(1)對于基于壓電效應產(chǎn)生超聲波的硬件系統(tǒng),傳感器與被測固件之間的耦合劑受溫度或應力影響,耦合效果變差,夾層氧化皮變厚,灰塵增加,使得超聲波路徑發(fā)生改變,超聲波入射與出射過程中出現(xiàn)能量損失和頻率變化[7-8],回波波形幅值減小,信噪比降低。
(2)對于基于電磁效應產(chǎn)生超聲波的硬件系統(tǒng),在長期使用的過程中,傳感器與硬件電路的屏蔽能力下降,信號受到電磁干擾,產(chǎn)生噪聲,同樣也會導致信噪比降低[9]。
(3)受溫度變化的影響,半導體元件的工作狀態(tài)發(fā)生改變,導致硬件電路中出現(xiàn)零點漂移的現(xiàn)象,在回波信號中也會產(chǎn)生噪聲信號[10]。
綜合上述3個問題可以看出,長期應用中的超聲波在線監(jiān)測設備采集到的回波信號,其回波幅值、信噪比均會受到環(huán)境影響而降低,因而計算出與幅值、信噪比相關的波形特征點的變化,即可反映硬件系統(tǒng)的健康狀況。
超聲波測厚原理示意如圖1所示,通過超聲波換能器裝置向被測試件發(fā)射超聲波,超聲波在經(jīng)過第一表面和第二表面時分別發(fā)生反射,被超聲波換能器接收[11-12]。
圖1 超聲波測厚原理示意
圖2是以壓電超聲波設備為例,檢測工件時,超聲波換能器接收的回波信號。圖2中的第一表面第一次回波簡稱第一回波,指的是換能器下表面反射的超聲波,第二表面第一次回波簡稱第二回波,指的是被測管道下表面反射的超聲波,理想情況下,采集到的回波信號都可以看成是由激勵信號在不同表面反射得到的。
圖2 超聲波換能器接收回波
通過軟件找到兩次回波的波峰位置之間的時間差T,再根據(jù)超聲波的聲速υ,即可計算出被測件的厚度d為
(1)
在超聲波測厚技術中,最核心的問題是對于超聲波飛行時間(TOF)的測量[13-15]。
需要注意的是,對于壓電超聲設備,這里的第一回波通常取換能器下表面反射的回波,超聲波透射效果變差,在超聲波聲能總量保持不變的情況下,反射產(chǎn)生的第一回波信號將有所增強,第二回波信號將有所減弱。
另外,對于電磁超聲設備中電磁干擾的問題與兩種設備中硬件電路產(chǎn)生漂移誤差的問題,整個回波信號均會受到影響,噪聲增加。
最后,在第一回波與第二回波之間的飛行時間段內(nèi),理想情況下,信號為零,該段波形的大小即可反映噪聲信號的大小。
回波信號在不同接收段內(nèi)的波形受到環(huán)境的影響不盡相同,因此,提出將回波信號在時間域上的分布大致分為3個階段:第一回波段、第二回波段和二者之間的超聲波飛行時間段。
分區(qū)后的回波信號如圖3所示。
圖3 分區(qū)后的回波信號
具體分區(qū)由超聲波設備的激勵波周期數(shù)n,周期T,采樣頻率f確定。此時一個回波段內(nèi)的采樣點數(shù)N可由式(2)計算得到。
N=nTf
(2)
針對環(huán)境因素帶來的影響,對應在波形上體現(xiàn)出不同的特征。
(1)超聲波能量消耗增加或者激勵電壓降低引起的回波信號能量減弱,體現(xiàn)在回波信號上的特征即是兩次回波波峰的降低。
(2)壓電超聲設備中耦合效果變差引起的第二回波能量減少,第一回波能量增加,體現(xiàn)在回波信號上的特征即是第二回波波峰降低,第一回波波峰升高,因此可以使用二者的幅值之比,來反映耦合效果是否變差。
(3)壓電超聲設備中耦合效果變差引入的頻率噪聲、電磁超聲設備中電磁干擾引入的雜散噪聲以及電路元件漂移帶來的漂移噪聲,體現(xiàn)在整個回波信號上最易區(qū)分的即是TOF階段產(chǎn)生的雜散波。
對于不同的測厚系統(tǒng)而言,其本身的噪聲水平不同,此時所取的閾值也不同。對于壓電超聲測厚系統(tǒng),其干擾少,噪聲信號幾乎在波形上體現(xiàn)不出,所以取其幅值閾值為2 mV;而對于電磁超聲測厚系統(tǒng),由于其高壓電路與電路本身電磁干擾較大,噪聲信號明顯,所以取其幅值閾值為20 mV。
(4)最后,回波信號的信噪比即可由第一回波波峰與TOF階段波峰之比得到。由于單個區(qū)段內(nèi)存在若干尖峰,若只選用單個尖峰求取峰值容易產(chǎn)生偏差,因此取依峰值大小排序的前若干尖峰點的峰值平均值作為該區(qū)段的峰值。
對于實際應用的超聲波設備,由于選取的是周期數(shù)為2T的激勵[16],得到的回波信號中普遍存在有5個尖峰點,第一回波主尖峰點為3個,第二回波主尖峰點為2個,所以第一回波段和第二回波段分別取前3個和前2個尖峰值,再求取平均。對于TOF段,取該段內(nèi)至多兩個尖峰峰值的平均值。
在最簡化的情況下,回波信號提取4個特征值,由特征值引申計算出兩個引申值(見表1)。
表1 特征值及引申值列表
得到波形信號的特征值與引申值后,通過與系統(tǒng)初始工作條件下的特征值與引申值進行對比,求取變化量,來判斷當下系統(tǒng)的硬件健康狀況。
方法應用初期,提示維修保養(yǎng)的變化率閾值由經(jīng)驗值50%確定,同樣的,TOF段極值點個數(shù)變化量閾值定為5個,信噪比變化量閾值對于壓電與電磁系統(tǒng)分別定為15 dB與4 dB。在后期的持續(xù)應用中,該值應當根據(jù)維護設備的反饋作出調(diào)整。
下面采用具體的超聲波回波數(shù)據(jù)進行計算和處理過程的詳細說明。
(1)首先提取出四川東北某大型氣田壓電超聲設備3號探頭于2017年6月25日測量的序號為32的初始數(shù)據(jù),以及2018年4月8日測量的序號為678的長時間監(jiān)測后的數(shù)據(jù),并對兩個波形進行上述的分區(qū)操作。超聲波測厚儀安裝現(xiàn)場如圖4所示。
圖4 超聲波測厚儀安裝現(xiàn)場
圖5 軟件信號識別效果截圖
(2)進行幅值特征點的提取。首先,將整個信號進行絕對值處理,依次取出波形信號中的極值點,將所有極值點的位置信息放入數(shù)組pos()中;然后根據(jù)操作人員輸入的幅值門限,將pos()中的幅值低于該門限值的點的位置信息刪除,得到pos1(),得到如圖5所示的軟件自動生成的幅值特征點識別界面。時間序列單位為10 ns,幅值單位為mV。
圖5中,左側數(shù)據(jù)區(qū)域為軟件操作員輸入的想要調(diào)用的數(shù)據(jù)信息以及該臺設備在提取特征時的閾值參數(shù)。
(3)在每個分區(qū)內(nèi)選擇定量的幅值最大的紅色特征點,對其幅值絕對值求平均,得到3個分區(qū)中的幅值特征V1,VTOF,V2。并記錄TOF段高于閾值條件的特征點數(shù)目,得到數(shù)目特征NTOF。
(4)通過這4個特征值,計算出回波幅值比μ和信噪比SNR。
(5)對比兩次波形的特征值與引申值,根據(jù)變化量與變化率閾值,作出系統(tǒng)健康狀況診斷。通過信號的特征提取,得到軟件分析結果為:① 第一回波幅值減弱了20.9%,第二回波幅值減弱了33.3%;② 識別的干擾特征點數(shù)增加了1個,干擾信號幅值不變,抗干擾性能幾乎保持不變;③ 回波幅值比減小了16%,雖然有所下降,但基本與初始狀態(tài)保持一致;④ 信號的信噪比降低了2.042 dB;⑤ 總體來說,需要進行持續(xù)的觀察。
通過以上軟件分析結果可知:在初始安裝時,第二回波的幅值就非常小,因此其回波幅值比也非常小。相對于設備剛安裝時的初始狀態(tài),序號為678的數(shù)據(jù)無論是第一回波還是第二回波,幅值都減小了,這證明激勵電壓降低,但信噪比幾乎還保持在一個非常高的狀態(tài);相比于初始狀態(tài),回波幅值比減小了16%,證明第二回波的幅值減弱效果更為明顯,因此可以看出,耦合效果相對于初始時也出現(xiàn)了下降。綜合來說,信號依舊有非常好的信噪比,但是傳感器性能和耦合效果都有所下降。
圖6 軟件信號識別效果截圖
為了驗證系統(tǒng)軟件識別和提取的準確性,又提取了一組典型的含有干擾波的超聲波信號,圖6為在四川東北某油氣凈化廠的電磁超聲設備的2號探頭所測得的初始數(shù)據(jù)(2017年10月12日)和長時間監(jiān)測后的數(shù)據(jù)(2018年4月3日),其時間序列單位為10 ns,幅值單位為mV。軟件的特征分析結果為:① 第一回波幅值增加了9.7%,第二回波幅值減弱了7.4%;② 識別的TOF段干擾特征點數(shù)增加了8個,干擾幅值增加了1.75倍,抗干擾性能下降,需要考慮進行維修;③回波幅值比減弱了15.6%,雖然有所下降,但基本與初始狀態(tài)保持一致;④ 信號的信噪比降低了7.995 dB,需要考慮進行維修;⑤ 總體來說,設備噪聲過大,需要考慮進行維修。
基于電磁原理的超聲波信號,雖然第二回波幅值非常高,但是自始至終都有干擾信號存在,相對于初始時期,序號為10663的信號的兩個回波的幅值比雖然有所降低,但是依然大于0.7,減小了15.6%,證明信號回波的幅值并沒有減小多少,而相對于初始狀態(tài)信噪比卻嚴重降低,且TOF區(qū)識別的特征點數(shù)也增加了不少,因此可以看出,通過長時間的監(jiān)測,信號受到的干擾變強,且干擾信號的個數(shù)變多,因此從一定程度上也能說明系統(tǒng)的抗干擾性能變?nèi)酢?/p>
使用軟件提取采集波形的特征值的方法,識別與判斷硬件設備的健康狀況,節(jié)約了工程人員的勞動力,提高了監(jiān)管效率。對于長期在線監(jiān)測的超聲波設備,該方法在實際應用中,會定時采集管道波形,在此波形基礎上進行高效地采樣與反饋,便于現(xiàn)場操作員及時維修或者更換設備硬件,提高了設備在線監(jiān)測的準確度,并且系統(tǒng)改進成本低,實用性強。但是該方法在長期應用中也需要進一步研究與完善,即目前判定為噪聲的電壓閾值為2 mV與20 mV,信噪比變化量閾值為15 dB與4 dB,TOF段極值點個數(shù)變化量閾值為5個,以及峰值變化率閾值為50%,這些都是憑借工程人員經(jīng)驗而定的,在系統(tǒng)持續(xù)使用中,這些經(jīng)驗確定的閾值需要通過數(shù)據(jù)訓練以達到最優(yōu)的區(qū)分效果。
基于特征點的系統(tǒng)性能識別與判斷方法是人工智能在腐蝕監(jiān)測方面的應用,也是未來監(jiān)測維護課題的發(fā)展方向。理論分析與實踐表明,基于特征點的系統(tǒng)性能識別與判斷方法可以有效應用于超聲波的在線監(jiān)測中,具有較強的實用意義。