基于蒙特卡洛模擬的井眼軌跡測斜計(jì)算方法對比分析研究

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

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

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

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

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

圖1 兩測點(diǎn)間滑動距離不同對測段的影響

一、井眼軌跡計(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í)值。在本研究中，井眼軌跡測斜計(jì)算方法即概率統(tǒng)計(jì)計(jì)算模型，模型中的誤差變量為井斜角誤差、方位角誤差和井深誤差。本研究主要圍繞最小曲率法、曲率半徑法和樣條曲線法展開蒙特卡洛模擬分析。

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

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

(1)

(2)

(3)

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

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

1999年威廉姆森對于各種測量工具和測量類型，在特定井眼條件下給定誤差的加權(quán)函數(shù)，測量誤差得到了很好的理解和建模。加權(quán)函數(shù)考慮了誤差源對方位角和井斜角的影響，測量工具通常根據(jù)一個標(biāo)準(zhǔn)偏差來反應(yīng)測量精度，加權(quán)函數(shù)對每個工具都是唯一的，對于這一分析，測得的深度、井斜角和方位角的在一西格瑪精度下的標(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)分布得到的測深、井斜角和方位角進(jìn)行模擬即可得井眼軌跡不確定性范圍散點(diǎn)圖。

(4)

(5)

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

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

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

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

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

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

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

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

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

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

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

三、結(jié)論與建議

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

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

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