王宏安，陳國(guó)明

(1.中國(guó)石油大學(xué)(華東)海洋油氣裝備與安全技術(shù)研究中心，青島 266580；2.中石化勝利石油工程有限公司鉆井工藝研究院，東營(yíng) 257017)

油氣輸送管道作為石油、天然氣和成品油最主要、安全有效的運(yùn)輸方式，在服役過(guò)程中受到介質(zhì)腐蝕、地質(zhì)環(huán)境、人為破壞等因素影響產(chǎn)生各種各樣缺陷，嚴(yán)重時(shí)可能發(fā)生管道泄漏，造成經(jīng)濟(jì)損失和環(huán)境破壞。目前常用的管道檢測(cè)方法包括漏磁[1-2]、超聲[3]、渦流[4]等，漏磁檢測(cè)由于具有自動(dòng)化程度高、不需要耦合劑、管道適應(yīng)能力強(qiáng)等特點(diǎn)，成為油氣管道檢測(cè)最常用的方法之一，得到廣泛應(yīng)用[5-6]。

漏磁內(nèi)檢測(cè)完成后，數(shù)據(jù)判讀流程主要包括數(shù)據(jù)預(yù)處理、管段劃分、缺陷區(qū)域標(biāo)記、缺陷重構(gòu)和缺陷評(píng)估等部分[7]。工程管道漏磁檢測(cè)數(shù)據(jù)量巨大，需研究快速?gòu)穆┐艛?shù)據(jù)中實(shí)現(xiàn)缺陷標(biāo)記，尋找出感興趣的區(qū)域，為缺陷重構(gòu)與評(píng)估提供基礎(chǔ)數(shù)據(jù)。目前缺陷區(qū)域標(biāo)記方法主要有閾值法、平均值幅值法、分段識(shí)別法、多尺度窗缺陷檢測(cè)法等[8-11]。

目前的方法都需要根據(jù)檢測(cè)器、待檢測(cè)管道和管道磁化情況設(shè)置不同的閾值，算法根據(jù)閾值對(duì)缺陷區(qū)域進(jìn)行判定。閾值選擇往往需要較多工程經(jīng)驗(yàn)，而且對(duì)于多種壁厚的工程管道，需要多次設(shè)置閾值，總起來(lái)說(shuō)合理選擇閾值對(duì)缺陷區(qū)域標(biāo)記效果影響較大。

針對(duì)漏磁數(shù)據(jù)缺陷區(qū)域標(biāo)記的問(wèn)題，提出了一種基于聚類法的油氣管道缺陷區(qū)域標(biāo)記方法，由于該方法不需要閾值選擇，可以對(duì)不同檢測(cè)器和不同壁厚數(shù)據(jù)進(jìn)行適應(yīng)，具有較廣泛的適應(yīng)性。

1 算法簡(jiǎn)介

1.1 K均值聚類方法

聚類是將數(shù)據(jù)集中的樣本劃分為若干個(gè)通常是不相交的子集，每個(gè)子集稱之為簇，通過(guò)這樣的劃分，每一個(gè)簇內(nèi)數(shù)據(jù)十分相似，不同簇內(nèi)的數(shù)據(jù)又非常不同，每個(gè)簇可能對(duì)應(yīng)于一些潛在的概念[12]。

K均值聚類是一種基于距離的聚類算法，是最常用的聚類算法之一，具有簡(jiǎn)單、快速的特點(diǎn)，被廣泛使用。對(duì)于給定樣本集X={x1,x2,…,xn}，其中n為樣本個(gè)數(shù)。根據(jù)樣本的某種屬性特征將樣本集劃分為k個(gè)簇，然后隨機(jī)選擇k個(gè)數(shù)據(jù)點(diǎn)作為的初始聚類中心，計(jì)算樣本到初始聚類中心的距離將其分配給最近的類別，重新劃分k個(gè)聚類中心，不斷迭代分配整個(gè)過(guò)程，直到聚類中心基本不變。K均值聚類方法的目標(biāo)是簇間樣本距離最大化，簇內(nèi)樣本距離最小化，最終以簇內(nèi)方差作為聚類準(zhǔn)則函數(shù)，方差越小，聚類效果越好[13]。

1.2 基于K均值聚類的缺陷區(qū)域標(biāo)記方法

油氣管道漏磁內(nèi)檢測(cè)的數(shù)據(jù)通常不能直接進(jìn)行缺陷區(qū)域標(biāo)記，一般需要進(jìn)行數(shù)據(jù)預(yù)處理、管段劃分后才進(jìn)行缺陷區(qū)域的標(biāo)記。

本文方法主要針對(duì)管段劃分完成后，將數(shù)據(jù)分塊，分別提取數(shù)據(jù)塊的特征，根據(jù)這些特征采用K均值聚類的方法，將數(shù)據(jù)塊分類為正常數(shù)據(jù)塊和缺陷數(shù)據(jù)塊兩個(gè)類別，對(duì)缺陷數(shù)據(jù)塊進(jìn)行標(biāo)記，實(shí)現(xiàn)管段內(nèi)缺陷區(qū)域標(biāo)記。

基于K均值聚類的缺陷區(qū)域標(biāo)記方法主要步驟如下。

步驟1 讀取原始管段數(shù)據(jù)，如圖1(a)所示。

步驟2 按照一定長(zhǎng)寬將原始數(shù)據(jù)分塊。

漏磁檢測(cè)的原理決定了對(duì)于小缺陷檢測(cè)能力有限，以長(zhǎng)、寬、深為5 mm×5 mm×1 mm缺陷為檢測(cè)最小缺陷。根據(jù)漏磁檢測(cè)原理，5 mm×5 mm×1 mm缺陷產(chǎn)生的漏磁信號(hào)在長(zhǎng)度和寬度信號(hào)影響范圍大于10 mm。因此，將數(shù)據(jù)塊大小設(shè)置為10 mm×10 mm，分塊后的結(jié)果如圖1(b)所示。

步驟3 計(jì)算數(shù)據(jù)塊特征，并采用K均值聚類將數(shù)據(jù)分成2個(gè)類別。

對(duì)于缺陷區(qū)域在邊界處漏磁信號(hào)會(huì)出現(xiàn)較大的變化，最終反映在數(shù)據(jù)梯度上發(fā)生波動(dòng)，因此分別提取每個(gè)數(shù)據(jù)塊的軸向(檢測(cè)器前進(jìn)方向)梯度最大值、周向(傳感器分布方向)梯度的最大值、軸向梯度最大值與最小值的差值、周向梯度的最大值與最小值的差值4個(gè)值作為特征，并將結(jié)果賦值到數(shù)據(jù)塊矩陣中。

對(duì)所有數(shù)據(jù)塊設(shè)置兩個(gè)聚類中心，分別計(jì)算不同數(shù)據(jù)塊特征的歐式距離，迭代計(jì)算直到聚類中心不再發(fā)生變化，將所有數(shù)據(jù)塊分成了兩個(gè)類別，輸出數(shù)據(jù)塊的分類結(jié)果，如圖1(c)所示。

步驟4 根據(jù)所占比例判斷數(shù)據(jù)塊的是否是缺陷數(shù)據(jù)。

對(duì)于一個(gè)存在缺陷的管段，缺陷漏磁數(shù)據(jù)在整體數(shù)據(jù)中所占比例是較小的，通過(guò)對(duì)大量已有檢測(cè)器檢測(cè)數(shù)據(jù)分析發(fā)現(xiàn)，整體比例都低于10%。因此在這里假定屬于缺陷類別的數(shù)據(jù)塊不超過(guò)總數(shù)據(jù)塊的10%。于是，根據(jù)步驟3得到分類結(jié)果，判斷數(shù)據(jù)塊數(shù)量較少的類別是否超過(guò)了總數(shù)據(jù)塊的10%，如果超過(guò)了則判斷兩個(gè)類別均為無(wú)缺陷區(qū)域，如果未超過(guò)，則判斷數(shù)據(jù)塊數(shù)量小的類別為缺陷類別，并在數(shù)據(jù)塊矩陣中進(jìn)行標(biāo)記。

步驟5 對(duì)矩陣中缺陷區(qū)域擴(kuò)展和合并，如圖1(d)所示。

步驟6 提取缺陷區(qū)域的邊界值，形成矩形框，如圖1(e)所示，記錄矩形框位置信息，為后續(xù)缺陷重構(gòu)提供基礎(chǔ)數(shù)據(jù)。最終的缺陷框在原始的數(shù)據(jù)中的標(biāo)記結(jié)果如圖1(f)所示。

圖1 算法示例

2 實(shí)驗(yàn)結(jié)果與討論

2.1 實(shí)驗(yàn)數(shù)據(jù)

在試驗(yàn)管段設(shè)置人工缺陷，采用管道內(nèi)漏磁檢測(cè)器牽拉試驗(yàn)和水驅(qū)試驗(yàn)得到實(shí)驗(yàn)數(shù)據(jù)，牽拉試驗(yàn)如圖2所示。

圖2 牽拉試驗(yàn)

試驗(yàn)管段包括Φ273 mm×12.7 mm、Φ273 mm×7.1 mm、Φ457 mm×14.3 mm、Φ457 mm×7.1 mm,共4種規(guī)格。所有數(shù)據(jù)已完成數(shù)據(jù)的預(yù)處理、管段劃分，采用的數(shù)據(jù)為軸向分量數(shù)據(jù)。

2.2 實(shí)驗(yàn)結(jié)果

2.2.1Φ273 mm牽拉管道實(shí)驗(yàn)結(jié)果

利用管道漏磁檢測(cè)器在Φ273 mm管道上進(jìn)行試驗(yàn)，提取數(shù)據(jù)并進(jìn)行預(yù)處理和管段劃分。在12.7、7.1 mm壁厚的管道識(shí)別結(jié)果如圖3所示?？梢钥闯鲈摲椒ㄔ趦煞N不同壁厚的管道所有缺陷均完全識(shí)別出。

圖3 Φ273 mm牽拉管道缺陷識(shí)別結(jié)果

2.2.2Φ457 mm牽拉管道實(shí)驗(yàn)結(jié)果

為了進(jìn)一步驗(yàn)證該方法在不同口徑檢測(cè)器的效果，在Φ457 mm管道上進(jìn)行試驗(yàn)。14.3、7.1 mm壁厚的管道識(shí)別結(jié)果如圖4所示。從圖4可以看出，所有缺陷均完全識(shí)別出。

圖4 Φ457 mm牽拉管道缺陷識(shí)別結(jié)果

2.2.3Φ273 mm水驅(qū)管道實(shí)驗(yàn)結(jié)果

通常在牽拉管道中管道運(yùn)行較為平穩(wěn)，數(shù)據(jù)質(zhì)量較高。為更真實(shí)的模擬實(shí)際的檢測(cè)工況，選取Φ273 mm水驅(qū)管道，壁厚12.7 mm，并設(shè)置了208個(gè)人工缺陷進(jìn)行測(cè)試，正確識(shí)別缺陷合計(jì)205個(gè)，檢出率98.6%，如圖5所示。

圖5 Φ273 mm水驅(qū)管道缺陷識(shí)別結(jié)果

2.3 結(jié)果分析

在牽拉管道和水驅(qū)管道的實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)，該方法在Φ273 mm和Φ457 mm的牽拉管道缺陷檢出率達(dá)到了100%，在Φ273 mm水驅(qū)管道檢出率為98.6%，達(dá)到了工程檢測(cè)要求。

水驅(qū)管道人工缺陷識(shí)別率低于牽拉管道主要原因是：牽拉管道試驗(yàn)過(guò)程中檢測(cè)器運(yùn)行速度穩(wěn)定，波動(dòng)小，管道壁厚較為均勻，檢測(cè)過(guò)程中引起的噪聲和干擾較少，數(shù)據(jù)質(zhì)量較高。但在水驅(qū)管道檢測(cè)過(guò)程中受介質(zhì)、管道彎頭等因素影響，運(yùn)行速度不穩(wěn)定，同時(shí)由于采用水作為驅(qū)動(dòng)介質(zhì)，在管道內(nèi)壁會(huì)出現(xiàn)不同程度的輕微腐蝕，整體的數(shù)據(jù)質(zhì)量相對(duì)于牽拉數(shù)據(jù)較低。在進(jìn)行聚類的過(guò)程中，易產(chǎn)生干擾，出現(xiàn)識(shí)別錯(cuò)誤。

整體來(lái)說(shuō)，本文方法在牽拉管道和水驅(qū)管道識(shí)別結(jié)果均符合管道內(nèi)檢測(cè)檢出率的要求，達(dá)到了漏磁檢測(cè)缺陷量化與重構(gòu)的要求。

3 結(jié)論

針對(duì)漏磁數(shù)據(jù)缺陷區(qū)域標(biāo)記的問(wèn)題，提出了一種基于K均值聚類的油氣管道缺陷區(qū)域標(biāo)記方法，并通過(guò)管道內(nèi)檢測(cè)數(shù)據(jù)進(jìn)行了驗(yàn)證及分析。結(jié)果表明，該方法可有效識(shí)別出漏磁數(shù)據(jù)中的缺陷區(qū)域，識(shí)別準(zhǔn)確度滿足工程要求。該方法由于不需要根據(jù)檢測(cè)器和管道情況單獨(dú)設(shè)置閾值，通過(guò)在不同管徑、不同壁厚的檢測(cè)數(shù)據(jù)試驗(yàn)表明，該方法具有較廣泛的適應(yīng)性。