使用最優化理論推導全力扭方位計算公式

白鳳俠

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

使用最優化理論推導全力扭方位計算公式

白鳳俠

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

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

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

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

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

1 扭方位計算的基本公式

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

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

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

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

2 全力扭方位的兩種算法

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

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

2．1 劉修善方法

劉修善［4，5］給出：

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

文獻［4，5］指出：根據微分學原理，裝置角ω應滿足下面的條件：

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

然后將式（6）代入式（8），經過復雜的計算之后得：

2．2 韓志勇方法

韓志勇［6］給出：

將式（10）代入下式：

經過復雜的計算之后得出的裝置角解也是式（9）。

2．3 注解

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

第一，式（7）與式（8）或式（11）的等價性并沒有經過嚴格的數學證明。

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

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

3 全力扭方位的新算法

劉修善算法和韓志勇算法都是直接將裝置角ε作為優化變量來求解一個間接的優化問題，由于優化的目標函數比較復雜，數學推導過程也比較復雜。

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

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

由于目標函數｜Δφ｜的導數具有不連續點，定義一個等價的、連續可導的新目標函數：f（α，Δφ）＝Δφ）2。

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

3．1 拉格朗日乘數法求解

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

則有：

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

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

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

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

故：

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

3．2 裝置角的計算

從式（3），得：

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

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

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

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

4 結論

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

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

［1］向軍文．定向對接連通井軌跡設計［J］．探礦工程（巖土鉆掘工程），2011，38（5）：11－14．

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

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

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

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

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

［7］趙鳳治，尉繼英．約束最優化計算方法［M］．北京：科學出版社，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年畢業于遼河石油學校采油工程專業，長期從事采油工藝技術研發等工作。

［責任編輯 郭華玉］