井眼轨道参数的插值计算
- 格式:doc
- 大小:440.50 KB
- 文档页数:17
井眼轨道的插值法
刘修善;王新清
【期刊名称】《石油钻采工艺》
【年(卷),期】1997(019)002
【摘要】井眼轨道值在石油工程中应用得十分广泛,而且不同于一般数学插值问题。
根据方式和井眼轨道的形成规律,提出了圆柱螺线,空间圆弧和自然参数三种典型的插值模型。
理论手算例表明,这些插值公式和对内插和外推都是适用的,可以满足各种井眼轨道插值计算的需要。
【总页数】5页(P11-14,25)
【作者】刘修善;王新清
【作者单位】大庆石油学院;大庆石油管理局
【正文语种】中文
【中图分类】TE243
【相关文献】
1.基于拉格朗日插值法的GPS卫星轨道位置拟合 [J], 刘金健
2.基于拉格朗日插值法的GPS卫星轨道位置拟合 [J], 刘金健;
3.石油工程科技名词的规范化使用——井眼轨迹与井眼轨道 [J], 陈会年
4.利用井眼偏移量模式预测井眼轨道 [J], 齐林;周大千;王新清;邸百英
5.实钻井眼轨迹坐标计算的样条插值法的讨论 [J], 陈铁铮;鲁港
因版权原因,仅展示原文概要,查看原文内容请购买。
井眼轨道参数的插值计算由于实钻井眼轨道的测点与钻柱单元体的划分可能并不一致,因此钻柱单元体边界点对应的井眼轨道参数必须靠插值计算获得。
插值结果的准确与否,对钻柱单元体的受力计算有着直接的影响。
因此,提高插值计算的精度具有重要意义。
由于测点是离散的,无法知道各测段内井眼轨道的实际形态,所以测段内某点几何参数的计算方法都是建立在一定假设的基础上的。
这些计算方法多数是将测段内的井眼轨道假设为直线、折线和曲线等,早期,由于计算机能力的限制,以平均角法和平衡正切法为代表直线或折线假设,因其计算简单快速,曾经被广泛应用,但随着钻井技术的发展,弯曲的井眼轨迹增多,如果仍采用直线或折线假设,则计算精度相对较低。
由于计算技术的高速发展,直线或折线假设,目前几乎淘汰,取而代之的是以圆柱螺线和空间圆弧曲线等为代表的曲线假设,大行其道。
在进行插值计算时,各插值点的坐标增量可以采用不同的计算方法,但坐标值的累加形式是相同的,即(X 为东向位移,Y 为北向位移, Z 为垂直向位移,S 为水平位移)⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧∆+=∆+=∆+=∆+=∆+=∆+=φφφa a αS S S Z Z Z YY Y X X X 121212121212所以,在以下的计算方法中将只给出坐标增量的计算式。
典型轨迹模型插值1、正切法:正切法又称下切点法,或下点切线法。
此法假定两相邻测点之间的孔段为一条直线,长度等于测距,该直线的井斜角和井斜方位角等于下测点的井斜角和井斜方位角,整个钻孔轨迹是直线与直线相连接的空间折线。
正切法井身轨迹计算图如图1所示,1、2 是孔身轨迹上相邻的两个测点,1′、2′是 1、2 两个测点的水平投影。
该测段的井斜角和井斜方位角等于下测点 2 的井斜角和井斜方位角。
对于切线法,上下两个相邻测点间各参数的计算公式如下:222222cos sin sin sin sin cos φαφαααL Y L X L S L Z ∆=∆∆=∆∆=∆∆=∆式中:Z ∆——测段上下测点间垂直深度的分量(增量)(以下同); L ∆——测段上下测点间沿钻孔轴线的距离(以下同);Y ∆∆X ——分别为测段上下测点间水平位移在 X 轴(西东方向)的分量(增量);水平位移在 Y 轴(南北方向)的分量(增量)(以下同);22 φα——分别为测段下测点的井斜角和井斜方位角。
石油地质与工程2011年3月PETROLEUM GEOLOGY AND ENGINEERING第25卷第2期文章编号:1673-8217(2011)02-0085-04定向井方位漂移轨道设计问题的数值算法刘开伦1,鲁港2(1.中国石油长城钻探工程公司钻井一公司,辽宁盘锦124010;2.中国石油辽河油田公司勘探开发研究院)摘要:方位漂移是实际钻井过程中普遍存在的现象,常规的轨道设计方法没有考虑方位漂移的影响,增加了钻井施工过程中的井眼轨迹控制的难度。
考虑方位漂移影响的轨道设计方法采用井段的自然曲线模型,能够较好地反映地层自然造斜和方位漂移规律。
所产生的约束方程组为三元非线性方程组,通常使用的数值迭代方法不易给定迭代初始值、迭代收敛性能差、收敛速度慢。
推导出了用稳斜井段井斜角来表示的初始方位角和稳斜段长的解析公式,将三元方程组的求解问题归结为以稳斜井段井斜角为未知数的一元非线性方程方程的求解问题,并给出了求解该一元非线性方程的数值迭代算法。
理论分析和大量实际算例表明,文中新算法具有收敛速度快、迭代初值相关性小、全局范围收敛等良好性质。
关键词:钻井设计;井眼轨道;方位漂移;非线性方程组;数值算法中图分类号:TE22文献标识码:A在钻井过程中,实钻井眼轨迹普遍存在着方位漂移现象[1-7]。
由于方位漂移,实钻轨迹与设计轨道产生较大偏差,给施工造成很大风险。
解决这个问题的一种方法是,在井眼轨道设计阶段就将方位漂移考虑在内,为此提出了一些新的设计方法[8-16]。
在这些方位漂移轨道设计方法中,以井段自然曲线模型为基础的设计方法较为实用。
在方位漂移轨道设计问题中,在给定的设计条件下,一般需要求解一个三元非线性方程组。
在使用数值迭代方法求解时,由于该方程组的特殊性,一般的迭代格式往往都是发散的(即使给定一个很好的迭代初值时也是如此)。
刘修善等[9]借鉴罚函数原理的基本思想,构造出了一种求解该方程组的数值算法。
卢毓周等[16]使用化简消元技术给出了直-增-稳-降-稳型漂移轨道设计的一种方法,将三元非线性方程组化简为一元非线性方程,并给出了具体的数值算法。
井眼轨迹计算及可视化研究现状陈刚;杨雪;潘保芝;蒋必辞【摘要】对现有的井眼轨迹计算方法进行对比分析,得出圆柱螺线法和最小曲率法计算精度较高,具有较好适用性.通过总结井眼轨迹可视化研究相关算法及可视化软件发展现状认为,可视化研究趋于普遍适用性、功能多元性和更加智能性方向发展.同时,对利用方位伽马测井曲线识别储层界面方法进行探索性研究.【期刊名称】《世界地质》【年(卷),期】2015(034)003【总页数】12页(P830-841)【关键词】水平井;井眼轨迹;可视化;插值计算;界面识别【作者】陈刚;杨雪;潘保芝;蒋必辞【作者单位】中煤科工集团西安研究院,西安710077;吉林大学地球探测科学与信息学院,长春130026;吉林大学地球探测科学与信息学院,长春130026;中煤科工集团西安研究院,西安710077【正文语种】中文【中图分类】P634.6近年来国内水平井、分支井逐年增多,直观了解水平井在地下的钻进情况,需对地下钻井进行透视观测,轨迹的可视化与控制就成为石油工业,尤其是钻井工业的热点。
通过三维动态显示井眼在地下钻进状态,及时调整钻进方向,使钻进轨迹保持在目的层中,防止邻井间相碰、大幅节约开发成本。
此外,可视化研究与随钻测井相结合,在油藏描述、分析地下地质构造、进行储层评价等方面具有重要意义[1,2]。
可视化图形技术是利用计算机图形学和现代计算机的OPP、DDE、OLE技术,把空间三维物体的基本特征数据以一种有效的方式输入计算机,利用近十几年来在井眼轨迹的可视化研究上,国内外在方法和技术上都有很大进步,技术由单一到复杂再到逼真。
由于国外在该领域起步较早,现有的技术相对于国内较成熟,软件方面功能也相对齐全和稳定。
国内在这方面的研究正处于发展阶段,经历从无到有,从简单到多元在不断地完善。
国内在可视化软件开发方面,在功能的多元化、适用的普遍性等方面还存在着不足。
因此,笔者对井眼轨迹可视化研究中的算法,包括井眼轨迹计算,插值计算进行了分析。
3钻柱摩阻计算与分析 (1)3.1 摩阻研究的意义、现状及存在问题 (1)3.1.1 研究目的和意义 (1)3.1.2 研究现状和发展趋势 (2)3.1.3 摩阻分析存在的主要问题 (3)3.2 钻柱动力学基本方程 (5)3.3 钻柱摩阻预测与分析 (11)3.3.1 摩阻分析方法与对象 (11)3.3.2 三维刚杆摩阻模型 (13)3.3.3 三维软杆摩阻模型 (22)3.3.4 摩阻分析结果的表述 (29)3.3.5 摩阻预测与分析 (31)3.3.6 摩擦系数反演 (35)3.3.7 摩阻影响因素分析 (35)3.3.8 计算结果验证与分析 (37)3.4 钻柱的屈曲 (39)3.4.1 斜直井中杆管柱屈曲的微分方程 (39)3.4.2 水平井段钻柱几何线性屈曲的数学模型 (44)3.4.3 水平井段钻柱几何非线性屈曲的数学模型 (45)3.4.4 斜直井段钻柱正弦屈曲和螺旋屈曲的临界载荷 (46)3.4.5 无重钻柱的几何线性螺旋屈曲 (50)3.4.6 无重钻柱的几何非线性螺旋屈曲 (52)3.5 实钻井眼轨道参数的插值计算 (55)3.5.1 概述 (55)3.5.2 典型轨迹模型插值 (55)3.5.3 典型曲线的曲率和挠率 (58)4 钻柱强度分析与校核 (64)4.1 钻柱应力状况分析 (64)4.1.1 轴向拉力产生的拉应力 (64)4.1.2 扭矩产生的剪应力 (64)4.1.3 钻柱弯曲产生的弯曲应力 (65)4.1.4 屈曲产生的应力和接触压力 (65)4.1.5 内压力产生的拉应力 (66)4.2 钻柱强度分析模型 (67)4.2.1 钻柱抗拉强度 (67)4.2.2 钻柱抗弯强度 (68)4.2.3 钻柱抗扭强度 (69)4.2.4 钻柱三轴应力强度 (71)4.3 钻柱强度校核 (72)3钻柱摩阻计算与分析3.1 摩阻研究的意义、现状及存在问题3.1.1 研究目的和意义随着油气田勘探开发的进行,钻井重点向深部、西部和海上发展。
◄钻井完井►doi:10.11911/syztjs.2023004引用格式:来建强,鲁港,周超,等. 井眼轨道模型中的数值积分计算[J]. 石油钻探技术,2023, 51(3):45-50.LAI Jianqiang, LU Gang, ZHOU Chao, et al. Numerical integral calculation in borehole trajectory model [J]. Petroleum Drilling Techniques ,2023,51(3):45-50.井眼轨道模型中的数值积分计算来建强1, 鲁 港2, 周 超3, 鲁天骐4(1. 北京四利通控制技术股份有限公司, 北京 101102;2. 中国石油辽河油田分公司勘探开发研究院, 辽宁盘锦 124010;3. 中国石油集团长城钻探工程有限公司工程技术研究院 辽宁盘锦 124010;4. 开姆尼茨工业大学测量与传感器工程系, 德国开姆尼茨 09126)摘 要: 为了提高钻井设计软件中所使用定积分的计算速度,从井眼轨道模型的基本公式出发,经过积分参数变换,将恒工具面法北东坐标定积分中的被积函数需要7次三角函数和1次对数函数的运算减少为3次三角函数和1次对数函数的运算,将空间圆弧法水平投影长度定积分的被积函数中的三角函数运算简化掉,从而提高定积分的计算速度;使用变步长Simpson 数值积分公式,计算恒工具面法北东坐标定积分和空间圆弧法水平投影长度定积分,可将计算精度控制到任意精度。
算例计算表明,本文算法除了提高了计算速度,计算精度还与公开数据完全相同。
本文给出的算法是对井眼轨道计算的进一步改进和完善,可应用于钻井软件开发中。
关键词: 定向钻井;井眼轨道;恒工具面模型;空间圆弧模型;数值积分中图分类号: TE21 文献标志码: A 文章编号: 1001–0890(2023)03–0045–06Numerical Integral Calculation in Borehole Trajectory ModelLAI Jianqiang 1, LU Gang 2, ZHOU Chao 3, LU Tianqi4(1. Beijing Sleton Control Technology Co., LTD, Beijing, 101102, China ; 2. Exploration and Development Research Institute,PetroChina Liaohe Oilfield Company, Panjin, Liaoning, 124010, China ; 3. Engineering Technology Research Institute, CNPC Greatwall Drilling Company, Panjin, Liaoning, 124010, China ; 4. Professorship of Measurement and Sensor Technology, Chemnitz University of Technology, Chemnitz, 09126, Germany )Abstract: In order to improve calculation speed of the definite integrals in the drilling design software, the basic formula of the borehole trajectory model and the integral parameter conversion were adopted. The calculation of seven trigonometric and one logarithmic functions was reduced to that of three trigonometric and one logarithmic functions required by the integrand of the definite integral in the north and east coordinates by the constant tool face method. The trigonometric functions in the integrand of the definite integral with a horizontal projection length by the space arc method were simplified, so as to improve the calculation speed of the definite integral. Using the variable step size Simpson numerical integration formula to calculate the above two kinds of definite integrals, the calculation accuracy can be controlled to arbitrary precision. The calculation shows that the new algorithm in this paper not only improves the calculation speed but also has the same calculation accuracy as the public data. The algorithm given in this paper is further improves and refines of the borehole trajectory calculation and can be applied in the software development of drilling.Key words: directional drilling; borehole trajectory; constant tool face model; spatial arc model; numerical integral在进行井眼轨道设计[1–7]、随钻测量数据处理时[8–15],需要计算北东坐标、垂深、水平投影长度等参数,如果井眼轨道上每一点的井斜角和方位角与井深呈线性关系,则这些参数可以通过解析计算公收稿日期: 2021-12-07;改回日期: 2023-02-23。
井眼轨道参数的插值计算由于实钻井眼轨道的测点与钻柱单元体的划分可能并不一致,因此钻柱单元体边界点对应的井眼轨道参数必须靠插值计算获得。
插值结果的准确与否,对钻柱单元体的受力计算有着直接的影响。
因此,提高插值计算的精度具有重要意义。
由于测点是离散的,无法知道各测段内井眼轨道的实际形态,所以测段内某点几何参数的计算方法都是建立在一定假设的基础上的。
这些计算方法多数是将测段内的井眼轨道假设为直线、折线和曲线等,早期,由于计算机能力的限制,以平均角法和平衡正切法为代表直线或折线假设,因其计算简单快速,曾经被广泛应用,但随着钻井技术的发展,弯曲的井眼轨迹增多,如果仍采用直线或折线假设,则计算精度相对较低。
由于计算技术的高速发展,直线或折线假设,目前几乎淘汰,取而代之的是以圆柱螺线和空间圆弧曲线等为代表的曲线假设,大行其道。
在进行插值计算时,各插值点的坐标增量可以采用不同的计算方法,但坐标值的累加形式是相同的,即(X 为东向位移,Y 为北向位移, Z 为垂直向位移,S 为水平位移)⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧∆+=∆+=∆+=∆+=∆+=∆+=φφφa a αS S S Z Z Z YY Y X X X 121212121212所以,在以下的计算方法中将只给出坐标增量的计算式。
典型轨迹模型插值1、正切法:正切法又称下切点法,或下点切线法。
此法假定两相邻测点之间的孔段为一条直线,长度等于测距,该直线的井斜角和井斜方位角等于下测点的井斜角和井斜方位角,整个钻孔轨迹是直线与直线相连接的空间折线。
正切法井身轨迹计算图如图1所示,1、2 是孔身轨迹上相邻的两个测点,1′、2′是 1、2 两个测点的水平投影。
该测段的井斜角和井斜方位角等于下测点 2 的井斜角和井斜方位角。
对于切线法,上下两个相邻测点间各参数的计算公式如下:222222cos sin sin sin sin cos φαφαααL Y L X L S L Z ∆=∆∆=∆∆=∆∆=∆式中:Z ∆——测段上下测点间垂直深度的分量(增量)(以下同); L ∆——测段上下测点间沿钻孔轴线的距离(以下同);Y ∆∆X ——分别为测段上下测点间水平位移在 X 轴(西东方向)的分量(增量);水平位移在 Y 轴(南北方向)的分量(增量)(以下同);22 φα——分别为测段下测点的井斜角和井斜方位角。
2、平均角法平均角法井眼轨迹计算图如图所示,1、2 是孔身轨迹上相邻的两个测点,1′、2′是 1、2 两个测点的水平投影,该测段的井斜角和井斜方位角等于上下两个测点的井斜角和井斜方位角的平均值。
假设测段内的井眼轨道为一条直线,其方向是上、下两侧点井眼方向的平均值,则有⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧==∆=∆∆=∆∆=∆∆=∆cc c cc c c c L S L Z L Y L X φφααααφαφαsin cos sin sin cos sin 其中⎪⎪⎩⎪⎪⎨⎧+=+=222121φφφαααc c 需要注意的是,当某测点的井斜角等于零时,是没有井斜方位角的。
“没有井斜方位角”,并不等于“井斜方位角等于零”。
这种情况下的平均井斜方位角可作如下处理:当01=α时,2φφ=c ;02=α时,1φφ=c ;3、平衡正切法这种方法认为井斜角和方位角,在测段的开始和末尾全部切线补偿。
根据测量的井斜角和方位角,用三角函数平均值确定钻孔轴线坐标值。
它也是把小段钻孔轴线作为折线来处理,并与以下要叙述的最小曲率法有相似的计算公式,仅少了一项修正系数。
该法得出一平滑曲线,较接近两测点间实际的钻孔轴线。
直观地看,两测点间距离越大,可能产生的误差越大。
平衡正切法井眼轨迹计算图如图所示,1、2 是孔身轨迹上相邻的两个测点,1′、2′是 1、2 两个测点的水平投影。
上段直线用上测段的井斜角和井斜方位角,下段直线用下测点的井斜角和井斜方位角。
这个钻孔轨迹是一条折点更多的空间折线。
假设测段内的井眼轨道为折线,两个线段的长度均等于测段长度的一半,其方向分别与上、下侧点的井眼方向相同,则有()()()()B L B L B S B L B Z B L B Y B L B X BL L S L Z L Y L X ≥∆⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧==-∆+=∆-∆+=∆-∆+=∆-∆+=∆<∆⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧==∆=∆∆=∆∆=∆∆=∆2221212211221111111111sin sin cos cos sin sin sin sin cos sin cos sin sin cos sin sin cos sin φφααααααφαφαφαφαφφααααφαφα其中2)(12L L B -= . 当测点的φ = 0时,处理方法同平均角法。
4、曲率半径法这种方法是角(其切线)在每一测段开始与末尾经常被描述为一曲线,它将代表钻孔的真实轴线。
此曲线具有球面圆弧形状平滑,可用圆周或球面的一部分表示。
圆弧的精确确定由两个方向的矢量和已知的两测点间的距离所给定。
由于这种假设,因测点间大距离所造成的误差小于其他计算方法,从而该法成为钻孔空间坐标计算最精确的方法之一。
曲率半径法井眼轨迹计算图假设测段内的井眼轨道在垂直剖面图和水平投影图上均为圆弧,其有下面几种计算方法: 第一种表达形式:()()()()⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧∆+=∆+=-=∆-=∆-=∆-=∆r SR L R S R Z r Y r X 111111cos cos sin sin cos cos sin sin φφααααααφφφφ 式中()12211212cos cos φφααR r ααL L R --=--=第二种表达形式:)sin (sin )cos (cos )cos (cos )sin (sin 12212112φφφφαααα-=∆-=∆-∆=∆-=∆r Y r X L S R Z式中:1212221180)cos (cos 180φφφαααπφαααπα-=∆-=∆⎪⎭⎫ ⎝⎛⋅∆⋅∆-∆=⋅∆∆=L r L R第三种表达形式:21221221212112180)sin )(sin cos (cos 180)cos )(cos cos (cos 180)cos (cos 180)sin (sin ⎪⎭⎫⎝⎛⋅∆⋅∆--∆=∆⎪⎭⎫⎝⎛⋅∆⋅∆--∆=∆⋅∆-∆=∆⋅∆-∆=∆πφαφφααπφαφφααπαααπαααL Y L X L S L Z这是最常用的公式。
第四种表达形式:22180cos sin 2sin 2sin 4180sin sin 2sin 2sin 4180sin 2sin 2180cos 2sin2⎪⎭⎫ ⎝⎛⋅∆⋅⋅∆⋅∆⋅∆=∆⎪⎭⎫ ⎝⎛⋅∆⋅⋅∆⋅∆⋅∆=∆⋅∆⋅∆⋅∆=∆⋅∆⋅∆⋅∆=∆παφαφαπαφαφαπαααπαααcc cc ccL Y L X L S L Z式中:221ααα+=c ,221φφφ+=c 。
上述公式中,在分母中都有α∆ 和φ∆,只要其中一个为零,都会使公式无法计算,所以在实际应用中要考虑到以下几种特殊情况:a) 第一种特殊情况,21αα=,21φφ≠,即0,0≠∆=∆φα。
此时测段计算公式如下:πφφφαπφφφααα180sin sin sin 180cos cos sin sin cos 21221222⋅∆-⋅∆=∆⋅∆-⋅∆=∆∆=∆∆=∆L Y L X L S L Z b)第二种特殊情况,21αα≠,21φφ=,即0,0=∆≠∆φα。
此时测段计算公式如下:πφαααπφααααααααα180cos cos cos 180sin cos cos cos cos cos sin 2212212112⋅⋅∆-⋅∆=∆⋅⋅∆-⋅∆=∆∆-⋅∆=∆∆-⋅∆=∆L Y L X L S L Z c)第三种特殊情况,21αα=,21φφ=,即0,0=∆=∆φα。
此时实际上是按着正切法的公式进行计算:222222cos sin sin sin sin cos φαφαααL Y L X L S L Z ∆=∆∆=∆∆=∆∆=∆d )第四种特殊情况,21αα≠,且其中之一等于零,则为零的该测点的方位角是不存在的。
此时,可按二测点方位角相等来处理,然后代入第二种特殊情况的公式中计算。
5、校正平均角法校正平均角法是在曲率半径法公式的基础上,简化处理而导出的一种新方法。
其简化思路如下: 将 sin x 展开成麦克劳林无穷级数的形式:-+-+-=!9!7!5!3sin 9753x x x x x x级数收敛很快,可近似取前两项,即3361!3sin x x x x x -=-=。
把曲率半径法计算公式的第四种表达式中的2sinα∆和2sin φ∆作上述近似处理,得到: )2411(24822sin )2411(24822sin 2323φφφφφααααα∆-∆=∆-∆=∆∆-∆=∆-∆=∆。
将此二式代入曲率半径法计算公式的第四种表达式中,并忽略去高次微量,可使公式大为简化,即可得校正平均角法的计算公式:cc cc ccL Y L X L S L Z φαφαφαφαααααcos sin )241(sin sin )241(sin )241(cos )241(222222∆∆+∆-=∆∆∆+∆-=∆∆∆-=∆∆∆-=∆。
令2412α∆-=H f ,24122φα∆+∆-=A f 可将上式简化。
H f 和 A f 是两个小于 1 而接近于 1 的数。
当α∆ 和 φ∆ 足够小时, H f 和 A f 可近似看作等于 1,则公式完全变成了平均角法的公式。
所以,可把 H f 和 A f 看作是一个校正系数。
从公式的形式上看,它是在平均角法的基础上乘以校正系数H f 和A f ,因此取名叫校正平均角法。
平均角法是直线法,而校正平均角法在实质上是曲线法,是从曲线法的曲率半径法推演出来的。
校正平均角法形式上是直线法,实质上是曲线法,其假设更接近真实井身,因而较直线法更为精确。
6、最小曲率法最小曲率法假设两测点间的井段是一段平面上的圆弧,圆弧在两端点处与上下二测点处的井身方向线相切,即在上、下二测点的井身方向一定的情况下,把测段看成圆弧曲线,乃是所有曲线中曲率最小的曲线,所以被定名为最小曲率法。
与曲率半径法相同,圆弧的精确确定也是由测段开始与末尾两个方向的 矢量和已知的两测点间的距离所给定。
假设测段内的井眼轨道为空间圆弧,则有()()()()[]⎥⎦⎤⎢⎣⎡∆⋅-+∆⋅++=∆-∆=-=∆∆∆⋅+=∆+=∆+=∆+=∆)/tan()sin sin cos cos (cos cos sin )/tan()sin cos cos sin (cos sin sin arctan )/sin(cos sin )/cos(cos arccos2/tg 2sin sin cos cos sin sin sin sin cos sin cos sin 1111111111111111111R L ωφωφαφαR L ωφωφαφαφR L ωαR L ααφφφφφααλS ααλZ φαφαλY φαφαλX (3.108)其中()12212112cos sin sin cos cos cos 2t gφφααααεεL L R RL R λ-+=-=∆⋅= ()()⎥⎦⎤⎢⎣⎡---=2112112tg tg cos cos sin tg ααφφαφφω7、自然参数法(同下面的自然参数模型, 可以参考)假设测段内的井眼轨道的井斜和方位随段长均匀变化,则有()()[]()()[]()()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧∆+=∆+=-=∆-=∆∆+--∆+-=∆∆+-+∆+-=∆Lk φφL k ααk ααS k ααZ k L k ααk L k ααY k L k ααk L k ααX φαααB B B B A A A A B B B B A A A A 1111cos cos sin sin 2)sin(sin )sin(sin 2)cos(cos )cos(cos (3.109) 其中 φαB φαA B A φαk k k k k k φααφααL L φφk L L ααk +=-=+=-=--=--=111112121212 8、弦步法弦步法也假设相邻两测点间的井身轴线为一空间平面上的圆弧曲线。