楊 杰 李紅梅 趙春田, 楊洪禮

基于相似波形的漏磁檢測(cè)缺陷信號(hào)小波分解及提取方法

楊 杰1李紅梅2趙春田1,2楊洪禮3

（1. 南方科技大學(xué)，廣東深圳 518055；2. 四川大學(xué)，成都 610207；3. 山東科技大學(xué)，山東 青島 266590）

漏磁檢測(cè)技術(shù)被廣泛應(yīng)用于鐵磁構(gòu)件缺陷的檢測(cè)和量化評(píng)估，其對(duì)缺陷形狀、尺寸的三維反演成像、量化能力，尤其是對(duì)常見的自然（腐蝕、疲勞）復(fù)雜缺陷的三維成像能力，是漏磁檢測(cè)技術(shù)水平的核心標(biāo)志之一。由于實(shí)際缺陷形狀不規(guī)則、常包含多個(gè)坑點(diǎn)，致使其漏磁檢測(cè)信號(hào)彼此影響。為提升復(fù)雜缺陷的成像精度和數(shù)據(jù)處理速度，從漏磁信號(hào)里自動(dòng)分解并剝離出對(duì)應(yīng)缺陷的信號(hào)是很重要的前置工作。本文提出一種基于相似波形的漏磁檢測(cè)缺陷信號(hào)小波分解及提取方法，用于自動(dòng)分解并拾取每個(gè)缺陷的漏磁檢測(cè)信號(hào)，并自動(dòng)確定每個(gè)缺陷的位置及邊界，以備后續(xù)利用其完成對(duì)缺陷的三維反演成像。該方法基于實(shí)際檢測(cè)到的缺陷漏磁信號(hào)的形態(tài)特征，匹配與其波形最相似的小波基；用選取的小波基對(duì)每個(gè)通道的漏磁檢測(cè)信號(hào)進(jìn)行多尺度分解，得到每個(gè)通道信號(hào)的小波高頻系數(shù)；再根據(jù)小波系數(shù)特征確定每個(gè)通道信號(hào)中存在的缺陷數(shù)量和缺陷所在的位置，并提取每個(gè)通道各處缺陷的漏磁信號(hào)；最終通過判定缺陷信號(hào)相鄰?fù)ǖ赖倪B通性，將所在連通區(qū)域內(nèi)的缺陷信號(hào)進(jìn)行集成，得到各處缺陷的信息集，從而實(shí)現(xiàn)對(duì)所有缺陷漏磁檢測(cè)信號(hào)的提取。本文提出的方法可用于對(duì)復(fù)雜缺陷漏磁檢測(cè)信號(hào)的準(zhǔn)確和快速提取、分離和集成，為后續(xù)缺陷的三維反演成像奠定基礎(chǔ)。

漏磁檢測(cè)；相似小波基；小波分解；缺陷信號(hào)提取

0 引言

近年來，管道運(yùn)輸以其安全、便捷、高效等優(yōu)點(diǎn)在石油、天然氣等能源運(yùn)輸行業(yè)被廣泛應(yīng)用。截止到2020年底，我國(guó)油氣長(zhǎng)輸管道里程累計(jì)達(dá)到16.9萬km。相關(guān)數(shù)據(jù)表明，我國(guó)超過半數(shù)的管道已經(jīng)進(jìn)入服役的中后期。為避免發(fā)生管道泄漏事故和造成重大的生命財(cái)產(chǎn)損失，對(duì)管道進(jìn)行有效的高精度完整性檢測(cè)是必要的。管道漏磁內(nèi)檢測(cè)法具有操作簡(jiǎn)便、可靠性好、檢測(cè)效果好等優(yōu)勢(shì)，被廣泛用于對(duì)油氣管道中缺陷的檢測(cè)和量化。

漏磁檢測(cè)[1-3]可實(shí)現(xiàn)對(duì)油氣管道內(nèi)、外表面的金屬損傷、焊縫等缺陷的有效無損檢測(cè)，其原理示意圖如圖1所示。在利用磁體對(duì)管壁進(jìn)行飽和磁化的狀態(tài)下，如圖1（a）所示，當(dāng)管壁中不存在缺陷時(shí)，磁通回路沒有擾動(dòng)，空氣中的磁敏傳感器拾取到的漏磁場(chǎng)不會(huì)發(fā)生畸變；如圖1（b）所示，當(dāng)管壁中存在缺陷時(shí)，管體中的磁力線在缺陷界面處發(fā)生折射，導(dǎo)致管體外空氣中漏磁場(chǎng)信號(hào)的畸變。通過分析檢測(cè)漏磁場(chǎng)信號(hào)是否有畸變，可實(shí)現(xiàn)對(duì)管體缺陷的識(shí)別；通過分析畸變漏磁場(chǎng)信號(hào)的形態(tài)、幅值等特征可實(shí)現(xiàn)對(duì)缺陷的粗略量化。在實(shí)際管道檢測(cè)工程中，經(jīng)常將沿管道長(zhǎng)度方向（軸向）設(shè)定為軸，管道周向方向設(shè)定成軸，管道徑向方向設(shè)置成軸（見圖1）。磁敏傳感器拾取的漏磁信號(hào)幅值為電壓值，單位為mV。

圖1 漏磁檢測(cè)原理示意圖

近年來，隨著我國(guó)裝備制造能力的大幅提升，國(guó)產(chǎn)漏磁檢測(cè)機(jī)器人的硬件裝備水平與國(guó)外的差距在快速縮小，但在缺陷的量化評(píng)定精度、速度和可靠性等方面，與國(guó)外相比仍存在很多不足。僅就缺陷的量化評(píng)定能力而言，目前國(guó)內(nèi)外的漏磁檢測(cè)技術(shù)仍處于對(duì)缺陷的定位查找及對(duì)其尺寸的粗略量化階段，對(duì)管道中缺陷形狀的三維量化反演成像能力仍然不足。本研究團(tuán)隊(duì)在利用實(shí)際檢測(cè)信號(hào)反演缺陷處磁荷分布規(guī)律的基礎(chǔ)上[4]，提出缺陷深度場(chǎng)概念，并建立缺陷三維輪廓反演的方法及算法，可實(shí)現(xiàn)對(duì)任意復(fù)雜缺陷的三維輪廓的反演成像[5]。以上理論創(chuàng)新技術(shù)已實(shí)現(xiàn)轉(zhuǎn)化，具有廣闊的工業(yè)應(yīng)用前景。進(jìn)一步地，為提升對(duì)缺陷的量化評(píng)定能力和處理速度，從漏磁檢測(cè)信號(hào)中自動(dòng)識(shí)別、分離和提取缺陷的漏磁信號(hào)是很重要的前置工作。由于實(shí)際工程缺陷形狀不規(guī)則，多個(gè)坑點(diǎn)信號(hào)常近距離分布、或一個(gè)區(qū)域包含多個(gè)坑點(diǎn)缺陷，致使其漏磁檢測(cè)信號(hào)彼此關(guān)聯(lián)、相互影響，降低了后期缺陷識(shí)別及成像量化精度。因此，為提升漏磁檢測(cè)對(duì)復(fù)雜缺陷的識(shí)別和量化成像精度，需要對(duì)上述近距離分布的多缺陷進(jìn)行分解剝離和單獨(dú)提取。目前應(yīng)用較廣泛的缺陷識(shí)別提取方法主要有幅值差分閾值法[6-7]，但其對(duì)近距離分布復(fù)雜缺陷的識(shí)別提取精度不高、信號(hào)提取自動(dòng)化程度低、對(duì)數(shù)據(jù)樣本和人工經(jīng)驗(yàn)的依賴性較高。即便是同一組數(shù)據(jù)，如果由不同的人員處理，其對(duì)缺陷信號(hào)的識(shí)別提取結(jié)果可能不同。因缺陷數(shù)據(jù)分析依賴人工、費(fèi)時(shí)費(fèi)力，導(dǎo)致數(shù)據(jù)分析結(jié)果常在檢測(cè)完成后2～4月才能給出，嚴(yán)重降低了檢測(cè)結(jié)果的時(shí)效性。

漏磁檢測(cè)缺陷信號(hào)的提取，尤其是對(duì)相互關(guān)聯(lián)的、跨越多通道的復(fù)雜缺陷信號(hào)的提取，目前還沒有可適用的高效、成熟算法，不利于后續(xù)的缺陷三維反演成像及量化評(píng)定。針對(duì)以上問題，本文給出一種基于相似波形的漏磁檢測(cè)缺陷信號(hào)小波分解及自動(dòng)提取方法，從缺陷漏磁檢測(cè)信號(hào)的波形特征出發(fā)，匹配與其波形相似度高的小波基，用其對(duì)漏磁檢測(cè)信號(hào)進(jìn)行多尺度分解，探究小波高頻系數(shù)與缺陷信號(hào)的對(duì)應(yīng)規(guī)律，從而實(shí)現(xiàn)對(duì)缺陷信號(hào)的自動(dòng)提取。該方法不僅能提高缺陷的提取精度和處理速度，也適用于對(duì)復(fù)雜缺陷信號(hào)的剝離和提取，為后續(xù)的缺陷三維反演及量化評(píng)定奠定基礎(chǔ)。

1 相似小波分解提取方法

1.1 小波多分辨率分析

小波分析[8-9]自20世紀(jì)70年代被提出，在信號(hào)分析及數(shù)學(xué)領(lǐng)域得到快速發(fā)展，很好地彌補(bǔ)了許多其他信號(hào)分析方法（如傅里葉分析方法）的不足，在理論和實(shí)際中都得到了廣泛應(yīng)用。相較于傅里葉分析方法，小波分析方法的一個(gè)顯著優(yōu)勢(shì)是可以對(duì)信號(hào)進(jìn)行局部分析。1986年，Meyer給出了一組小波，對(duì)其進(jìn)行伸縮和平移可構(gòu)成一組空間序列的標(biāo)準(zhǔn)化正交基，顯著推動(dòng)了小波的實(shí)際應(yīng)用。1988年，Mallat在做小波構(gòu)造研究時(shí)提出多分辨率分析[10-13]的理念，從數(shù)學(xué)定義和空間概念對(duì)正交小波[11]都進(jìn)行了形象的解釋和描述，并提出目前在信號(hào)去噪領(lǐng)域中廣泛采用的算法——Mallat算法[14]，使小波分析方法在理論和應(yīng)用層面都實(shí)現(xiàn)了跨越。在小波的計(jì)算可行性問題被解決后，它的獨(dú)特優(yōu)勢(shì)，例如可檢測(cè)信號(hào)的突變點(diǎn)和奇異性等局部信息的能力，便顯著地發(fā)揮了出來，解決了信號(hào)分析領(lǐng)域的諸多難題。

信號(hào)經(jīng)小波分解后映射到不同尺度的函數(shù)空間，并得到不同尺度下相互正交的時(shí)頻系數(shù)。不同的尺度對(duì)應(yīng)不同的頻率范圍，即可實(shí)現(xiàn)在不同頻率范圍內(nèi)觀察信號(hào)，并獲取其多頻域的時(shí)頻信息。另外，小波的可構(gòu)造性極大提高了小波分析在實(shí)際應(yīng)用中的靈活性和便捷性。

1.2 漏磁檢測(cè)缺陷信號(hào)的小波提取方法

基于相似波形小波分解的漏磁檢測(cè)缺陷信號(hào)提取技術(shù)路線如圖2所示，本文對(duì)漏磁檢測(cè)缺陷信號(hào)的提取可分解為五步：①對(duì)檢測(cè)信號(hào)去噪；②選取最優(yōu)小波基；③對(duì)檢測(cè)信號(hào)進(jìn)行小波分解；④初步鎖定缺陷區(qū)域；⑤精確提取各缺陷信號(hào)。

圖2 基于相似波形小波分解的漏磁檢測(cè)缺陷信號(hào)提取技術(shù)路線

1）小波基的選擇

由于皮爾遜系數(shù)是對(duì)稱的，兩個(gè)向量的平移或尺度變化都不會(huì)影響皮爾遜系數(shù)值的變化。兩個(gè)向量間的皮爾遜系數(shù)及其所對(duì)應(yīng)的相關(guān)性關(guān)系見表1。

表1 Pearson系數(shù)與兩個(gè)向量的相關(guān)性

由此，計(jì)算不同的小波基與漏磁檢測(cè)缺陷信號(hào)的Pearson系數(shù)，選取Pearson系數(shù)值高的小波基作為對(duì)漏磁檢測(cè)缺陷信號(hào)分解用的小波函數(shù)，可提升對(duì)單個(gè)缺陷識(shí)別拾取的精準(zhǔn)度。以工程中實(shí)際獲取的漏磁檢測(cè)缺陷信號(hào)的徑向分量為基準(zhǔn)，本文計(jì)算了其與典型的幾種小波基函數(shù)的Pearson系數(shù)，見表2。Rbio3.1小波與漏磁檢測(cè)缺陷信號(hào)徑向分量的Pearson系數(shù)為0.964 9，相似度較高，兩者的波形相似性如圖3所示。

2）缺陷信號(hào)提取

根據(jù)波形相似原則，采用Rbio3.1小波基對(duì)漏磁檢測(cè)信號(hào)進(jìn)行小波多尺度分解。

表2 不同小波基與漏磁檢測(cè)缺陷信號(hào)的Pearson系數(shù)

圖3 Rbio3.1小波與缺陷信號(hào)徑向分量的波形相似性

（1）缺陷信號(hào)的單通道分解及提取

用Rbio3.1小波對(duì)如圖4（a）所示的含有單個(gè)缺陷的管道漏磁檢測(cè)單通道信號(hào)做4級(jí)分解，結(jié)果如圖4（b）和圖4（c）所示。圖4中，CA和CL（=1, 2, 3, 4）分別表示小波低頻和高頻系數(shù)值。由圖4結(jié)果可見，高頻系數(shù)CL可較好地拾取并表征缺陷信號(hào)。在無缺陷區(qū)域，CL接近于0；在有缺陷區(qū)域，CL存在1個(gè)極小負(fù)峰、2個(gè)極大正峰，單個(gè)缺陷信號(hào)經(jīng)小波分解后得到的高頻系數(shù)的分布規(guī)律具有惟一性，沒有衍生出多余旁瓣峰，可據(jù)此對(duì)單個(gè)缺陷進(jìn)行識(shí)別，有利于后續(xù)對(duì)復(fù)雜缺陷的分解和拾取。本文中，選取CL2對(duì)缺陷進(jìn)行識(shí)別和提取，尺度選取原則為：

①該尺度分解得到的小波低頻系數(shù)與原信號(hào)波形相似。

②該尺度分解得到的小波高頻系數(shù)，在有缺陷部分能量集中、無旁瓣，在無缺陷部分的值接近于0。

③在滿足①②條件下選低尺度。依據(jù)小波多尺度分解的模極大值傳遞規(guī)律[21]，小波分解尺度越低，極大值點(diǎn)傳遞時(shí)誤差越小。

依據(jù)小波變換基本理論，具有顯著特征波形的信號(hào)，經(jīng)小波分解后得到的系數(shù)的模極大值點(diǎn)，對(duì)應(yīng)原信號(hào)的特征點(diǎn)。CL2與缺陷特征點(diǎn)對(duì)應(yīng)關(guān)系如圖5所示，用Rbio3.1小波對(duì)圖5（a）所示的管道漏磁檢測(cè)缺陷信號(hào)進(jìn)行2級(jí)分解，得到的CL2如圖5（b）所示。圖5（b）中標(biāo)注的模極大值點(diǎn)與圖5（a）中標(biāo)注的特征點(diǎn)從左到右存在一一對(duì)應(yīng)的關(guān)系。因此，可先通過對(duì)CL2模極大值的拾取實(shí)現(xiàn)對(duì)缺陷及缺陷所在邊界位置的確定，再對(duì)應(yīng)出其在漏磁檢測(cè)信號(hào)中的位置，從而可從原漏磁檢測(cè)信號(hào)中截取出缺陷信號(hào)。

圖4 對(duì)信號(hào)的小波多尺度分解結(jié)果

（a）管道漏磁檢測(cè)缺陷信號(hào)

（b）信號(hào)經(jīng)分解后的二尺度小波高頻系數(shù)

圖5 CL2與缺陷特征點(diǎn)對(duì)應(yīng)關(guān)系

綜上所述，一個(gè)具有明確特征波形的信號(hào)，選取與其特征波形相似度大的小波基，對(duì)其進(jìn)行小波分解，得到與原信號(hào)對(duì)應(yīng)的小波低頻系數(shù)和高頻系數(shù)。其中，低頻系數(shù)波形與原信號(hào)波形相似，高頻系數(shù)放大表征原信號(hào)中的局部細(xì)節(jié)特征信息，且高頻系數(shù)的模大值點(diǎn)與原信號(hào)的特征點(diǎn)具有惟一映射關(guān)系。由此，先在小波高頻系數(shù)中根據(jù)特征提取出缺陷信號(hào)，后根據(jù)其映射關(guān)系提取出其在原信號(hào)中的位置信息，從而完成缺陷漏磁信號(hào)的提取。該方法的優(yōu)越性表現(xiàn)在，基于相似波形的小波多分辨率分析對(duì)特征信號(hào)的辨認(rèn)有惟一性，可提高復(fù)雜缺陷區(qū)域中對(duì)單個(gè)缺陷的識(shí)別率，進(jìn)而提高缺陷的提取精度。

漏磁檢測(cè)單通道缺陷信號(hào)提取的方法步驟如下：

①對(duì)檢測(cè)信號(hào)去噪。用基于Mallet算法的小波分解與重構(gòu)法實(shí)現(xiàn)。

②小波分解并確定閾值。用Rbio3.1小波對(duì)漏磁檢測(cè)信號(hào)進(jìn)行二級(jí)分解，根據(jù)分解后得到的二尺度高頻系數(shù)（CL2）確定閾值。經(jīng)驗(yàn)證，非缺陷信號(hào)對(duì)應(yīng)的CL2模值不超過缺陷信號(hào)對(duì)應(yīng)的CL2模值的30%。由此，找到分解后得到的CL2模的最大值max，將max（∈[0.1, 0.3]）作為閾值去篩取缺陷信息。

③拾取缺陷數(shù)據(jù)。選定步驟②中的，在步驟②得到的CL2中，將模值大于max的數(shù)據(jù)提取出來。

④缺陷數(shù)據(jù)分段。在步驟③得到的缺陷信息中，找到其中的局部極小值點(diǎn)，以極小值點(diǎn)為基準(zhǔn)對(duì)缺陷信息分段。

⑤確定每個(gè)缺陷信號(hào)段的邊界并截取。以每個(gè)局部極小值點(diǎn)所在位置為基準(zhǔn)，向左擴(kuò)散找到與其相鄰的局部極大值點(diǎn)，將其所在位置作為缺陷左邊界，將局部極小值點(diǎn)的位置作為缺陷右邊界。為減小誤差，左右邊界分別向左和向右各擴(kuò)散兩個(gè)單位，根據(jù)模極大值傳遞規(guī)律，將這些位置坐標(biāo)乘以4對(duì)應(yīng)到原信號(hào)中的位置，截取出原信號(hào)中邊界范圍內(nèi)的信號(hào)數(shù)據(jù)。

單通道缺陷信號(hào)提取流程如圖6所示。

圖6 單通道缺陷信號(hào)提取流程

（2）平面缺陷信號(hào)的分解和提取

基于一維單通道信號(hào)（缺陷長(zhǎng)度方向）的缺陷提取結(jié)果，缺陷在長(zhǎng)、寬平面內(nèi)的二維完整信號(hào)的提取，主要通過對(duì)缺陷邊界位置的連通性判斷來實(shí)現(xiàn)。長(zhǎng)、寬平面缺陷信號(hào)提取流程如圖7所示，包含以下主要步驟：

圖7 長(zhǎng)、寬平面缺陷信號(hào)提取流程

①劃分缺陷區(qū)域。按照?qǐng)D6所示缺陷數(shù)據(jù)拾取流程，對(duì)平面漏磁信號(hào)逐通道進(jìn)行缺陷識(shí)別，獲取所有存在缺陷的通道。將其中具有連通性的通道劃定為一個(gè)初級(jí)缺陷區(qū)域。

②在每個(gè)初級(jí)缺陷區(qū)域進(jìn)行單個(gè)缺陷的剝離提取，主要步驟如下：

a）單個(gè)缺陷的提取。按照?qǐng)D6所示流程依次在步驟①劃定的初級(jí)缺陷區(qū)域內(nèi)進(jìn)行每個(gè)通道內(nèi)含獨(dú)立缺陷的劃分和提取，并記錄每個(gè)獨(dú)立缺陷的位置和長(zhǎng)度。

b）獨(dú)立缺陷長(zhǎng)度的確定。在含缺陷數(shù)目最多的通道中，取各缺陷在不同通道中長(zhǎng)度的最大值作為其長(zhǎng)度，并記錄其所在的邊界位置。

c）獨(dú)立缺陷寬度的確定。在步驟b）中確定的各缺陷的長(zhǎng)度邊界范圍內(nèi)，獲取存在缺陷的通道，將其中具有連通性的通道集成為一個(gè)缺陷，其所包含的通道邊界表征了缺陷的寬度，記錄其位置。

d）從原漏磁檢測(cè)信號(hào)中，截取出由以上步驟b）和c）中確定的缺陷長(zhǎng)度和寬度邊界位置范圍內(nèi)的漏磁檢測(cè)數(shù)據(jù)，存儲(chǔ)起來，完成單個(gè)缺陷的提取。

其中，本文將相鄰間距小于等于2個(gè)采樣單位的通道認(rèn)為是具有連通性的通道。

2 缺陷提取結(jié)果驗(yàn)證

2.1 一維漏磁檢測(cè)缺陷信號(hào)提取

為了驗(yàn)證上述漏磁檢測(cè)缺陷信號(hào)的單通道分解與提取方法的可行性和準(zhǔn)確度，用Matlab對(duì)實(shí)際管道漏磁檢測(cè)的工程數(shù)據(jù)進(jìn)行了缺陷提取。選取其中存在缺陷的300條單通道信號(hào)進(jìn)行了缺陷提取驗(yàn)證，提取閾值系數(shù)設(shè)定為0.1。運(yùn)行結(jié)果顯示，本方法對(duì)其中缺陷信號(hào)的有效識(shí)別和提取率為98%，圖8為其中的兩組缺陷提取效果。圖8（a）為內(nèi)含單個(gè)缺陷，圖8（b）為內(nèi)含3個(gè)相鄰缺陷。圖8中，origin標(biāo)識(shí)線為原漏磁檢測(cè)信號(hào)，01、02、03標(biāo)識(shí)線分別為所提取的缺陷編號(hào)。由圖可見，本方法可對(duì)相鄰的多個(gè)缺陷進(jìn)行有效剝離和單獨(dú)提取，且被提取的缺陷信號(hào)完整、無冗余信息。

圖8 單通道缺陷信號(hào)的提取結(jié)果

2.2 二維平面缺陷信號(hào)的提取

為了驗(yàn)證上述漏磁檢測(cè)缺陷平面信號(hào)提取方法的有效性，用Matlab對(duì)實(shí)際工程的管道漏磁檢測(cè)信號(hào)中的缺陷二維信號(hào)進(jìn)行了提取，將提取閾值系數(shù)設(shè)定為0.1。第一批測(cè)試數(shù)據(jù)內(nèi)含20個(gè)缺陷，本方法自動(dòng)識(shí)別和有效提取了其中的19個(gè)缺陷信號(hào)，各缺陷信號(hào)都被單獨(dú)且完整地提取和保存。第二批測(cè)試數(shù)據(jù)內(nèi)含50個(gè)缺陷，本方法有效識(shí)別和提取了其中的48個(gè)缺陷信號(hào)，有效提取率為96%。

圖9（a）為內(nèi)含3個(gè)密集缺陷的仿真漏磁檢測(cè)信號(hào)。缺陷信號(hào)彼此影響，所跨越的通道范圍相互覆蓋。缺陷提取結(jié)果如圖9（b）～圖9（d）所示，3個(gè)密集缺陷都能夠被單獨(dú)識(shí)別并剝離提取出來。

圖9 多個(gè)缺陷的二維信號(hào)提取結(jié)果

圖10（a）為實(shí)際工程管道的漏磁檢測(cè)信號(hào)1，內(nèi)含一個(gè)長(zhǎng)軸向缺陷。圖10（b）為缺陷信號(hào)的提取結(jié)果，可見缺陷信息被完整提取且無冗余信息。圖11（a）為實(shí)際工程管道的漏磁檢測(cè)信號(hào)2，其內(nèi)含3個(gè)形狀不規(guī)則、相互影響的密集缺陷，圖11（b）～圖11（d）為其缺陷提取結(jié)果，可見3個(gè)缺陷的單獨(dú)識(shí)別和提取效果良好。

圖10 工程管道漏磁檢測(cè)信號(hào)1的二維缺陷提取結(jié)果

3 結(jié)論

本文基于小波多分辨率分析方法，給出了一種漏磁檢測(cè)缺陷信號(hào)提取方法。本方法基于相似波形小波分解思想，經(jīng)過波形相似度計(jì)算，選取了與缺陷漏磁信號(hào)相似度高的Rbio3.1小波，并應(yīng)用其對(duì)實(shí)際漏磁檢測(cè)信號(hào)進(jìn)行了分解和提取。測(cè)試結(jié)果表明，本文方法可有效從內(nèi)含多個(gè)密集缺陷的信號(hào)區(qū)域，分離和提取出單個(gè)獨(dú)立的缺陷信號(hào)，并可精確地確定各缺陷所在的位置和邊界，實(shí)現(xiàn)對(duì)相互影響的復(fù)雜缺陷進(jìn)行分別識(shí)別和單獨(dú)提取。本文方法的創(chuàng)新點(diǎn)及優(yōu)勢(shì)如下：

圖11 工程管道漏磁檢測(cè)信號(hào)2的二維缺陷提取結(jié)果

1）通過選取與缺陷信號(hào)波形相似度高的小波基函數(shù)對(duì)漏磁信號(hào)進(jìn)行小波分解，建立了小波高頻系數(shù)與單個(gè)缺陷信號(hào)的映射關(guān)系，并基于此進(jìn)行缺陷識(shí)別和提取，解決了漏磁檢測(cè)對(duì)密集型復(fù)雜缺陷的單獨(dú)提取難題。

2）開發(fā)了漏磁檢測(cè)缺陷提取算法，編寫了相應(yīng)的計(jì)算程序，從實(shí)際工程管道的漏磁檢測(cè)數(shù)據(jù)里，提取出了96%以上的二維缺陷信號(hào)，實(shí)現(xiàn)了缺陷提取過程的全程自動(dòng)化。

相對(duì)于現(xiàn)有方法，本文方法提高了對(duì)密集型復(fù)雜缺陷的分解和單獨(dú)識(shí)別能力，提升了缺陷提取的準(zhǔn)確率和速度，為隨后的缺陷輪廓反演成像計(jì)算奠定了基礎(chǔ)，在實(shí)際管道漏磁檢測(cè)工程中有廣泛的應(yīng)用前景。

[1] 孟祥吉, 宋兵臣, 劉健, 等. 管道部件及典型缺陷漏磁內(nèi)檢測(cè)圖像化顯示研究[J]. 管道技術(shù)與設(shè)備, 2021(1): 26-32.

[2] 戴光, 吳忠義, 朱祥軍, 等. 管道內(nèi)外壁缺陷的漏磁檢測(cè)[J]. 無損檢測(cè), 2018, 40(3): 19-23, 28.

[3] 楊理踐, 耿浩, 高松巍. 長(zhǎng)輸油氣管道漏磁內(nèi)檢測(cè)技術(shù)[J]. 儀器儀表學(xué)報(bào), 2016, 37(8): 1736-1746.

[4] LI Hongmei, ZHANG Fuchen, YANG Bin, et al. Dis- tribution characteristics of calculated magnetic charges around discontinuous structures in magnetic memory testing[J]. International Journal of Applied Electro- magnetics and Mechanics, 2019, 59(4): 1321-1329.

[5] LI Hongmei, HUANG Ranran, ZHAO Chuntian, et al. 3D reconstructing of arbitrary defects with magnetic flux leakage testing signals[C]//2020 IEEE Far East NDT New Technology and Application Forum (FENDT), Kunming, 2020: 51-55.

[6] 王婷婷. 金屬表面缺陷特征智能提取及特征分析的方法研究[D]. 沈陽: 東北大學(xué), 2017.

[7] LI Min, LI Xue, GAO Chenxing, et al. Acoustic microscopy signal processing method for detecting near-surface defects in metal materials[J]. NDT and E International, 2019, 103: 133-140.

[8] MALLAT S G. Multiresolution approximations and wavelet orthonormal bases of L2(R)[J]. Transactions of the American Mathematical Society, 1989, 315(1): 69-87.

[9] CHUI C K. An introduction to wavelets[M]. American: Academic Press Inc., 1992.

[10] 周志昊. 基于小波變換的電力變壓器直流偏磁時(shí)勵(lì)磁電流研究[J]. 電氣技術(shù), 2020, 21(6): 69-72.

[11] 胡昌華, 張軍波, 夏軍, 等. 基于MATLAB的系統(tǒng)分析與設(shè)計(jì)—小波分析[M]. 西安: 西安電子科技大學(xué)出版社, 2000.

[12] 楊福生. 小波變換的工程分析與應(yīng)用[M]. 北京: 北京科學(xué)出版社, 2000.

[13] 喬蘇朋, 楊艷, 陳世群, 等. 光伏陣列故障檢測(cè)方法綜述[J]. 電氣技術(shù), 2021, 22(7): 1-6.

[14] 米正英, 王瑜, 王立東, 等. 電弧信號(hào)小波去噪分析[J]. 電氣技術(shù), 2020, 21(3): 103-107, 116.

[15] 肖迎群, 馮良貴, 何怡剛. 基于小波分形和核判別分析的模擬電路故障診斷[J]. 電工技術(shù)學(xué)報(bào), 2012, 27(8): 230-238.

[16] 武佳卉, 邵振國(guó), 楊少華, 等. 數(shù)據(jù)清洗在新能源功率預(yù)測(cè)中的研究綜述和展望[J]. 電氣技術(shù), 2020, 21(11): 1-6.

[17] 崔芮華, 曹歡. 基于相空間重構(gòu)的航空電弧故障識(shí)別方法[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(增刊1): 243- 250.

[18] 鐘建林, 何友, 任獻(xiàn)彬. 基于波形相似度的容差模擬電路軟故障診斷[J]. 電工技術(shù)學(xué)報(bào), 2012, 27(8): 222-229.

[19] 彭海. 皮爾遜相關(guān)系數(shù)應(yīng)用于醫(yī)學(xué)信號(hào)相關(guān)度測(cè)量[J]. 電子世界, 2017(7): 163.

[20] 陳世群, 高偉, 陳孝琪, 等. 一種基于極限學(xué)習(xí)機(jī)和皮爾遜相關(guān)系數(shù)的光伏陣列故障快速診斷方法[J].電氣技術(shù), 2021, 22(10): 57-64.

[21] 朱延功, 高學(xué)山, 劉嵩, 等. 使用小波分析方法提取焊縫位置信息[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2001, 33(3): 389-392.

Wavelet decomposition and extraction method of defect magnetic flux leakage testing signals based on waveform similarity

YANG Jie1LI Hongmei2ZHAO Chuntian1,2YANG Hongli3

(1. Southern University of Science and Technology, Shenzhen, Guangdong 518055; 2. Sichuan University, Chengdu 610207;3. Shandong University of Science and Technology, Qingdao, Shandong 266590)

Magnetic flux leakage testing (MFLT) technology has been widely used in engineering to identify and measure defects in ferromagnetic components. Its capability in reconstructing defects in the three-dimensional (3D) form, especially for general complex defects formed naturally (such as by corrosion and/or fatigue), is one of the primary indexes representing the technological advances. Because the actual defect shape is irregular and often contains multiple pits, its magnetic flux leakage detection signals affect each other. To improve the accuracy of reconstructing defects in 3D and speed up the data processing, the prerequisite is to automatically strip off and extract the valid defect data from the MFLT signals. To effectively locate, separate, and integrate information for reconstructing each defect, a wavelet decomposition and extraction method is proposed. Based on the morphological characteristics of the MFLT signals, a wavelet basis is selected that matches the characteristics of the defect signals closely. This wavelet basis is used in the multi-scale decomposing of the MFLT signals and in the calculating of the associated wavelet high frequency coefficients. Then the wavelet coefficient is utilized in identifying and locating the defects in each channel of the signals, as well as in extracting the MFLT data corresponding to each defect. Finally, by examining the signals of each defect in a channel and integrating all the associated data in the adjacent channels, a complete set of data related to each defect is obtained. The proposed method can be applied to quickly identify, extract and integrate the MFLT data for all defects, which provides the basis for the 3D defect reconstruction.

magnetic flux leakage testing; similar wavelet basis; wavelet decomposition; defect signal extract

2022-01-06

2022-02-26

楊 杰（1997—），女，山西省大同市人，碩士研究生，主要從事數(shù)據(jù)分析與缺陷檢測(cè)方面的研究工作。

深圳市科技研發(fā)資金（JCYJ20190809144809345）

四川大學(xué)人才項(xiàng)目（YJ202091）

中央高?；究蒲谢?/p>