陳庚緒，樊 銳，李紅星，王新剛，白雪龍，王越支

(1長(zhǎng)江大學(xué)石油工程學(xué)院 2中國(guó)石油集團(tuán)渤海鉆探工程有限公司第二鉆井工程分公司 3中海油服油技塘沽作業(yè)公司)

井眼軌跡測(cè)斜計(jì)算方法很多，常見(jiàn)的有平均角法、校正平均角法、平衡正切法、最小曲率法、曲率半徑法和弦步法等，不論哪種方法都是將測(cè)點(diǎn)之間的測(cè)段假定為某種特定的線段，但是測(cè)段間真實(shí)客觀存在的井眼軌跡形狀是無(wú)法預(yù)知的，故測(cè)斜計(jì)算得到的井眼軌跡與真實(shí)客觀存在的井眼軌跡必然不相同。針對(duì)井眼軌跡測(cè)斜計(jì)算誤差問(wèn)題，陳煒卿、管志川等學(xué)者通過(guò)構(gòu)造理論參考曲線，對(duì)不同測(cè)斜計(jì)算方法計(jì)算的平均井眼位置進(jìn)行了大樣本空間對(duì)比分析，得出所有測(cè)斜計(jì)算方法表現(xiàn)出來(lái)的誤差特性都很相似的結(jié)論，并用統(tǒng)一的解析式表示了所有的測(cè)斜計(jì)算方法所產(chǎn)生的誤差概率密度函數(shù)，但是誤差和不確定度是兩個(gè)不同的概念，他們沒(méi)有從不確定度的角度去進(jìn)一步分析[1-2]。

最小曲率法是建立在如下假設(shè)基礎(chǔ)上的：①最小曲率法是基于分段近似，假設(shè)測(cè)點(diǎn)之間的曲率為定值;②最小曲率法是一階連續(xù)的，但它在二階導(dǎo)數(shù)和三階導(dǎo)數(shù)中表現(xiàn)出不連續(xù)性。對(duì)于給定的一組測(cè)點(diǎn)，實(shí)際井眼軌跡可能是曲線段和直線段的混合曲線，而不是連續(xù)的曲線，就算是先進(jìn)的旋轉(zhuǎn)導(dǎo)向工具(RSS)也不能鉆出連續(xù)曲線。在眾多井眼軌跡計(jì)算方法中，樣條曲線法基于線性代數(shù)方法，對(duì)于一階，二階和三階都是連續(xù)的。樣條曲線法能有效使得給定測(cè)段的曲率和曲率變化率連續(xù)(二階和三階)[3-6]。

圖1顯示了用最小曲率法計(jì)算兩測(cè)點(diǎn)之間的三種可能的井眼軌跡，由與測(cè)斜數(shù)據(jù)相匹配的圓弧表示。中間的圓弧為最小曲率法計(jì)算的井眼軌跡，因測(cè)點(diǎn)間可能是復(fù)合鉆進(jìn)和滑動(dòng)鉆進(jìn)交叉進(jìn)行且兩者之比無(wú)法確定，故另外兩條井眼軌跡不同于用最小曲率法計(jì)算的井眼軌跡，從而導(dǎo)致井眼曲率(二階)不連續(xù)。這種不連續(xù)性導(dǎo)致由每個(gè)圓弧的半徑表示的恒定曲率半徑產(chǎn)生了差異。最小曲率法與樣條曲線法的區(qū)別是一個(gè)只能保證函數(shù)值連續(xù)，另一個(gè)除了保證函數(shù)值連續(xù)還保證它的導(dǎo)數(shù)也連續(xù)。

圖1 兩測(cè)點(diǎn)間滑動(dòng)距離不同對(duì)測(cè)段的影響

一、井眼軌跡計(jì)算方法的蒙特卡洛模擬

蒙特卡洛模擬是基于概率統(tǒng)計(jì)的一種算法[7]，若需模擬井眼軌跡誤差不確定度，需要先建立概率統(tǒng)計(jì)模型，其次確定模型中誤差變量數(shù)據(jù)的分布，然后根據(jù)誤差變量的精度要求確定模擬抽樣次數(shù)，最后在隨機(jī)抽樣產(chǎn)生結(jié)果的基礎(chǔ)上實(shí)現(xiàn)概率統(tǒng)計(jì)分析，隨機(jī)抽樣的次數(shù)越多，井眼軌跡位置越逼近真實(shí)值。在本研究中，井眼軌跡測(cè)斜計(jì)算方法即概率統(tǒng)計(jì)計(jì)算模型，模型中的誤差變量為井斜角誤差、方位角誤差和井深誤差。本研究主要圍繞最小曲率法、曲率半徑法和樣條曲線法展開(kāi)蒙特卡洛模擬分析。

1.井眼軌跡測(cè)斜計(jì)算模型

最小曲率法和曲率半徑法是常見(jiàn)的井眼軌跡測(cè)斜計(jì)算方法，這兩種計(jì)算方法不再累述。樣條曲線法假設(shè)井眼軌跡是一個(gè)由三次多項(xiàng)式組成的函數(shù)，且函數(shù)一階導(dǎo)、二階導(dǎo)在所有測(cè)點(diǎn)都是連續(xù)的，其測(cè)斜計(jì)算模型為[7]：

(1)

(2)

(3)

式(1)～式(3)中：ΔV、ΔN和ΔE—分別為垂深、北坐標(biāo)和東坐標(biāo)增量，m；I和A—井斜角和方位角，°；ΔI和ΔA—井斜角和方位角增量，°；ΔL—測(cè)段長(zhǎng)，m。

2.井眼軌跡誤差變量分布及精度的確定

1999年威廉姆森對(duì)于各種測(cè)量工具和測(cè)量類型，在特定井眼條件下給定誤差的加權(quán)函數(shù)，測(cè)量誤差得到了很好的理解和建模。加權(quán)函數(shù)考慮了誤差源對(duì)方位角和井斜角的影響，測(cè)量工具通常根據(jù)一個(gè)標(biāo)準(zhǔn)偏差來(lái)反應(yīng)測(cè)量精度，加權(quán)函數(shù)對(duì)每個(gè)工具都是唯一的，對(duì)于這一分析，測(cè)得的深度、井斜角和方位角的在一西格瑪精度下的標(biāo)準(zhǔn)偏差分別為0.3 m、0.25°和0.36°(即σH=0.3，σI=0.25，σA=0.36)，整口井的標(biāo)準(zhǔn)偏差為常量。圖2為標(biāo)準(zhǔn)正太分布圖。

圖2 標(biāo)準(zhǔn)正態(tài)分布

井眼軌跡空間位置用隨機(jī)正態(tài)分布得到的測(cè)深、井斜角和方位角進(jìn)行模擬即可得井眼軌跡不確定性范圍散點(diǎn)圖。

(4)

(5)

式中：I—井斜角，°；A—方位角，°。

3.三種井眼軌跡測(cè)斜計(jì)算方法的蒙特卡洛模擬結(jié)果

本研究是比較不同井眼軌跡測(cè)斜計(jì)算方法在假設(shè)分布條件下所產(chǎn)生的井眼位置不確定性，井眼軌跡測(cè)斜計(jì)算誤差指的就是由某種測(cè)斜計(jì)算方法得到的井眼位置與真實(shí)井眼位置在相同測(cè)深處的井眼位置之間的距離。因此，首先需要定義真實(shí)井眼位置或者一個(gè)可以作為真實(shí)值的測(cè)量值。自動(dòng)尋北高分辨率連續(xù)速率陀螺儀不受磁場(chǎng)干擾，而最小曲率法是目前井眼軌跡計(jì)算的行業(yè)標(biāo)準(zhǔn)，因此，根據(jù)陀螺測(cè)量數(shù)據(jù)用最小曲率法計(jì)算的井眼位置被認(rèn)為是比較真實(shí)的井眼軌跡，可將其作為擬定真值[8]。

由于測(cè)點(diǎn)之間為曲線的假設(shè)相對(duì)來(lái)說(shuō)更加符合井眼軌跡的實(shí)際情況，故本文選取中海油遼東灣某井的隨鉆測(cè)量和陀螺測(cè)量數(shù)據(jù)，就最小曲率法、曲率半徑法和樣條曲線法三種井眼軌跡計(jì)算方法展開(kāi)了蒙特卡洛模擬分析研究。圖3為三種方法的對(duì)比圖。

圖3 三種方法井眼軌跡位置的蒙特卡洛模擬對(duì)比

二、三種井眼軌跡計(jì)算方法蒙特卡洛模擬統(tǒng)計(jì)對(duì)比分析

在上述計(jì)算實(shí)例中，樣條曲線法較最小曲率法和曲率半徑法在井眼軌跡不確定誤差橢圓/球方面有了微小的變化。使用樣條曲線法計(jì)算得到的井眼軌跡不確定誤差橢圓/球的大小略大于最小曲率法和螺旋圓柱法，但是樣條曲線法計(jì)算得到的不確定性誤差橢圓/球的中心位置卻更集中于擬定真值附近。若忽略置信區(qū)間的影響，則樣條曲線法計(jì)算得到的井眼軌跡位置的不確定度因接近擬定真值而較最小曲率法和曲率半徑法分別降低了約78%和91%。具體結(jié)果見(jiàn)表1和表2。

表1 蒙特卡洛模擬在井底測(cè)點(diǎn)的統(tǒng)計(jì)分析

表2 蒙特卡洛模擬在井底測(cè)點(diǎn)的不確定度分析

與1 m間隔的高分辨率測(cè)量得到的擬定真值相比，無(wú)論采用何種井眼軌跡計(jì)算方法得到的不確定度誤差橢球都會(huì)在橢球中心處發(fā)生偏移。原因是系統(tǒng)誤差未經(jīng)過(guò)充分修正且測(cè)點(diǎn)之間的高狗腿度很難精確的去確定。圖4是井眼坐標(biāo)系下的井眼軌跡誤差橢球各半軸的概率密度圖?？梢钥闯觯?0 m測(cè)量間隔計(jì)算的井眼軌跡誤差橢球中心之間的差異顯示了置信區(qū)間和概率密度之間的區(qū)別。

圖4 井眼坐標(biāo)系下井底測(cè)點(diǎn)誤差橢球各半軸的概率密度圖

三、結(jié)論與建議

(1)井眼軌跡測(cè)斜計(jì)算方法的蒙特卡羅模擬表明，各種井眼軌跡測(cè)斜計(jì)算方法在誤差橢球不確定度表征上的差異不大是可以忽略不計(jì)的。忽略置信區(qū)間的影響，針對(duì)30 m的標(biāo)準(zhǔn)測(cè)量間隔距離，樣條曲線法計(jì)算的井眼軌跡位置更準(zhǔn)確，但與最小曲率法和曲率半徑法一樣受到系統(tǒng)誤差的影響。

(2)計(jì)算井眼軌跡位置不確定度的大小不僅受各種誤差因素(標(biāo)準(zhǔn)偏差)的影響，還取決于誤差模型(本文為井眼軌跡測(cè)斜計(jì)算方法)，蒙特卡洛計(jì)算機(jī)模擬簡(jiǎn)化了繁瑣的計(jì)算過(guò)程，使得井眼軌跡測(cè)斜計(jì)算方法對(duì)比起來(lái)更加簡(jiǎn)單明了。

(3)針對(duì)最小曲率法計(jì)算的井眼軌跡，可以采取相應(yīng)措施來(lái)大大提高測(cè)量的準(zhǔn)確性，比如在滑動(dòng)鉆進(jìn)之間的中點(diǎn)處進(jìn)行測(cè)量或在已知的狗腿角較大的井段加密測(cè)量，這都將對(duì)井眼軌跡位置和井眼位置產(chǎn)生顯著影響。