李翔宇，高憲文，李琨，侯延彬（東北大學信息科學與工程學院，遼寧 沈陽 089；渤海大學工學院，遼寧 錦州 03）

?

基于多源信息特征融合的抽油井動液面集成軟測量建模

李翔宇1，高憲文1，李琨2，侯延彬1
（1東北大學信息科學與工程學院，遼寧 沈陽 110819；2渤海大學工學院，遼寧 錦州 121013）

摘要：針對傳統(tǒng)抽油井動液面（DLL）檢測只能依靠人工操作回聲儀測試，無法實時在線檢測的問題，提出基于多源信息特征融合的抽油井動液面集成軟測量新方法。采用快速傅里葉變換（FFT）將抽油機懸點載荷及振動時域信號轉(zhuǎn)換成頻域信號；采用核主元分析（KPCA）提取懸點載荷及振動頻譜和電功率、井口油、套壓時域信號非線性特征；利用改進的模糊交互式自組織數(shù)據(jù)分析聚類（ISODATA）和高斯過程回歸（GPR）融合時頻信息特征,建立多個動態(tài)子模型；利用權(quán)重優(yōu)化證據(jù)理論（D-S）構(gòu)造的概率分配函數(shù)作為權(quán)值因子，對子模型輸出進行集成以得到最終的DLL預測值。油田現(xiàn)場應用證明了該方法的有效性。

關(guān)鍵詞：信息融合；動液面；高斯過程回歸；預測；石油；動態(tài)建模

2015-11-06收到初稿，2016-03-14收到修改稿。

聯(lián)系人：高憲文。第一作者：李翔宇（1982—），男，博士研究生。

Received date: 2015-11-06.

Foundation item: supported by the National Natural Science Foundation of China (61573088, 61403040, 61433004).

引　言

游梁式抽油機井有桿泵采油（簡稱有桿抽油系統(tǒng)或有桿泵抽油井）是國內(nèi)外油田普遍采用的機械采油方式，是油田開采工藝流程的主要能耗、物耗單元，如何實現(xiàn)其優(yōu)化運行，提高機采系統(tǒng)效率，一直是該領(lǐng)域研究重點[1-3]。實時在線檢測井下動液面位置（dynamic liquid level，DLL）對實現(xiàn)采油過程優(yōu)化運行和油田節(jié)能降耗減排意義重大[4-5]。

在實際生產(chǎn)中，一般是在停機狀態(tài)下人工操作回聲儀探測井下動液面，該測試方法工作強度大、成本高、效率低，無法實時在線檢測，并存在一定安全隱患[6-7]。文獻[6-7]以地面示功圖、油井產(chǎn)量、井口油、套管壓力等易測關(guān)鍵生產(chǎn)參數(shù)為輔助變量，分別采用最小二乘支持向量機（LSSVM）[8]和高斯過程回歸（GPR）[9-11]建立了井下動液面單一軟測量模型。但采油過程作為典型的連續(xù)生產(chǎn)工業(yè)過程，過程對象具有機理復雜、多變量、強耦合、非線性、時變性、大時滯和隨機干擾大等特性[12-13]，且采集到的數(shù)據(jù)存在多工況，其擾動幅度和對象特性均不同，導致單一軟測量模型存在模型學習時間過長、過程特性匹配不佳、信息融合不充分、泛化性差、精度低等問題[13-15]。研究表明，多模型集成軟測量建模方法能夠有效提高模型精度、泛化能力、有效性及可信度[16-19]。在多模型軟測量建模中，通常采用k-均值、模糊C均值等聚類方法對樣本集進行聚類，但上述聚類方法需要事先確定聚類中心和數(shù)目，算法對初始值的選取較敏感，計算過程易陷入局部最優(yōu)[19-21]。針對回歸問題，常用的多模型集成方法有簡單平均集成、基于切換策略的集成、基于多元線性回歸的集成、基于加權(quán)策略的集成等方法，樣本劃分不精確及聚類結(jié)果不佳等因素都會導致多模型集成輸出精度降低、融合能力和預測效果不理想[17-18]。此外，文獻[22]指出，動液面等井下工況參數(shù)與抽油機懸點載荷、振動信號頻譜強相關(guān)，且頻域特征明顯，但文獻[6-7]并未提取載荷、振動的頻域信息。然而，頻譜數(shù)據(jù)卻存在超高維和共線性問題，不利于構(gòu)建有效的軟測量模型[18-19]。

工業(yè)過程一般都是動態(tài)的，但大多數(shù)軟測量建模方法是基于過程的靜態(tài)假設，即認為生產(chǎn)過程運行比較平穩(wěn)，建模使用的樣本是與時間無關(guān)的離散數(shù)據(jù)，但當系統(tǒng)的動態(tài)變化比較明顯時，靜態(tài)軟測量模型無法反映工業(yè)過程的動態(tài)信息和全局特性，造成模型的適應性差，預測精度大大下降，甚至造成模型失效[17,23]。文獻[7]采用具有外部輸入自回歸滑動平均模型（ARMA）結(jié)構(gòu)的GPR軟測量模型來提高對井下動液面動態(tài)響應特性的預測能力，但單一模型結(jié)構(gòu)仍然導致動液面預測效果不佳。文獻[17]采用差分自回歸滑動平均（ARIMA）模型對靜態(tài)條件下LSSVM多模型輸出進行動態(tài)校正，以提高模型的動態(tài)響應特性，但由于該軟測量模型的訓練樣本和輸入變量中不包含過去時刻的信息，因而本質(zhì)上仍屬靜態(tài)模型。

針對上述問題，本文提出基于多源信息特征融合和多模型集成的抽油井動液面軟測量建模新方法。首先利用核主元分析（KPCA）方法[23-28]非線性提取懸點載荷及振動頻譜、電動機輸入功率、井口油、套管壓力時域信號特征，消除數(shù)據(jù)相關(guān)性，降低樣本的維數(shù)，從而降低模型的復雜性，提高其泛化能力。在不事先確定聚類數(shù)目和中心的前提下，采用改進的模糊ISODATA算法[29]將特征提取后的樣本集快速、有效地自適應融合拆分為多個與對象特征更加吻合的樣本子集。利用GPR處理高維數(shù)、小樣本、非線性等復雜問題具有很好的適應性，且泛化能力強，與神經(jīng)網(wǎng)絡、支持向量機相比，具有容易實現(xiàn)、超參數(shù)自適應獲取、非參數(shù)推斷靈活以及輸出具有概率意義的優(yōu)點[10]，進行多模型動態(tài)軟測量建模。最后結(jié)合權(quán)重優(yōu)化D-S證據(jù)融合方法[30]，針對證據(jù)間的沖突問題，采用基于加權(quán)平均修正法的證據(jù)體修改方法，將修正后的概率分配函數(shù)作為權(quán)值因子對子模型的輸出進行有效集成，以提高模型的預測精度和泛化能力。

1　有桿泵采油過程機理分析

有桿泵采油是一個典型的具有非線性、機電液強耦合等綜合復雜特性的物理過程，其原理如圖1所示。油田常用的CYJ10-3-53HB型抽油機沖程3 m，沖次5 min?1時，懸點運動規(guī)律如圖2所示。上沖程時，掛在驢頭上的懸繩器通過井下數(shù)百至數(shù)千米長、重達數(shù)噸的抽油桿柱帶動深井泵柱塞上行，柱塞上的游動閥受油管內(nèi)液柱壓力而關(guān)閉，泵內(nèi)壓力降低，固定閥在油、套管環(huán)型空間內(nèi)氣、液柱壓力（通常稱為沉沒壓力或泵口壓力，主要取決于動液面位置）的作用下被打開，泵內(nèi)吸入液體。此時，考慮沉沒壓力影響后的液柱載荷逐漸作用于抽油桿柱，完成加載過程。下沖程時，抽油桿柱和柱塞在其自重和懸繩器的約束下下行，泵內(nèi)壓力升高，固定閥關(guān)閉，游動閥被頂開，泵向油管內(nèi)排液。此時，液柱載荷逐漸作用于油管柱，完成卸載過程。此外，考慮到驢頭變速運動和周期性作用于抽油桿柱的液柱載荷的激勵，抽油桿柱產(chǎn)生的彈性振動載荷、抽油桿柱和液柱的慣性載荷以及井下摩擦載荷亦作用于懸點，因此，抽油機懸點載荷存在周期性大幅度波動。油田實測載荷、扭矩波形如圖3所示。

圖1　有桿泵采油工藝及動液面軟測量策略Fig.1　Rod pumping process and soft sensor strategy for DLL

圖2　抽油機懸點運動規(guī)律Fig.2　Polished rod motion law of beam pumping units

圖3　有桿抽油系統(tǒng)負荷/扭矩曲線Fig.3　Load and torque curves of rod pumping systems

抽油井動液面是指油井正常生產(chǎn)時的油套環(huán)空液面，是采油過程中油氣從地層向井筒滲流與深井泵向地面排液相互作用的宏觀表征。動液面深度越大，油井生產(chǎn)壓差越大，油井產(chǎn)量越高，井內(nèi)液體的有效舉升高度越大，抽油機做功越多，電動機平均輸入功率也越大，反之亦然。在微觀層面上，動液面位置直接影響流體進泵運動規(guī)律，從而影響懸點載荷、泵的充滿程度和泵效，并在示功圖形狀特征上有所反映。上沖程中，在沉沒壓力的作用下，井內(nèi)液體克服固定閥阻力進入泵內(nèi)，作用在柱塞底部而產(chǎn)生向上的載荷，其大小與動液面位置強相關(guān)，是上沖程懸點靜載荷的重要組成部分，直接影響懸點振動特性和動載荷。當動液面深度過大時，泵的沉沒度過小，沉沒壓力過小，加之由于井底流動壓力降低造成原油中所溶解的天然氣大量析出，井液中游離氣含量增大，導致泵的充滿程度不足，泵效低，懸點示功圖通常呈如圖4所示“刀把”形狀[7]，其中，圖4(a)、(b)分別為深井泵處于輕度和嚴重充滿程度不足工況時的實測示功圖。此外，動液面位置影響油套環(huán)空的物質(zhì)平衡，直接反映于井口套壓，并最終影響井口油壓。

圖4　有桿抽油系統(tǒng)懸點示功圖Fig.4　Surface dynamometer cards of rod pumping systems

2　動液面軟測量建模策略

結(jié)合油田普遍采用的有桿泵抽油井，提出由時頻轉(zhuǎn)換、非線性特征提取、樣本聚類、軟測量建模、多模型集成這5個部分組成的井下動液面軟測量策略如圖1所示。時頻轉(zhuǎn)換模塊求得懸點載荷及振動信號的頻譜，非線性特征提取模塊提取懸點載荷/振動頻域信號和電功率、井口油/套管壓力時域信號的非線性特征，樣本聚類模塊將樣本集融合分成多個樣本子簇，軟測量模型建立特征變量與動液面間的函數(shù)關(guān)系，多模型集成模塊則融合子模型預測輸出，獲得最終動液面估計值。

動液面參數(shù)、抽油機懸點載荷、振動、電動機輸入功率以及井口油、套管壓力信號間的關(guān)系為

式中，Ld表示動液面深度；AP、fP、Av、fv分別表示懸點載荷、振動頻域信號的幅值和頻率；Pm表示電動機瞬時輸入功率時域信號；pc表示井口套壓時域信號；pt表示井口油壓時域信號；?V表示為建模動態(tài)項。

3　KPCA非線性特征提取

KPCA算法對非線性系統(tǒng)具有更好的特征提取能力[24]，通過求解核矩陣的特征值與特征向量，避免了特征空間中求解特征向量的復雜問題[31]。采用非線性映射φ(?)，將數(shù)量為N的輸入樣本集X={ xi,i = 1,2, L ,N }映射到一個高維特征空間，φ( xi)協(xié)方差矩陣C由式（2）計算[24-28,31]

式中，xi為X的第i行向量。通過解特征值實現(xiàn)主成分計算，找到滿足式（3）的特征值與特征向量。

v可由映射到特征空間的所有樣本線性表示為

式中，iα為方程系數(shù)。

定義一個N×N的矩陣Kij

求解式（4）特征向量v的問題可轉(zhuǎn)變?yōu)榍缶仃嘖ij的特征向量α的問題。

式中，α為iα所組成的列向量。

歸一化特征向量v，此時樣本φ( X )在v上的投影為

式中，h(X)為非線性主成分分量；特征值λi小成分hi，通常情況下可視為是噪聲引起的；比值表示分量hi在整體方差中的貢獻程度。主成分按式（8）選取

式中，Z為選取的主成分數(shù)量，E為選取閾值。

以上推導假設∑φ(xi)=0，而實際情況不一定成立，因此，將式（5）中Kij中 心化處理為

式中，L為N×N單位矩陣。

4　改進的模糊ISODATA聚類

模糊ISODATA（interactive self-organizing data analysis technique algorithm）是一種無監(jiān)督聚類算法，能自動根據(jù)樣本的特征進行聚類中心的修正，但仍存在高維數(shù)據(jù)距離計算不準確、容易陷入局部最優(yōu)解以及需要設定聚類數(shù)目等缺點，針對上述問題，采用基于Hsim函數(shù)的距離度量、模擬退火算法（SA）的全局尋優(yōu)，引入ISODATA算法中的“合并”和“分裂”機制實現(xiàn)聚類過程中聚類數(shù)目的自動修正，改進后算法具體步驟如下[29-30]。

（1）設定初始參數(shù)，包括期望得到的聚類數(shù)c、每一類中允許的最小類別數(shù)目Mm、類內(nèi)樣本分布標準差SD、兩分類間最小距離Dm、每次迭代允許合并的最大聚類對數(shù)L、允許的最大迭代次數(shù)Ml、初始溫度和每個溫度下迭代的次數(shù)R。

隨機產(chǎn)生初始分類矩陣（2≤c≤N），有

采用Chebyshev距離計算初始聚類中心

給定一個足夠大的初始溫度T0，初始最優(yōu)解表示為：計算初始解對應的目標函數(shù)值J(S(0))，在每個溫度值T需要的迭代次數(shù)為R。令s =0。

（2）進行聚類的“分裂”與“合并”，它們的執(zhí)行次數(shù)不超過允許的最大迭代次數(shù)。

其中，1≤h≤c，1≤i≤N。

隨機更新兩分類間最小距離（Dm）值，搜索范圍為(0,1]。此時新解為：計算新解所對應的目標函數(shù)值f(S(s))。

本文采用Xie-Beni聚類有效性指標，定義如下

式中，vi是第i類的中心，uik是數(shù)據(jù)xk屬于類i的模糊隸屬度，XB(U, c)越小說明數(shù)據(jù)的聚類效果越好。

（4）計算差值?XB= XB(S(s+1))?XB(S(s))。若?XB≤0，則接受新解為當前解；若?XB＞0，則以概率P接受新解。

當e= random[0, 1]＜P，則接受新解S(s)作為當前解，否則采用前解。

（5）對當前溫度T降溫，將步驟（2）～步驟（4）迭代L次。

（6）當滿足循環(huán)終止條件時，算法輸出當前解為最優(yōu)解，結(jié)束算法；當循環(huán)不滿足終止條件時，按照退火方式降低溫度，轉(zhuǎn)到步驟（2），令s= s +1，繼續(xù)迭代。本文采用的退火方式為

式中，t為步長，α為常數(shù)，取值0.5～0.99，T0為初始溫度。

5　GPR動態(tài)建模

式中，y=[y1, … , yn]T為訓練樣本輸出構(gòu)成的向量；k(x?)=[C(xi,x?)]N×1為測試輸入和訓練樣本輸入間的協(xié)方差向量；為訓練樣本輸入間的N×N協(xié)方差矩陣；C(x?, x?)為測試輸入與其自身的協(xié)方差。式（17）和式（18）表明，GPR利用指定的協(xié)方差函數(shù)和訓練樣本信息實現(xiàn)對測試樣本的預測，給出包含均值和方差的預測統(tǒng)計學特性，這是GPR區(qū)別于ANN、SVM等方法非常突出的優(yōu)點[9]。

協(xié)方差函數(shù)的確定是建立高斯過程模型的關(guān)鍵步驟，它包含了對期望函數(shù)特性的假設。對任意一組輸入，協(xié)方差函數(shù)應滿足其產(chǎn)生的協(xié)方差矩陣為對稱半正定?？紤]到系統(tǒng)的平穩(wěn)性，下列徑向基函數(shù)是最常用的一類協(xié)方差函數(shù)

式中，υ0表示服從高斯分布的噪聲方差；υ1表示局部相關(guān)性的程度；ωd為模型的測度參數(shù)；δij是Kronecker算子。

當上述協(xié)方差函數(shù)類型確定后，通常使用極大似然、交叉驗證和馬爾可夫鏈蒙特卡羅3種方法對其超參數(shù)進行調(diào)整，本文通過極大化對數(shù)似然函數(shù)的方法求得

優(yōu)化過程中，要計算對數(shù)似然函數(shù)對各參數(shù)的導數(shù)

GPR除了可用于建立靜態(tài)非線性映射，基于具有外部輸入自回歸模型結(jié)構(gòu)的GPR還可用于動態(tài)系統(tǒng)建模[32]，系統(tǒng)輸出y(k)與過去n個輸出和過去n個輸入服從如下非線性關(guān)系

式中，f()?為非線性函數(shù)，本文利用GPR對其進行擬合；ε(k )為白噪聲。

6　權(quán)重優(yōu)化D-S多模型集成

6.1權(quán)重優(yōu)化D-S證據(jù)融合方法

將聚類所得c個子模型作為D-S證據(jù)理論中的辨識框架，則可將任一子模型視為焦元Aj( j= 1,2, L ,c )。對于N個樣本數(shù)據(jù)X={ xi}分別求出其對應于第j個子模型（也即第1個焦元A1）的模糊隸屬度，根據(jù)D-S證據(jù)合成規(guī)則對上述N個概率分配函數(shù)進行融合，將融合后的結(jié)果作為第j個子模型的基本概率分配函數(shù)[17,30]

D-S證據(jù)理論的主要缺陷是不能有效處理沖突證據(jù)之間的組合問題，為了避免這種情況出現(xiàn)，本文采用基于加權(quán)平均修正法的證據(jù)體修改方法。為了進一步降低沖突，本文采用一種證據(jù)間支持度的計算方法，并由其計算證據(jù)的可信度，但是還存在沖突度高的問題；采用證據(jù)體不確定度函數(shù)AM衡量，并利用負指數(shù)函數(shù)來抑制不確定度，使證據(jù)變得清晰。由證據(jù)的可信度和不確定度計算各條證據(jù)的權(quán)重，根據(jù)Dempster組合規(guī)則進行合成。為了避免產(chǎn)生可信度小而證據(jù)權(quán)重大的情況，采用PSO算法選擇合理的支持度維數(shù)和不確定度抑制因子，削弱證據(jù)合成的沖突度，保證證據(jù)合成的合理性。具體修正步驟如下[30]。

定義設有辨識框架Θ={A1, A2, L , Ac}，X為2Θ上的一個子集，有N個證據(jù)體mi( i= 1,2, L ,N )，m(·)為Θ上給定的基本概率分配函數(shù)，給出單個焦元Y∈Θ的pignistic概率分布如下

其中，|X|為子集X中焦元的個數(shù)。那么，兩個證據(jù)的支持度可以定義為

（1）計算證據(jù)體的可信度Credp(mi)為

（2）采用負指數(shù)函數(shù)形式計算證據(jù)體mi的不確定度Unceα(mi)

式中，AM(mi)表示證據(jù)體不確定度函數(shù)，其定義為

（3）由證據(jù)的可信度和不確定度求取每條證據(jù)所對應的權(quán)重為

（4）以各證據(jù)間沖突程度最小為目標，采用PSO算法優(yōu)化支持度維數(shù)p和不確定度抑制因子α。

（5）修正后的概率分配函數(shù)值可表示為

6.2多模型集成輸出

將式（30）得到的c個證據(jù)概率分配函數(shù)作為子模型的權(quán)值因子，對子模型進行集成，得到的多模型輸出為

7　工業(yè)應用

以如圖5所示的某油田J13-21井生產(chǎn)過程為例，按實驗設計方案逐漸調(diào)整抽油機沖次，記錄整個過程的關(guān)鍵生產(chǎn)參數(shù)信號，某一樣本的懸點載荷/位移/振動、井口油/套壓和電機電流/功率信號波形如圖6～圖8所示。具體的建模步驟如下。

圖5　智能油井網(wǎng)絡診斷控制系統(tǒng)Fig.5　Diagnosis and control system of internet-based intelligent oil wells

（1）數(shù)據(jù)處理。對現(xiàn)場采集到的數(shù)據(jù)，采用基于穩(wěn)健位置估計的方法[27]剔除離群點，將其中360組數(shù)據(jù)樣本隨機分成兩組，260組數(shù)據(jù)作為訓練樣本，100組用于檢驗模型的有效性。將懸點載荷、加速度時域信號通過FFT轉(zhuǎn)化為頻域信號，如圖9所示，標準化處理后進行主元分析。

圖6　懸點載荷/位移/振動信號波形Fig.6　Load, position and vibration signals of polished rod

圖7　井口油/套壓Fig.7　Tubing/casing pressure of wellhead

圖8　電機瞬時電流/功率波形Fig.8　Instantaneous input current/power of motor

圖9　懸點載荷/振動信號頻域波形Fig.9　Frequency domain curves of load and vibration signals of polished rod

（2）采用KPCA對時/頻信號進行非線性特征提取。本文選擇徑向基函數(shù)K=exp[? (x? xi)2/2 σ2]作為KPCA核函數(shù)，E =85%，分析結(jié)果如表1所示。

表1　主元分析結(jié)果Table 1　Results of principal component analysis

（3）利用改進的模糊ISODATA對訓練樣本進行聚類。設定初始參數(shù)，期望得到的聚類數(shù)目c =4，每一類中允許的最小類別數(shù)目Mm=5，類內(nèi)樣本分布標準差Ms= 0.1×10?3，兩分類間最小距離Md= 0.3，每次迭代允許合并的最大聚類對數(shù)L =1，允許的最大迭代次數(shù)Ml= 1000。最終聚類數(shù)目c =6。

（4）動態(tài)GPR子模型建立。利用GPR方法對6個樣本子集進行訓練學習，建立6個動態(tài)子模型。協(xié)方差函數(shù)超參數(shù)θ中w1～w9的初始值均為1，1υ 和0υ分別取1和0.001，即假設每個輔助變量對輸出預測的貢獻相同[9]，使用共軛梯度法，通過式（21）搜索得到與輸入對應的6個子模型超參數(shù)最優(yōu)值。根據(jù)對采油過程中動液面變化規(guī)律的觀察和油田采油工藝專家經(jīng)驗，提取油井前一天關(guān)鍵生產(chǎn)參數(shù)特征作為軟測量模型的輸入，即式（22）中取n =1。

（5）權(quán)重優(yōu)化D-S子模型集成輸出。按照式（24）～式（30）計算得到各個子模型所對應的權(quán)重，然后利用式（31）對各子模型的輸出進行證據(jù)融合，得到多模型集成輸出。

為驗證本文方法的有效性，與另外5種不同軟測量方法作比較，模型訓練和測試結(jié)果如圖10和圖11所示，誤差分析如表2和表3所示。GPR、KPCA-GPR方法采用單一靜態(tài)模型，GPR中輔助變量采用機理分析結(jié)合人工經(jīng)驗獲取，KPCA-GPR中采用多源信息融合KPCA特征提取獲取輔助變量，上述兩種方法由于模型過于單一導致預測精度不高。KPCA-FKCM-GPR和KPCA-ISODATA-DGPR方法采用多模型結(jié)構(gòu)，使用加權(quán)策略集成多模型輸出，后者使用動態(tài)建模方法，但聚類隸屬度的準確性對權(quán)值因子影響較大，使得上述兩種模型擬合效果仍不十分理想。ILDGPR是文獻[7]所建增量學習動態(tài)GPR動液面軟測量方法。圖12為現(xiàn)場應用過程中15組動液面明顯變化時軟測量值與實際值的比較，表4為預測誤差分析。

圖10　軟測量模型訓練結(jié)果比較Fig.10　Comparison of model training results

圖11　軟測量模型測試結(jié)果比較Fig.11　Comparison of model test results

圖12　軟測量模型預測結(jié)果比較Fig.12　Comparison of model prediction results

表2　訓練誤差分析Table 2　Training error analysis

表3　測試誤差分析Table 3　Test error analysis

表4　預測誤差分析Table 4　Prediction error analysis

由上述結(jié)果可知，在使用KPCA對采油過程關(guān)鍵生產(chǎn)參數(shù)時/頻特征進行非線性提取的基礎(chǔ)上，基于改進的模糊ISODATA和權(quán)重優(yōu)化D-S的GPR動態(tài)集成建模與傳統(tǒng)的單模型及多模型方法相比，預測精度和泛化能力有了較大改善。本文方法不但充分考慮了懸點載荷、振動、電機瞬時功率等過程參數(shù)與井下動液面間存在相關(guān)性、互補性和冗余性的現(xiàn)象，而且采用改進的模糊ISODATA動態(tài)聚類算法，在不預先確定聚類數(shù)目的前提下，實現(xiàn)對訓練樣本的自動聚類劃分，利用權(quán)重優(yōu)化D-S證據(jù)融合方法，加權(quán)平均修正相應焦元下各證據(jù)體的基本概率分配值，對多模型輸出進行有效集成，動態(tài)GPR建模改善系統(tǒng)的動態(tài)響應特性，因而，在對抽油井動液面的軟測量建模中取得了較好的擬合效果。

8　結(jié)　論

本文針對采油過程機電液耦合非線性機理復雜、井下工況的復雜性和不可預見性影響因素多等特點以及井下動液面實時在線檢測困難，提出了基于多源信息特征融合的抽油井動液面集成軟測量新方法。該方法解決了抽油機懸點測試信號特征難以提取、模型輸入變量維數(shù)過多造成訓練速度慢和過擬合、樣本聚類精度不高、多模型輸出合理集成以及靜態(tài)模型動態(tài)響應特性不佳等問題。該方法能夠比較準確地檢測抽油井動液面位置，具有廣泛的應用前景。進一步的研究方向是通過井下動液面、產(chǎn)液量等關(guān)鍵生產(chǎn)參數(shù)對有桿抽油系統(tǒng)綜合運行狀態(tài)進行推理判別。

References

[1]李兆敏, 林日億, 付路長, 等. 有桿抽油系統(tǒng)效率分析及抽汲參數(shù)優(yōu)化設計 [J]. 石油學報, 2005, 26 (5): 102-106. LI Z M, LIN R Y, FU L C, et al. Efficiency analysis of sucker-rod pumping system and optimal design of swabbing parameter [J]. Acta Petrolei Sinica, 2005, 26 (5): 102-106.

[2]姚春冬. 提高抽油機井系統(tǒng)效率的計算機仿真分析 [J]. 石油學報, 2005, 26 (4): 106-110. YAO C D. Computer simulation for enhancing system efficiency of rod pumping well [J]. Acta Petrolei Sinica, 2005, 26 (4): 106-110.

[3]董世民, 張喜順, 吳長杰, 等. 抽油機井區(qū)塊整體節(jié)能參數(shù)仿真優(yōu)化方法 [J]. 石油學報, 2010, 31 (3): 475-479. DONG S M, ZHANG X S, WU C J, et al. Simulation and optimization method for the integral energy-saving suction parameters of rod pumping wells in an oil field [J]. Acta Petrolei Sinica, 2010, 31 (3): 475-479.

[4]林日億, 孫茂盛, 張邵東, 等. 有桿抽油泵沉沒度的優(yōu)化設計方法[J]. 石油大學學報, 2005, 29 (4): 87-90. LIN R Y, SUN M S, ZHANG S D, et al. Optimization design method to determine submergence depth of sucker rod pump [J]. Journal of the University of Petroleum, 2005, 29 (4): 87-90.

[5]董世民, 邢明明, 張紅, 等. 基于系統(tǒng)效率最高的抽油機合理流壓研究 [J]. 燕山大學學報, 2011, 35 (4): 296-301. DONG S M, XING M M, ZHANG H, et al. Study of reasonable flowing pressure of rod pumping wells based on maximum system efficiency [J]. Journal of Yanshan University, 2011, 35 (4): 296-301.

[6]王通, 高憲文, 劉文芳. 自適應軟測量方法在動液面預測中的研究與應用 [J]. 化工學報, 2014, 65 (12): 4898-4904. WANG T, GAO X W, LIU W F. Adaptive soft sensor method and application in determination of dynamic fluid levels [J]. CIESC Journal, 2014, 65 (12): 4898-4904.

[7]李翔宇, 高憲文, 侯延彬. 基于在線動態(tài)高斯過程回歸抽油井動液面軟測量建模 [J]. 化工學報, 2015, 66 (6): 2150-2158. LI X Y, GAO X W, HOU Y B. Online dynamic Gaussian process regression for dynamic liquid level soft sensing of sucker-rod pumping well [J]. CIESC Journal, 2015, 66 (6): 2150-2158.

[8]SUYKENS J A K, VANDEWALLE J. Least squares support vector machine classifiers [J]. Neural Processing Letters, 1999, 9 (3): 293-300.

[9]王華忠. 高斯過程及其在軟測量建模中的應用 [J]. 化工學報, 2007, 58 (11): 2840-2845. WANG H Z. Gaussian process and its application to soft-sensor modeling [J]. Journal of Chemical Industry and Engineering (China), 2007, 58 (11): 2840-2845.

[10]何志昆, 劉光斌, 趙曦晶, 等. 高斯過程回歸方法綜述 [J]. 控制與決策, 2013, 28 (8): 1121-1129. HE Z K, LIU G B, ZHAO X J, et al. Overview of Gaussian process regression [J]. Control and Decision, 2013, 28 (8): 1121-1129.

[11]于濤, 王建林, 何坤, 等. 基于MPCA-GP的發(fā)酵過程分階段軟測量建模方法 [J]. 儀器儀表學報, 2013, 34 (12): 2703-2708. YU T, WANG J L, HE K, et al. Staged soft-sensor modeling method for fermentation process based on MPCA-GP [J]. Chinese Journal of Scientific Instrument, 2013, 34 (12): 2703-2708.

[12]湯健, 鄭秀萍, 趙立杰, 等. 基于頻域特征提取與信息融合的磨機負荷軟測量 [J]. 儀器儀表學報, 2010, 31 (10): 2161-2167. TANG J, ZHENG X P, ZHAO L J, et al. Soft sensing of mill load based on frequency domain feature extraction and information fusion [J]. Chinese Journal of Scientific Instrument, 2010, 31 (10): 2161-2167.

[13]雷瑜, 楊慧中. 基于高斯過程和貝葉斯決策的組合模型軟測量 [J].化工學報, 2013, 64 (12): 4434-4438. LEI Y, YANG H Z. Combination mode soft sensor based on Gaussian process and Bayesian committee machine [J]. CIESC Journal, 2013, 64 (12): 4434-4438.

[14]湯健, 柴天佑, 叢秋梅, 等. 基于EMD和選擇性集成學習算法的磨機負荷參數(shù)軟測量 [J]. 自動化學報, 2014, 40 (9): 1853-1866. TANG J, CHAI T Y, CONG Q M, et al. Soft sensor approach for modeling mill load parameters based on EMD and selective ensemble learning algorithm [J]. Acta Automatica Sinica, 2014, 40 (9): 1853-1866.

[15]李修亮，蘇宏業(yè)，褚健. 基于在線聚類的多模型軟測量建模方法[J]. 化工學報, 2007, 58 (11): 2834-2839. LI X L, SU H Y, CHU J. Multiple models soft-sensing technique based on online clustering arithmetic [J]. Journal of Chemical Industry and Engineering (China), 2007, 58 (11): 2834-2839.

[16]吳敏, 丁雷, 曹衛(wèi)華, 等. 鉛鋅燒結(jié)過程燒穿點的集成預測模型[J]. 控制理論與應用, 2009, 26 (7): 739-744. WU M, DING L, CAO W H, et al. An integrated prediction model for burn-through-point in lead-zinc sintering process [J]. Control Theory & Applications, 2009, 26 (7): 739-744.

[17]王振雷, 唐苦, 王昕. 一種基于D-S和ARIMA的多模型軟測量方法 [J]. 控制與決策, 2014, 29 (7): 1160-1166. WANG Z L, TANG K, WANG X. A multi-model soft sensing method based on D-S and ARIMA model [J]. Control and Decision, 2014, 29 (7): 1160-1166.

[18]湯健, 柴天佑, 趙立杰, 等. 基于振動頻譜的磨礦過程球磨機負荷參數(shù)集成建模方法 [J]. 控制理論與應用, 2012, 29 (2): 183-191. TANG J, CHAI T Y, ZHAO L J, et al. Ensemble modeling for parameters of ball-mill load in grinding process based on frequency spectrum of shell vibration [J]. Control Theory & Applications, 2012, 29 (2): 183-191.

[19]湯健, 柴天佑, 余文, 等. 在線KPLS建模方法及在磨機負荷參數(shù)集成建模中的應用 [J]. 自動化學報, 2013, 39 (5): 471-486. TANG J, CHAI T Y, YU W, et al. On-line KPLS algorithm with application to ensemble modeling parameters of mill load [J]. Acta Automatica Sinica, 2013, 39 (5): 471-486.

[20]李麗娟, 宋坤, 趙英凱. 基于仿射傳播聚類的ARA發(fā)酵過程建模[J]. 化工學報, 2011, 62 (8): 2116-2121. LI L J, SONG K, ZHAO Y K. Modeling of ARA fermentation based on affinity propagation clustering [J]. CIESC Journal, 2011, 62 (8): 2116-2121.

[21]張文清, 傅雨佳, 楊慧中. 基于改進聚類和加權(quán)bagging的多模型軟測量建模 [J]. 化工學報, 2012, 63 (9): 2697-2702. ZHANG W Q, FU Y J, YANG H Z. Multi-model soft-sensor modeling based on improved clustering and weighted bagging [J]. CIESC Journal, 2012, 63 (9): 21697-2702.

[22]何巖峰, 吳曉東, 韓國慶, 等. 示功圖頻譜分析新方法 [J]. 石油學報, 2008, 29 (4): 619-624. HE Y F, WU X D, HAN G Q, et al. Frequency spectrum analysis method for recognition of dynamometer card [J]. Acta Petrolei Sinica, 2008, 29 (4): 619-624.

[23]曹鵬飛, 羅雄麟. 化工過程軟測量建模方法研究進展 [J]. 化工學報, 2013, 64 (3): 788-800. CAO P F, LUO X L. Modeling of soft sensor for chemical process [J]. CIESC Journal, 2013, 64 (3): 788-800.

[24]湯健, 趙立杰, 岳恒, 等. 基于多源數(shù)據(jù)特征融合的球磨機負荷軟測量 [J]. 浙江大學學報, 2010, 44 (7): 1406-1413. TANG J, ZHAO L J, YUE H, et al. Soft sensor for ball mill load based on multi-source data feature fusion [J]. Journal of Zhejiang University, 2010, 44 (7): 1406-1413.

[25]MOGHADDAM B. Principal manifolds and probabilistic subspaces for visual recognition [J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2002, 24 (6): 780-788.

[26]TWINING C J, TAYLOR C J. The use of kernel principal component analysis to model data distributions [J]. Pattern Recognition, 2003, 36 (1): 217-227.

[27]李海波, 柴天佑, 岳恒. 浮選工藝指標KPCA-ELM軟測量模型及應用 [J]. 化工學報, 2012, 63 (9): 2892-2898. LI H B, CHAI T Y, YUE H. Soft sensor of technical indices based on KPCA-ELM and application for flotation process [J]. CIESC Journal, 2012, 63 (9): 2892-2898.

[28]王強, 田學民. 基于KPCA-LSSVM的軟測量建模方法 [J]. 化工學報, 2011, 62 (10): 2813-2817. WANG Q, TIAN X M. Soft sensing based on KPCA and LSSVM [J]. CIESC Journal, 2011, 62 (10): 2813-2817.

[29]LI K, GAO X W, ZHOU H B, et al. Fault diagnosis for down-hole conditions of beam pumping unit based on improved fuzzy ISODATA dynamic clustering algorithm [C]//Proc. 2013 10th International Conference on Fuzzy Systems and Knowledge Discovery. Shenyang, China: IEEE Computer Society, 2013: 279-284.

[30]李琨. 基于示功圖的游梁式抽油機井下故障診斷方法研究[D]. 沈陽: 東北大學, 2013. LI K. Study on diagnostic method for down-hole faults of the beam pumping unit based on dynamometer card [D]. Shenyang: Northeastern University, 2013.

[31]田中大, 高憲文, 李琨. 基于KPCA與LSSVM的網(wǎng)絡控制系統(tǒng)時延預測方法 [J]. 系統(tǒng)工程與電子技術(shù), 2013, 35 (6): 1281-1285. TIAN Z D, GAO X W, LI K. Networked control system time-delay prediction method based on KPCA and LSSVM [J]. Systems Engineering and Electronics, 2013, 35 (6): 1281-1285.

[32]PETELIN D, GRANCHAROVA A, KOCIJAN J. Evolving Gaussian process models for prediction of ozone concentration in the air [J]. Simulation Modelling Practice and Theory, 2013, 24 (6): 68-80.

Ensemble soft sensor modeling for dynamic liquid level of oil well based on multi-source information feature fusion

LI Xiangyu1, GAO Xianwen1, LI Kun2, HOU Yanbin1
(1College of Information Science & Engineering, Northeastern University, Shenyang 110819, Liaoning, China;2College of Engineering, Bohai University, Jinzhou 121013, Liaoning, China)

Abstract:The dynamic liquid level (DLL) of an oil well is traditionally measured onsite by using the acoustic method. This method, however, has its limitation in determining real-time dynamic liquid level. A new ensemble soft-sensor approach of DLL based on the multi-source information feature fusion was proposed. The polish rod load and vibration signal in the time domain was transformed into the frequency domain using fast Fourier transform (FFT). The kernel principal component analysis (KPCA) was used to extract the nonlinear feature of the load and vibration spectral signal and the power, casing head pressure andtubing head pressure time signal. The improved fuzzy interactive self-organizing data analysis technique algorithm (ISODATA) and Gaussian process regression (GPR) were used to fuse time/frequency information feature and establish multiple sub-models. Then, the final DLL prediction model was obtained through the ensemble of the sub-models based on the weight factor calculated by optimized-weighted Dempster-Shafer (D-S) theory. The oil field application showed the validity of the proposed method.

Key words:information fusion; dynamic liquid level; Gaussian process regression; prediction; petroleum; dynamic modeling

中圖分類號：TP 273

文獻標志碼：A

文章編號：0438—1157（2016）06—2469—11

DOI：10.11949/j.issn.0438-1157.20151673

基金項目：國家自然科學基金項目（61573088，61403040，61433004）。

Corresponding author:Prof. GAO Xianwen, gaoxianwen@ise.neu.edu.cn