使用最優(yōu)化理論推導(dǎo)全力扭方位計算公式

白鳳俠

（中國石油遼河油田公司鉆采工藝研究院，遼寧 盤錦124010）

使用最優(yōu)化理論推導(dǎo)全力扭方位計算公式

白鳳俠

（中國石油遼河油田公司鉆采工藝研究院，遼寧 盤錦124010）

全力扭方位的數(shù)學(xué)本質(zhì)是在給定的約束條件下求裝置角，使得方位變化達到最大。文獻中已有的幾種求解方法在求解過程中沒有考慮約束條件，這使得有關(guān)計算公式缺乏嚴(yán)格的數(shù)學(xué)基礎(chǔ)。使用有約束極值問題的標(biāo)準(zhǔn)求解方法——拉格朗日乘數(shù)法對全力扭方位問題重新進行了數(shù)學(xué)表述、并推導(dǎo)出了最優(yōu)解的計算公式。

鉆井工程；全力扭方位；裝置角；拉格朗日法

扭方位計算是鉆井施工現(xiàn)場常用的一種工程計算［1－3］，在劉修善［4，5］、韓志 勇［6］的著作 中有非常 詳細的討論。全力扭方位是一種特殊形式的扭方位，它要求在給定條件下使得方位變化達到最大。全力扭方位計算實際上是求解一個極值問題，在文獻中［4－6］提供了幾種不同的求解方法，但是細讀之下感覺在數(shù)學(xué)推導(dǎo)上不是很嚴(yán)格，都沒有考慮極值問題中的約束條件；從數(shù)學(xué)觀點上來看，有約束極值問題與無約束極值問題是有本質(zhì)區(qū)別的，盡管在某些特定條件下、有約束極值問題與無約束極值問題的解是相同的。本文遵循有約束極值問題的標(biāo)準(zhǔn)求解方法——拉格朗日乘數(shù)法，給出了全力扭方位這一有約束極值問題的標(biāo)準(zhǔn)求解過程，可以看做是對文獻中各種方法的一種補充和完善。

約定：除非特別指明，本文公式中的角度參數(shù)的物理單位為rad，長度參數(shù)的物理單位為m，井眼曲率的物理單位為m－1。

1 扭方位計算的基本公式

假設(shè)當(dāng)前井底的井斜角為α1，方位角為φ1，初始裝置角為ω，扭方位井段終點的井斜角為α，方位角為φ，井段為圓弧井段，彎曲角為ε，則有［4］：

式（4）中：Δφ＝φ－φ1。

另外，當(dāng)井段長度ΔL和井眼曲率K給定時，彎曲角可以由下式計算：

在四個獨立參數(shù)α、φ、ε、ω中，只要給定其中的兩個參數(shù)，就能從方程組（3）～（4）求出另外的兩個參數(shù)。

2 全力扭方位的兩種算法

全力扭方位是一種特殊形式的扭方位，就是在給定井段長度和井眼曲率（即是相當(dāng)于彎曲角ε為已知參數(shù)）的情況下，求初始裝置角，使得井段末端方位變化最大。

因為既可以是方位增加、也可以是方位減少，故全力扭方位就是要求｜Δφ｜為最大的一個極值問題。目前，求這個極值問題有兩種方法，盡管數(shù)學(xué)求解過程有些不同，但最終結(jié)果相差不多。

2．1 劉修善方法

劉修善［4，5］給出：

式（6）實際上是式（2）除以式（1）再化簡的結(jié)果。

文獻［4，5］指出：根據(jù)微分學(xué)原理，裝置角ω應(yīng)滿足下面的條件：

由于上式比較復(fù)雜，可以使用下面的等價方程：

然后將式（6）代入式（8），經(jīng)過復(fù)雜的計算之后得：

2．2 韓志勇方法

韓志勇［6］給出：

將式（10）代入下式：

經(jīng)過復(fù)雜的計算之后得出的裝置角解也是式（9）。

2．3 注解

上述兩種方法中都有三個地方需要加以注意：

第一，式（7）與式（8）或式（11）的等價性并沒有經(jīng)過嚴(yán)格的數(shù)學(xué)證明。

第二，式（7）只是這個極值問題的必要條件、而不是充分條件，所以所得到的解（9）是不是極值點需要加以驗證（也有可能是駐點，即既不是極大值點、也不是極小值點）。

第三，由于式（6）或式（10）本身是有理三角函數(shù)，推導(dǎo)出式（8）或式（11）左邊的表達式也是一個比較復(fù)雜的求導(dǎo)和化簡過程。

3 全力扭方位的新算法

劉修善算法和韓志勇算法都是直接將裝置角ε作為優(yōu)化變量來求解一個間接的優(yōu)化問題，由于優(yōu)化的目標(biāo)函數(shù)比較復(fù)雜，數(shù)學(xué)推導(dǎo)過程也比較復(fù)雜。

從空間圓弧井段的幾何學(xué)意義上來看，當(dāng)井段彎曲角給定之后，井段的形狀是確定的，當(dāng)井段末端的井斜角和方位角變化時，所有這樣的井段形成一個以井段起點處井眼方向線為對稱軸的旋轉(zhuǎn)曲面——曲圓錐，井段的末端在曲圓錐的底圓上。全力扭方位的目標(biāo)就是在曲圓錐的底圓上選擇某個點或多個點，使得｜Δφ｜達到最大。

如果這樣來認(rèn)識全力扭方位問題，可以形成下面的數(shù)學(xué)問題（二元約束極值問題）：求α、Δφ，使得在滿足約束方程（4）的條件下，求目標(biāo)函數(shù)｜Δφ｜的極值。

由于目標(biāo)函數(shù)｜Δφ｜的導(dǎo)數(shù)具有不連續(xù)點，定義一個等價的、連續(xù)可導(dǎo)的新目標(biāo)函數(shù)：f（α，Δφ）＝Δφ）2。

如果該極值問題有解，再根據(jù)式（1）～（3）的任一式反求裝置角ω。

3．1 拉格朗日乘數(shù)法求解

拉格朗日乘數(shù)法［7］是求解等式約束極值問題的一種有效的方法。定義拉格朗日函數(shù)如下：

則有：

在由式（16）和式（18）得：

這個方程只包含未知數(shù)α，化簡后求得：

在扭方位時，一般來說井段初始井斜角α1＜π／2，故：

將式（20）代入式（16）得：

故：

上面的結(jié)果表明，全力扭方位問題有兩個解：一個解是方位角增加方向上的，另一個解是方位角減少方向上的，但是這兩個解的井斜角和方位角增量絕對值是相同的，可以看成是幾何意義上的對稱解。還應(yīng)該注意到，上述推導(dǎo)過程完全沒有涉及裝置角的概念，這說明全力扭方位問題的數(shù)學(xué)本身是與裝置角無關(guān)的。

3．2 裝置角的計算

從式（3），得：

再將式（20）代入式（24），得：

這個結(jié)果與劉修善方法和韓志勇方法的結(jié)果相同，見式（9）。

由于三角函數(shù)的周期性，從式（25）反求裝置角時有兩個解，令：

是［0，π］上的解，則另一個解為：ω2＝2π－ω1。根據(jù)韓志勇［6］的分析，ω1對應(yīng)于增方位的解Δφ，ω2對應(yīng)于減方位的解－Δφ。

4 結(jié)論

1）在全力扭方位問題的求解方法中，都是以裝置角為參數(shù)求解一個一元極值問題。但是在井段彎曲角固定的情況下，方位角的變化值是裝置角和井段末端井斜角的二元函數(shù)，故以一元極值問題來求解全力扭方位問題不具有嚴(yán)謹(jǐn)性；而且全力扭方位問題的各個參數(shù)之間具有一定的約束關(guān)系，但是原來的解法沒有考慮這些約束條件。

2）本文提供了一個新的解法。將全力扭方位問題描述為二元約束極值問題，并且使用標(biāo)準(zhǔn)的拉格朗日乘子法求出了該問題的解。與原來解法相比較，本文解法在數(shù)學(xué)上是嚴(yán)謹(jǐn)?shù)?，所得到的解具有簡單明了的幾何直觀性。

［1］向軍文．定向?qū)舆B通井軌跡設(shè)計［J］．探礦工程（巖土鉆掘工程），2011，38（5）：11－14．

［2］李景東．哈達門溝金礦區(qū)復(fù)雜地層中深孔鉆進實踐［J］．探礦工程（巖土鉆掘工程），2011，37（6）：20－23．

［3］隆東，張新剛，岳剛，等．H024U井施工工藝及精確中靶技術(shù)措施［J］．探礦工程（巖土鉆掘工程），2011，38（3）：5－12．

［4］劉修善，王珊，賈仲宣，等．井眼軌道設(shè)計理論與描述方法［M］．哈爾濱：黑龍江科學(xué)技術(shù)出版社，1993：257－279．

［5］周大千，劉修善，齊林，等．井眼軌道實用理論基礎(chǔ)［M］．北京：石油工業(yè)出版社，1993：59－72．

［6］韓志勇．定向鉆井設(shè)計與計算［M］．東營：中國石油大學(xué)出版社，2007：285－300．

［7］趙鳳治，尉繼英．約束最優(yōu)化計算方法［M］．北京：科學(xué)出版社，1991．

On Deducing Calculation Formula for Full－Twist Azimuth Based on Optimization Theory

BAI Feng－xia

（Drilling ＆Production Technology Research Institute of Liaohe Oilfield Company，CNPC，Panjin，Liaoning，124010，China）

Mathematics essence of full－twist azimuth is calculating setting angle with given constrain condition to obtain maximum azimuth change．Some of existed calculation methods in the literature neglect constrain condition in the course of proof，thus making the related calculation formula lack of strict foundation of mathematics．The author has given anew mathematics formulation of full－twist azimuth by means of standard method to calculate constrained extreme－value problems（Lagrange Multiplier method），and deduced the calculation formula for optimal solution．

Drilling Engineering；Full－Twist Azimuth；Setting Angle；Lagrangian Method

TE22

A

1009—301X（2012）04—0016—03

2012－03－08

白鳳俠（1968－），女，助理工程師，1988年畢業(yè)于遼河石油學(xué)校采油工程專業(yè)，長期從事采油工藝技術(shù)研發(fā)等工作。

［責(zé)任編輯 郭華玉］