曲漢武 王智明 張 崢 邵天宇 郭心宇

（中海油田服務(wù)股份有限公司，北京 101149）

隨鉆測井是指測井儀器在鉆進時對井下的工程參數(shù)和地質(zhì)參數(shù)進行測量并上傳。在鉆進過程中，井下測量傳感器測得工程參數(shù)和地層參數(shù)。這些測得的參數(shù)（通常為模擬信號）通過數(shù)據(jù)編碼器轉(zhuǎn)換為數(shù)字信號，經(jīng)過控制電路調(diào)制后傳遞給驅(qū)動電路，驅(qū)動控制電動機運動。電動機按照控制電路給定的控制信號運動，帶動擺動閥泥漿脈沖發(fā)生器轉(zhuǎn)子按照相應(yīng)的軌跡旋轉(zhuǎn)或擺動。擺動閥脈沖發(fā)生器的定轉(zhuǎn)子剪切流經(jīng)的流體，產(chǎn)生泥漿壓力波信號。這些泥漿壓力波信號，經(jīng)過鉆桿內(nèi)泥漿傳輸?shù)降孛媪⒐?。?shù)據(jù)采集系統(tǒng)采集地面立管上壓力傳感器的壓力信號，通過解調(diào)系統(tǒng)解析井下的壓力信號，并將傳輸?shù)哪酀{脈沖信號轉(zhuǎn)換為井下工程參數(shù)和地層參數(shù)[1-5]。

對于擺動閥脈沖器而言，軸系轉(zhuǎn)動慣量和粘滯摩擦系數(shù)可以用于分析和求解電動機的機械運動方程。這對于優(yōu)化擺動閥轉(zhuǎn)子運動控制軌跡和擺動閥脈沖器本身機械結(jié)構(gòu)參數(shù)具有十分重要的意義[6-10]。

目前，國內(nèi)外對電動機轉(zhuǎn)子軸系轉(zhuǎn)動慣量的測量普遍采用基于電機轉(zhuǎn)子連續(xù)旋轉(zhuǎn)的方式[11-12]，而擺動閥脈沖器由于連接電動機轉(zhuǎn)子和擺動閥閥片的扭桿尾端被固定鎖死（如圖1 所示），且扭桿本身存在機械限位（如圖2 所示），導(dǎo)致電動機只能在有限的角度內(nèi)擺動，無法進行連續(xù)旋轉(zhuǎn)，無法用常規(guī)的方法測量軸系的轉(zhuǎn)動慣量和粘滯摩擦系數(shù)。

圖1 擺動閥脈沖器結(jié)構(gòu)

圖2 擺動閥脈沖器扭桿機械限位結(jié)構(gòu)

文章提出了一種通過電動機在有限角度空間內(nèi)做周期性擺動來測量轉(zhuǎn)子轉(zhuǎn)動慣量和粘滯摩擦系數(shù)的方法，突破了傳統(tǒng)的電動機需要在連續(xù)旋轉(zhuǎn)狀態(tài)下才能測量轉(zhuǎn)動慣量的局限。

1 擺動閥脈沖器的機械平衡方程

擺動閥脈沖器在空氣中擺動（此時水力轉(zhuǎn)矩為0）的機械運動方程為

它可以描述為輸出的電磁轉(zhuǎn)矩等于慣性力、扭軸的彈性力、粘滯摩擦力之和。其中：擺動閥位于中心位置為t=0 的初始時刻，J為轉(zhuǎn)子的轉(zhuǎn)動慣量，ω為擺動閥脈沖器的機械角速度，ks為扭軸的彈性系數(shù)，θ為擺動閥片所在位置的機械角度，kb為粘滯摩擦系數(shù)，Te為永磁同步電機輸出的電磁轉(zhuǎn)矩。

擺動閥脈沖器選擇的電動機一般為表貼式永磁同步電機，控制方法采用最大轉(zhuǎn)矩電流比控制，即控制過程中保持直軸電流Id為0，因此Te可以表示為

式中：km為電動機的電磁轉(zhuǎn)矩系數(shù)，與電動機自身的結(jié)構(gòu)參數(shù)有關(guān)；Iq為交軸電流，可以由三相電流Ia、Ib、Ic經(jīng)過Clark 變換（式3）和Park 變換（式4）計算后得到。

式中：Iα和Iβ為兩相靜止坐標系下的電流；Id和Iq分別為旋轉(zhuǎn)坐標系下的直軸電流和交軸電流；φe為電機A 相繞組方向轉(zhuǎn)動至轉(zhuǎn)子N 極方向所經(jīng)過的電角度。

2 軸系轉(zhuǎn)動慣量和粘滯摩擦系數(shù)的測量原理

由擺動閥脈沖器在空氣中擺動的機械平衡方程式（1）可知，擺動閥轉(zhuǎn)子軸系的轉(zhuǎn)動慣量J和粘滯摩擦系數(shù)kb可以通過測量多組位置角度θ、角速度ω和交軸電流Iq求解方程計算得到。然而，擺動閥脈沖器本身由于存在機械限位，且扭軸尾端被鎖死，電機本身不能通過一般調(diào)速的控制方法連續(xù)旋轉(zhuǎn)，只能在有限的角度空間內(nèi)做反復(fù)式擺動，導(dǎo)致實際測量過程中無法獲得穩(wěn)態(tài)的轉(zhuǎn)速ω和交軸電流Iq。通過測量得到瞬態(tài)角度θ、角速度ω以及交軸電流Iq均存在較大的誤差，且數(shù)據(jù)采集過程中可能存在數(shù)據(jù)不同步問題，因此測量得到的轉(zhuǎn)動慣量J和粘滯摩擦系數(shù)kb均存在較大誤差，無法指導(dǎo)工程應(yīng)用。

文章利用擺動閥脈沖器在空氣中擺動的機械平衡方程，通過控制脈沖器在有限空間內(nèi)做正弦擺動，采集多組瞬態(tài)角度θ、角速度ω以及交軸電流Iq數(shù)據(jù)，通過對其在整數(shù)周期內(nèi)進行加權(quán)計算的方法來消除采集誤差，并且利用不同頻率下多組數(shù)據(jù)線性擬合的方法來消除控制誤差，最終達到精確測量轉(zhuǎn)動慣量J和粘滯摩擦系數(shù)kb的目的。

如圖3 所示，對驅(qū)動擺動閥脈沖器擺動的永磁同步電機采用位置閉環(huán)的控制方式進行控制時，通過調(diào)節(jié)控制器的參數(shù)，可以確保擺動閥片的實際位置角度較為準確地跟蹤輸入的位置角度。

圖3 擺動閥脈沖器的位置閉環(huán)控制方式

當控制誤差為0 時，輸入角度隨時間變化的公式為

式中：θm為擺動幅度；f為擺動頻率，設(shè)擺動閥位于中心位置為t=0 的初始時刻。

設(shè)擺動周期T=1/f，將式（1）兩邊同時乘以θ，在時間T內(nèi)求積分，即可約掉角速度ω，得

式中：km、ks、θm、f等均為已知量；θ為測量值；Iq可以通過測量三相電流后計算，可由式（3）、式（4）得到。

同理，將式（1）兩邊同時乘以ω后，在時間T內(nèi)求積分，可以約掉等號右邊的第1 項和第2 項，得

式中：km為已知量；w、Iq可以通過測量和計算得到。

實際應(yīng)用中，測量的Iq、θ為離散量。設(shè)定采樣周期Ts=T/N，則式（6）經(jīng)離散化變形后可以得到轉(zhuǎn)動慣量的計算公式為

式（7）經(jīng)過離散化變形后，可以得到粘滯摩擦系數(shù)的計算公式為

式（8）和式（9）是在控制誤差為0 的假設(shè)條件下得到的。然而，實際上受控制精度的限制，控制誤差不可能完全為0，如圖4 所示。

為進一步減小誤差，可以采集不同頻率下的多組數(shù)據(jù)，然后利用式（8）和式（9），通過線性擬合的方式得到轉(zhuǎn)動慣量和粘滯摩擦系數(shù)。因此，引入x、y、z、w這4 個量，令

依然控制擺動閥輸入的位置波形θ*=θmsin2πft，選擇M種擺動頻率f∈[f1,f2,…,fm]，在每種頻率下擺動相同的時間（保證為M種頻率下所有擺動周期的整數(shù)倍），將采樣值分別代入式（8）～式（13），計算可得

式中：xj(j=1,2,…,M)、yj(j=1,2,…,M)、zj(j=1,2,…,M)、wj(j=1,2,…,M)為對應(yīng)M種擺動頻率下計算的x、y、z、w的值。對式（14）進行最小二乘擬合，可以得到擺動閥轉(zhuǎn)子的軸系轉(zhuǎn)動慣量值J。對式（15）進行最小二乘擬合，可以得到擺動閥擺動過程中的粘滯摩擦系數(shù)kb。

3 測量流程

步驟1：對擺動閥脈沖器的永磁同步電機采用位置閉環(huán)和最大轉(zhuǎn)矩比電流的控制方式進行控制。

步驟2：設(shè)置輸入位置為θ*=θmsin2πft的運動軌跡進行運動，以Ts為采樣周期記錄時間Ta（Ta為擺動周期的整數(shù)倍，且Ta=NTs）內(nèi)位置傳感器采樣的擺動位置角度θi(i=1,2,…,N)和三相電流Iai、Ibi、Ici(i=1,2,…,N)，進而計算出對應(yīng)的ωi和Iqi(i=1,2,…,N)。

步驟3：利用式（10）、式（11）、式（12）、式（13），計算得到x1、y1、z1、w1。

步驟4：變化擺動頻率fi∈[f1,f2,…,fm]，重復(fù)步驟2 和步驟3，得到每種頻率下的xj、yj、zj、wj(j=1,2,…,M)。

步驟5：根據(jù)式（14）和式（15），利用最小二乘擬合的方法，計算得到擺動閥轉(zhuǎn)子的軸系轉(zhuǎn)動慣量值J和擺動閥擺動過程中的粘滯摩擦系數(shù)kb。

4 實例驗證

選用中海油服自主研發(fā)的A 型擺動閥脈沖器進行實驗測試，如圖5 所示。擺動閥脈沖器的相關(guān)參數(shù)，如表1 所示。利用最大轉(zhuǎn)矩比電流的控制方式控制擺動閥轉(zhuǎn)子，按照式（5）的正弦運動軌跡進行擺動，分別選擇擺動頻率范圍為fi∈[6,7,8,…,25]的20 種頻率進行采樣，采樣時記錄每種頻率fi下每個采樣點j處的位置θ(i,j)、角速度ω(i,j)和q軸電流Iq(i,j)，利用式（10）～式（13）分別計算相應(yīng)的xi、yi、zi、wi(i=1,2,…,20)，分別對(xi,yi)和(zi,wi)做最小二乘擬合，擬合結(jié)果如圖6 和圖7 所示。

表1 實驗使用低的脈沖器參數(shù)表

圖5 擺動閥脈沖器測試試驗

圖6 轉(zhuǎn)動慣量的擬合曲線

圖7 粘滯摩擦系數(shù)的擬合曲線

從圖6 和圖7 可以看出，不同頻率下計算的x、y數(shù)據(jù)線性相關(guān)度較高，符合式（14）中兩者的正比關(guān)系。如表2 所示，擬合得到的軸系轉(zhuǎn)動慣量值為J為0.002 206 kg·m2，同SolidWorks 軟件建模計算得到轉(zhuǎn)動慣量值十分接近，驗證了方法計算的準確性。

表2 測量值與仿真值對比 單位：kg·m2

此外，從圖7 的擬合結(jié)果可知，計算的粘滯摩擦系數(shù)kb為0.002 510 1，不同頻率下計算的z、w兩組數(shù)據(jù)線性相關(guān)度同樣非常高，符合式（15）兩者的正比關(guān)系。

5 結(jié)語

通過電動機在有限角度空間內(nèi)做周期性擺動來測量轉(zhuǎn)子轉(zhuǎn)動慣量和粘滯摩擦系數(shù)的方法，突破了傳統(tǒng)電動機需要在連續(xù)旋轉(zhuǎn)狀態(tài)下才能測量轉(zhuǎn)動慣量的局限。通過擺動閥脈沖器實際測量案例驗證可知，測量方法測量計算得到的擺動閥脈沖器軸系轉(zhuǎn)子轉(zhuǎn)動慣量和粘滯摩擦系數(shù)具有較高的精度，可以廣泛應(yīng)用于擺動閥轉(zhuǎn)子運動控制軌跡的優(yōu)化及擺動閥脈沖器本身機械結(jié)構(gòu)參數(shù)的優(yōu)化。