， ，

(中國石油大學(xué)(北京)地球物理與信息工程學(xué)院，北京 102249)

0 引言

石油在能源結(jié)構(gòu)中占據(jù)至關(guān)重要的地位，采油是石油生產(chǎn)過程中最基礎(chǔ)和重要的環(huán)節(jié)。目前常用的采油方法中，有桿泵采油法是應(yīng)用最為廣泛的[1]。因此，有效地監(jiān)測抽油機(jī)井的工作狀態(tài)，快速準(zhǔn)確地識別油井故障并采取合理有效的措施，對提高采油效率，增進(jìn)油田效益有著舉足輕重的作用。

油井實時示功圖是油井工況的一個重要表征，它能夠直觀地反映出抽油機(jī)在采油過程中發(fā)生的各種異常狀況，同時聯(lián)系地質(zhì)情況及井下技術(shù)狀況等因素，可以識別出井下泵發(fā)生的故障類型。根據(jù)示功圖進(jìn)行工況分析的主要方法有“五指式動力儀”分析法、地面示功圖分析法和井下泵示功圖診斷法。相比于其它方法，地面示功圖分析法具有更直觀便捷，簡單易行的優(yōu)點。因此，根據(jù)油井示功圖判斷油井的工況是目前被廣泛應(yīng)用的方法[2]。

在使用示功圖進(jìn)行工況診斷的過程中，如何準(zhǔn)確地對示功圖進(jìn)行特征提取，直接關(guān)系到能否精確地對油井存在的故障做出判斷，因此，必須要選取恰當(dāng)?shù)奶卣魈崛》椒╗3]。目前，用于油井工況分析的示功圖特征提取的方法主要有基于幾何不變矩的矩特征向量分析法、基于網(wǎng)格法的灰度矩陣統(tǒng)計量法和基于DFT的傅里葉描述子法[4-6]。相比于其它方法，傅里葉描述子法具有計算方法簡潔高效、計算量小、冗余信息少的優(yōu)勢，可以比較準(zhǔn)確地提取出示功圖的圖像特征。

在傳統(tǒng)的傅里葉描述子法中，只提取了示功圖數(shù)據(jù)經(jīng)DFT得到的各次諧波的幅度譜，并沒有對各次諧波的相位譜進(jìn)行整合，相當(dāng)于只提取了圖像波動的幅度情況，并沒有提取圖像波動的位置信息，從而丟失了示功圖曲線形狀的凸起或凹陷這類代表位置的有用信息，對圖像識別的準(zhǔn)確性有較大的影響。本文提出了一種融合信號DFT的幅值與相位信息的特征提取方法，更加完整地提取了示功圖的特征信息，并構(gòu)造出能在極大程度上反映示功圖圖形特征的特征向量。此外，相對于現(xiàn)在被廣泛應(yīng)用的二分類支持向量機(jī)，本文采用一種基于多分類向量機(jī)的分類識別方法[7-8]，并優(yōu)選了核函數(shù)，進(jìn)一步提高了示功圖診斷的準(zhǔn)確率。

1 示功圖特征分析及特征向量的構(gòu)造

1.1 示功圖特征理論分析

示功圖是指在一個采油沖程過程中，以抽油機(jī)相對下死點的位移S為橫軸，以抽油機(jī)的懸點載荷P為縱軸，繪制出一條關(guān)于P-S的曲線，它是一個封閉曲線[1]。

理論示功圖是指只考慮懸點所承受的靜載荷及由其引起抽油機(jī)桿柱及有關(guān)柱的彈性變形，而不考慮其它因素影響時，所得到的示功圖曲線。典型示功圖是在理論示功圖的基礎(chǔ)上，只考慮某單一因素影響的載荷隨位移的變化關(guān)系曲線。

圖1為抽油機(jī)井常見工況的典型示功圖圖例[9]。由圖可見，正常示功圖近似于平行四邊形(理論示功圖)，上下邊出現(xiàn)輕微的震蕩波浪線；各類典型故障的示功圖圖形均存在一定的凸起或凹陷，且凸起或凹陷的位置不同。

圖1 典型示功圖圖例

在油田現(xiàn)場，由示功儀所采集的示功圖一般是由下死點(圖1(a)中A點)開始，經(jīng)過一段時間達(dá)到加載終止點(圖1(a)中B點)，該過程中得到的曲線稱為增載線(圖1(a)中AB段)；保持一定時間后又由上死點(圖1(a)中C點)開始，經(jīng)過一段時間達(dá)到卸載終止點(圖1(a)中D點)，該過程中得到的曲線稱為卸載線(圖1(a)中CD段)。

結(jié)合各類故障示功圖的形成因素，對圖1中各圖進(jìn)行分析可知，各類典型故障的示功圖圖形的特征主要表現(xiàn)為：

1)氣鎖：增載線和卸載線均有一定弧度。

2)氣體影響：增載線、卸載線均都比較平緩，但增載線相對于卸載線更陡一些。

3)供液不足：增載線和卸載線相互平行，且均比較陡直。

4)固定閥漏失：卸載線比增載線平緩，示功圖的左下角為圓弧。

5)游動閥漏失：增載線比卸載線平緩，示功圖的右上角為圓弧。

6)雙凡爾漏失：增載線和卸載線都比較平緩，形狀呈水平橢圓形。

7)游動閥關(guān)閉遲緩：加載線延長且形狀較陡直，示功圖的左上方有缺失。

8)柱塞脫出工作筒：增載線正常，卸載線不太明顯，圖形右下角有波浪曲線。

9)上碰泵：圖形右上角有一個環(huán)狀凸出。

10)下碰泵：圖形左下角有一個環(huán)狀凸出。

1.2 示功圖特征DFT分析

在油田現(xiàn)場，在一個沖程內(nèi)，由示功儀所采集的實測示功圖數(shù)據(jù)通常為200或250個點的位移與載荷數(shù)據(jù)，從而可以得到一組序列值(sn,pn),n=0,1,…,N-1，其中sn為位移，pn為載荷。將位移sn和載荷pn組合成如下一維復(fù)數(shù)序列：

zn=sn+jpn

n=(0,1,…,N-1)

(1)

對于封閉的曲線，該序列是以點數(shù)N為周期的周期序列。對其做DFT，得到：

k=0,1,…,N-1;n=0,1,…,N-1

(2)

Zk的模|Zk|稱為DFT的幅值，反映的是示功圖序列波動的各次諧波的大小，其相角h(k)稱為DFT的相位，包含著其波動的位置信息。

圖2為對各個典型示功圖曲線作250點DFT得到的幅值譜|Zk|的曲線。由于幅度譜反映的是示功圖波動的大小情況，因此幅度譜分析結(jié)果與前文所述的各類典型故障的示功圖圖形均存在一定的凸起或凹陷相對應(yīng)。

根據(jù)圖2可知，某些圖形特征相似的故障類型的示功圖，如“氣體影響”與“供液量不足”、“游動閥關(guān)閉遲緩”與“柱塞脫出工作筒”均存在較大的面積缺失現(xiàn)象；再如“上碰泵””與“下碰泵”均存在環(huán)狀凸出現(xiàn)象等，它們兩兩之間的DFT幅度譜差異并不大，因此，難以對這些故障種類的示功圖進(jìn)行準(zhǔn)確地區(qū)分。

圖2 典型故障示功圖的DFT幅值信息

圖3為分別對各個典型示功圖的上行數(shù)據(jù)(增載線)和下行數(shù)據(jù)(卸載線)作250點DFT得到的相位譜h(k)的曲線。由于相位譜反映的是示功圖波動的位置情況，因此相位譜分析結(jié)果與前文所述的各類典型故障的示功圖圖形的凸起或凹陷位置不同相對應(yīng)。

根據(jù)圖3可知，對于前文所述的DFT幅度譜差異不大的故障類型的示功圖，如“氣體影響”與“供液量不足”、“游動閥關(guān)閉遲緩”與“柱塞脫出工作筒”、“上碰泵”與“下碰泵”，它們兩兩之間的相位譜則有明顯的差異，可以較好地分析出示功圖特征的位置，從而較好地對這些故障種類的示功圖進(jìn)行準(zhǔn)確地區(qū)分。

圖3 典型故障示功圖的DFT相位信息

1.3 特征向量構(gòu)造

為了克服傳統(tǒng)傅里葉描述子方法在用于示功圖識別的缺陷，本文將示功圖序列的DFT幅度譜與相位譜進(jìn)行融合，構(gòu)造了能夠充分反映出示功圖圖形的形狀輪廓特征、細(xì)節(jié)波動與凹凸位置信息的特征向量。同時，根據(jù)DFT的定義可知，此特征向量的各個分量之間均線性無關(guān)并且相互正交，從而最大程度地獲取了示功圖在其特征空間的信息，因而有效地提高示功圖類型識別的準(zhǔn)確率。由歐拉變換可知，

Zk=|Zk|ejh(k)=|Zk|cos(h(k))+j|Zk|sin(h(k))

(3)

公式(3)中的實部和虛部均同時包含了DFT的幅度信息和相位信息，因此，為了簡化計算，僅選取實部來構(gòu)造特征向量。又由于信號的主要能量集中在其低次諧波分量上，因此，在兼顧精度及計算量的前提下，選取DFT的直流項及前9次諧波來構(gòu)造特征向量，則最終構(gòu)造的融合了幅度與相位信息的示功圖特征向量如(4)式所示：

(4)

這樣，對每一個示功圖都可以構(gòu)造出一個線性無關(guān)的特征向量，以此特征向量作為該示功圖的特征，輸入到由訓(xùn)練樣本訓(xùn)練得到的分類判別器中，進(jìn)行人工智能判別，從而對示功圖進(jìn)行正確的分類識別，得到油井的實際工況。

2 支持向量機(jī)算法

2.1 經(jīng)典支持向量機(jī)算法

支持向量機(jī)(SVM)是Corinna Cortes和Vapnik等人于20世紀(jì)90年代提出的一種機(jī)器學(xué)習(xí)方法，它主要應(yīng)用于模式識別領(lǐng)域，是一種結(jié)合了統(tǒng)計學(xué)理論的VC維理論以及結(jié)構(gòu)風(fēng)險最小化原則的模式識別方法，在有限的訓(xùn)練樣本中尋求學(xué)習(xí)能力與模型復(fù)雜度融合后的最佳折中結(jié)果，從而得到最小誤差分類器。

支持向量機(jī)的基本思想是：首先，判斷樣本數(shù)據(jù)是否線性可分：1)若樣本為線性可分，則直接尋找最優(yōu)分類面進(jìn)行分類；2)若樣本為線性不可分，則引入兩個變量(松弛系數(shù)和懲罰分量)進(jìn)行輔助分析，通過非線性映射將樣本映射到高維特征空間中，在高維空間中尋找最優(yōu)分類面對樣本進(jìn)行分類[10]；同時，它通過使用結(jié)構(gòu)風(fēng)險最小化原理，也在一定程度上對分類器實現(xiàn)全局的最優(yōu)效果進(jìn)行了保障。

SVM只用于兩類問題的分類，解決多分類問題時通常是利用若干個兩分類決策器解決多分類問題，具體實現(xiàn)方法大致可分為以下幾種：

1)一對一分類方法(one-versus-one)：該方法將多分類問題轉(zhuǎn)化為若干個子分類問題，任意兩類構(gòu)造一個兩分類分類器，對所有的兩分類分類器都進(jìn)行判斷，最終的結(jié)果為所占比重最大的類，如圖4所示。當(dāng)子分類的數(shù)量過多時，計算量很大，會大幅度的降低訓(xùn)練與測試的速度。分類結(jié)果可能會出現(xiàn)分類重疊現(xiàn)象。

圖4 一對一分類方式示意圖

2)一對多分類方法(one-versus-rest)：該方法將問題轉(zhuǎn)化為多個兩分類問題，將某類看作兩分類問題中的一類，將其余的所有類看作另一類進(jìn)行判別分類，依次循環(huán)，最終結(jié)果為分類函數(shù)值最大的類，如圖5所示。若訓(xùn)練樣本數(shù)目過大，會導(dǎo)致訓(xùn)練過程耗時很長。分類結(jié)果可能會出現(xiàn)分類重疊現(xiàn)象或不可分類現(xiàn)象。

圖5 一對多分類方法示意圖

3)二叉樹分類方法(Binary Tree)：該方法是將樣本數(shù)據(jù)分為兩大類，然后繼續(xù)劃分每個大類為兩個類，直至所有類不能再劃分為止，如圖6所示(圖中示例n=4)。在訓(xùn)練學(xué)習(xí)的過程中，存在錯誤向下累積現(xiàn)象。

圖6 二叉樹方法示意圖

4)DDAG算法(Decision Directed Acyclic Graph)：該算法以有向無環(huán)圖為理論基礎(chǔ)，將所有分類器構(gòu)成有向無環(huán)圖，每次參與分類的兩個類別可無重復(fù)的任意選取，如圖7所示(圖中示例n=4)。在訓(xùn)練的過程中，同樣存在錯誤向下累積現(xiàn)象，而且訓(xùn)練過程的時間通常較長。

圖7 有向無環(huán)方法示意圖

2.2 LIBSVM算法

針對經(jīng)典SVM的不足之處，許多研究人員對SVM進(jìn)行了更為深入地研究，他們在將公式變形后研究出應(yīng)用范圍更具體且具有某些特定優(yōu)勢的變形算法。其中，由臺灣大學(xué)林智仁博士等[11-12]于2001年開發(fā)設(shè)計的LIBSVM算法融合了收縮和緩存技術(shù)，

具有程序小，運用靈活；輸入?yún)?shù)少，操作簡單；是開源的，方便改進(jìn)、易于擴(kuò)展等特點，同時，具有較好性能又適用于Windows、Unix等多種操作系統(tǒng)，因此，目前被許多國際著名研究機(jī)構(gòu)所采用作為訓(xùn)練算法。

LIBSVM作為通用的SVM軟件包，該軟件對SVM所涉及的參數(shù)調(diào)節(jié)相對比較少，提供了很多默認(rèn)的參數(shù)設(shè)置，利用這些默認(rèn)參數(shù)可以有效地解決包括分類問題在內(nèi)的多種類型的問題。同時，它還有了交互檢驗的功能，利用所提供的多種常用核函數(shù)可以將交叉驗證的參數(shù)選擇變得更為精準(zhǔn)。該軟件可以解決C-SVM、ν-SVM、ε-SVR和ν-SVR等問題，包括基于一對一算法的多類模式識別問題。

所以，本文選用了LIBSVM支持向量機(jī)分類器，它能夠?qū)⒅С窒蛄繖C(jī)、分布估計與回歸完整結(jié)合在一起，從而對不同類型的示功圖進(jìn)行判別與診斷。LIBSVM分類器的設(shè)計與應(yīng)用步驟如圖8所示。

圖8 LIBSVM分類器的設(shè)計與應(yīng)用步驟

3 應(yīng)用實例

3.1 實驗過程

本論文研究的是一種對油井實測示功圖進(jìn)行分類識別的方法，該方法的實驗過程分為圖像預(yù)處理、訓(xùn)練和測試3個階段。整體實驗過程的流程如圖9所示。

圖9 實驗過程流程圖

圖像預(yù)處理階段主要是消除圖像中的無關(guān)信息，增強有用信息的可檢測性，同時最大限度地簡化數(shù)據(jù)，從而提高特征提取及匹配識別的準(zhǔn)確性。主要包括歸一化處理、形態(tài)學(xué)處理和坐標(biāo)數(shù)字化處理。

訓(xùn)練階段是根據(jù)已知故障類型的示功圖的面積特征，將故障示功圖粗略分為若干個集合，然后分別在各集合中對由訓(xùn)練樣本構(gòu)造的特征向量進(jìn)行訓(xùn)練，最終得到多分類判別器。

測試階段是根據(jù)待識別故障類型的示功圖的面積特征來確定其所屬集合，然后將由其構(gòu)造的特征向量輸入到相應(yīng)的多分類判別器中，從而實現(xiàn)示功圖的分類判別。

訓(xùn)練階段和測試階段在構(gòu)造特征向量時均采用利用前文介紹的方法，在得到融合DFT幅度譜與相位譜各特征成分的特征向量后，均采用匹配高斯徑向基核函數(shù)的支持向量機(jī)對示功圖進(jìn)行訓(xùn)練和測試。

3.2 實驗參數(shù)選擇

本實驗在Matlab仿真平臺上進(jìn)行仿真驗證。由于LIBSVM算法的準(zhǔn)確率主要與訓(xùn)練和測試時所選取的核函數(shù)種類有關(guān)，而核函數(shù)的正確選取依賴于產(chǎn)生分類問題的實際問題的特點，不同的實際問題對相似程度有著不同的度量，因此，選擇正確的核函數(shù)有助于提高分類準(zhǔn)確率。

核函數(shù)的種類由參數(shù)t決定。t=0時為線性核；t=1為多項式核；t=2時為高斯徑向基函數(shù)；t=3時為Sigmod核函數(shù)；t=4為自己定制的任意一種的核函數(shù)。表1為采用不同核函數(shù)時識別準(zhǔn)確率的對比。對表1中數(shù)據(jù)進(jìn)行對比分析可知：

1)對于相同的數(shù)據(jù)，不同核函數(shù)的識別準(zhǔn)確率不同；

2)當(dāng)樣本數(shù)據(jù)量較少時，高斯徑向基核函數(shù)識別準(zhǔn)確率更高。

表1 不同核函數(shù)時識別準(zhǔn)確率對比 %

因此，考慮到實際采集到的示功圖樣本數(shù)量較少的情況，本實驗在訓(xùn)練階段和測試階段支持向量機(jī)的核函數(shù)均選取高斯徑向基核函數(shù)。

3.3 實驗結(jié)果及分析

本實驗選取了氣鎖、氣體影響、供液不足、固定閥漏失、游動閥漏失、雙爾凡漏失游動閥關(guān)閉遲緩、柱塞脫出工作筒、上碰泵和下碰泵，共10 種油田實際生產(chǎn)中易出現(xiàn)的故障類型的共275組示功圖數(shù)據(jù)進(jìn)行分析研究。利用現(xiàn)場采集的實測示功圖建立樣本庫，并根據(jù)面積特征劃分為若干個集合，在各集合中根據(jù)每類典型故障示功圖樣本的數(shù)量情況，適當(dāng)選取其中一部分作為訓(xùn)練樣本，其余部分劃歸為測試樣本。在構(gòu)造特征向量時，分別采取僅使用幅值信息和融合幅值信息與相位信息兩種方法，對相同的樣本分別進(jìn)行訓(xùn)練和測試，測試結(jié)果如表2所示。

表2 不同特征向量構(gòu)造方法的識別準(zhǔn)確率對比

其中，方法一為僅利用DFT幅度譜構(gòu)造特征向量；方法二為融合DFT幅度譜與相位譜構(gòu)造特征向量。

結(jié)合前文對各類典型故障的示功圖圖形的理論分析和DFT分析結(jié)果，對表2中數(shù)據(jù)進(jìn)行對比分析可知：

1)本文提出的方法二的識別準(zhǔn)確率的平均值較方法一的識別準(zhǔn)確率的平均值有明顯提高；

2)方法一識別準(zhǔn)確率較高的均為示功圖圖形特征明顯且與其他故障類型的圖形相似度較低的故障類型，如氣鎖、雙凡爾漏失等，對于這些故障類型的示功圖，本文提出的方法二可將方法一的識別準(zhǔn)確率繼續(xù)保持；

3)方法一識別準(zhǔn)確率較低的故障類型多為圖形特征不明顯或與其他故障類型的圖形相似度較高的故障類型，如氣體影響、供液不足、固定閥漏失、游動閥關(guān)閉遲緩、柱塞脫出工作筒、下碰泵等，對于這些故障類型的示功圖，本文提出的方法二的識別準(zhǔn)確率相較于方法一的識別準(zhǔn)確率則有大幅提高。

4 結(jié)論

1)提出了一種融合DFT幅值與相位信息的示功圖特征提取方法，相較于傳統(tǒng)的傅里葉方法，本文方法不僅可以反映示功圖曲線的波動特征，而且可以反映示功圖曲線凹凸的位置特征，能夠更加全面和準(zhǔn)確地表達(dá)出示功圖所包含的特征信息。

2)相對于現(xiàn)在被廣泛應(yīng)用的二分類支持向量機(jī)，本文采用一種基于多分類向量機(jī)的分類識別方法，并選取高斯徑向基核函數(shù)，進(jìn)一步提高了示功圖診斷的準(zhǔn)確率。

3)實測表明，本文提出的示功圖識別方法對油井示功圖故障類型的診斷更加準(zhǔn)確有效，可為油井生產(chǎn)的實時分析與優(yōu)化控制提供技術(shù)支撐。