CN107153738A - 一种基于恒工具面角的井眼轨道建模方法 - Google Patents

一种基于恒工具面角的井眼轨道建模方法 Download PDF

Info

Publication number
CN107153738A
CN107153738A CN201710342223.XA CN201710342223A CN107153738A CN 107153738 A CN107153738 A CN 107153738A CN 201710342223 A CN201710342223 A CN 201710342223A CN 107153738 A CN107153738 A CN 107153738A
Authority
CN
China
Prior art keywords
mrow
msub
msubsup
msup
alpha
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710342223.XA
Other languages
English (en)
Inventor
唐雪平
洪迪峰
盛利民
窦修荣
王家进
高文凯
王鹏
彭烈新
滕鑫淼
吕海川
蓝蓉
范锦辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
CNPC Engineering Technology R&D Co Ltd
Original Assignee
China National Petroleum Corp
CNPC Drilling Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, CNPC Drilling Research Institute Co Ltd filed Critical China National Petroleum Corp
Priority to CN201710342223.XA priority Critical patent/CN107153738A/zh
Publication of CN107153738A publication Critical patent/CN107153738A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种轨道建模方法,属于井下钻探技术领域,具体是涉及一种基于恒工具面角的井眼轨道建模方法。本发明建立限定目标点位置,以及同时限定目标点位置和入靶方向两种典型三维恒工具面角井眼轨道设计模型,能根据不同的设计要求灵活、快速求解设计约束方程,使设计具普遍性和灵活性,为井眼轨道设计和控制提供理论依据。

Description

一种基于恒工具面角的井眼轨道建模方法
技术领域
本发明涉及一种轨道建模方法,属于井下钻探技术领域,具体是 涉及一种基于恒工具面角的井眼轨道建模方法。
背景技术
随着定向钻井技术的发展,对井眼轨道设计和控制提出了更高的 要求。在定向钻井作业中,涉及的井眼轨道多为三维井眼轨道,如侧 钻井水平井、地质导向钻井和待钻井眼轨道,要求井眼轨道控制尽量 符合设计要求,且易于定向控制,确保有效、快速地钻达目的层位。
恒工具面角曲线井眼轨道的突出特点为井眼曲率、井斜角变化率 和工具面角均为常数,便于定向操作,可降低钻井作业费用,作为一 种典型的井眼轨道模型用于定向钻井中。由于恒工具面角曲线的北、 东坐标涉及到含三角函数的数值积分计算,至使井眼轨道模型求解困 难,需要研究可行、有效的数值计算方法来实现。
空间斜平面圆弧模型广泛应用于定向钻井中,三维圆弧轨道需不 断变化工具面角才能实现,实际施工难以做到,不能精确控制钻进轨 迹,且作业效率低。现有恒工具面角模型通常针对某一特定问题提出, 缺乏普遍适用性,且参数求解方式单一。为了满足实际工程应用的需 要,井眼轨道模型要有通用性,求解参数方式应有多样性,然而目前 这方面的研究还不完善和系统,难以满足定向钻井轨道设计及控制的 需要。
在三维井眼轨道设计和控制时,起始点的位置、井眼方向和目标 点的位置已确定,而目标点的井眼方向可分为不限定和限定两种类 型,分别对应定向井和水平井要求。为此建立两种典型的恒工具面角 井眼轨道模型,能演化出多种剖面形式,可灵活、快速求解,以满足 多种井眼轨道设计和控制的需要。
发明内容
本发明的目的在于建立限定目标点位置,以及同时限定目标点位 置和入靶方向两种典型三维恒工具面角井眼轨道设计模型,能根据不 同的设计要求灵活、快速求解设计约束方程,使设计具普遍性和灵活 性,为井眼轨道设计和控制提供理论依据。
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种基于恒工具面角的井眼轨道建模方法,基于以下约束建立轨 道模型:
其中,ΔX=XT-XS,ΔY=YT-YS,ΔZ=ZT-ZS,并且,
α=αi+Kαixi(曲线段i)
(曲线段i)
式中,Li为直线段长度,单位m;Si为曲线段长度,单位m;αi为井斜角,单位rad;φi为方位角,单位rad;Ki为井眼曲率,单位 rad/m;Kαi为井斜角变化率,单位rad/m;Kαi=Kicosωi;ωi为工具面角, rad;XS、YS、ZS分别为钻井起始点的北坐标、东坐标和垂深,m; XT、YT、ZT分别为钻井目标点的北坐标、东坐标和垂深,单位m; ΔX、ΔY、ΔZ分别为钻井目标点和起始点的北坐标、东坐标和垂深差, 单位m。
优化的,上述的一种基于恒工具面角的井眼轨道建模方法,所述 轨道为三段制,并且基于以下公式求解空间斜平面:
式中,L为起始点和目标点的直线距离,单位m;Ts为L在起始点 切线矢量上的投影长度,单位m;L1、L2为圆弧段前后两直线段长度, 单位m。
优化的,上述的一种基于恒工具面角的井眼轨道建模方法,给定 目标点的井斜角α2,基于下式求解圆弧段两切线交点到起始点的长 度:
式中,a=cos2α1-cos2α2,b=-2[ΔZcosα1-Tscos2α2],c=ΔZ2-L2cos2α2。当 α2>90°时,取“+”;当α2≤90°时,取“-”。
优化的,上述的一种基于恒工具面角的井眼轨道建模方法,给定 目标点的方位角φ2,基于下式求解圆弧段两切线交点到起始点的长 度:
基于下式求解圆弧段两切线交点到目标始点的长度:
优化的,上述的一种基于恒工具面角的井眼轨道建模方法,基于 下式求解圆弧段弯曲角:
θ=arccos[(ΔXsinα1cosφ1+ΔYsinα1sinφ1+ΔZcosα11)/λ2]
基于下式求解直线段L1的长度:
优化的,上述的一种基于恒工具面角的井眼轨道建模方法,所述 轨道为五段制,并且基于给定变量L1、L2、L3、K1和K2中的任意4个 参数求解另一参数;或者,当K1=K2时,基于变量L1、L2、L3,求解K1和K2
优化的,上述的一种基于恒工具面角的井眼轨道建模方法,给定 L1、L3、K1和K2,求解L2的具体步骤为:
步骤1,以L2为基本求解变量,再选择α2和φ2变量为未知数, 作为方程组(1)的3个求解参数;
步骤2,根据设计起始点和目标点的坐标及方向,以及给定L1、 L3、K1和K2已知参数,求对应空间圆弧轨道模型的解,并将其结果L2、 α2和φ2作为对应恒工具面角模型解的初值;
步骤3,求Kα1、Kα2、S1和S2,确定两曲线段方位角的增减符号;
式中,增斜取“+”,降斜取“-”,稳斜时,Kα1=0,当α1=0时,Kα1=K1
同理,可求得Kα2和S2
步骤4,用拟牛顿法求解方程组(1);
步骤5,若方程组有解,则计算轨道节点参数,若无解,则提示无 解,可调整已知参数或重新确定参数求解方式进行求解。
一种基于恒工具面角的井眼轨道规划方法,采用上述任一种井眼 轨道建模方法建立轨道模型。
因此,本发明具有如下优点:(1)三段制实现限定目标点位置及 五段制实现同时限定目标点位置和入靶方向两种典型三维恒工具面 角井眼轨道设计,可满足定向井和水平井设计一般需求,能演化出一 到五段制多种剖面形式,还可扩充到五段制以上剖面,具有一般性; (2)参数求解方式组合多,并给出了对应初值求解方法,能根据不 同的设计要求灵活、快速求解设计约束方程,使设计具普遍性和灵活 性,为井眼轨道设计和控制提供理论依据。
附图说明
图1单恒工具面角段轨道(三段制轨道)示意图;
图2双恒工具面角段轨道(五段制轨道)示意图;
图3三段制轨道求解设计流程示意图;
图4五段制轨道求解设计流程示意图;
图5恒工具面角轨道设计方法流程图。
图6井眼轨道控制方法。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具 体的说明。
实施例:
1、恒工具面角轨道设计模型
井眼轨道设计采用北东地坐标系OXYZ,O为井口,X轴指向正北, Y轴指向正东,Z轴指向铅垂方向。为了满足定向井和水平井轨道设 计需要,建立“直线段-恒工具面角段-直线段”和“直线段-恒工具 面角段-直线段-恒工具面角段-直线段”两类设计模型,其中直线段 长度可为零,能够满足限定目标位置及同时限定目标位置和方向的两 种典型轨道设计要求。
两种恒工具面角轨道设计模型如图1和2所示,其约束方程为:
其中,ΔX=XT-XS,ΔY=YT-YS,ΔZ=ZT-ZS,并且,
式中,Li为直线段长度,m;Si为曲线段长度,m;αi为井斜 角,rad;φi为方位角,rad;Ki为井眼曲率,rad/m;Kαi为井斜角变 化率,rad/m;Kαi=Kicosωi;ωi为工具面角,rad。
当恒工具面角段初始井斜角为0时,则该段退化二维圆弧井段。 由于直线段长度可为零,能演化出多种井眼轨道设计剖面形式,以满 足井眼轨道设计和控制的需要。
2、恒工具面角轨道设计模型求解方法
轨道设计模型含有3个独立约束方程式,可求解3个未知数。为 了满足设计灵活求解的需要,以直线段长L1、L2、L3或曲率K1、K2为1个基本求解变量,选择另外2个变量为未知数,其余参数为已知 数进行求解。由于约束方程中涉及含三角函数的数值积分计算,采用拟牛顿数值法进行求解,可避免复杂的求导运算。应用拟牛顿法求解 时,求解变量初值选择极为重要。选择可行解附近值为求解变量初值, 可快速收敛到正确结果。在确定3个求解变量组合后,先按对应的圆 弧轨道模型求解,其结果作为恒工具面角轨道设计模型的初值,再用 拟牛顿法求解,从而使复杂的非线性方程组求解问题得到解决。
2.1、单恒工具面角轨道设计模型求解
对单恒工具面角轨道设计模型,给定变量L1、L2和K1中的任意两 个参数,可求解另一参数,以及在限定目标点的井斜角或方位角情况 下,给定3个变量中的任意一个,可求解另外两个变量,共计9种求 解组合,如图3所示。
空间斜平面三段制求解公式如下:
其中,L为起始点和目标点的直线距离,m;Ts为L在起始点切线 矢量上的投影长度,m。
当给定目标点的井斜角α2时,可求圆弧段两切线交点到起始点的 长度为:
其中,a=cos2α1-cos2α2,b=-2[ΔZcosα1-Tscos2α2],c=ΔZ2-L2cos2α2。当 α2>90°时,取“+”;当α2≤90°时,取“-”。
当给定目标点的方位角φ2时,可求圆弧段两切线交点到起始点的 长度为:
圆弧段两切线交点到目标始点的长度为:
当限定目标点的井斜角或目标点的方位角时,由(5)或(6)及 (7)式,可求出λ1,λ2
圆弧段弯曲角为:
θ=arccos[(ΔXsinα1cosφ1+ΔYsinα1sinφ1+ΔZcosα11)/λ2] (9)
进而求得以下关系式:
由(2)~(9)式,给定设计变量L1、L2、K1中的任意一个可解 析求得另外两个变量。
2.2、双恒工具面角轨道设计模型求解
对双恒工具面角轨道设计模型,给定变量L1、L2、L3、K1和K2中 的任意4个参数,可求解另一参数,以及给定变量L1、L2、L3,求解 K1=K2,共计6种求解组合,如图4所示。
在确定3个求解变量组合和求解恒工具面角轨道设计模型初值 后,再用拟牛顿法求解,可快速收敛到恒工具面角轨道设计非线性方 程组解的正确结果。
以双恒工具面角轨道设计模型,以给定L1、L3、K1和K2,求解L2为例,介绍设计模型求解步骤:
步骤1,以L2为基本求解变量,再选择α2和φ2变量为未知数, 作为方程组(1)的3个求解参数;
步骤2,根据设计起始点和目标点的坐标及方向,以及给定L1、 L3、K1和K2已知参数,求对应空间圆弧轨道模型的解,并将其结果L2、 α2和φ2作为对应恒工具面角模型解的初值;
步骤3,求Kα1、Kα2、S1和S2,确定两曲线段方位角的增减符号;
式(10)中,增斜取“+”,降斜取“-”,稳斜时,Kα1=0,当α1=0时, Kα1=K1
同理,可求得Kα2和S2
步骤4,用拟牛顿法求解方程组(1);
步骤5,若方程组有解,则计算轨道节点参数。
根据模型求解方法,可编制两种典型设计模型的15种求解组合 计算程序,求解灵活、快速,可满足定向设计、水平井、侧钻井及待 钻轨道设计和控制需求。
3、井眼轨道设计方法
恒工具面角井眼轨道设计方法如图5所示,具体包括以下步骤。
步骤1,选择轨道设计模型
根据井眼轨道设计和控制需要,选择轨道设计模型。当需要限定 井眼轨道入靶方向时,则需选择双恒工具面角轨道设计模型。
步骤2,确定设计基础数据
给定井眼轨道设计的起始点和目标点的已知数据,即起始点坐标 (XS,YS,ZS)及方向(α1,φ1),目标点坐标(XT,YT,ZT)。对双恒工具 面角轨道设计模型,给出目标点方向(α3,φ3)。对单恒工具面角轨道 设计模型,若限定目标点井斜角或方位角时,给出目标点的井斜角α2或方位角φ2。对于侧钻井或井眼轨道调整设计,通过测斜数据处理方 法得到设计起始点的数据。
步骤3,计算分析原始数据
根据步骤2中的数据,可计算出起始点和目标点间的空间距离L 和起始点到目标点的最小井眼曲率Kmin
若三段制轨道目标点在起始点切线上或五段制轨道起始点和目 标点切线共线时,则稳斜钻到目标点即可,转步骤8;否则到步骤4。
步骤4,确定参数求解方式
对单恒工具面角轨道设计模型,在3个设计参数L1、L2、K1中, 可任选2个为已知参数,求解另一参数。若限定目标点的井斜角α2或方位角φ2时,可任选L1、L2、K1中一个为已知参数,求解另外两个 参数。共计9种求解组合,如图3所示。
对双恒工具面角轨道设计模型,在5个设计参数L1、L2、L3、K1、 K2中,可任选4个为已知参数,求解另一参数,以及给定设计参数L1、 L2、L3,求K1=K2。共计6种求解组合,如图4所示。
设计参数的取值范围参考步骤3中的计算结果L、Kmin和工程上给 定许可的最大井眼曲率Kmax
步骤5,给定已知参数
根据轨道设计和控制需求,给定步骤4参数L1、L2、L3、K1、K2中选择的已知变量。
步骤6,计算求解变量初值
按对应空间斜平面轨道设计模型进行求解,将其结果作为恒工具 面角轨道设计模型求解变量的初值。
步骤7,求解恒工具面角轨道设计模型
应用拟牛顿法求解恒工具面角轨道设计模型。若无解,转步骤5 调整已知参数,或转步骤4重新选择参数求解方式。
步骤8,计算轨道节点参数
求出未知参数后,按恒工具面角轨道模型计算节点的井眼轨道参 数,包括井深,井斜角、方位角和坐标数据等。
步骤9,计算轨道明细数据
按一定步长计算曲线段和稳斜段轨道参数。
步骤10,输出轨道设计数据。
4、井眼轨道控制方法
井眼轨道控制方法根据设计轨道剖面形式、设计井眼曲率大小、 井眼轨道控制工具和要求来确定。控制方法如图6所示,基本步骤如 下:
步骤1,确定控制目标及要求
根据设计轨道数据,确定控制目标及允许偏差,满足油气勘探开 发及钻井工艺要求。
步骤2,设计BHA及钻柱结构
根据轨道控制目标及要求,确定钻井轨道控制方式,设计下部钻 具组合BHA,包括钻头、BHA类型、稳定器个数和外径、钻具结构尺 寸、螺杆钻具的弯角大小、旋转导向工具类型及规格;根据轨道设计、 控制工艺、机械钻速和钻头使用时间等,预计钻井进尺,进行摩阻、 扭矩分析和钻柱强度校核,优化钻柱结构。
步骤3,工具配置、测试及下钻
按设计BHA和钻柱结构组合钻具,对测量控制工具进行配置和测 试,测试合格后,组合钻具下钻。
步骤4,按设计参数导向钻进
钻具下钻到底后,按设计井眼轨道和钻井参数进行导向钻进。在 钻进过程中,随钻测量各种参数,并上传到地面,实时监测导向情况 和地层变化等。
步骤5,暂停钻进作业测斜
当钻达预定深度时或需要测斜时,暂停钻进作业,静态测量井斜 数据,确保测斜数据准确。
步骤6,轨迹跟踪及偏差分析
根据测斜数据,进行轨迹跟踪及偏差分析。实钻轨迹和设计轨道 的偏差包括距离偏差、方向偏差、井斜角和方位角偏差等。若偏差在 允许范围内,继续按原设计轨道钻进,直至钻达预定目标;若偏差超 出允许范围,且不在钻具可控范围内时,则结束钻进作业。
步骤7,轨道调整设计
若偏差超出允许范围或者地质目标发生变化,且在钻具可控范围 内时,则进行轨道调整设计。根据当前井底位置、井斜角、方位角和 调整目标点数据,合理选择剖面类型进行待钻轨道设计,并按待钻设 计进行轨道控制。
采用本发明提出的恒工具面角轨道模型进行控制时,易于常规导 向工具的定向施工操作,也可减少旋转导向系统的工具面角下传指令 数目,节约定向作业时间,提高工作效率。
步骤8,结束钻进作业
当钻达预定目标或需更换钻头等时,则结束钻进作业。
5、效果评估
本实施例的方法可以满足限定目标点位置以及同时限定目标点 位置和方向的恒工具面角井眼轨道设计和控制要求。模型可实现多种 剖面组合形式,具有普遍适用性;给出了恒工具面角井眼轨道约束方 程求解方法,设计和控制方法简单、实用,易于井眼控制,可提高设 计效率,降低定向作业时间和成本,实现安全、快速、优质钻完井作 业的目的。
下面采用两个具体实施例评估本实施例的效果。
例1:某定向井当前点的坐标为:XS=100m,YS=―100m,ZS=2000m; 井斜角和方位角分别为α1=10°,φ1=315°。目标点的坐标为: XT=200m,YT=―1100m,ZT=4000m。若继续按当前方向钻进L1=200m, 再用造斜能力为K1=3°/30m的工具来调整井眼轨道。设计结果为以 恒工具面角ω=-33.28°钻进S1=235.54m后再用稳斜钻具钻进 L2=1821.75m可准确达到设计目标点。其井眼轨道数据见表1。
表1定向井轨道设计数据(三段制轨道)
例2:某三维水平井终点坐标XT=31.70m,YT=-1288.80m, ZT=869.42m,水平段井斜角α3=90.36°、方位角φ3=270°,段长 L3=941.72m。从垂直井段造斜,造斜点深L1=532m,曲率K1=K2=6°/30m, 用恒工具面角法设计井眼轨道。
应用双恒工具面角轨道设计模型,求出对应圆弧轨道稳斜段长 度、井斜角及方位角,作为拟牛顿法求解恒工具面角轨道约束方程的 初值,经2次数值迭代计算,求得结果为L2=80.706m,α2=46.164°, φ2=277.686°,井眼轨道数据见表2。
表2水平井轨道设计数据(五段制轨道)
由表2可知,第二个曲线井段为造斜扭方位井段,按恒工具面进 行轨道控制,随着井深的变化,其工具面角恒为-8.87°,易于定向 施工作业;当采用旋转导向系统时,可显著减少地面下传工具面角指 令数,可提高作业效率和控制精度,达到快速、优质的钻井目的。
例3某水平井设计的着陆点井斜和方位为α3=88°,φ3=50°。 当前点井斜角和方位角为α1=76°,φ1=54°,距着陆点的坐标分别 为:ΔX=52.70m,ΔY=57.51m,ΔZ=12m,如何设计轨道使其准确着 陆。
因当前点距着陆点的垂深仅有12m,因此设计一种最简单的双恒 工具面角轨道,即恒工具面角段+恒工具面角段形式(K1=K2),以便一 套常规弯螺杆导向工具能使其准确着陆。设计结果为K1=K2=8.327°/ 30m,轨道数据见表3。
表3水平井轨道设计数据(同曲率双恒工具面角段)
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明 所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或 补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权 利要求书所定义的范围。

Claims (7)

1.一种基于恒工具面角的井眼轨道建模方法,其特征在于,基于以下约束建立轨道模型:
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>L</mi> <mi>i</mi> </msub> <msub> <mi>sin&amp;alpha;</mi> <mi>i</mi> </msub> <msub> <mi>cos&amp;phi;</mi> <mi>i</mi> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <msub> <mi>S</mi> <mi>i</mi> </msub> </msubsup> <msub> <mi>sin&amp;alpha;cos&amp;phi;dx</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>&amp;Delta;</mi> <mi>X</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>L</mi> <mi>i</mi> </msub> <msub> <mi>sin&amp;alpha;</mi> <mi>i</mi> </msub> <msub> <mi>sin&amp;phi;</mi> <mi>i</mi> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <msub> <mi>S</mi> <mi>i</mi> </msub> </msubsup> <msub> <mi>sin&amp;alpha;sin&amp;phi;dx</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>&amp;Delta;</mi> <mi>Y</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>L</mi> <mi>i</mi> </msub> <msub> <mi>cos&amp;alpha;</mi> <mi>i</mi> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mfrac> <mrow> <msub> <mi>sin&amp;alpha;</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>sin&amp;alpha;</mi> <mi>i</mi> </msub> </mrow> <msub> <mi>K</mi> <mrow> <mi>&amp;alpha;</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>=</mo> <mi>&amp;Delta;</mi> <mi>Z</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mrow> <mo>(</mo> <mi>n</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mi>n</mi> <mo>=</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中,ΔX=XT-XS,ΔY=YT-YS,ΔZ=ZT-ZS,并且,
α=αi+Kαixi(曲线段i)
式中,Li为直线段长度,单位m;Si为曲线段长度,单位m;αi为井斜角,单位rad;φi为方位角,单位rad;Ki为井眼曲率,单位rad/m;Kαi为井斜角变化率,单位rad/m;Kαi=Kicosωi;ωi为工具面角,rad;XS、YS、ZS分别为钻井起始点的北坐标、东坐标和垂深,m;XT、YT、ZT分别为钻井目标点的北坐标、东坐标和垂深,单位m;ΔX、ΔY、ΔZ分别为钻井目标点和起始点的北坐标、东坐标和垂深差,单位m。
2.根据权利要求1所述的一种基于恒工具面角的井眼轨道建模方法,其特征在于,所述轨道为三段制,并且基于以下公式求解空间斜平面:
<mrow> <msub> <mi>K</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msqrt> <mrow> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> <mrow> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>L</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>L</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mn>2</mn> <msub> <mi>L</mi> <mn>1</mn> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>L</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>L</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>+</mo> <mn>2</mn> <msubsup> <mi>K</mi> <mn>1</mn> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msup> <mrow> <mo>(</mo> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mrow> </msqrt> </mrow>
<mrow> <msub> <mi>L</mi> <mn>2</mn> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>+</mo> <msubsup> <mi>L</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mn>2</mn> <msub> <mi>L</mi> <mn>1</mn> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <mn>2</mn> <msubsup> <mi>K</mi> <mn>1</mn> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msup> <mrow> <mo>(</mo> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mrow> </msqrt> </mrow>
式中,L为起始点和目标点的直线距离,单位m;Ts为L在起始点切线矢量上的投影长度,单位m;L1、L2为圆弧段前后两直线段长度,单位m。
3.根据权利要求2所述的一种基于恒工具面角的井眼轨道建模方法,其特征在于,给定目标点的井斜角α2,基于下式求解圆弧段两切线交点到起始点的长度:
<mrow> <msub> <mi>&amp;lambda;</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>b</mi> <mo>&amp;PlusMinus;</mo> <msqrt> <mrow> <msup> <mi>b</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>a</mi> <mi>c</mi> </mrow> </msqrt> </mrow> <mrow> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> </mrow>
式中,a=cos2α1-cos2α2,b=-2[ΔZcosα1-Tscos2α2],c=ΔZ2-L2cos2α2;当α2>90°时,分子中的符号取“+”;当α2≤90°时,分子中的符号取“-”。
4.根据权利要求2所述的一种基于恒工具面角的井眼轨道建模方法,其特征在于,给定目标点的方位角φ2,基于下式求解圆弧段两切线交点到起始点的长度:
<mrow> <msub> <mi>&amp;lambda;</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>X</mi> <mi> </mi> <msub> <mi>tan&amp;phi;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>Y</mi> </mrow> <mrow> <msub> <mi>sin&amp;alpha;</mi> <mn>1</mn> </msub> <msub> <mi>cos&amp;phi;</mi> <mn>1</mn> </msub> <msub> <mi>tan&amp;phi;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>sin&amp;alpha;</mi> <mn>1</mn> </msub> <msub> <mi>sin&amp;phi;</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mrow>
基于下式求解圆弧段两切线交点到目标始点的长度:
<mrow> <msub> <mi>&amp;lambda;</mi> <mn>2</mn> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msub> <mi>T</mi> <mi>s</mi> </msub> <msub> <mi>&amp;lambda;</mi> <mn>1</mn> </msub> <mo>+</mo> <msup> <msub> <mi>&amp;lambda;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </msqrt> <mo>.</mo> </mrow>
5.根据权利要求4所述的一种基于恒工具面角的井眼轨道建模方法,其特征在于,基于下式求解圆弧段弯曲角:
θ=arccos[(ΔXsinα1cosφ1+ΔYsinα1sinφ1+ΔZcosα11)/λ2]
基于下式求解直线段L1的长度:
<mrow> <msub> <mi>L</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>&amp;lambda;</mi> <mn>1</mn> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <msub> <mi>K</mi> <mn>1</mn> </msub> </mfrac> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mfrac> <mi>&amp;theta;</mi> <mn>2</mn> </mfrac> <mo>.</mo> </mrow>
6.根据权利要求1所述的一种基于恒工具面角的井眼轨道建模方法,其特征在于,所述轨道为五段制,并且基于给定变量L1、L2、L3、K1和K2中的任意4个参数求解另一参数;或者,当K1=K2时,基于变量L1、L2、L3,求解K1和K2
7.根据权利要求6所述的一种基于恒工具面角的井眼轨道建模方法,其特征在于,给定L1、L3、K1和K2,求解L2的具体步骤为:
步骤1,以L2为基本求解变量,再选择α2和φ2变量为未知数,作为方程组(1)的3个求解参数;
步骤2,根据设计起始点和目标点的坐标及方向,以及给定L1、L3、K1和K2已知参数,求对应空间圆弧轨道模型的解,并将其结果L2、α2和φ2作为对应恒工具面角模型解的初值;
步骤3,求Kα1、Kα2、S1和S2,确定两曲线段方位角的增减符号;
<mrow> <msub> <mi>K</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mo>&amp;PlusMinus;</mo> <mfrac> <msub> <mi>K</mi> <mn>1</mn> </msub> <mrow> <mn>1</mn> <mo>+</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;phi;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>l</mi> <mi>n</mi> <mfrac> <mrow> <msub> <mi>tan&amp;alpha;</mi> <mn>2</mn> </msub> </mrow> <mrow> <msub> <mi>tan&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>,</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>&gt;</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow>
式中,增斜取“+”,降斜取“-”,稳斜时,Kα1=0,当α1=0时,Kα1=K1
<mrow> <msub> <mi>S</mi> <mn>1</mn> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mo>(</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>/</mo> <msub> <mi>K</mi> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </msub> </mrow> </mtd> <mtd> <mrow> <mo>(</mo> <msub> <mi>K</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> <mo>&amp;NotEqual;</mo> <mn>0</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>|</mo> <msub> <mi>&amp;phi;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;phi;</mi> <mn>1</mn> </msub> <mo>|</mo> <mo>/</mo> <mi>K</mi> </mrow> </mtd> <mtd> <mrow> <mo>(</mo> <msub> <mi>K</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mn>0</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
同理,可求得Kα2和S2
步骤4,用拟牛顿法求解方程组(1);
步骤5,若方程组有解,则计算轨道节点参数,若无解,则提示无解,可调整已知参数或重新确定参数求解方式进行求解。
CN201710342223.XA 2017-05-16 2017-05-16 一种基于恒工具面角的井眼轨道建模方法 Pending CN107153738A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710342223.XA CN107153738A (zh) 2017-05-16 2017-05-16 一种基于恒工具面角的井眼轨道建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710342223.XA CN107153738A (zh) 2017-05-16 2017-05-16 一种基于恒工具面角的井眼轨道建模方法

Publications (1)

Publication Number Publication Date
CN107153738A true CN107153738A (zh) 2017-09-12

Family

ID=59793362

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710342223.XA Pending CN107153738A (zh) 2017-05-16 2017-05-16 一种基于恒工具面角的井眼轨道建模方法

Country Status (1)

Country Link
CN (1) CN107153738A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107762411A (zh) * 2017-12-05 2018-03-06 重庆科技学院 连续管钻井轨道纠偏方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1936264A (zh) * 2005-09-20 2007-03-28 中国石油化工股份有限公司 一种钻井井眼轨道设计方法
CN102393533A (zh) * 2011-07-29 2012-03-28 中国石油化工股份有限公司 一种利用钻井轨迹校正vsp上行转换波的处理方法
US20130140037A1 (en) * 2010-08-24 2013-06-06 Jose J. Sequeira, JR. System and method for planning a well path
CN103699738A (zh) * 2013-12-20 2014-04-02 中国石油大学(华东) 一种控制点约束条件下的侧钻水平井轨道设计方法
CN102425374B (zh) * 2011-11-21 2014-07-23 中国石油集团钻井工程技术研究院 一种井眼轨道控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1936264A (zh) * 2005-09-20 2007-03-28 中国石油化工股份有限公司 一种钻井井眼轨道设计方法
US20130140037A1 (en) * 2010-08-24 2013-06-06 Jose J. Sequeira, JR. System and method for planning a well path
CN102393533A (zh) * 2011-07-29 2012-03-28 中国石油化工股份有限公司 一种利用钻井轨迹校正vsp上行转换波的处理方法
CN102425374B (zh) * 2011-11-21 2014-07-23 中国石油集团钻井工程技术研究院 一种井眼轨道控制方法
CN103699738A (zh) * 2013-12-20 2014-04-02 中国石油大学(华东) 一种控制点约束条件下的侧钻水平井轨道设计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HARRIS J 等: "Designing snakey: a tangible user interface supporting well path planning", 《CONFERENCE ON HUMAN-COMPUTER INTERACTION. SPRINGER, BERLIN, HEIDELBERG》 *
唐雪平 等: "三维井眼轨道设计模型及应用", 《数学的实践与认识》 *
方敏 等: "恒装置角曲线的半解析计算", 《同济大学学报(自然科学版)》 *
李国栋: "恒装置角模式井眼轨迹控制计算方法研究", 《中外能源》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107762411A (zh) * 2017-12-05 2018-03-06 重庆科技学院 连续管钻井轨道纠偏方法
CN107762411B (zh) * 2017-12-05 2019-03-01 重庆科技学院 连续管钻井轨道纠偏方法

Similar Documents

Publication Publication Date Title
CN102425374B (zh) 一种井眼轨道控制方法
CN104453713B (zh) 一种利用旋转导向工具控制井眼轨迹的方法
US8612084B2 (en) System and method for autonomous navigation of a tracked or skid-steer vehicle
CN103883312B (zh) 一种导向钻井入靶形势预测的通用方法
US9932819B2 (en) Method of orienting a second borehole relative to a first borehole
CN105114059A (zh) 一种三维定向井轨迹设计及控制参数判别方法
CN106050223A (zh) 一种全向静磁场救援井随钻定位方法及定位系统
CN103883249A (zh) 一种基于旋转导向钻井的水平井着陆控制方法
CN110147588A (zh) 一种基于钻进时间最短的三维水平井轨道优化设计方法
CN105631753A (zh) 一种基于地层产状的水平井油藏剖面建模方法
CN103699738A (zh) 一种控制点约束条件下的侧钻水平井轨道设计方法
CN103883307B (zh) 一种中靶优先着陆控制的通用导向钻井方法
CN107191177A (zh) 恒工具面角钻进轨迹控制方法
CN107153738A (zh) 一种基于恒工具面角的井眼轨道建模方法
CN106761685B (zh) 利用工具面角识别井眼轨迹模式的方法
CN107869309A (zh) 小井斜情况下的井眼轨迹控制方法及装置
CN107133413A (zh) 恒工具面角曲线井眼轨道设计方法
CN109711023A (zh) 处理三维绕障井的双二维轨迹优化方法
CN111414705A (zh) 一种丛式井井眼轨迹防碰设计的三维空间可穿越评估方法
CN103883306B (zh) 一种基于复合导向钻井的水平井方向优先着陆控制方法
CN106567705B (zh) 一种用工具面角控制井眼轨迹的方法
CN114517687A (zh) 一种非开挖施工智能导向辅助系统
CN103883252A (zh) 一种基于滑动导向钻井的水平井着陆控制方法
CN103883251A (zh) 一种基于旋转导向钻井的水平井方向优先着陆控制方法
Elshafei et al. Optimization of rotary steerable drilling

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170912