，， ，， ，

(1.中國石油大學(xué)(北京) 機(jī)械與儲運工程學(xué)院，北京 102249；2.中國石油天然氣股份有限公司 塔里木油田分公司，庫爾勒 841000)

由于磁記憶檢測技術(shù)對應(yīng)力集中敏感而成為鐵磁性構(gòu)件早期損傷檢測的重要方法。鐵磁性材料應(yīng)力和變形的集中區(qū)域會在地磁場和工作載荷下發(fā)生不可逆的磁疇組織重新排列，導(dǎo)致應(yīng)力集中區(qū)域產(chǎn)生磁場畸變，表面出現(xiàn)漏磁場，再通過對漏磁場的檢測來進(jìn)行應(yīng)力集中區(qū)域的診斷[1]。磁記憶現(xiàn)象的基本理論表明，采用求梯度等信號處理方法可找出法向分量零點或切向分量最大點[2-4]。隨著近幾年的發(fā)展，磁記憶檢測已經(jīng)由單通道檢測開始邁入多通道檢測，但多通道信號梯度的處理上依然延續(xù)以往的在檢測方向上求梯度的方法，易造成與檢測方向同向的應(yīng)力集中區(qū)域信號反應(yīng)不靈敏或漏檢的問題[5-7]。已有學(xué)者研究了采用切向和法向分量兩種信號聯(lián)合進(jìn)行檢測的方法，但這會增加對硬件系統(tǒng)的要求。筆者探討了一種將切向分量方向與法向分量方向聯(lián)合梯度求解的方法，某點的梯度由檢測方向上的點與垂直方向上等距的點共同求得。對鉆桿的磁記憶檢測信號進(jìn)行兩種梯度聯(lián)合求解的處理方法，比只從檢測方向求梯度的方法對畸變磁場更為靈敏。

1 現(xiàn)有的多通道磁記憶檢測信號梯度算法

現(xiàn)有多通道磁記憶檢測信號處理的過程一般為均值降噪、數(shù)據(jù)歸一化及求取檢測方向上的梯度。均值降噪是為了剔除信號曲線上的毛刺，即信號中夾雜的高頻環(huán)境噪聲，方便識別磁場畸變信號的突變點；數(shù)據(jù)歸一化是先求取所有傳感器整體信號均值，然后將每個傳感器信號序列都減去自身均值，即偏移到零軸附近，最后將各傳感器信號加上整體均值，以偏移回整體基線附近，達(dá)到消除傳感器差異的目的；梯度表征磁記憶信號的突變情況，梯度突變點可表征漏磁場的位置，梯度最大值位置可表征漏磁場變化最劇烈的位置。梯度的求解過程為兩點間磁場強(qiáng)度H的變化量與兩點間距L的比值，梯度的表達(dá)式如式(1)所示[8]。

(1)

式中：KN為N點的梯度；HN,HN+i分別為第N點和第N+i點的磁場強(qiáng)度；Lx，N,Lx，N+i分別為第N點和第N+i點的磁記憶信號的位置坐標(biāo)。

該梯度算法是在檢測方向上進(jìn)行單一的梯度求解(見圖1)，假設(shè)進(jìn)行n通道檢測，各通道每隔1 mm均采集1個數(shù)據(jù)點，間隔3 mm進(jìn)行梯度求解，則第1個通道的第1個點的梯度如式(2)所示，此后窗口依次滑動遍歷求取其他位置的梯度。

(2)

可見該梯度的求解過程只與本通道數(shù)據(jù)有關(guān)，與相鄰?fù)ǖ罒o任何關(guān)系。當(dāng)應(yīng)力集中區(qū)域沿檢測方向呈細(xì)條狀且足夠長時，該通道漏磁場信號均異常大，但只在該區(qū)域的端部產(chǎn)生梯度突變，無法對應(yīng)力集中全區(qū)域產(chǎn)生梯度異常響應(yīng)，不能很好地反映應(yīng)力集中區(qū)域。應(yīng)力集中區(qū)域與檢測方向同向的磁記憶檢測示意如圖2所示。

圖1 沿檢測方向求梯度算法示意

圖2 應(yīng)力集中區(qū)域與檢測方向同向的磁記憶檢測示意

2 多通道磁記憶檢測信號的多方向梯度算法

2.1 多方向梯度算法的原理

在多方向上進(jìn)行梯度求解，指的是每點的梯度求解與4周等距的4個相鄰點有關(guān)，而非只與檢測方向上的本通道信號相關(guān)。從檢測原理上分析，多方向的梯度值表征該點的漏磁場相對于相鄰等距位置的漏磁場的突變情況，不再僅是相對于一點的突變情況，因此從理論上避免了因掃描方向帶來的漏檢，此時的梯度如式(3)所示。

(3)

式中：K(i,j)為第j通道第i點的梯度值；H(i,j)為第j通道第i點的磁記憶信號；H(i,j-1)和H(i,j+1)分別為第j-1和第j+1通道第i點的磁記憶信號；H(i-k,j)和H(i+k,j) 分別為第j通道第i-k點和第i+k點的磁記憶信號，且這兩點與H(i,j)的距離等于H(i,j)與H(i,j-1)的距離，即4點與待求點距離均相等；Lx(i,j)為(i,j)點的位置坐標(biāo);Lx(i,j-1)為(i,j-1)點的位置坐標(biāo)。

從數(shù)學(xué)計算上分析，即通過與待求點上下左右等距的特定梯度的卷積核對信號進(jìn)行線性卷積。卷積公式[9-10]為

(4)

式中：f為待梯度處理信號；g為梯度卷積核;y為信號梯度卷積值。

2.2 多方向梯度算法的實現(xiàn)

對于n通道，數(shù)據(jù)長度為m，通道物理間隔為L，單通道采集頻率L/k的信號進(jìn)行數(shù)值遍歷的算法實現(xiàn)如下所述。

循環(huán)i從2到n-1，循環(huán)j從k到m。

(5)

循環(huán)j結(jié)束，循環(huán)i結(jié)束，循環(huán)j從k到m。

(6)

(7)

循環(huán)j結(jié)束。

對于圓柱形構(gòu)件，計算第1通道時，應(yīng)將第n通道和第2通道作為相鄰?fù)ǖ捞幚?計算第n通道時，應(yīng)將第n-1通道和第1通道作為相鄰?fù)ǖ肋M(jìn)行處理。對于板型構(gòu)件，計算第1通道和第n通道時，只需將相鄰空區(qū)域補(bǔ)0即可。

構(gòu)造梯度卷積核的方法如下，假設(shè)第n通道檢測，相鄰?fù)ǖ篱g隔距離為1個距離單位(這里不設(shè)具體單位)，各通道每0.5個距離單位均采集1個數(shù)據(jù)點，為與通道間隔等距，需選定檢測方向間隔為1進(jìn)行聯(lián)合求解，則梯度卷積核表達(dá)式為

(8)

此時的計算過程為梯度處理前矩陣和梯度卷積核矩陣對應(yīng)位相乘后求和，再除以相鄰?fù)ǖ赖奈锢砭嚯x，即得梯度處理前矩陣與梯度卷積核交叉部分中心的梯度，多方向求梯度算法示意如圖3所示。

圖3 磁記憶檢測信號的多方向求梯度算法示意

3 兩種梯度算法應(yīng)用對比

檢測裝置采用中國石油大學(xué)(北京)安全監(jiān)測與智能診斷實驗室研發(fā)的16通道鉆桿磁記憶檢測儀，沿周向均布的磁記憶傳感器環(huán)在電機(jī)驅(qū)動下運轉(zhuǎn)，傳感器貼合鉆桿表面進(jìn)行軸向勻速掃描，并將鉆桿表面磁記憶信號通過無線傳送給上位機(jī)，多通道鉆桿磁記憶檢測原理示意如圖4所示。

圖4 多通道鉆桿磁記憶檢測原理示意

檢測對象為加厚的表面存在由腐蝕坑組成軸向條形腐蝕區(qū)域的φ76.2 mm鉆桿，將傳感器環(huán)的第15號通道對準(zhǔn)該損傷區(qū)域進(jìn)行檢測(見圖5)。

圖5 鉆桿損傷區(qū)域?qū)嵨?/p>

將檢測信號進(jìn)行兩種梯度算法的處理，每種算法均對梯度進(jìn)行曲線繪制，結(jié)果如圖6，7所示。

圖6 磁記憶檢測信號的沿檢測單方向求梯度處理結(jié)果

圖7 磁記憶檢測信號的多方向求梯度處理結(jié)果

從圖6可以看出，對同一個檢測信號數(shù)據(jù)計算梯度，沿檢測單方向求梯度時，梯度只受單數(shù)據(jù)點的影響，在梯度圖上異常部位的曲線顯示并不明顯；從圖7可明顯看出，第15號通道信號因為腐蝕坑而劇烈跳動，這是因為該通道梯度求解過程中與相鄰的第14，16號通道正常信號相關(guān)。對比按照通道進(jìn)行排列的梯度瀑布圖可以看出，單方向求梯度時其他通道均存在低頻高賦值跳動，而多方向梯度瀑布圖只存在高頻低賦值跳動。從梯度瀑布圖第15號通道的信號對比可以看出，多方向梯度信號曲線跳動更劇烈，這是因為多方向多點求梯度的過程本身也是對噪聲信號的抑制以及對異常信號的放大。

4 結(jié)語

對多通道磁記憶信號梯度求解的算法原理和算法過程進(jìn)行分析，在多通道磁記憶梯度求解過程中采用多方向梯度算法，并通過構(gòu)造梯度卷積核進(jìn)行了鉆桿磁記憶信號切向分量多方向梯度求解的應(yīng)用，從理論和實際應(yīng)用上證明了多方向梯度算法比沿檢測單方向求梯度效果更好。給鉆具磁記憶檢測及圓柱形、板型構(gòu)件的多通道磁記憶檢測信號處理提供了一種參考。

[1] LI C, DONG L, WANG H. Metal magnetic memory technique used to predict the fatigue crack propagation behavior of 0.45%C steel[J]. Journal of Magnetism and Magnetic Materials, 2016, 405:150-157.

[2] DOUBOV A A. The method of metal magnetic memory-the new trend in engineering diagnostic[J]. Welding in the World, 2005, 49(9): 314-319.

[3] 任吉林, 林俊明, 任文堅,等. 金屬磁記憶檢測技術(shù)研究現(xiàn)狀與發(fā)展前景[J]. 無損檢測, 2012, 34(4): 3-11.

[4] 任吉林, 劉海超， 宋凱. 金屬磁記憶檢測技術(shù)的興起與發(fā)展[J]. 無損檢測, 2016,38(11) :7-15.

[5] 周成森. 基于梯度法的金屬磁記憶檢測儀[D].南京：南京航空航天大學(xué),2012.

[6] 欒明, 樊建春, 張來斌, 等. 一種新的多探頭金屬磁記憶檢測儀信號處理方法[J]. 無損檢測, 2011, 33(1): 54-57.

[7] 任吉林, 王進(jìn), 范振中, 等. 一種磁記憶檢測定量分析的新方法[J].儀器儀表學(xué)報,2010,31(2):431-436.

[8] 孫玉江,樊建春,祖強(qiáng),等. 基于磁記憶技術(shù)的鉆桿檢測系統(tǒng)研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù),2014,10(9):129-133.

[9] 胡廣書.數(shù)字信號處理——理論、算法與實現(xiàn)[M].北京:清華大學(xué)出版社,2003.

[10] 陳佳益,黃楠,熊剛強(qiáng)，等. 基于置信區(qū)間的自適應(yīng)加權(quán)均值濾波算法[J]. 南京理工大學(xué)學(xué)報,2017,41(3):307-312.