CN103883249B - 一种基于旋转导向钻井的水平井着陆控制方法 - Google Patents
一种基于旋转导向钻井的水平井着陆控制方法 Download PDFInfo
- Publication number
- CN103883249B CN103883249B CN201310146371.6A CN201310146371A CN103883249B CN 103883249 B CN103883249 B CN 103883249B CN 201310146371 A CN201310146371 A CN 201310146371A CN 103883249 B CN103883249 B CN 103883249B
- Authority
- CN
- China
- Prior art keywords
- target
- point
- delta
- landing
- coordinate
- 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.)
- Active
Links
- 238000005553 drilling Methods 0.000 title claims abstract description 83
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000013461 design Methods 0.000 claims abstract description 20
- 230000008569 process Effects 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 238000013213 extrapolation Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000008685 targeting Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004179 indigotine Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
Landscapes
- Earth Drilling (AREA)
Abstract
本发明公开了一种基于旋转导向钻井的水平井着陆控制方法,其包括以下步骤:计算井底点的轨迹参数;选取入靶点位置,并计算从井底点到入靶点的坐标增量;设计着陆轨迹的水平投影图,得到着陆轨迹在其上的曲率、入靶方位角以及水平投影长度;设计着陆轨迹的垂直剖面图,得到着陆轨迹在其上的曲率、入靶井斜角和井段长度;校核算得的入靶井斜角和方位角是否满足工程要求;继续优化着陆控制方案从而确定出最优的着陆控制方案;计算着陆轨迹的分点参数,并以图表形式输出设计结果。本发明在优先满足水平井中靶要求的前提下,通过校核入靶方向,采用单一的钻井工艺和技术参数实现着陆入靶,从而用最简单的工艺和最少的工序满足了水平井着陆入靶的轨迹控制要求。
Description
技术领域
本发明涉及石油钻井工程领域,尤其涉及一种水平井旋转导向钻井的着陆控制方法。
背景技术
井眼轨迹控制是一个复杂的多扰动控制过程,要使实钻轨迹与设计轨道完全吻合是不可能的,工程上允许二者之间存在一定的误差。当二者误差较大时,则需要修正设计从当前井底到靶点的井眼轨道。这种修正轨道(也称待钻井眼)设计主要有两种方案:一是中靶控制方案,中靶控制方案只要求击中给定的靶区,而对进入靶区的井斜角和方位角没有严格的限制。该方案的典型井身剖面为“直线段—曲线段—直线段”剖面,最简化的井身剖面为“曲线段—直线段”剖面,如图1所示;二是软着陆控制方案,软着陆控制方案既给定了入靶点的空间位置,也给定了入靶点的井眼方向。该方案的典型井身剖面为“直线段—曲线段—直线段—曲线段—直线段”剖面,最简化的井身剖面为“曲线段—直线段—曲线段”剖面,如图2所示。
现有的井眼轨迹控制技术无论是中靶控制方案还是软着陆控制方案,至少需要2个井段甚至多达5个井段。而每个井段会采用不同的导向钻井方式和工艺技术参数,并且涉及到若干次的起下钻(起下钻次数=井段数-1)。在钻井施工过程中,钻头距离靶区窗口越近,其轨迹控制要求越高。水平井着陆控制的关键阶段往往位于距靶区窗口数十米的范围内,此时不仅要满足着陆入靶要求,还应尽量采用最简单的工艺及工序,减少施工难度,提高井身质量。另外,现有的着陆控制也不涉及靶区窗口(靶平面)、入靶井斜角和方位角的校核等问题,也没有控制方案优化的技术方法。
综上,现有的着陆控制技术存在以下缺点:(1)工艺复杂,需要多个井段才能实现中靶着陆控制;(2)中靶控制方案不涉及靶平面、入靶井斜角和方位角的校核等问题;(3)没有中靶控制方案的优化方法。
发明内容
本发明针对现有技术中的上述问题,提出了一种基于旋转导向钻井的水平井着陆控制方法。该方法包括以下步骤:
S101、根据随钻测量仪获取的实钻轨迹测斜数据,按实际使用的导向钻井工艺,采用外推法计算井底点的轨迹参数,所述轨迹参数包括所述井底点的井斜角、方位角和空间坐标;
S102、在靶平面上选定入靶点的位置并基于所述井底点的轨迹参数计算从所述井底点到所述入靶点的坐标增量,所述入靶点的位置用所述入靶点在所述靶平面的坐标来表示,所述坐标增量为空间坐标增量;
S103、在旋转导向钻井条件下,根据所计算出的井底点与入靶点之间的坐标增量及所述轨迹参数,按入靶位置要求及圆柱螺线模型设计着陆轨迹的水平投影图,得到着陆轨迹在水平投影图上的曲率、入靶方位角以及水平投影长度,其计算公式如下:
其中,
Δh=ΔNb,esinφb-ΔEb,ecosφb
Ch=ΔNb,ecosφb+ΔEb,esinφb
式中,κh为着陆轨迹在水平投影图上的曲率;ΔNb,e、ΔEb,e分别为井底点到入靶点的北坐标和东坐标的坐标增量;φb为井底点的方位角;φe为入靶方位角;ΔSb,e为水平投影长度;
S104、在旋转导向钻井条件下,根据井底点与入靶点之间的坐标增量、所述轨迹参数和步骤S103中所计算出的水平投影长度,按入靶位置要求及圆柱螺线模型设计着陆轨迹的垂直剖面图,得到着陆轨迹在垂直剖面图上的曲率、入靶井斜角和井段长度,其计算公式如下:
其中,
Δv=ΔHb,esinαb-ΔSb,ecosαb
Cv=ΔHb,ecosαb+ΔSb,esinαb
式中,κv为着陆轨迹在垂直剖面图上的曲率;ΔHb,e为井底点到入靶点的垂深坐标的坐标增量;αb为井底点的井斜角;αe为入靶井斜角;ΔLb,e为井段长度;
S105、校核由步骤S103和S104算得的入靶井斜角和方位角是否满足工程要求,如果满足要求,则着陆控制方案可行,执行下面的步骤,否则,重新选取入靶点位置返回到步骤S102,重复执行以获得满足工程要求的入靶方向,所述入靶方向包括入靶井斜角和方位角;
S106、继续优化着陆控制方案,将靶区窗口划分成多个网格单元,分别将每个纵横网格线的交点作为一个入靶点位置,然后采用步骤S102到S104的方法计算各入靶位置对应的入靶方向,按照工程要求从一系列入靶位置、入靶方向中选出最优值,从而确定出最优的着陆控制方案;
S107、根据所述最优的着陆控制方案及工艺技术参数,按圆柱螺线模型计算着陆轨迹的分点参数,并以图表形式输出设计结果,作为水平井着陆控制施工的依据。
根据本发明的一个实施例,在所述步骤S101中,按照以下步骤计算所述井底点的井斜角、方位角及空间坐标:
S201、利用随钻测量仪器获取一系列测点Mi(i=1,2,…,n)的测斜数据,所述测斜数据包括井深、井斜角和方位角;
S202、根据实际钻井工艺条件选用相应的井眼轨迹模型,在滑动导向钻井、旋转导向钻井和复合导向钻井条件下,宜分别选用空间圆弧模型、圆柱螺线模型和自然曲线模型作为井眼轨迹模型;
S203、根据最后两测点的测斜数据计算最后测段的轨迹特征参数,如果实钻轨迹采用滑动导向钻井方式,所述轨迹特征参数为最后测段的曲率半径和测点Mn-1处的工具面角,按以下公式进行计算:
其中,
ΔLn-1,n=Ln-Ln-1
Δφn-1,n=φn-φn-1
cosεn-1,n=cosαn-1cosαn+sinαn-1sinαncosΔφn-1,n
式中,Ln-1和Ln分别是最后两测点Mn-1和Mn的井深,αn-1和αn分别是最后两测点Mn-1和Mn的井斜角,φn-1和φn分别是最后两测点Mn-1和Mn的方位角,εn-1,n是最后测段的弯曲角;如果εn-1,n为零,表明最后测段为直线段,则不需要计算曲率半径和工具面角;
S204、计算井底点的井斜角、方位角及空间坐标,
在滑动导向钻井条件下,按照以下公式计算井底点的井斜角、方位角及空间坐标:
当εn-1,n=0时
当εn-1,n≠0时
其中,
式中,αb和φb分别表示井底点的井斜角和方位角,εn-1,b表示实钻轨迹从测点Mn-1到井底点的弯曲角,Nb、Eb、Hb分别表示井底点的北坐标、东坐标和垂深坐标,Nn-1、En-1、Hn-1分别为最后测点之上的一个测点Mn-1的北坐标、东坐标和垂深坐标,ΔLn,b为最后测点Mn到井底点的井段长度,ΔLn-1,b为最后测点之上的一个测点Mn-1到井底点的井段长度。
根据本发明的另一个实施例,在所述步骤S102中,按照以下步骤选定入靶点位置并计算着陆轨迹的坐标增量:
S301、建立以首靶点为原点的坐标系t-xyz,其中,x轴铅垂向上,y轴水平向右,z轴为靶平面的法线方向;
S302、在靶平面上选取入靶点位置的坐标(xe,ye)并计算其对应的空间坐标,其计算公式如下:
式中,Ne、Ee、He分别为入靶点的北坐标、东坐标和垂深坐标,Nt、Et、Ht分别为设定的首靶点的北坐标、东坐标和垂深坐标,φz为靶平面的法线方位,xe和ye为入靶点在靶平面上的坐标;
S303、根据所计算的井底点和入靶点的空间坐标,计算从井底点到入靶点的空间坐标增量,其公式为:
本发明带来了以下有益效果:
(1)结合旋转导向钻井的工艺技术特点,在优先满足水平井中靶要求的前提下,通过校核入靶方向,采用单一的钻井工艺技术参数实现着陆中靶,从而用最简单的工艺和最少的工序(最少的起下钻次数)来满足水平井着陆中靶的轨迹控制要求,技术方案简明、实用性强。
(2)给出了滑动导向钻井条件下当前钻头位置和井眼方向的计算方法,弥补了实钻轨迹监测计算与着陆控制方案设计之间的重要环节,提高了科学性和实用性。
(3)通过建立靶平面方程,将着陆控制方案与靶区有机地结合起来,提出了包括入靶位置网格细化、入靶方向校核等内容的优化方法,从而能更好地设计着陆控制方案。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
图1是现有技术中靶控制方案的典型井身剖面;
图2是现有技术软着陆控制方案的典型井身剖面;
图3是本发明的技术原理示意图;
图4是本发明的着陆轨迹控制方法流程图;
图5是本发明的计算井底点轨迹参数的流程图;
图6是本发明的计算着陆轨迹坐标增量的流程图;
图7是本发明的优化着陆控制方案的网格划分示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
另外,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图3显示了本发明的技术原理示意图。钻井过程中,设计轨道往往要求通过首靶点t,实钻轨迹已钻达井底点b(当前的钻头位置)。而着陆轨迹是从井底点b开始钻达入靶点e的待钻轨迹,因此,着陆控制方案就是要设计出着陆轨迹以及钻井工艺技术参数。
着陆轨迹既可以继续采用当前的导向钻井方式而只改变工艺技术参数,也可以更换导向钻井方式并重新设计工艺技术参数。换句话说,井底点b以上的实钻轨迹与井底点b以下的着陆轨迹既可以采用相同的导向钻井方式,也可以采用不同的导向钻井方式。为叙述方便且不失一般性,在本具体实施方式中假设实钻轨迹采用了滑动导向钻井技术,其井眼轨迹符合空间圆弧模型;着陆轨迹将采用旋转导向钻井技术,其井眼轨迹符合圆柱螺线模型。因此,着陆轨迹是一条等变螺旋角的圆柱螺线,其垂直剖面图和水平投影图都是圆弧,如图3所示。对于实钻轨迹和着陆轨迹都采用旋转导向钻井技术等其它情况,在本发明的技术原理和方法基础上,本领域内的技术人员不难做出相应的改进或变型,因此本发明并非局限于以下具体的实施内容。
实施例一:
图4显示了本发明的着陆轨迹控制方法流程图。
在步骤S101中,计算井底点b的轨迹参数。图5显示了本发明的计算井底点b的轨迹参数的流程图。在一个实施例中,可以按照以下步骤计算井底点b的轨迹参数:
步骤S201,利用随钻测量仪器获取一系列测点Mi(i=1,2,…,n)的测斜数据,该测斜数据包括井深Li、井斜角αi和方位角φi。这里,随钻测量仪器可以选择MWD等仪器。
步骤S202,根据实际钻井工艺条件选用相应的井眼轨迹模型。
计算井底点b的轨迹参数应依据实钻轨迹所采用的导向钻井方式来选用井眼轨迹模型。在滑动导向钻井、旋转导向钻井和复合导向钻井条件下,宜分别选用空间圆弧模型、圆柱螺线模型和自然曲线模型作为井眼轨迹模型。本实施例给出了滑动导向钻井条件下的井底点轨迹参数计算方法,对于旋转导向钻井、复合导向钻井等其它钻井方式,在本发明的技术原理和方法基础上,本领域内的技术人员不难做出相应的改进或变型,所以本发明的保护范围并不限于滑动导向钻井方式。
步骤S203,根据最后两测点Mn-1和Mn的测斜数据计算最后测段[Ln-1、Ln]的轨迹特征参数。
在滑动导向钻井条件下,该轨迹特征参数为最后测段[Ln-1、Ln]的曲率半径R和测点Mn-1处的工具面角ω,可按以下公式进行计算:
其中,ΔLn-1,n=Ln-Ln-1(2)
Δφn-1,n=φn-φn-1(3)
cosεn-1,n=cosαn-1cosαn+sinαn-1sinαncosΔφn-1,n(4)
式中,Ln-1和Ln分别是最后两测点Mn-1和Mn的井深,单位均为m;αn-1和αn分别是最后两测点Mn-1和Mn的井斜角,单位均为(°);φn-1和φn分别是最后两测点Mn-1和Mn的方位角,单位均为(°);εn-1,n是最后测段的弯曲角。如果εn-1,n为零,表明最后测段为直线段,则不需要计算曲率半径和工具面角。
步骤S204、计算井底点b的井斜角、方位角及空间坐标等轨迹参数。
在滑动导向钻井条件下,按照以下公式计算井底点b的井斜角、方位角及空间坐标:
当εn-1,n=0时
当εn-1,n≠0时
其中
式中,αb和φb分别表示井底点b的井斜角和方位角,单位均为(°);εn-1,b表示实钻轨迹从测点Mn-1到井底点b的弯曲角;Nb、Eb、Hb分别表示井底点b的北坐标、东坐标和垂深坐标,单位均为m;Nn-1、En-1、Hn-1分别为最后测点之上的一个测点Mn-1的北坐标、东坐标和垂深坐标,单位均为m;ΔLn,b为最后测点Mn到井底点b的井段长度,单位为m,ΔLn-1,b为最后测点之上的一个测点Mn-1到井底点b的井段长度。
如图4所示,在步骤S102中,选定入靶点e的位置并计算着陆轨迹的坐标增量。该坐标增量为从井底点b到入靶点e的空间坐标增量。图6为计算从井底点b到入靶点e坐标增量的流程图。在一个实施例中,按照以下步骤和公式计算着陆轨迹的坐标增量:
步骤S301,建立以首靶点t为原点的坐标系t-xyz。其中,x轴铅垂向上,y轴水平向右,z轴为靶平面的法线方向。
步骤S302,在靶平面上选取入靶点e的位置坐标并计算其对应的空间坐标。在靶平面上,选定入靶点坐标(xe,ye),计算入靶点e的空间坐标的公式如下:
式中,Ne、Ee、He分别为入靶点e的北坐标、东坐标和垂深坐标,单位均为m;Nt、Et、Ht分别为设定的首靶点t的北坐标、东坐标和垂深坐标,单位均为m;φz为靶平面的法线方位,单位为(°)。
步骤S303,计算着陆轨迹的坐标增量。根据所计算的井底点b和入靶点e的空间坐标,计算从井底点b到入靶点e的空间坐标增量,其公式为:
如图4所示,步骤S103中,按入靶位置要求设计着陆轨迹的水平投影图。在旋转导向钻井条件下,根据所计算出的井底点b与入靶点e之间的坐标增量及井底点b的轨迹参数,按圆柱螺线模型设计着陆轨迹的水平投影图,得到着陆轨迹在水平投影图上的曲率、水平投影长度以及入靶方位角等参数,其计算公式如下:
着陆轨迹在水平投影图上的曲率:
其中
Δh=ΔNb,esinφb-ΔEb,ecosφb(15)
式中,ΔNb,e、ΔEb,e为井底点b到入靶点e的北坐标和东坐标的坐标增量,单位均为m;φb表示井底点b的方位角,单位为(°);κh为井眼轨迹在水平投影图上的曲率,单位是(°)/m。
着陆轨迹的入靶方位角:
其中,
Ch=ΔNb,ecosφb+ΔEb,esinφb(17)
式中,φe为着陆轨迹的入靶方位角,单位是(°)。
着陆轨迹的水平投影长度:
式中,ΔSb,e为着陆轨迹的水平投影长度。
步骤S104中,按入靶位置要求设计着陆轨迹的垂直剖面图。在旋转导向钻井条件下,根据井底点b与入靶点e之间的坐标增量、井底点b的轨迹参数和步骤S103中所计算出的着陆轨迹参数,按圆柱螺线模型设计着陆轨迹的垂直剖面图,得到着陆轨迹在垂直剖面图上的曲率、入靶井斜角、井段长度等参数,其计算公式如下:
着陆轨迹在垂直剖面图上的曲率:
其中
Δv=ΔHb,esinαb-ΔSb,ecosαb(20)
式中,ΔHb,e为井底点b到入靶点e的垂深坐标的坐标增量,单位为m;αb为井底点b的井斜角,单位是(°);κv为井眼轨迹在垂直剖面图上的曲率,单位是(°)/m。
着陆轨迹的入靶井斜角:
其中
Cv=ΔHb,ecosαb+ΔSb,esinαb(22)
式中,αe为着陆轨迹的入靶井斜角,单位是(°)。
着陆轨迹的井段长度:
式中,ΔLb,e为着陆轨迹的井段长度。
步骤S105中,校核入靶方向是否满足要求。校核由公式(21)和(16)算得的入靶井斜角αe和入靶方位角φe是否满足工程要求。如果满足要求,则着陆控制方案可行,执行后续的设计工作;否则,重新选取入靶点位置返回到步骤S102,执行以获得满足工程要求的入靶方向。其中,入靶方向包括入靶井斜角和方位角。
多数情况下着陆控制方案的入靶位置和方向与设计的入靶位置和方向越接近越好。然而,其衡量标准是一种综合指标,至于什么样的方案最优应根据实际的工程情况来确定。例如,当选取的入靶位置与设计轨道相同时,可能井眼方向就会相差很大,此时就未必是个好方案。再如,若入靶位置偏左,但入靶方向偏右,即使与设计轨道的入靶方向相差较大,也可能是一种好方案,等等。正是由于这些原因,提出了下述优化着陆控制方案的步骤S106。
步骤S106,继续优化着陆控制方案。当完成步骤S105之后,就得到了一个满足入靶位置和入靶方向要求的着陆控制方案,但是不一定是最优方案。为了得到最优的着陆控制方案,可以将靶区窗口(靶平面)划分成多个网格单元,将每个纵横网格线的交点作为一个入靶位置,如图7所示。然后,用上述方法求得相应的入靶井斜角、入靶方位角以及着陆轨迹的其它参数,进而可从中优选出最优的着陆控制方案。为了减少计算量,可先用间距较大的纵横网格,然后选取较优方案的区域,进一步细划网格线,不断优化着陆控制方案,直到满足优选控制方案的间距要求为止。通过这种循环往复的优化过程,保证能设计出最优的着陆控制方案。
步骤S107中,输出设计结果。当完成步骤S106之后,就确定了最优的着陆控制方案及工艺技术参数,包括着陆轨迹在垂直剖面图和水平投影图上的曲率、入靶方向、水平投影长度和井段长度等参数。为了便于具体实施该控制方案,还需要按圆柱螺线模型计算出着陆轨迹上各分点的轨迹参数,并以图表等形式来输出着陆控制方案的设计结果,作为水平井着陆控制施工的依据。
由以上实施例得,在本发明的实施过程中,在优先满足水平井中靶要求的前提下,通过校核入靶方向,采用单一的钻井工艺技术参数实现着陆入靶,从而用最简单的工艺和最少的工序(最少的起下钻次数)满足了水平井着陆入靶的轨迹控制要求。
本发明给出了滑动导向钻井条件下当前钻头位置和井眼方向的计算方法,弥补了实钻轨迹监测计算与着陆控制方案设计之间的重要环节,提高了科学性和实用性。
另外,本发明通过建立靶平面方程,将着陆控制方案与靶区有机地结合起来,通过将靶区窗口划分成若干个网格单元并逐步细化的方法,可得到最优的着陆控制方案。
实施例二
下面以某实际水平井为例来具体说明按照本发明的技术原理和步骤如何设计着陆控制方案。
某水平井设计轨道的节点数据见表1,其中首靶点t的坐标及靶区参数为:首靶点t的垂深Ht=1500m、水平位移At=280m、平移方位及靶平面法线方位角靶窗宽度wt=20m、宽度ht=6m。当进入着陆井段后,采用滑动导向钻井工艺,钻至井深L132=1557m时(测点编号132),井斜角α132=65.5°、方位角φ132=63.2°、北坐标N132=94.36m、东坐标E132=172.72m、垂深H132=1480.53m。继续钻进至L133=1567m时(测点编号133),测得α133=67.86°、φ133=60.75°,且钻头距测点ΔLn,b=16m。现改用旋转导向钻井工艺继续钻进,试设计着陆控制方案。
表1某水平井设计轨道节点数据
根据本发明的技术方案,设计着陆控制方案包含以下步骤:
滑动导向钻井所钻出的井眼轨迹更符合空间圆弧模型,其轨迹特征参数为曲率半径(或井眼曲率)和工具面角。对于最后测段[1557m、1567m],首先,用公式(1)~(4)计算曲率半径R和测点M132处的工具面角ω:
ΔL132,133=1567-1557=10m
Δφ132,133=60.75-63.20=-2.45°
接着,用公式(5)~(11),计算井底点b的井斜角、方位角及空间坐标:
然后,根据选定的入靶点位置,由公式(12)和(13)计算着陆轨迹的坐标增量。由表1知:首靶点t坐标为(140.00,242.49,1500.00)。在靶平面上,若选取入靶点e坐标为(0.5,3.0),则其空间坐标为
那么,从井底点b到入靶点e的坐标增量为
旋转导向钻井所钻出的井眼轨迹更符合圆柱螺线模型。按照圆柱螺线模型,由公式(14)和(15)计算着陆轨迹的在水平投影图上的曲率:
Δh=30.96×sin56.99°-50.32×cos56.99°≠0
接着,由公式(16)和(17)计算着陆轨迹的入靶方位角:
Ch=30.96×cos56.99°+50.32×sin56.99°=59.07m
Ch+ΔNb,e=59.06+30.96≠0
由公式(18)计算着陆轨迹的水平投影长度:
按照圆柱螺线模型,由公式(19)和(20)计算着陆轨迹的在垂直剖面图上的曲率:
Δv=9.48×sin71.71°-59.09×cos71.71°≠0
由公式(21)和(22)计算着陆轨迹的入靶点井斜角:
Cv=9.48×cos71.71°+59.09×sin71.71°=59.08m
Cv+ΔHb,e=59.08+9.48≠0
最后,由公式(23)计算着陆轨迹的井段长度:
由上述结果知:αe-αt=0.06°,φe-φt=-0.20°,即入靶井斜角和方位角与设计轨道的首靶井斜角和方位角符合很好,因此该着陆控制方案可行。
因此,在该实施例中,若选取入靶点e的靶点坐标系下坐标为(0.5,3.0),则在旋转导向钻井条件下着陆控制方案的主要工艺参数为:着陆轨迹在垂直剖面图上的曲率为9.16°/30m,在水平投影图上的曲率为1.42°/30m。该着陆控制方案的节点数据见表2。
表2着陆控制方案的节点数据
在该实施例的实施过程中,为了得到最优的控制方案,可以用纵横网格线将靶区划分成若干个网格单元,每个纵横网格线的交点作为一个入靶点位置,重复上述方法和步骤,可得到:在-3≤xe≤3、-10≤ye≤10的靶区(即整个靶区窗口)范围内,取纵横网格线间距均为1m,则入靶井斜角和方位角的设计结果见表3和表4。
表3整个靶区窗口的入靶井斜角数据
表4整个靶区窗口的入靶方位角数据
由上述结果知:若要求|φe-φt|≤1.5°,则2≤ye≤4(见表4阴影部分);进而再要求|αe-αt|≤1.0°,则-1≤xe≤1(见表3阴影部分)。在-1≤xe≤1、2≤ye≤4的靶区(即部分靶区窗口)范围内,若横向网格线间距取为0.25m,纵向网格线间距取为0.20m,则入靶井斜角和方位角的设计结果见表5和表6。
表5网格细化后的部分靶区窗口入靶井斜角数据
表6网格细化后的部分靶区窗口入靶方位角数据
按上述方法和步骤,可逐步细化,最后得到最优入靶点e。如果只要求入靶点e位置位于靶区内,而期望入靶井眼方向(αe,φe)与设计入靶井眼方向(αt,φt)相同,则应xe=0.47m、ye=3.10m。此时,着陆轨迹在垂直剖面图和水平投影图上的曲率分别为9.13°/30m和1.52°/30m,入靶点参数为:Le=1643.10m,αe=90.00°,φe=59.99°,Ne=137.32m,Ee=244.04m,He=1499.53m。
显然,入靶点位置和入靶井眼方向的最优点往往并不重合。换句话说,当xe=ye=0时,一般难以满足αe=αt、φe=φt。然而,在一定的入靶点位置(xe,ye)和入靶方向(αe,φe)允许范围内,本发明可以设计出满足要求的着陆控制方案,并且可逐步优化方案。
虽然本发明所揭露的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (3)
1.一种基于旋转导向钻井的水平井着陆控制方法,其特征在于,包括以下步骤:
S101、根据随钻测量仪获取的实钻轨迹测斜数据,按实际使用的导向钻井工艺,采用外推法计算井底点(b)的轨迹参数,所述轨迹参数包括所述井底点(b)的井斜角、方位角和空间坐标;
S102、在靶平面上选定入靶点(e)的位置并基于所述井底点(b)的轨迹参数计算从所述井底点(b)到所述入靶点(e)的坐标增量,所述入靶点(e)的位置用所述入靶点(e)在所述靶平面的坐标来表示,所述坐标增量为空间坐标增量;
S103、在旋转导向钻井条件下,根据所计算出的井底点(b)与入靶点(e)之间的坐标增量及所述轨迹参数,按入靶位置要求及圆柱螺线模型设计着陆轨迹的水平投影图,得到着陆轨迹在水平投影图上的曲率、入靶方位角以及水平投影长度,其计算公式如下:
其中,
Δh=ΔNb,esinφb-ΔEb,ecosφb
Ch=ΔNb,ecosφb+ΔEb,esinφb
式中,κh为着陆轨迹在水平投影图上的曲率;ΔNb,e、ΔEb,e分别为井底点(b)到入靶点(e)的北坐标和东坐标的坐标增量;φb为井底点(b)的方位角;φe为入靶方位角;ΔSb,e为水平投影长度;
S104、在旋转导向钻井条件下,根据井底点(b)与入靶点(e)之间的坐标增量、所述轨迹参数和步骤S103中所计算出的水平投影长度,按入靶位置要求及圆柱螺线模型设计着陆轨迹的垂直剖面图,得到着陆轨迹在垂直剖面图上的曲率、入靶井斜角和井段长度,其计算公式如下:
其中,
Δv=ΔHb,esinαb-ΔSb,ecosαb
Cv=ΔHb,ecosαb+ΔSb,esinαb
式中,κv为着陆轨迹在垂直剖面图上的曲率;ΔHb,e为井底点(b)到入靶点(e)的垂深坐标的坐标增量;αb为井底点(b)的井斜角;αe为入靶井斜角;ΔLb,e为井段长度;
S105、校核由步骤S103和S104算得的入靶井斜角和方位角是否满足工程要求,如果满足要求,则着陆控制方案可行,执行下面的步骤,否则,重新选取入靶点位置返回到步骤S102,重复执行以获得满足工程要求的入靶方向,所述入靶方向包括入靶井斜角和方位角;
S106、继续优化着陆控制方案,将靶区窗口划分成多个网格单元,分别将每个纵横网格线的交点作为一个入靶点位置,然后采用步骤S102到S104的方法计算各入靶位置对应的入靶方向,按照工程要求从一系列入靶位置、入靶方向中选出最优值,从而确定出最优的着陆控制方案;
S107、根据所述最优的着陆控制方案及工艺技术参数,按圆柱螺线模型计算着陆轨迹的分点参数,并以图表形式输出设计结果,作为水平井着陆控制施工的依据。
2.如权利要求1所述的方法,其特征在于,在所述步骤S101中,按照以下步骤计算所述井底点(b)的井斜角、方位角及空间坐标:
S201、利用随钻测量仪器获取一系列测点Mi(i=1,2,…,n)的测斜数据,所述测斜数据包括井深、井斜角和方位角;
S202、根据实际钻井工艺条件选用相应的井眼轨迹模型,在滑动导向钻井、旋转导向钻井和复合导向钻井条件下,宜分别选用空间圆弧模型、圆柱螺线模型和自然曲线模型作为井眼轨迹模型;
S203、根据最后两测点(Mn-1)和(Mn)的测斜数据计算最后测段的轨迹特征参数,如果实钻轨迹采用滑动导向钻井方式,所述轨迹特征参数为最后测段的曲率半径(R)和最后测点之上的一个测点(Mn-1)处的工具面角(ω),按以下公式进行计算:
其中,
ΔLn-1,n=Ln-Ln-1
Δφn-1,n=φn-φn-1
cosεn-1,n=cosαn-1cosαn+sinαn-1sinαncosΔφn-1,n
式中,Ln-1和Ln分别是最后两测点(Mn-1)和(Mn)的井深,αn-1和αn分别是最后两测点(Mn-1)和(Mn)的井斜角,φn-1和φn分别是最后两测点(Mn-1)和(Mn)的方位角,εn-1,n是最后测段的弯曲角;如果εn-1,n为零,表明最后测段为直线段,则不需要计算曲率半径和工具面角;
S204、计算井底点(b)的井斜角、方位角及空间坐标,
在滑动导向钻井条件下,按照以下公式计算井底点(b)的井斜角、方位角及空间坐标:
当εn-1,n=0时
当εn-1,n≠0时
其中,
式中,αb和φb分别表示井底点(b)的井斜角和方位角,εn-1,b表示实钻轨迹从最后测点之上的一个测点Mn-1到井底点(b)的弯曲角,Nb、Eb、Hb分别表示井底点(b)的北坐标、东坐标和垂深坐标,Nn-1、En-1、Hn-1分别为最后测点之上的一个测点Mn-1的北坐标、东坐标和垂深坐标,ΔLn,b为最后测点Mn到井底点(b)的井段长度,ΔLn-1,b为最后测点之上的一个测点Mn-1到井底点(b)的井段长度。
3.如权利要求1所述的方法,其特征在于,在所述步骤S102中,按照以下步骤选定入靶点(e)位置并计算着陆轨迹的坐标增量:
S301、建立以首靶点(t)为原点的坐标系t-xyz,其中,x轴铅垂向上,y轴水平向右,z轴为靶平面的法线方向;
S302、在靶平面上选取入靶点(e)位置的坐标(xe,ye)并计算其对应的空间坐标,其计算公式如下:
式中,Ne、Ee、He分别为入靶点(e)的北坐标、东坐标和垂深坐标,Nt、Et、Ht分别为设定的首靶点(t)的北坐标、东坐标和垂深坐标,φz为靶平面的法线方位,xe和ye为入靶点(e)在靶平面上的坐标;
S303、根据所计算的井底点(b)和入靶点(e)的空间坐标,计算从井底点(b)到入靶点(e)的空间坐标增量,其公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310146371.6A CN103883249B (zh) | 2013-04-24 | 2013-04-24 | 一种基于旋转导向钻井的水平井着陆控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310146371.6A CN103883249B (zh) | 2013-04-24 | 2013-04-24 | 一种基于旋转导向钻井的水平井着陆控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103883249A CN103883249A (zh) | 2014-06-25 |
CN103883249B true CN103883249B (zh) | 2016-03-02 |
Family
ID=50952344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310146371.6A Active CN103883249B (zh) | 2013-04-24 | 2013-04-24 | 一种基于旋转导向钻井的水平井着陆控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103883249B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104453714B (zh) * | 2014-11-02 | 2016-04-13 | 中国石油集团钻井工程技术研究院 | 一种旋转导向工具的控制方法 |
CN106050215A (zh) * | 2015-04-03 | 2016-10-26 | 普拉德研究及开发股份有限公司 | 目标属性的直接控制 |
US11319780B2 (en) | 2015-04-03 | 2022-05-03 | Schlumberger Technology Corporation | Wellsite system services |
CN105114059B (zh) * | 2015-08-03 | 2017-10-13 | 中国地质大学(武汉) | 一种三维定向井轨迹设计及控制参数判别方法 |
CN106761685B (zh) * | 2015-11-23 | 2022-05-06 | 中国石油化工股份有限公司 | 利用工具面角识别井眼轨迹模式的方法 |
CN106869790B (zh) * | 2017-02-24 | 2018-12-11 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种页岩气水平井快速精细地质导向方法 |
CN108798644A (zh) * | 2018-06-04 | 2018-11-13 | 北京六合伟业科技股份有限公司 | 一种钻测测斜数据qc校核方法 |
CN108894768B (zh) * | 2018-06-25 | 2021-05-14 | 中国地质大学(武汉) | 一种基于蝙蝠算法和井壁稳定的钻进轨迹设计方法与系统 |
CN114746841A (zh) | 2019-10-28 | 2022-07-12 | 吉奥奎斯特系统公司 | 钻井活动建议系统和方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5141060A (en) * | 1989-11-01 | 1992-08-25 | Teleco Oilfield Services Inc. | Method for optimizing of stabilizer positioning in a bottomhole assembly to eliminate the effects of borehole inclination |
EP0816627A2 (en) * | 1996-07-03 | 1998-01-07 | Kubota Corporation | Underground drilling method |
US6523623B1 (en) * | 2001-05-30 | 2003-02-25 | Validus International Company, Llc | Method and apparatus for determining drilling paths to directional targets |
CN1508380A (zh) * | 2002-12-18 | 2004-06-30 | 杜晓瑞 | 超短曲率水平井钻井和完井的方法及其装置 |
RU2270907C1 (ru) * | 2004-07-26 | 2006-02-27 | ОАО НПО "Буровая техника" | Способ проводки направленной скважины по плавной траектории |
CN101929310A (zh) * | 2010-07-21 | 2010-12-29 | 北京航空航天大学 | 一种基于有源交变磁场信息的钻井轨道引导定位方法 |
-
2013
- 2013-04-24 CN CN201310146371.6A patent/CN103883249B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5141060A (en) * | 1989-11-01 | 1992-08-25 | Teleco Oilfield Services Inc. | Method for optimizing of stabilizer positioning in a bottomhole assembly to eliminate the effects of borehole inclination |
EP0816627A2 (en) * | 1996-07-03 | 1998-01-07 | Kubota Corporation | Underground drilling method |
US6523623B1 (en) * | 2001-05-30 | 2003-02-25 | Validus International Company, Llc | Method and apparatus for determining drilling paths to directional targets |
CN1508380A (zh) * | 2002-12-18 | 2004-06-30 | 杜晓瑞 | 超短曲率水平井钻井和完井的方法及其装置 |
RU2270907C1 (ru) * | 2004-07-26 | 2006-02-27 | ОАО НПО "Буровая техника" | Способ проводки направленной скважины по плавной траектории |
CN101929310A (zh) * | 2010-07-21 | 2010-12-29 | 北京航空航天大学 | 一种基于有源交变磁场信息的钻井轨道引导定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103883249A (zh) | 2014-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103883249B (zh) | 一种基于旋转导向钻井的水平井着陆控制方法 | |
CN103883250B (zh) | 一种基于滑动导向钻井的水平井方向优先着陆控制方法 | |
CN103883255B (zh) | 一种基于连续导向钻井的水平井着陆轨迹控制方法 | |
CN105114059B (zh) | 一种三维定向井轨迹设计及控制参数判别方法 | |
CN104481398B (zh) | 一种三维水平井钻井施工方法 | |
CN103883312B (zh) | 一种导向钻井入靶形势预测的通用方法 | |
CN108571287B (zh) | 基于过程控制的井眼轨迹控制系统 | |
CN102900366A (zh) | 一种水平定向钻自由轨迹规划及纠偏方法 | |
CN103758455B (zh) | 一种利用造斜工具钻井的方法及装置 | |
CN107869346B (zh) | 基于导向钻具工作特性的井眼轨迹预测方法 | |
CN103883307B (zh) | 一种中靶优先着陆控制的通用导向钻井方法 | |
CN111119855A (zh) | 基于地层和钻具特性的井眼轨迹预测方法 | |
CN106567675B (zh) | 旋转导向钻井的井眼轨迹控制方法 | |
CN103883251B (zh) | 一种基于旋转导向钻井的水平井方向优先着陆控制方法 | |
CN103883252B (zh) | 一种基于滑动导向钻井的水平井着陆控制方法 | |
CN103883306B (zh) | 一种基于复合导向钻井的水平井方向优先着陆控制方法 | |
CN106761685B (zh) | 利用工具面角识别井眼轨迹模式的方法 | |
CN103883253B (zh) | 一种基于复合导向钻井的水平井着陆控制方法 | |
CN109711023B (zh) | 处理三维绕障井的双二维轨迹优化方法 | |
CN107882505B (zh) | 基于导向钻具工作特性的井眼轨迹控制方法 | |
CN103883254B (zh) | 一种基于导向钻井方向优先着陆控制的通用方法 | |
CN105089643B (zh) | 一种邻井距离扫描的解析方法 | |
CN111119735A (zh) | 基于地层特性和井眼轨迹的导向钻井控制方法 | |
CN105332693A (zh) | 一种钻头水平偏移轨迹获取方法 | |
CN104153764A (zh) | 一种气体钻井裸眼段井筒形状确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |