張 濤 王凱凱 李玉梅 周文彬 鄧靖宇
(1.北京信息科技大學(xué) 高動態(tài)導(dǎo)航技術(shù)北京市重點實驗室 2.北京信息科技大學(xué) 現(xiàn)代測控技術(shù)教育部重點實驗室 3.中國石油井控應(yīng)急救援響應(yīng)中心)
鉆井作業(yè)是油氣開采的重要過程,如何處理好、分析好、利用好鉆井過程中的隨鉆測量數(shù)據(jù)是安全、高效進(jìn)行鉆井作業(yè)的重中之重[1-4]。井下環(huán)境復(fù)雜,底部鉆具組合在作業(yè)過程中會受到各種激勵作用,產(chǎn)生復(fù)雜的振動現(xiàn)象。通過準(zhǔn)確的隨鉆測量數(shù)據(jù)對井下鉆頭運動狀態(tài)進(jìn)行分析、反演和控制,能夠減輕鉆頭異常振動的危害,提高石油天然氣的開采效率。
對于振動數(shù)據(jù)的位移重建,許自順等[5]研究了時域積分在振動信號中的應(yīng)用;王劍等[6]研究了振動加速度數(shù)值積分的Lagrange多項式擬合方法,發(fā)現(xiàn)在低采樣頻率下此方法優(yōu)于辛普森方法;成勛等[7]利用小波去噪的方法對振動數(shù)據(jù)進(jìn)行去噪,得到良好的位移重建數(shù)據(jù);許曉紅等[8]使用組合辛普森的方式對振動的時域積分進(jìn)行了研究。
對于近鉆頭的三軸振動數(shù)據(jù)處理分析,國內(nèi)外學(xué)者進(jìn)行了大量的研究。A.HOHL等[9-11]對井下高頻扭轉(zhuǎn)振動進(jìn)行了數(shù)學(xué)建模分析,提出了優(yōu)化井下鉆具組合和鉆井操作的建議,以改進(jìn)鉆井性能;CHEN W.等[12]利用深度學(xué)習(xí)的方法識別底部鉆具的工作狀態(tài),加快了數(shù)據(jù)解釋的工作效率;E.CAYEUX等[13]根據(jù)井下高頻數(shù)據(jù)重建當(dāng)前測量工具所處位置的三維管道運動,以告知鉆井作業(yè)團(tuán)隊當(dāng)前鉆頭狀態(tài);蘇曉眉等[14]通過K-Means聚類算法對井下三軸振動數(shù)據(jù)進(jìn)行分析,成功實現(xiàn)了對沉砂卡鉆的預(yù)測;任凱[15]根據(jù)跳鉆、渦動和黏滑振動的動力學(xué)模型和井下三軸振動數(shù)據(jù)的分析,提出了振動工況的分級和減振措施。
伴隨著虛擬現(xiàn)實技術(shù)的日漸成熟,虛擬現(xiàn)實技術(shù)的高沉浸性、高交互性等優(yōu)點為數(shù)字虛擬井場帶來了新的發(fā)展。筆者依據(jù)隨鉆測量過程中振動數(shù)據(jù)的特點,提出一種滑動窗口的積分方法,并建立單自由度振動模型進(jìn)行分析驗證。利用該滑動窗口積分方法將三軸振動數(shù)據(jù)轉(zhuǎn)換為位移數(shù)據(jù),建立虛擬地下鉆井場景,利用位移數(shù)據(jù)和歷史數(shù)據(jù)驅(qū)動井下鉆頭運動,形成了鉆頭反演運動和數(shù)據(jù)解釋新方式,對虛擬現(xiàn)實井場的完善具有一定的指導(dǎo)意義。
加速度、速度和位移是衡量振動強度的3個指標(biāo),但由于作業(yè)場景的限制,常常不能同時測量這3種參數(shù)。為了獲得其他參數(shù),一般需要微分或積分操作。井下底部鉆具組合的振動由三軸加速度計測得,為了重建鉆頭的位移,由加速度信號積分得到位移信號是比較可行的方案。
直接時域積分和直接頻域積分是2種常用的積分方法,前一種方法是在時間域中通過累積梯形求積公式或者其他求積公式計算:
(1)
式中:a、b是積分區(qū)域的上下限;n為將[a,b]等分為多少個間距;f(x)是被積信號。
還有一種積分方法則是先將數(shù)據(jù)進(jìn)行頻域變換,在頻域中進(jìn)行積分計算,然后再進(jìn)行反傅里葉變換得到結(jié)果[16],即:
(2)
式中:S(t)是時間段內(nèi)位移數(shù)據(jù)信號;wk=2πkFs/N;Fs是采樣頻率;N是采樣數(shù)據(jù)量;A(k)是加速度信號的頻譜;k是非負(fù)整數(shù)。
根據(jù)井下三軸振動信號的特點:隨著鉆進(jìn)過程推進(jìn),會在不同時間段產(chǎn)生新噪聲,即便經(jīng)過濾波處理也無法完全去除噪聲的影響。為了減少其他時間段噪聲對當(dāng)前時間段積分的影響,本文提出一種滑動窗口積分方法:選取合適的窗口大小,在時間域內(nèi)以相同時間步長滑動,依次對窗口內(nèi)的加速度信號進(jìn)行2次頻域積分得到該段的位移重建數(shù)據(jù)。
設(shè)a(t)+b1(t)+b2(t)為帶噪加速度信號。其中,a(t)為加速度信號,t0到t1窗口內(nèi)沒有噪聲;b1(t)為t1到t2窗口的噪聲;b2(t)是t2到t3窗口的噪聲。則其直接頻域積分得到的位移重建信號為:
(3)
滑動窗口積分得到的位移重建信號為:
(4)
式中:A1(k)、A2(k)、A3(k)分別為3個時間段的加速度信號;S1(t)、S2(t)、S3(t)分別為3個時間段的位移重建信號;B1(k)和B2(k)分別為b1(t)與b2(t)的頻譜。
對于三軸振動數(shù)據(jù)來說,直接頻域積分需要在整段頻域內(nèi)進(jìn)行2次積分,會擴(kuò)大噪聲的影響,得到的位移絕對值偏大,而滑動窗口積分則可以在當(dāng)前窗口內(nèi)積分,不受窗口外噪聲影響,相當(dāng)于又做了一次濾波,積分過程受噪聲影響較小,所以位移值能更好地貼近真實位移。
為了能夠更好地分析滑動窗口積分、時域積分和頻域積分的誤差情況,引入平均峰值誤差、平均最大值誤差和均方誤差這3個指標(biāo)[17]。
平均峰值誤差可以看出積分信號和實際信號的整體最大峰值和最小峰值的情況:
(5)
式中:A(t)是重建的位移數(shù)據(jù);S(t)是實際的位移數(shù)據(jù)。
平均最大值誤差可以看出重建數(shù)據(jù)和實際數(shù)據(jù)的局部最大峰值和最小峰值的情況:
(6)
由均方誤差可以看出重建數(shù)據(jù)和實際數(shù)據(jù)之間偏離的情況:
(7)
利用Simulink軟件建立單自由度振動模型(見圖1),以驗證算法的可行性。其中作用力是頻率為10 Hz的正弦波信號,施加時間為12 s,模型采樣頻率為1 000 Hz。使用單位門信號控制噪聲的作用時間,從第1 s開始到第12 s每秒加入一個噪聲,每次增加5 Hz,頻率從10~60 Hz分布。通過傳感器模塊獲取振動模型的加速度、速度、位移數(shù)據(jù),并將模擬信號轉(zhuǎn)換成數(shù)字信號,同時輸出保存加速度、速度、位移數(shù)據(jù)。
圖1 單自由度振動模型Fig.1 Single-degree-of-freedom vibration model
以單自由度振動模型為基礎(chǔ),仿真時間12 s,采樣頻率為1 000 Hz,得到120 000個實際加速度和位移采樣數(shù)據(jù)點。
滑動窗口積分相當(dāng)于以整體時域內(nèi)的若干段分別代替整體進(jìn)行傅里葉變換,然后再分別進(jìn)行頻域積分,所以窗口的選擇決定了每次積分段內(nèi)的噪聲數(shù)量。鑒于窗口太小會導(dǎo)致頻譜能量泄露增加,造成積分曲線窗口處不穩(wěn)定,所以以采樣頻率為滑動窗口的初始大小L,設(shè)定滑動窗口積分的窗口大小分別為iL(i=1,2,3,4,5,6),確定窗口大小和均方誤差的關(guān)系,以確定最好的積分效果。同時計算該窗口內(nèi)主頻率幅值和該段窗口總幅值的比,作為窗口大小的依據(jù),得到不同窗口大小下位移重建的均方誤差關(guān)系曲線,如圖2所示。由圖2可以發(fā)現(xiàn),當(dāng)主頻率幅值比降低,即當(dāng)前窗口內(nèi)的噪聲數(shù)量增多時,滑動窗口積分的均方誤差逐漸增大,滑動窗口積分對實際位移數(shù)據(jù)的擬合效果在逐步下降。這是由于滑動窗口內(nèi)的噪聲增多,導(dǎo)致正常加速度信號在2次頻域積分內(nèi)的積分區(qū)間擴(kuò)大,造成了其他頻率的累積誤差,降低了數(shù)據(jù)擬合的能力。
圖2 幅值比與均方誤差關(guān)系圖Fig.2 Relation between amplitude ratio and mean square error
為此,本文取滑動窗口大小為采樣頻率的大小,對同一段加速度信號進(jìn)行位移重建積分。對3種積分方法使用同樣的積分計算流程:首先對加速度信號進(jìn)行多項式擬合消除趨勢項,一次積分得到速度信號;然后再使用多項式擬合消除一次積分產(chǎn)生的趨勢項,二次積分得到位移信號。為了方便結(jié)果曲線的觀察,分別放大觀察1~5 s、5~8 s、8~12 s的曲線,如圖3~圖5所示。
圖3 1~5 s秒位移重建曲線Fig.3 1~5 s displacement reconstruction curve
圖4 5~8 s位移重建曲線Fig.4 5~8 s displacement reconstruction curve
從圖3~圖5可以看出:直接時頻域積分都明顯受到噪聲信號的影響,發(fā)生“擺偏”現(xiàn)象,而滑動窗口積分的位移重建相對較好,能夠擬合實際位移曲線。
為了量化誤差,引入平均峰值誤差、平均最大值誤差和均方誤差這3個指標(biāo),得到結(jié)果如表1所示。
表1 3種積分位移重建誤差Table 1 Displacement reconstruction error of 3 kinds of integrations
由于多個噪聲的存在,即便做了趨勢項消除,也無法完全根除積分期間產(chǎn)生的非線性趨勢項的影響,這使直接時域積分重建的位移數(shù)據(jù)發(fā)生“擺偏”。對于直接頻域積分來說,不進(jìn)行頻率截斷,即不將主頻率之外的頻率置零處理,直接頻域積分會受到頻率相近噪聲的影響,和時域積分一樣產(chǎn)生“擺偏”。由于計算值和實際值的差值過大,部分區(qū)域誤差達(dá)到幾毫米,所以導(dǎo)致時頻域的平均最大值誤差非常大。
而滑動窗口積分的各個窗口之間的積分過程相對獨立,對每個積分窗口來說相當(dāng)于做了頻率截斷處理,降低了當(dāng)前積分窗口內(nèi)的噪聲數(shù)量,因而能夠相對較好地擬合實際的位移信號。由均方誤差也可以看出,直接時頻域積分的均方誤差是滑動窗口積分的20倍左右。這充分證明了滑動窗口積分運用于振動位移重建的可行性。但是由于矩形窗口的加入會導(dǎo)致頻譜能量泄露,所以滑動窗口積分在滑動窗口交界處和實際位移信號相似度較小。
本文采用的三軸振動數(shù)據(jù)為新疆某油田鉆井過程中的實測數(shù)據(jù),數(shù)據(jù)測量點在近鉆頭測量短節(jié)。該測量數(shù)據(jù)記錄了井下鉆具組合在鉆井過程中的鉆速、鉆壓、扭矩和三軸加速度數(shù)據(jù),信號采樣頻率為400 Hz。
為了重建底部鉆具組合的振動位移信號,首先對井下測量的三軸振動電壓數(shù)據(jù)進(jìn)行數(shù)據(jù)的預(yù)處理操作,以得到加速度數(shù)據(jù);然后對三軸加速度數(shù)據(jù)以采樣頻率為窗口大小分別進(jìn)行2次滑動窗口積分得到位移數(shù)據(jù),如圖6所示。從圖6可以發(fā)現(xiàn),滑動窗口積分較好地由振動加速度數(shù)據(jù)恢復(fù)了振動的位移數(shù)據(jù),位移數(shù)據(jù)大致在毫米級波動,處于正常鉆進(jìn)過程。
圖6 三軸位移重建Fig.6 Triaxial displacement reconstruction
為了完整與友好地展示底部鉆具組合的運動狀態(tài),建立基于虛擬現(xiàn)實技術(shù)的地下鉆井場景。該地下鉆井場景與虛擬井場中的鉆井平臺聯(lián)動,屬于地下演示模塊。
首先,需要明確的是地下鉆井場景中法向、切向和軸向振動的正方向。以底部鉆具組合模型的前后方向為法向振動方向,其中正方向向前;以底部鉆具組合模型的左右方向為切向振動方向,其中正方向向右;以底部鉆具組合模型的上下方向為軸向振動方向,其中正方向向下。三軸振動正方向如圖7所示。
圖7 三軸振動正方向Fig.7 Positive direction of triaxial vibration
其次,虛擬地下鉆井場景應(yīng)當(dāng)包含以下幾個模塊:數(shù)據(jù)流展示模塊、底部鉆具組合演示模塊、地層模塊及用戶交互模塊。
數(shù)據(jù)流展示模塊:以實際數(shù)據(jù)為基礎(chǔ),通過實時動態(tài)折線圖的方式進(jìn)行鉆井過程中實測的鉆速、三軸振動等信號的數(shù)據(jù)回放。為了實現(xiàn)折線圖的動態(tài)移動方式,以相同的時間幀傳入數(shù)據(jù)流,將數(shù)據(jù)歸一化處理,并依次繪制數(shù)據(jù),實現(xiàn)舊數(shù)據(jù)的移除與新數(shù)據(jù)的注入,從而實現(xiàn)動態(tài)折線圖的更新繪制。
底部鉆具組合演示模塊:以實際振動數(shù)據(jù)為基礎(chǔ),選取底部鉆具組合模型的本地坐標(biāo)中心為振動中心,根據(jù)三軸振動數(shù)據(jù)重建的三維位移數(shù)據(jù)實時驅(qū)動模型的振動,反演鉆具組合的中心偏移情況;以實際鉆速數(shù)據(jù)驅(qū)動鉆具組合的轉(zhuǎn)動速度。同時,考慮到虛擬地下場景運行環(huán)境的顯示穩(wěn)定性和運行設(shè)備的性能,將底部鉆具組合的驅(qū)動頻率設(shè)為200幀,即在虛擬地下場景中每2 s才能運行完1 s實際數(shù)據(jù)。圖8為底部鉆具組合模型。
圖8 底部鉆具組合模型Fig.8 Bottomhole assembly model
地層模塊:為了減少地層空間建模體積,提高虛擬鉆井場景的性能,通過地層向上移動的方式反襯鉆頭鉆進(jìn)的過程。構(gòu)建動態(tài)地層著色程序(shader),根據(jù)鉆速大小改變地層上移速度,鉆速快,地層上移快,鉆速降低,降低地層上移速度,以動態(tài)模擬底部鉆具組合鉆進(jìn)的畫面,優(yōu)化視覺效果。圖9是地層空間的俯視圖。
圖9 地層空間俯視圖Fig.9 Top view of formation space
用戶交互模塊:用戶通過虛擬手柄選擇是否啟動底部鉆具組合運動反演。選擇是,則暫時切斷地下鉆井場景與地上鉆井平臺的聯(lián)動,轉(zhuǎn)換為鉆頭運動反演情景,確保虛擬地下環(huán)境的運行效率。
地下鉆井場景的虛擬場景如圖10所示。圖10中:中間是底部鉆具組合模型,實時通過實際數(shù)據(jù)驅(qū)動鉆具組合的運動反演;左右兩邊的UI是實際數(shù)據(jù)流折線圖,左邊依次是法向振動加速度和切向振動加速度,右邊依次是軸向振動加速度和鉆速,展現(xiàn)當(dāng)前歷史數(shù)據(jù)的大小。
圖10 虛擬地下鉆井場景Fig.10 Virtual underground drilling scene
基于以上工作,本文建立了基于虛擬現(xiàn)實技術(shù)的底部鉆具的運動反演場景,用戶可以佩戴虛擬現(xiàn)實設(shè)備進(jìn)入地下鉆井場景,通過相應(yīng)的手柄交互,實現(xiàn)底部鉆具組合的運動反演觀測。
本文針對振動數(shù)據(jù)濾波后仍存在噪聲干擾的問題,提出滑動窗口積分的位移重建方法,首先建立單自由度振動模型進(jìn)行算法的驗證;然后將滑動窗口積分方法用于井下振動數(shù)據(jù)的位移重建,建立了基于虛擬現(xiàn)實技術(shù)的底部鉆具組合的運動反演觀測場景,用戶可以通過虛擬現(xiàn)實設(shè)備進(jìn)入場景。研究結(jié)果表明:
(1)通過原理分析與仿真對比的研究手段可以發(fā)現(xiàn),直接時頻域積分在加速度信號的積分過程中受噪聲干擾較大,而本文提出的滑動窗口積分相當(dāng)于在直接頻域積分的基礎(chǔ)上做了一次頻率截斷,降低了當(dāng)前窗口內(nèi)的噪聲數(shù)量,能夠更好地擬合實際位移曲線。
(2)利用滑動窗口積分的方法對井下振動信號進(jìn)行位移重建,構(gòu)建基于虛擬現(xiàn)實技術(shù)的地下鉆井場景,根據(jù)鉆速數(shù)據(jù)和位移數(shù)據(jù)驅(qū)動底部鉆具組合運動反演,利用虛擬現(xiàn)實設(shè)備即可進(jìn)行運動觀測和數(shù)據(jù)分析。表明了通過實際鉆井?dāng)?shù)據(jù)重現(xiàn)鉆井過程的可行性,拓展了虛擬現(xiàn)實技術(shù)在石油鉆井井場的應(yīng)用范圍。