抽油機(jī)光桿位移的無(wú)線測(cè)試方法解析

付 濤

(中石化勝利油田分公司生產(chǎn)管理部，山東 東營(yíng) 257000)

0 引言

示功圖是判斷抽油機(jī)井工況的基本依據(jù)，也是采油廠日常生產(chǎn)管理的常規(guī)測(cè)試項(xiàng)目［1－2］。示功圖測(cè)量?jī)x是載荷與位移參數(shù)的綜合測(cè)量?jī)x器。對(duì)于位移測(cè)量，目前較常用的方法是拉繩式的直接位置測(cè)量法。

拉繩式直接位置測(cè)量方法以光柵、光電編碼器等作為傳感器［3］，通過(guò)計(jì)算拉繩的行程圈數(shù)來(lái)獲得油桿的位置。這種位置測(cè)量方式結(jié)構(gòu)簡(jiǎn)單且容易實(shí)現(xiàn)，但必須將傳感器的不同部件分別安裝在抽油機(jī)的活動(dòng)部位和固定部位，同時(shí)要求這些部件在很長(zhǎng)的運(yùn)動(dòng)軌道上嚙合良好，這在抽油機(jī)井污染嚴(yán)重、震動(dòng)和沖擊很大的環(huán)境中難以滿足。因此，許多功圖測(cè)量?jī)x改為利用加速度測(cè)量來(lái)間接進(jìn)行位置測(cè)量。

基于加速度的位移測(cè)量屬于間接測(cè)量方法，能夠?qū)崿F(xiàn)位移的無(wú)線測(cè)試，簡(jiǎn)化了測(cè)試程序，提高了功圖測(cè)量?jī)x的可靠性和位移測(cè)試數(shù)據(jù)的準(zhǔn)確性［4－6］。然而，由于存在加速度的零點(diǎn)校正和死點(diǎn)加速度振蕩等問(wèn)題，基于加速度的光桿位移測(cè)試方法也有其自身的不足之處。

1 系統(tǒng)組成與數(shù)據(jù)特點(diǎn)

無(wú)線位移測(cè)試系統(tǒng)的基本思路是沿用加速度值推算位移的方法。因此，為了研究抽油機(jī)光桿位移的無(wú)線測(cè)試方法，搭建了一個(gè)基于PC的光桿加速度無(wú)線數(shù)據(jù)采集試驗(yàn)系統(tǒng)。測(cè)試系統(tǒng)結(jié)構(gòu)如圖1所示。測(cè)試系統(tǒng)采用的加速度傳感器為ADXL202。

無(wú)線位移測(cè)試系統(tǒng)由無(wú)線傳感和數(shù)據(jù)接收處理兩部分組成。加速度傳感器、單片機(jī)系統(tǒng)和無(wú)線發(fā)送模塊整體封裝于殼體中組成無(wú)線傳感部分，通過(guò)卡口快速加裝的方式安裝于光桿上。無(wú)線接收模塊與PC機(jī)相連組成數(shù)據(jù)接收處理部分，在PC機(jī)上可通過(guò)軟件實(shí)現(xiàn)加速度波形的實(shí)時(shí)顯示與存儲(chǔ)，加速度的采樣頻率為50 Hz。

圖1中，ADXL202是AD公司推出的一種低成本、低功耗的單片雙軸加速度傳感器。它既可測(cè)量動(dòng)態(tài)加速度，又能測(cè)量靜態(tài)加速度［7］。

某抽油機(jī)光桿所對(duì)應(yīng)的1 min內(nèi)加速度曲線如圖2所示。

圖2 實(shí)測(cè)光桿加速度曲線Fig.2 The actual acceleration of polished rod

從圖2可以看出，受抽油桿自身震動(dòng)以及電磁干擾的影響，傳感器所測(cè)得的加速度波形包含大量的雜波分量。因此，在利用加速度值推算位移時(shí)，必須對(duì)加速度數(shù)據(jù)進(jìn)行濾波處理。

2 位移測(cè)算方法解析

2.1 理論計(jì)算方法

當(dāng)抽油桿上下運(yùn)動(dòng)時(shí)，加速度傳感器輸出A(t)的表達(dá)式為:

式中:a(t)為抽油桿運(yùn)動(dòng)的加速度;g(t)為重力加速度。

由于抽油桿上下作周期性運(yùn)動(dòng)，其速度呈周期性變化，即:

式中:v(t)為抽油桿運(yùn)動(dòng)的速度;T為一個(gè)沖程周期。

對(duì)a(t)積分，可得:

考慮到抽油桿上下作周期性運(yùn)動(dòng)的特點(diǎn)，若將v(0)看成下死點(diǎn)的初始速度，則v(0)=0。

于是，有:

位移同樣可以通過(guò)對(duì)速度積分得到，即:

若將下死點(diǎn)的位移S(0)作為參考點(diǎn)，則抽油桿上下作周期性運(yùn)動(dòng)的相對(duì)位移為:

由于重力加速度的變化很小，可看作常數(shù)。令其為G，它就是加速度的零點(diǎn)值。作零點(diǎn)校正時(shí)，一般采用測(cè)量前先存儲(chǔ)零點(diǎn)值、計(jì)算時(shí)將測(cè)量值減去零點(diǎn)值的方法。零點(diǎn)值由于各種原因會(huì)有微小變化，因此，用減去預(yù)先存儲(chǔ)的零點(diǎn)值的方式作零點(diǎn)校正后可能還存在殘差，且經(jīng)兩次積分后會(huì)產(chǎn)生很大的誤差。

在離散采樣的情況下，由加速度求解位移的算法步驟具體如下。

第一步，設(shè)N為一個(gè)沖程周期內(nèi)的加速度采樣點(diǎn)數(shù)，ΔTn為相鄰兩次采樣的時(shí)間間隔，an為第n個(gè)時(shí)間間隔的加速度采樣瞬時(shí)值，取下死點(diǎn)時(shí)的加速度值為第一個(gè)點(diǎn)，則計(jì)算校正零點(diǎn)值為:

第二步，計(jì)算第n個(gè)時(shí)間間隔的加速度采樣修正值為:

第三步，計(jì)算第n個(gè)時(shí)間間隔的速度為:

式中:t為采樣時(shí)間間隔(對(duì)于本測(cè)試系統(tǒng)，由于采樣頻率為 50 Hz，則 t=0.02 s)。若 n=1，則 vn=0。

第四步，計(jì)算速度的校正零點(diǎn)值為:

計(jì)算第n個(gè)時(shí)間間隔的速度修正值為:

第五步，計(jì)算第n個(gè)時(shí)間間隔的位移為:

2.2 濾波方法

在得到加速度測(cè)量值后，要計(jì)算抽油桿運(yùn)動(dòng)的相對(duì)位移，還必須確定積分求速度時(shí)的邊界條件——即速度的零點(diǎn)，一般是選擇光桿運(yùn)動(dòng)的下死點(diǎn)。

由于加速波形包含有大量的雜波分量，直接由原始加速度測(cè)量數(shù)據(jù)推算下死點(diǎn)是很困難的，因此，需要對(duì)加速度波形進(jìn)行濾波，再利用其判斷下死點(diǎn)。考慮到實(shí)現(xiàn)算法的有效性和實(shí)用性，選擇滑動(dòng)值濾波法對(duì)加速度波形進(jìn)行濾波［8－9］。

最簡(jiǎn)單的滑動(dòng)平滑算法是矩形或無(wú)加權(quán)的滑動(dòng)平均平滑，其實(shí)現(xiàn)方法是將時(shí)間序列中的每一個(gè)點(diǎn)用其相鄰的m個(gè)點(diǎn)的均值來(lái)代替。對(duì)于3點(diǎn)平滑(m=3)，則有:

式中:j=2，…，n－1;Sj為平滑時(shí)間序列的第j個(gè)點(diǎn);Yj為原始時(shí)間序列的第j個(gè)點(diǎn);n為時(shí)間序列的總點(diǎn)數(shù)。

平滑寬度m通常為奇數(shù)。如果時(shí)間序列包含的噪聲為白噪聲(即噪聲在各個(gè)頻段均勻分布)，噪聲標(biāo)準(zhǔn)差為s，則經(jīng)過(guò)一次無(wú)加權(quán)的滑動(dòng)平均平滑濾波后，其白噪聲的標(biāo)準(zhǔn)差變?yōu)?/p>

三角平滑是在矩形平滑的基礎(chǔ)上使用了加權(quán)平滑函數(shù)。對(duì)于一個(gè)5點(diǎn)的三角平滑，則有:

式中:j=3，…，n－2。

從式(14)可以看出，5點(diǎn)的三角平滑相當(dāng)于將3點(diǎn)矩形平滑的過(guò)程重復(fù)兩遍。因此，與矩形平滑相比，三角平滑能夠更有效地濾除高頻噪聲。經(jīng)過(guò)一次三角平滑濾波后，其白噪聲的標(biāo)準(zhǔn)差變?yōu)?/p>

當(dāng)平滑寬度越大時(shí)，時(shí)間序列的噪聲抑制程度也越大，與此同時(shí)，平滑過(guò)程對(duì)信號(hào)的干擾程度也越大。最優(yōu)的平滑寬度取決于信號(hào)的形狀和采樣間隔。對(duì)于尖峰型的信號(hào)，平滑比是很重要的一項(xiàng)參數(shù)。平滑比即平滑寬度和尖峰半寬的比值，一般來(lái)說(shuō)，增加平滑比會(huì)提高信噪比，不會(huì)對(duì)尖峰的位置造成影響，但會(huì)造成尖峰幅值與寬度的減小。原則上，平滑比不應(yīng)超過(guò)1，否則平滑后的時(shí)間序列的尖峰將會(huì)失真。如果時(shí)間序列平滑的目的是測(cè)試尖峰的位置，則可將平滑比設(shè)得大一些(介于0.5和1之間);而如果時(shí)間序列平滑的目的是更好地測(cè)試尖峰高度和寬度的真實(shí)值，則平滑比不應(yīng)大于0.2。

基于上述原則，選取平滑比為0.1，設(shè)“尖峰”半寬為20點(diǎn)，則平滑寬度為200點(diǎn)［10］，得到的平滑加速度曲線如圖3所示。

圖3 平滑濾波后的光桿加速度曲線Fig.3 The acceleration of polished rod after smooth filtering

2.3 下死點(diǎn)與沖程周期的確定

加速度的峰值點(diǎn)即為光桿運(yùn)動(dòng)的下死點(diǎn)。考慮到平滑后的時(shí)間序列會(huì)帶來(lái)信息損失，故先通過(guò)尋找平滑一階導(dǎo)數(shù)的下行過(guò)零點(diǎn)來(lái)計(jì)算得到峰值的位置;再通過(guò)最小二乘曲線擬合法，對(duì)未經(jīng)平滑的原始時(shí)間序列中峰值附近的數(shù)據(jù)段進(jìn)行處理，進(jìn)而計(jì)算峰值的位置;從而計(jì)算加速度時(shí)間序列的峰值。

時(shí)間序列的一階導(dǎo)數(shù)采用如下公式計(jì)算:

式中:X'j和Y'j為第j個(gè)點(diǎn)導(dǎo)數(shù)的X坐標(biāo)和Y坐標(biāo)值;n為時(shí)間序列的點(diǎn)數(shù);ΔX為X坐標(biāo)點(diǎn)與其相鄰點(diǎn)的差值。對(duì)于等間隔的三個(gè)相鄰點(diǎn)，也可采用如下公式計(jì)算某點(diǎn)平均斜率，即:

兩個(gè)相鄰峰值點(diǎn)之間的時(shí)間差即為每個(gè)沖程的周期。

2.4 位移測(cè)試結(jié)果分析

利用上述算法計(jì)算得到的一個(gè)沖程的光桿位移曲線、下死點(diǎn)提前100點(diǎn)所得的位移曲線以及下死點(diǎn)滯后100點(diǎn)所得的位移曲線分別如圖4所示。

圖4 光桿位移曲線Fig.4 The curves of displacement of polished rod

由圖4可知，確定邊界條件對(duì)光桿位移計(jì)算是十分重要的。

3 結(jié)束語(yǔ)

利用筆者設(shè)計(jì)的光桿加速度無(wú)線數(shù)據(jù)采集試驗(yàn)系統(tǒng)，獲取了抽油機(jī)光桿的加速度數(shù)據(jù)。針對(duì)所獲取數(shù)據(jù)的特點(diǎn)，研究了數(shù)據(jù)的濾波方法，并給出相關(guān)參數(shù)的選擇依據(jù)，進(jìn)而通過(guò)尋找平滑一階導(dǎo)數(shù)的下行過(guò)零點(diǎn)來(lái)計(jì)算得到峰值的位置;同時(shí)，通過(guò)最小二乘曲線擬合法獲得下死點(diǎn)位置，最終利用考慮抽油桿運(yùn)動(dòng)周期性的加速度二次積分算法得到位移。通過(guò)對(duì)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)的計(jì)算與分析，證明了該測(cè)算方法的有效性。

［1］劉益江，張學(xué)臣，李偉，等.抽油井示功圖綜合解釋［J］.油氣田地面工程，2007，26(8):3－5.

［2］高鵬.抽油井示功圖測(cè)試方法改進(jìn)探討［J］.內(nèi)蒙古石油化工，2007(12):55－56.

［3］袁蓓蓓.KSZ_1光柵長(zhǎng)度變送器［J］.自動(dòng)化儀表，1995，16(1):23－25.

［4］沈承虎，容太平，蔡德鈞，等.加速度傳感器在示功儀中的應(yīng)用［J］.儀表技術(shù)與傳感器，2000(1):40－42.

［5］于云華，張家珍，時(shí)海濤，等.基于加速度傳感器的油井示功圖位移測(cè)量技術(shù)研究［J］.電子產(chǎn)品世界，2009(9):25－30.

［6］容太平，沈承虎，袁中平，等.用加速度傳感器測(cè)量位移的原理與誤差分析［J］.華中理工大學(xué)學(xué)報(bào):自然科學(xué)版，2000，25(5):58－60.

［7］鹿麟，林凌，李剛.ADXL203型雙軸加速計(jì)在傾斜度測(cè)量中的應(yīng)用［J］.國(guó)外電子元器件，2007(7):61－64.

［8］劉藝柱，楊瑞蘭.采用滑動(dòng)平均濾波法提高硬幣識(shí)別準(zhǔn)確率的研究［J］.制造業(yè)自動(dòng)化，2010，32(1):42－44.

［9］胡松，江小煒，楊光.滑動(dòng)平均濾波在微弱脈沖信號(hào)檢測(cè)中的應(yīng)用［J］.計(jì)算機(jī)與數(shù)字工程，2007，35(10):169－171.

［10］鹿麟，林凌，李剛.滑動(dòng)平均數(shù)字濾波參數(shù)研究［J］.集美大學(xué)學(xué)報(bào):自然科學(xué)版，2006，11(4):381－384.