龔大洪， 陳元鵬， 楊文娟， 楊金生， 王志玲， 張樹義， 馬鴻彥

(1中國石油集團(tuán)渤海鉆探工程有限公司定向井技術(shù)服務(wù)分公司 2中國石油渤海鉆探工程技術(shù)研究院 3中國石油集團(tuán)渤海鉆探工程有限公司第五鉆井工程分公司 )

在鉆井過程中井眼軌跡的延伸是一個動態(tài)的過程，如果想實時監(jiān)控井眼軌跡并且對其有個更直觀的認(rèn)識，則在實現(xiàn)井眼軌跡三維可視化的基礎(chǔ)上通過合理的誤差模型進(jìn)一步實現(xiàn)三維井眼軌跡誤差橢球可視化對叢式井井眼防碰或救援定向井連通等具有很好的指導(dǎo)意義[1-7]。

目前，運用最為廣泛的井眼軌跡不確定性誤差描述模型有系統(tǒng)誤差模型和ISCWSA誤差模型。系統(tǒng)誤差模型為Wollf[8]和 Wardt于1981年提出的井眼軌跡不確定性分析模型 ，簡稱WdW模型，該模型只提出了六個誤差源，更適用于研究基于陀螺類測斜工具測量垂直井和中等斜度井時系統(tǒng)誤差引起的井眼位置不確定性，具有一定的局限性。ISCWSA誤差模型為井眼測量精度工業(yè)導(dǎo)向委員會通過的一種誤差模型，它以Williamson[9]于2000年提出的基于磁性儀器的誤差分析模型為基礎(chǔ)并且還在不斷的完善，該模型提供了井眼軌跡不確定性分析的數(shù)學(xué)框架，能綜合考慮系統(tǒng)誤差、隨機(jī)誤差和全局誤差，被廣泛應(yīng)用到石油行業(yè)中。

一、 ISCWSA誤差模型及測量誤差源分析

1. ISCWSA模型

磁性測量工具測量的井眼軌跡受到很多誤差項的影響，這些誤差項可以歸集為幾個大的誤差源且這些誤差源是相互獨立的，ISCWSA模型假設(shè)誤差源在數(shù)理統(tǒng)計上服從正態(tài)分布，那么在磁性測量工具測量的井眼軌跡上的每個測點的誤差源確定該點的其中一個誤差。將誤差與井眼軌跡三維空間位置相結(jié)合，就能用向量把誤差表達(dá)出來。井眼軌跡上連續(xù)兩個測點的位置向量用計算位移來表示，定義磁性測量工具測量的井眼軌跡上第k-1個測點和第k個測點之間的位置增量為Δrk，則第l個測段的第k個測點上(觀察點以上的測點)由第i個誤差源引起的誤差為：

(1)

式中:ei,l,k—第i個誤差在第l個測段第k個測點NEV軸導(dǎo)致的位置矢量標(biāo)準(zhǔn)差；σi，l—第l組測量第i個誤差源的標(biāo)準(zhǔn)差;pk—第k個測點上的測量矢量,p=[DIA]T；?p/?εi—第i個誤差對測量矢量的影響(一個3×1向量)，即誤差的權(quán)重函數(shù)；dr/dp—測量向量p對要求點井眼位置矢量(NEV軸)r的影響(一個3×3矩陣)。

屬于第L個測段的觀察點K處第i個誤差源引起的誤差為：

(2)

2. 測量誤差源分析

井眼誤差大致可以分為三種類型：系統(tǒng)誤差、隨機(jī)誤差和全局誤差。

2.1 系統(tǒng)誤差及其對應(yīng)權(quán)重函數(shù)

在測量井眼軌跡的過程中系統(tǒng)誤差又可以再次劃分為已定系統(tǒng)誤差和未定系統(tǒng)誤差。已定系統(tǒng)誤差的方向、大小是已知的。未定系統(tǒng)誤差由于其方向未知，大小也未知，所以一般情況下只能計算范圍，故已定系統(tǒng)誤差可以進(jìn)行修正，而未定系統(tǒng)誤差則不能進(jìn)行修正。

2.1.1 MWD傳感器誤差權(quán)重函數(shù)

MWD共有12個傳感器誤差源，每個誤差源對應(yīng)一個權(quán)重函數(shù)向量，可通過對井斜角、方位角求誤差源偏導(dǎo)數(shù)獲得。以X軸重力加速度計的偏移誤差和X軸重力加速度計的刻度誤差為例，通過偏微分求導(dǎo)，從而得到修正后的權(quán)重函數(shù)ABIX：

(3)

2.1.2 儀器不居中誤差權(quán)重函數(shù)

儀器不居中誤差可以細(xì)致的歸為底部鉆具組合垂曲誤差項和徑向不對稱誤差項兩類。底部鉆具組合垂曲只影響井斜角的測量，其大小與垂直于鉆具組合的重力分量成正比，底部鉆具組合垂曲誤差項的權(quán)重函數(shù)：

(4)

徑向不居中誤差項認(rèn)為在X,Y軸上存在相互獨立且固定的徑向不對稱誤差，它們在高邊和側(cè)向方向上的分量分別導(dǎo)致井斜角誤差和方位角誤差為：

(5)

(6)

2.1.3 鉆具磁干擾誤差權(quán)重函數(shù)

磁方位角誤差項權(quán)重函數(shù)為:

(7)

軸向磁干擾誤差項權(quán)重函數(shù)為：

(8)

2.1.4 深度刻度系數(shù)誤差權(quán)重函數(shù)

刻度誤差權(quán)重函數(shù)為：

(9)

2.2 隨機(jī)誤差及其對應(yīng)權(quán)重函數(shù)

隨機(jī)誤差通過多次測量估算平均值，這個平均值又稱為隨機(jī)誤差限。參考誤差權(quán)重函數(shù)為：

(10)

2.3 全局誤差及其權(quán)重函數(shù)

磁偏角固定數(shù)值誤差項權(quán)重函數(shù)：

(11)

磁偏角水平磁場強(qiáng)度(變值)誤差項權(quán)重函數(shù)為：

(12)

磁傾角誤差項權(quán)重函數(shù)為：

(13)

總磁場強(qiáng)度誤差項權(quán)重函數(shù)為：

(14)

拉伸類型誤差項權(quán)重函數(shù)為：

(15)

二、井眼軌跡不確定性誤差橢球的計算

在每個誤差項權(quán)重函數(shù)已經(jīng)被建立的基礎(chǔ)上，就能夠用向量的形式將每個誤差項所產(chǎn)生的誤差表示出來，然后根據(jù)概率統(tǒng)計相關(guān)知識建立協(xié)方差矩陣并得出誤差橢球來表達(dá)井眼位置的不確定性。

將隨機(jī)誤差、系統(tǒng)誤差、全局誤差三種誤差類型相加，就能算出井眼軌跡上任何測點的總的位置不確定性的協(xié)方差矩陣為：

(16)

三、誤差橢球的三維可視化

用Python語言對ISCWSA誤差模型進(jìn)行描述。根據(jù)以上推導(dǎo)公式，把ISCWSA誤差模型轉(zhuǎn)化為算法并實現(xiàn)在計算機(jī)上，井眼軌跡計算采用最小曲率法。整個程序采用模塊化面向?qū)ο缶幊棠Ｊ?，程序的靈活性和可重復(fù)利用性強(qiáng)，如圖1。

圖1 ISCWSA求井眼軌跡不確定性算法流程

由于誤差橢球半軸相對于井深來說太小，為了能直觀的在三維空間下顯示誤差橢球，在誤差橢球三維可視化的過程中一般需要把橢球半軸長度按照比例放大，放大比例可以根據(jù)實際情況而定。針對某一個點的和針對連續(xù)多個點的誤差橢球可視化圖形如圖2和圖3所示。

圖2 單點誤差橢球

圖3 連續(xù)誤差橢球

四、應(yīng)用實例

目前，MWD隨鉆測量工具在鉆井中的普遍使用，基于ISCWSA模型的井眼軌跡不確定性誤差橢球計算和三維可視化對現(xiàn)場鉆井施工有著很好的指導(dǎo)作用。例如：為了防止正鉆井與鄰井相碰，鉆井人員可以通過觀察鉆頭所在位置的誤差范圍是否與鄰井交碰來調(diào)整鉆進(jìn)方向；在打救援井時為了盡量讓其與事故井相碰，鉆井人員可以通過井眼軌跡控制來使誤差橢球中心靠近事故井，從而最大概率的使救援井與事故井相碰。

為驗證實現(xiàn)方法的可行性，選擇渤海某平臺上的兩口相鄰的井B1井和B2井來繪制誤差橢球，用來驗證三維誤差橢球可視化在防碰預(yù)測和井眼軌跡控制領(lǐng)域的指導(dǎo)意義，兩口井的井口間距為25 m，各種計算參數(shù)選取文獻(xiàn)[4]中提供的參數(shù)，井眼軌跡分段數(shù)據(jù)及相應(yīng)誤差分析結(jié)果分別如表1和表2所示。

表1 B1井井眼軌跡分段數(shù)據(jù)及誤差分析結(jié)果

表2 B2井井眼軌跡分段數(shù)據(jù)及誤差分析結(jié)果

圖4 兩口井防碰預(yù)警或連通示例圖

如圖4所示，將B2井作為參考井，B1井在鉆井過程中，在約2 600 m處誤差橢球即將與B2井相碰，若B1井繼續(xù)按照原來的方向鉆進(jìn)，則與B2井碰撞的概率增大，為防止兩井相碰，應(yīng)及時調(diào)整B1井井斜或方位來調(diào)整其鉆進(jìn)方向。

五、結(jié)束語

該研究用Python編程語言實現(xiàn)了井眼軌跡的三維可視化和基于ISCWSA法誤差橢球的計算，并展示了隨井深變化的參數(shù)信息，取得了良好的效果。在此基礎(chǔ)上，利用NumPy, SciPy, Matplotlib,和Axes 3D等函數(shù)庫實現(xiàn)了對三維空間井眼軌跡誤差橢球的可視化，可以直觀的顯示井眼軌跡上各個測點的誤差范圍，對叢式井井眼防碰或井眼軌跡控制領(lǐng)域起到了較好的參考和指導(dǎo)意義。