趙明君

(中國(guó)石油長(zhǎng)城鉆探工程公司工程技術(shù)部遼河分部,遼寧盤錦 124010）

井眼軌道的軟著陸設(shè)計(jì)模型的改進(jìn)解法

趙明君

(中國(guó)石油長(zhǎng)城鉆探工程公司工程技術(shù)部遼河分部,遼寧盤錦 124010）

井眼軌道的軟著陸設(shè)計(jì)模型的求解可以歸結(jié)為一個(gè)七元非線性方程組的求解問(wèn)題。前人給出了數(shù)值迭代求解算法,然而并沒(méi)有證明該迭代算法的收斂性,并且該算法是否收斂嚴(yán)重依賴于用戶給出的迭代初始值。通過(guò)一系列的消元、化簡(jiǎn)的數(shù)學(xué)技巧,將七元非線性方程組化簡(jiǎn)為一元多項(xiàng)式方程,并在此基礎(chǔ)上給出了軟著陸設(shè)計(jì)模型的一個(gè)新算法。理論分析和實(shí)際算例表明,新算法的主要計(jì)算工作量是求多項(xiàng)式方程的非負(fù)實(shí)數(shù)根,其他未知數(shù)與實(shí)數(shù)根是簡(jiǎn)單的函數(shù)關(guān)系,計(jì)算量很小。新算法克服了迭代算法的初值依賴性以及迭代過(guò)程可能發(fā)散等缺陷,并且在設(shè)計(jì)模型有多個(gè)解的情況下,可以同時(shí)求出這些解。

軟著陸;井眼軌道;鉆井設(shè)計(jì);多項(xiàng)式

劉修善等[1]將井眼軌道進(jìn)入水平井靶體的過(guò)程類比于飛機(jī)降落著陸過(guò)程,提出了井眼軌道的軟著陸設(shè)計(jì)模型,其宗旨是對(duì)入靶方向進(jìn)行限制。在水平井設(shè)計(jì)中,入靶方向需與靶體軸線方向相一致(簡(jiǎn)單的長(zhǎng)方體靶的情況下）。在進(jìn)行修正軌道設(shè)計(jì)時(shí),為了給后續(xù)的施工創(chuàng)造有利條件,如果其靶點(diǎn)不是最終目標(biāo)點(diǎn),而是原設(shè)計(jì)軌道上的某個(gè)中間點(diǎn),也必須限定合理的入靶方向。文獻(xiàn)[1]使用“圓弧段-直線段-圓弧段”實(shí)現(xiàn)了上述軟著陸模型,并提出了一種迭代算法進(jìn)行設(shè)計(jì)模型的數(shù)值求解。然而,文獻(xiàn)[1]并沒(méi)有討論該迭代算法的收斂性,算法是否收斂、收斂速度怎樣、收斂約束條件是什么等問(wèn)題都沒(méi)有解決。

本文使用“消元-化簡(jiǎn)-多項(xiàng)式方程求根”的方法給出了軟著陸設(shè)計(jì)模型求解的一個(gè)全新算法。

約定:除非特別指明,文中具有長(zhǎng)度量綱的變量的物理單位均為m,角度變量的物理單位均為rad。

1 數(shù)學(xué)模型

沿用文獻(xiàn)[1]中的數(shù)學(xué)符號(hào)。井眼軌道有圓弧段AB、直線段BC、圓弧段CT構(gòu)成(見(jiàn)圖1）。

圖1 軟著陸軌道的數(shù)學(xué)模型

過(guò)始點(diǎn)A和末點(diǎn)T分別作圓弧的切線,交穩(wěn)斜段的延長(zhǎng)線與D點(diǎn)和E點(diǎn)。若線段AD、DE和ET的長(zhǎng)度分別用u1、uW、u2來(lái)表示,則井眼軌道設(shè)計(jì)問(wèn)題滿足下面的幾個(gè)方程:

式中:αA、φA——分別為始點(diǎn)A的井斜角和方位角; αW、φW——分別為穩(wěn)斜段的井斜角和方位角;αT、φT——分別為末點(diǎn)T的井斜角和方位角;ΔXA,T、ΔYA,T、ΔZA,T——分別為始末點(diǎn)間的北坐標(biāo)增量、東坐標(biāo)增量和垂深增量。

另外,兩個(gè)圓弧段所對(duì)應(yīng)的圓心角θ1和θ2分別滿足:

式中:R1、R2——分別為第一個(gè)圓弧段和第二個(gè)圓弧段的曲率半徑。

在方程組(1）～(7）中,未知數(shù)為u1、uW、u2、αW、φW、θ1和θ2,共7個(gè),與方程個(gè)數(shù)相等,故方程組為定解問(wèn)題。

2 劉修善算法

然而,迭代法的收斂是有限定條件的,并且與迭代初始值有很大關(guān)系[2]。劉修善算法是否收斂沒(méi)有理論上的保證,而且如何選擇迭代初始值也沒(méi)有具體的說(shuō)明。

3 消元

解多元方程組的一般策略是消元和化簡(jiǎn),將未知數(shù)的個(gè)數(shù)降到最少,將復(fù)雜的方程化簡(jiǎn)到已知可解的方程[3]。本文遵循這一策略對(duì)方程組(1）～(7）進(jìn)行求解。記:

顯然,式(41）和式(45）構(gòu)成一個(gè)以u(píng)1和u2為未知數(shù)的二元線性代數(shù)方程組,方程組的系數(shù)為未知數(shù)u的多項(xiàng)式函數(shù)。使用克萊默法則[4]求得該線性代數(shù)方程組的形式解如下:

顯然,F(u）是一個(gè)以u(píng)為未知數(shù)的多項(xiàng)式方程。

5 新算法

根據(jù)以上的理論準(zhǔn)備,給出方程組(1）～(7）的一個(gè)新算法。

(1）使用實(shí)根分離算法[5]求出多項(xiàng)式方程(55）的全部非負(fù)實(shí)數(shù)根,記為νi,i=1,Λ,N,其中N為非負(fù)實(shí)數(shù)根個(gè)數(shù)。

從新算法的計(jì)算過(guò)程來(lái)看,主要計(jì)算工作量集中在求多項(xiàng)式方程(55）的實(shí)數(shù)根上,其余計(jì)算步驟僅僅是簡(jiǎn)單的公式計(jì)算。新算法從根本上解決了常規(guī)迭代算法固有的初值依賴性、迭代可能發(fā)散等缺陷;另外,如果井眼軌道設(shè)計(jì)問(wèn)題有多個(gè)解的話,新算法還能同時(shí)求出這些解,這個(gè)特點(diǎn)是其他算法不具備的。

在進(jìn)行數(shù)值計(jì)算時(shí),為了避免數(shù)字反復(fù)相乘使多項(xiàng)式系數(shù)過(guò)大、從而增大計(jì)算誤差,令:

顯然,方程f(x）=0等價(jià)于方程F(u）=0,故在新算法中使用方程f(x）=0來(lái)代替方程(55）。

6 算例

在某多目標(biāo)井設(shè)計(jì)中,要求上靶點(diǎn)的井斜角αA=75°、方位角φA=310°,下靶點(diǎn)的井斜角αB=91°、方位角φB=340°。如果這兩個(gè)靶點(diǎn)間的垂深增量ΔZA,T=30 m、南北坐標(biāo)增量ΔXA,T=180 m、東西坐標(biāo)增量ΔYA,T=-100 m,試設(shè)計(jì)該靶段的井眼軌道。

根據(jù)本文新算法,首先可以算出:

求出上述解之后,再根據(jù)圓弧軌道計(jì)算方法計(jì)算軌道其他參數(shù)[6～8]。

7 結(jié)論

(1）井眼軌道的軟著陸設(shè)計(jì)模型歸結(jié)為求解一個(gè)7元非線性方程組,通過(guò)消元和化簡(jiǎn),將7元非線性方程組的求解問(wèn)題歸結(jié)為等價(jià)的1元多項(xiàng)式方程的求解問(wèn)題,并在此基礎(chǔ)上,給出了軟著陸設(shè)計(jì)模型的一個(gè)新算法。新算法克服了常規(guī)迭代算法對(duì)初值的依賴性和無(wú)法證明迭代過(guò)程收斂等固有缺陷。

(2）實(shí)例計(jì)算表明,新算法具有非?？斓挠?jì)算速度,并且能夠正確判斷軟著陸設(shè)計(jì)模型是否有解。

[1] 劉修善,何樹(shù)山.井眼軌道的軟著陸設(shè)計(jì)模型及其應(yīng)用[J].天然氣工業(yè),2002,22(2）:43-45.

[2] 李慶揚(yáng),莫孜中,祁力群.非線性方程組的數(shù)值解法[M].北京:科學(xué)出版社,1987.

[3] 吳文俊.吳文俊論數(shù)學(xué)機(jī)械化[M].濟(jì)南:山東教育出版社, 1996.

[4] 編寫(xiě)組.數(shù)學(xué)手冊(cè)[M].北京:人民教育出版社,1979.103-111.

[5] 陸征一,何碧,羅勇.多項(xiàng)式系統(tǒng)的實(shí)根分離算法及其應(yīng)用[M].北京:科學(xué)出版社,2004.

[6] 魯港,王剛,孫忠國(guó),等.定向井鉆井中空間圓弧軌道計(jì)算的兩個(gè)問(wèn)題[J].石油地質(zhì)與工程,2006,20(6）:53-55.

[7] 魯港,李曉光,單俊峰,等.平均井眼曲率的計(jì)算[J].鉆采工藝,2007,30(4）:149-150,160.

[8] 張積鎖,魯港.圓弧井段井斜變化率和方位變化率的計(jì)算[J].石油地質(zhì)與工程,2007,21(4）:68-70.

I mproved Solution to DesignM odel of SoftLanding in Borehole Trajectory

ZHAO M ing-jun(LiaoheBranch of Engi-neering&Technology Research Institute,GreatWallDrilling Corporation,PetroChina,Panjin Liaoning 124010,China）

The solution for design model of soft landing in borehole trajectory can come down to the solution with 7-element nonlinear equations.The algorithm convergence of known numerical iterative algorithm has not been proved,and the con-vergence of this iterative algorithm severely depends on initial iteration value from the user.7-element nonlinear equations are simplified to be one-variable polynomial equation by element elimination and s implification;and based onwhich,a new algorithm for design model of soft landingwas developed.Theoretical analysis and practical calculation examples show that the new algorithm ismainly for the calculation of nonnegative real rootsof polynomial equation,and calculation of functional relationship between other unknown number and real root is very less.The new algorithm overcomes the dependence on the initial value and probable iterative process divergence,and can get the solutions at the same timewhen there are several de-sign model solutions.

soft landing;borehole trajectory;drilling design;polynomial equation

TE243

:A

:1672-7428(2010）05-0010-04

2010-01-16;

2010-04-08

趙明君(1963-）,男(漢族）,遼寧本溪人,中國(guó)石油長(zhǎng)城鉆探工程公司工程技術(shù)部遼河分部工程師,石油工程專業(yè),從事石油鉆井領(lǐng)域技術(shù)服務(wù)和管理工作,遼寧省盤錦市興隆臺(tái)區(qū)。