用于进行井下测量的系统和方法
【专利摘要】进行井下测量的系统和方法利用：中子源(12)，其用于在地下地层(12)中形成标记(122)；核辐射检测器(D、124)，其用于检测所述标记(122)发射的伽玛和/或其他核辐射；以及井下钻井过程监测单元(20)，其基于所述核辐射和所述中子源(120)和所述核辐射检测器(D、124)之间的轴向距离L的检测，生成关于自动化或其他钻井组件(4)的增量深度和/或关联机械钻速(ROP)的实时信息，而不需要借助复杂的有线或无线数据传输系统将实时增量深度和/或ROP数据传输到地表设备。
【专利说明】
用于进行井下测量的系统和方法
技术领域
[0001] 本发明涉及在在井下钻孔中进行井下测量的系统和方法。钻孔可涉及勘探或生产 诸如原油和/或天然气的烃流体。
【背景技术】
[0002] -般，使用可旋转钻柱，钻出用于生产来自井下地层的原油和/或天然气的钻孔。 钻柱通常包括一系列互连的钻杆片段。钻柱的井下端可通常设置有底部钻具组合(BHA)， BHA包括提供钻压的较重钻环的片段、随钻测量(MWD)设备、用于破碎地层在井下端的钻头。 钻机的用于保持钻柱的表面设置有驱动系统，驱动系统用于旋转钻柱，通常包括顶部驱动 或其他旋转台。
[0003] 实现实时地下导航的关键障碍之一在于地表和井下之间的通信瓶颈。当前可用的 随钻测量工具和随钻测井工具可测量除了深度之外的所有至关重要的井下信息。深度估计 井下的实时可用性可提供钻井操作中的实时自动化位舵和优化的新可能性。
[0004] 为了进行地下导航，最重要的测量值之一是任何时间点的钻孔深度。伴随着方位 角和旋转位移的方向测量，它提供钻头相对于表面位置的定位。在下面优化的钻孔轨迹中， 存在多种深度测量的使用，诸如，用于在地层中进行地理特征定位、用于遵循优化的钻孔轨 迹、用于计算套管鞋深度、用于估计水泥数量等。
[0005] 在深度测量的当前系统中，地表系统可记录钻台下方的钻柱的时间和长度。可使 用钻柱长度作为标准深度测量。可使用方钻杆补心或钻台作为基于陆地的钻机的基准和海 上钻机的平均海平面。使用钻柱(BHA和钻杆片段的组合长度)至顶部驱动器(或游车）的距 离和顶部驱动器(或游车)在井架中的位置来确定钻头的深度和机械钻速(R0P)。
[0006] 通过从绞车放出的钻井钢丝绳来测量游车的移动，这是用绞车旋转进行校准的或 者在钻速记录仪的帮助下进行测量的。可使用海上的升沉补偿器来消除浮式海上设施中的 升沉效应。尽管为了准确进行深度测量费尽全力，但当前系统由于与热膨胀、钻杆伸展、压 力效应、和/或钻压误差、钻柱传感器校准和升沉校正而出现误差。由于这些因素导致的深 度总误差可以是3000m深度内高达10-12m。已经努力计算与以上提到的现象相关的误差并 且使用此作为确定校正深度的调整，但它们没有将误差准确定量并且极少实际运用。深度 的井下测量有可能消除这些误差。
[0007] 美国专利5,341，886公开了用于控制旋转钻机的前进方向的方法和设备。该系统 包括钻柱、携载在钻柱上的可旋转钻头、可旋转头和有助于钻井方向变化的顺从性子组件。 该系统包括磁性标记组件，包括地层磁化器。地层中的磁性标记是通过磁化器用对应的电 流脉冲形成的。该工具还包括与磁化器相隔距离L的磁性传感器。当磁性传感器检测到磁性 标记时，促使磁化器形成另一个磁性标记。新磁性标记均分隔开距离L，从而允许测量增量 深度。
[0008] 实际上，美国专利5,341，886的系统证实有致使其不适于应用的缺陷。地层常常是 不可磁化的，从而在深度测量时产生显著误差。美国专利5,341，886指示了使用高强度磁性 脉冲的必要性，大约杆面有几千奥斯特。除了大量功率需求之外，高强度磁性脉冲会与工具 杆中的可磁化材料发生干扰。
[0009]美国专利7,283,910公开了用于记录地层并且获取地下信息以得到所关注参数的 方法和设备，这些参数可包括密度、孔隙率、声学反射率、核磁共振性质、或电阻。用多个传 感器获取这些参数。确定来自各个传感器的信号之间的时间间隔值。使用已知的传感器空 间间隔和时间间隔值，确定钻孔速率并且定义地下特征的递增深度。
[0010]美国专利7,283,910的设备发出核能(更特别地，伽玛射线）。伽玛射线源与被屏蔽 以不受射线源影响的两个或更多个伽玛射线检测器组合。在探针的操作期间，从射线源发 射的伽玛射线(或光子)进入地层并且通过光电吸收，通过康普顿散射，或者通过电子偶的 产生与地层材料的原子电子相互作用。传感器可包括核磁共振传感器。
[0011]在钻井的同时，两个或更多个传感器根据它们的间隔和机械钻速在不同时间经过 地层中的同一位置。可使用算法比较这些传感器的实时输出，以将源自同一地下地层的信 号的特征相关。因为这些信号对应于同一地下位置，所以可使用传感器之间的已知距离L来 计算沿着井下的前进以及机械钻速。即，钻柱在时间差a t = (t2-tl)期间沿着井下前进距 离L。机械钻速ROP = L/(t2-t 1)。通过对R0P求积分，可计算增量深度。
[0012]美国专利7,283,910的设备面临几个挑战。例如，在水平片段中，其中，钻孔通常在 地层内延伸，并且在同一岩性中，特性岩石属性特征不会大幅度改变以致足以区分特征和 噪声。另外，日志决不精确重复。连续伽玛射线测量的微量变化经常是由于到达检测器或传 感器的放射性脉冲的随机性质而导致的统计波动。例如，伽玛射线工具的准确的一般是5% 左右并且精确度与记录速度的平方根成反比并且将受到瞬时速度变化的影响。例如由于冲 蚀导致的流体侵蚀和孔径变化还可使测得的值随时间推移而变化。因 R0P的积分，在计算出 的增量深度中引入积分误差，该误差将随深度增大而增大。误差的大小取决于例如R0P变 化、传感器和测量频率之间的距离。最后，记录工具的精确度随着温度升高而减小，这样限 制了尚温应用。
[0013] 美国专利7,999,220公开了一种脉动中子源和用于钻孔记录的伽玛射线检测器的 组件。该检测器组件包括溴化镧(LaBr3)闪烁晶体和与晶体协作的数字化谱仪和与晶体协 作以使脉冲处理吞吐量最大。该组件可应用于钻孔记录方法，钻孔记录方法在严酷的钻孔 条件下使用伽玛射线进行测量。该系统尤其可应用于碳/氧记录。
[0014] 欧洲专利申请EP2615477公开了用于测量地层压裂中的压裂液的方位分布。工具 具有中子源和绕着工具的外周分隔的多个检测器。检测器彼此屏蔽，使得各检测器检测来 自它最靠近的钻孔区域和地层的伽玛射线。为了用该工具采集记录，中子源将高能种子发 送到周围地层中。中子由于散射而快速失去能量，此后，它们被周围环境内的各种原子吸 附。被散射和吸收的中子发射具有特征能量的伽玛射线。可测量这些伽玛射线与特征能量 的关系并且可确定特定物质的存在与否。
[0015]国际专利申请W02006/004740公开了一种井下脉动中子发射和检测技术，该技术 用于确定钻孔中流体的流动速率。
[0016]美国专利申请U2008/251710公开了一种电缆式悬挂钻孔记录工具，该工具装配有 脉动中子发射源和核辐射传感器和处理器，用于确定地层的硅含量，处理器位于地面上，可 被构造成估计借助用于悬挂记录工具的电缆发送到地表处理器的、基于传感器进行的测量 的记录工具的速率。对于想要估计正钻出的钻孔的深度的熟练钻井专家而言，显而易见，找 到这个问题的解决方案要翻阅US2008/251710，US2008/251710涉及用于确定地层中的硅含 量的电缆式记录工具。US2008/251710没有教导作为将钻孔中子发射源和核辐射传感器进 行的发送工具速率测量值发送到地面的处理器的替代，可将这些速率测量值发送到井下工 具，从而不再需要借助数据传输电缆将大量数据发送到地面的处理器，旋转钻井组件牵涉 到借助数据传输电缆将大量数据发送到地面的处理器是不且实际的或者成本高。
[0017]现有技术的文献提供了电缆作业的或相对不准确的钻孔深度测量工具。然而，无 法在不需要借助电缆将数据发送到地表处理器的情况下测量深度和/或无法以对于自动化 钻孔而言将是充分的准确度来测量深度。
[0018]因此，需要改善如上所述的方法和系统中的一个或多个。
[0019]此外，需要一种准确的井下钻井过程监测单元，该井下钻井过程监测单元可在钻 井期间提供自动化或其他的钻井组件的准确实时增量深度，而不需要复杂和脆弱的从钻井 组件的底部钻具组合(BHA)延伸到地面的有线或无线数据传输设备。
【发明内容】
[0020]按照本发明，提供了一种在地下钻孔中进行井下测量的系统，所述系统包括：
[0021 ]-中子源，其用于发送中子脉冲，以激活包围所述钻孔的地下地层中的靶物质； [0022]-核辐射检测器，其被设置成与核源相距轴向距离L，用于检测被激活的所述靶物 质发射的核辐射;以及
[0023]-井下钻井过程监测单元，其被构造成计基于核辐射和轴向距离L的检测来计算钻 孔中的钻井组件的增量深度。
[0024]所述井下钻井过程监测单元可布置在钻井组件的底部钻具组合(BHA)中并且被构 造成向BHA中的自动化井下钻头导航系统提供关于BHA的增量深度、关联的机械钻速(R0P)、 方位角和/或倾斜度的实时数据并且所述井下钻井过程监测单元可被构造成借助所述BHA 中的信号传输组件将所述实时数据传输到所述自动化钻头导航系统，而不需要将所述实时 数据传输到位于地表的数据处理设备。
[0025]井下自动化钻头导航系统还可包括用于基于存储在所述系统中的关于以下的数 据将所述BHA向着地下钻井靶导航的装置：
[0026]-所述钻孔的井口和所述地下钻井靶的三维地理坐标；以及
[0027] -基于所述钻井过程监测单元收集的关于所述钻孔轨迹的整个长度内所述BHA的 增量深度、方位角和倾斜度的实时数据进行所述钻孔轨迹的估计。
[0028] 可选地，所述检测器可以是圆柱形，具有至少5cm的高度和/或被构造成检测包括 氯-37、钠 -23和/或硅-30的靶物质发射的核辐射并且所述中子源和核辐射检测器(D)之间 的轴向距离L可等于或小于最大值Lmax，其中，轴向距离L大于沿着所述钻孔的轴的所述靶物 质的被激活核的扩散和/或L max在3米和6米之间和/或轴向距离是至少20cm。
[0029] 按照本发明，此外提供了一种在地下钻孔中进行井下测量的方法，所述方法包括 以下步骤：
[0030] -通过中子源发送中子脉冲，以激活包围所述钻孔的地下地层中的靶物质；以及
[0031] -致使被设置成与所述中子源相距轴向距离L的核发射检测器检测被激活的所述 靶物质发射的核辐射;以及
[0032]-致使井下钻孔过程监测单元基于所述核辐射和所述轴向距离L的检测来计算所 述钻孔中的钻井组件的增量深度。
[0033]所述井下钻孔过程监测单元可布置在钻井组件的底部钻具组合(BHA)中并且借助 BHA中的信号传输组件将关于所述BHA的增量深度、关联的机械钻速(R0P)、方位角和/或倾 斜度的实时信息提供到BHA中的自动化井下钻头导航系统，而不需要借助在BHA和地表钻孔 导航和控制设备之间复杂、易碎且成本高的无线或有线数据传输系统将实时数据传输到位 于地表的数据处理设备。
[0034]所述自动化钻头导航系统可包括基于关于以下的所述系统中存储的数据将所述 BHA向着地下钻井靶导航的装置：
[0035]-所述钻孔的井口和所述地下钻井靶的三维地理坐标；以及
[0036] -基于所述钻井过程监测单元收集的关于所述钻孔轨迹的整个长度内所述BHA的 增量深度、方位角和倾斜度的实时数据进行所述钻孔轨迹的估计并且当所述BHA已经达到 或几乎达到所述地下钻井靶时，所述自动化钻井过程可终止，由此，可从所述钻孔取出所述 钻井组件，所述钻孔可完成并且被转换成烃流体生产井，可由所述生产井来产生烃流体。
【附图说明】
[0037] 下文中，将参照附图更详细地并且以举例方式描述本发明，其中：
[0038] 图1示出包括钻井系统的实施例的钻孔的剖视图；
[0039] 图2示出传统定向钻井系统的剖视图；
[0040]图3示出本发明的定向钻井系统的实施例的剖视图；
[0041]图4A和图4B分别示出在时间tl和t2的用于测量增量深度的标记-检测器系统的示 例性实施例的剖视图；
[0042]图5示出用于测量增量深度的示例性系统的剖视图；
[0043 ]图6至图8示出本发明的系统的流程方案的实施例；
[0044]图9A和图9B分别示出时间tl和t2的本发明的系统的实施例的剖视图；
[0045] 图10示出本发明的系统的实施例的透视图；
[0046] 图11示出本发明的系统的流程方案；
[0047] 图12示出本发明的系统的实施例的剖视图；以及
[0048] 图13示出本发明的系统的流程方案。
【具体实施方式】
[0049]在说明书和附图中，类似的参考标号涉及类似的组件。所描绘的组件和所描述的 特征可按许多方式进行组合。
[0050]图1示出延伸到地下地层2中的钻孔1，地下地层2通常包括多个层5、7,每一层具有 不同的特性。层7可以例如是储层，储层包括诸如原油和/或天然气的烃流体。用于钻出钻孔 1的钻井组件3通常包括管状钻柱4,钻柱4从地表8上的钻机6延伸到钻孔1的井下端9。钻柱4 可具有由多段钻杆形成的第一或上段10。钻柱的第二或井下段可被称为底部钻具组合或 BHA 12AHA 12的井下端设置有钻头13。环面14在钻柱4和井眼1的壁之间延伸。钻机6通常 设置有旋转驱动系统15，旋转驱动系统15用于将钻柱绕着其纵轴旋转。
[0051] 底部钻具组合或BHA 12还可包括额外的设备，诸如钻环和钻柱稳定器。钻环可包 括较重的管道片段以提供钻压。
[0052] 在钻井期间，旋转驱动系统15进行操作，以将钻柱4在井眼1中旋转。可借助钻柱4 的内部流体通道将钻井流体从地表栗送到钻头13,钻头13排出钻井流体。钻井流体的回流 通过环状空间向上流向地表，带有因钻头13的切割动作而导致的钻肩。
[0053]烃行业中所关注的新兴领域之一是开发薄和/或致密的储层。这些储层需要比传 统烃储层更大的储层接触面积。增大储层接触面积的选择包括钻相对长的水平片段和/或 钻多个侧钻以形成根状多侧面井眼。
[0054] 这些最新水平的钻井应用中的大部分使用实时地层评价技术和复杂的井下工具 来找寻地理靶地带并且在向前钻井的同时保持与储层接触例如达数千米。在钻井的同时改 变井眼轨迹来实现与地理靶地带更大的接触或者避开风险地带被称为地质导向。地质导向 会得到复杂的三维井眼几何形状。
[0055] 地质导向系统可包括构建和维修成本较高的复杂机电系统。这个成本因素限制了 它们的使用并且会致使使用它们是不经济的。
[0056] 当前地质导向技术的另一种限制是其依赖于人类专业技能理解地下数据和钻井 参数以进行实时决策。一直依赖人类判断为误差留有空间和不一致并且增加了风险承担。
[0057] 可用更简单的工具构造和处理的自动化来实现对当前系统的改善，以使成本更低 并且使性能一致。
[0058]为了实现自动化地下导航的构思，主要的决策之一是决定系统的控制单元20或大 脑的位置。控制单元20将取得来自所有可能测量的输入并且用获取的情报处理它以进行实 时决定。如果系统的控制单元位于地面上，则所有井下测量值22必须被发送到地面，而关于 工具设置方面的决策24必须被发送到井下，如图2中所示。
[0059] 开发基于地表的自动化的关键挑战之一是当前通信系统的能力瓶颈。用于底部孔 和地表之间的通信的当前最新水平的技术是泥浆脉冲遥测技术，泥浆脉冲遥测技术遭遇由 于高频信号的衰减较高而导致随着底部孔的深度增加，数据携载能力降低。商用系统声明， 对于浅井而言，数据速率是大约40比特/秒，对于大位移井而言，减小至3比特/秒，其中在下 行(down-1 inking)操作中有中断的额外开销。这也是信息的地下测量和地表再现之间的时 滞，伴随着有限的数据传递速率严重限制了控制系统实现钻头和井眼轨迹的闭环控制和获 得高频、高分辨率数据采集井下的能力。
[0060] 另一个选择是将"大脑"20(即，控制单元)布置成较靠近钻头13,如图3中所示。与 地层评价相关的井下测量值(诸如，伽玛射线、电阻率、密度、声速记录等)可容易得到，并且 可被用作进行地质导向决策的基础。在井下定位控制单元20将大大减小与地表进行数据传 输的带宽需要。通过定位在井下，控制单元可实时利用高分辨率数据集，从而基本上消除了 时滞问题和数据压缩需要，这将增强决策的准确度和质量。
[0061] 可以使用底部钻具组合12中包括的随钻测量(MWD)系统获取所有至关重要的钻井 参数(诸如，钻压、钻头每分钟的转数(RPM)、钻井流体的流速、温度、压力等）。MWD系统能够 改善监测和优化。
[0062] 并非所有参数是在地下测量的。实际上，目前仍然是在地表测量钻头的深度并且 将把深度发送到BHA。深度的井下测量可不再需要将深度数据从地表发送到BHA并且可增强 自动化地下导航系统的实时能力和准确度。
[0063] 尽管过去已经研究了井下深度估计使其有可能消除基于地表的深度测量的误差， 但目前还没有可用的商用系统。在各种专利和公开文献中已经提出了许多思路，例如，使用 机械轮(US-5341886)、钻具表(US-5896939，US-2008/0105423)、磁性标记物-检测器系统 (US-5341886)、地理标记的相关性(US-7283910)。
[0064] 在以上的旋转之中，使用钻具表是实现最不复杂的。然而，深度估计的准确度常常 不及基于地表的测量。得自系统的深度估计相对粗糙，不准确（以钻杆接头长度为间隔），因 此不适于实时自动化应用。
[0065] 提议基于压力或重力测量的深度测量系统(US-4475386，US-8113041)可仅仅估计 垂直深度并且在井下环境中需要鲁棒的重力或压力模型以及高准确度测量。
[0066] 几个作者已经提议使用类似于用于火箭和飞机的惰性导航系统的加速度计(1^_ 2008/0105423，US-2009/0084546，US-2002/0195276，US-5585726)。然而，这些系统遭遇积 分漂移并且需要辅助系统来校正它们的估计。
[0067] 使用波速以及源和位于钻柱端部的接收器之间的波传播时间可提议一系列新构 思。这些构思的线索来自随钻进行的垂直地震剖面测量。应用的波可以是钻柱中的声波、泥 浆柱中的压力波或地下的电磁波。实现这些构思的最大挑战是创建准确的速率模型和时钟 的漂移。最新水平的基于石英的井下时钟可具有将估计的深度的准确度限于几米的毫秒级 测量准确度。
[0068] 使用地理标记的相关性进行井下深度估计的固有优点是，它可构建在当前随钻记 录系统上，不需要开发新的测量技术。另一方面，可基于模式识别方法用得自信号处理范畴 的图像分析或交叉相关来构建相关部分，这样可进一步减少总开发工作量。这也是优于其 他系统的优异选择，因为可用用于任何地理设置的可用记录数据来测试其可用性并且可针 对特定需要容易地调节系统。
[0069] 参照图4A和图4B，使用地理标记进行增量深度估计的构思涉及用两个或更多个传 感器(S1、S2)沿着钻孔进行连续测量，传感器(S1、S2)布置在例如BHA中并且分开已知距离 L。计算在这些测量中出现同一记录特征的时间差。测得的时间差和已知的传感器S1和S2之 间的空间间隔L允许计算平均机械钻速(R0P)。可对计算出的R0P求积分，以计算沿着钻孔1 的增量深度。
[0070] 在图4A中示出布置的示意图，其中，S1和S2是位于底部钻具组合中并且分隔距离L 的相同传感器。在图4A中，第一传感器S1在时间tl经过具有高伽玛射线特征的区域。第二传 感器S2在时间t2经过相同的高伽玛射线特征gl。因为传感器相隔距离L，所以钻柱在时间 (t2-tl)内已经沿着钻孔前进距离L。
[0071] 因此，平均机械钻速(R0P)是
[0073] 通过对R0P求积分，可计算深度变化Δ D:
[0074] AD = R0P*At (2)
[0075] 而可通过将各深度变化相加来计算增量深度D
[0077] 此构思的要素之一是把多个记录的相关联以识别所有这些记录上的地下的每个 特有点的出现。地下测量的相关面临例如以下挑战：
[0078] 1)钻孔测量会具有错误值，需要从待相关的数据集中过滤掉这些错误值；
[0079] 2)在沿着钻孔的同一位置进行井下测量可能在各个传感器上提供不同测量值。测 量之间的差异需要引入容差带宽以考虑不同工具的测量精确度；
[0080] 3)在同一地层内，沿着钻孔的任何岩石属性的测量值的变化会非常小。因此，为相 关而选择的岩石属性应该具有比预期最小变化小的相关测量精确度。这样会需要在不同地 层或不同传感器设计中使用不同类型的测量。
[0081] 为了应对以上提到的挑战，指示这种相关运用的算法需要是灵活的，以使算法能 够适于地下化学和岩石属性。通常，可使用任何合适的记录组件或一组记录组件进行相关， 只要岩石属性沿着钻孔深度充分变化。
[0082] 系统对于传感器位置是不敏感的。传感器可靠近或远离钻头。优选地，传感器被包 括在BHA 12中。还可针对不同地下位置的多个数据点的相关进行方位测量。
[0083]在图5中示出针对该构思的基本布置。第一传感器S1是钻孔中的前测量仪器。第二 传感器S2是后测量仪器。
[0084] 第一传感器S1收集到包括N个测量值的数据集Y。第二传感器S2收集包括Μ个值的 第二数据集X或模式X。本文中，Μ〈Ν。控制单元20(图5中未示出；参见例如图3)尝试使用预定 算法将数据集Μ匹配在数据集Ν内的特定位置。数据集Ν足够大，以包括对应于数据集Μ的所 有可能匹配。
[0085] 图6示出用于说明构思的使用两个传感器布置进行深度计算的算法的示例。可针 对不止一种岩石属性或更多个传感器来扩展图6的算法。以下，说明算法的主要步骤。
[0086] 在第一步骤50中，开始算法。算法可例如以设定的时间间隔来周期性开始。可供选 择地，可通过显著事件(例如，前传感器S1测得的显著增大或减小的值)触发算法。
[0087] 在第二步骤52中，第一传感器S1和第二传感器S2提供分别包括Ν个数据点和Μ个数 据点的数据集Υ和Χ(另外参见图5)。可通过添加最新获取的数据点并且去除最老的数据点 来动态更新数据集Υ和X。可针对各特定情形和构造来确定Μ和Ν的值。用来自后传感器S2的 数据集X和来自前传感器S1的数据集，在任何时间t s:
[0090] 各数据集中的各个数据点按时间排序，以能够进行比较并且找寻各个数据集之间 的相关性。本文中的比较可考虑预定容差，以处理地层特定特性并且保持预选的准确度。
[0091] 在第三步骤54中，可使用任何相关性算法来识别数据集X与数据集Y的可接受匹 配，使得：
[0095]下一个步骤56检查是否发现可接受匹配。如果未发现，则机械钻速(R0P)将被设置 成零(步骤58)并且深度将被更新(步骤60)。
[0096]如果在步骤56中发现可接受匹配，则将在步骤62中计算来自传感器S1和S2的测量 值之间的时间差At。
[0097]使用该时间差At，计算机械钻速(R0P)(步骤64)。如果传感器之间的空间间隔是 L ：
[0100] 随后，在步骤60中，更新深度。可使用将(步骤64中计算的）平均R0P乘以时间差Δ t 来计算深度增大AD。本文中的深度是沿着钻孔的增量深度。
[0101] 地理标记的相关性中的构成之一是匹配来自多个传感器的响应。为了处于特定目 的找寻最有效的模式匹配算法，探查不同的领域(诸如，文本编辑器中的串匹配、信号处理 和图像处理）。这些方法中的每种表现出自身的优点，从而致使其均最适于不同情形。优选 地，最终选择的算法是最有效的并且能被用于钻孔中的应用。基于下述的特定需要和数据 集的属性，进行本发明中的针对相关性的可能匹配算法的选择。
[0102] 1.数据集中的所有数据是数字。这样允许使用数字运算来比较可能的匹配。
[0103] 2.由于模式X和数据集Y的已知大小，导致可根据需要将导致不匹配的容差带外的 测量值的数量上限进行定性和选择，以简化复杂性。
[0104] 3.不匹配仅限于具有与数据集X中的值对应的数据集Y中的不同值。不匹配不包括 缺失数据集X和数据集Y中的值的情况。
[0105] 4.给定容差内的两个值应该被视为匹配。
[0106] 5.在时常没有合适匹配的情况下，可针对步骤64中的R0P计算，忽略各个测量值或 数据点（例如，特定时间ts的数据集X，或所述数据集X的数据点）。替代地，可使用前一 R0P来 进行步骤60中的深度计算，因为R0P不会从一个测量急剧变化成下一个。
[0107] 6.相比于编辑器中的文本搜索问题，数据集X和数据集Y的大小可被视为是相对小 的。因此，预期匹配算法禁止出现任何计算能力相关问题，即使是使用强力算法。
[0108] 基于独特需要和问题的简化，以下在本文中考虑并且说明用于模式匹配的两个可 能的解决方案。
[0109] 1.归一化的交叉相关算法
[0110] 交叉相关是在信号处理中用于识别两个信号之间的时滞的方法。在地表接收到记 录之后，在确定钻孔末端的平均机械钻速时使用该方法。交叉相关方法对于所有可能的匹 配计算相关系数。相关系数的最高值识别最有可能的匹配。模式的大小仍然保持可变的选 择并且可被最佳化。这种方法就计算而言可比其他算法的成本更高，因为需要检查所有可 能的匹配和/或可涉及更复杂的计算。
[0111] 如果分别通过前传感器S1和后传感器S2进行的测量取得N数据点的数据集Y和Μ数 据点的数据集X，使得适用公式⑷、(5)、(6)和⑴，贝嗵过以下给出任何子集c: 的相关系数：
[0113] 公式（10)中的分母用于将相关系数归一化，使得-1< rdl，其中，界限(值-1或1) 指示最大相关并且"〇"指示不相关。高的负相关指示高相关，但是一系列之一的倒数。
[0114] 交叉相关性的关键优点之一是操纵传感器中的任何漂移的能力。传感器漂移噪声 的任何偏移将不影响交叉相关系数，因为偏移还将反映在分子和分母中的被从数据集Y和 数据集X二者中减去的平均值，从而使交叉相关性系数不变。
[0115] II.混合匹配算法
[0116] 混合匹配算法模式是用于应用记录相关性的匹配算法，并且是强力方法和平方距 离方法的组合。相关算法首先使用强力方法基于拟合和容差标准来找寻可能的匹配。如果 不止一个匹配达标，则选择具有最小平方距离的匹配。
[0117] 混合匹配算法灵活地用于不同大小的数据集X和Y(分别地，数据点Μ和N的数量）。 这允许根据地层属性和钻孔参数的变化进行匹配。可基于测量的类型来选择容差和拟合需 要的值。
[0118] 以下，在本文中，针对控制拟合选择和基于统计的容差的伽玛射线测量来说明算 法。参照图7说明算法的流程图。
[0119] 算法开始于步骤70所指示的预定时间。
[0120] 在第二步骤72中，设置匹配需要μ(用单位％表达）。
[0121] 随后，在步骤74中，第一传感器S1和第二传感器S2提供分别包括Ν个和Μ个数据点 的数据集Υ和X。
[0122] 在步骤76中，并且进一步地，针对τ = η，计算容差τ，其中，i = [ 1，Μ]。
[0123] 对于各i = l、2、…(Ν-Μ+1)，算法检查是否| xi-yi | < η。初始地，在步骤78中选择数 据集X的第一值并且在步骤80中进行比较。
[0124] 对于j = l、2'"M(步骤82)，算法检查是否 |xj-yi+j-1| < Tj。
[0125] 如果计数2 μ，则算法结束。
[0126] 在步骤84中，如果实现了匹配需要，则接受数据集X。本文中的匹配需要是例如步 骤72中设置的匹配需要μ(用单位％表达）。
[0127] 如果数据集X被作为匹配接受，则在步骤86中，计算平方距离。
[0128] 针对数据集Υ中的每个值，重复以上步骤。在步骤90中设置阶梯增量，并且从步骤 78起重复匹配处理。
[0129] 当达到数据集Υ中的第(Ν-Μ+1)值时，即，当数据集X与整个数据集Υ相匹配时，步骤 88参考通向步骤92的数据，以得到步骤86中提供的所有计算出的平方距离的最小平方距离 匹配。
[0130] 测量相同地层属性的两个传感器S1和S2之间的距离可对于计算出的深度关联的 误差产生影响。这个效果是由于任何时间计算出的R0P是时间间隔内的R0P值的平均值的事 实导致的，因为传感器S1和S2经过同一地下位置。本文中计算出的ROP涉及作为如上所述的 作为计算出的井下的R0P。误差引入在地表记录的R0P和井下计算出的平均R0P之间的差异。
[0131] 经常，深度的计算值比深度的实际测量值高。在测试中，误差几乎从大约1% (对于 大约lm的传感器空间间隔距离L)线性增大至大约3%至3.5% (对于大约10m的距离L)。
[0132] 计算出的R0P波动越大，平均和关联误差的效果越大。如果可消除R0P的波动，则可 避免与R0P平均关联的误差，并且能实现理论上的零误差。
[0133] 这些值与混合匹配算法相关。可通过高斯分布逼近伽玛射线测量，从统计学上讲， 任何位置处的测量的大约68%应该产生一个标准偏差内（± 〇)的计数。对于高斯分布，可通 过总计数的平方根来逼近标准偏差。使用误差的传播规则，可通过以下给出来自同一设计 的不同传感器的两个数据集之间的总标准偏差(〇t)
[0135]其中，~和％是同一位置测得的计数的标准偏差。
[0137] 可使用总标准偏差(〇t)作为找寻可接受匹配的容差(τ)。从统计学上讲，来自对应 于同一位置的两个传感器的时间值中的大约68%应该处于± 〇t的容差带内，这样将提供针 对± 〇t容差的68 %的拟合需要。
[0138] 因为y(tP)是未知的并且需要找寻，所以总标准偏差可被逼近为：
[0140]在漂移的情况下，以上计算出的容差值不可保持有效。如果漂移即时增大，则可随 时间推移添加和更新引起漂移的额外项。可通过模式（即，数据集X)的平均值和可能匹配之 间的差异来逼近由于漂移导致的偏移，使得：
[0142]并且可通过以下给出新容差值τ
[0144] 将在其中搜索模式（即，数据集）的匹配的数据集Υ的大小优选地长得足以包含所 有情形下的模式。因需要将找寻故障匹配的任何可能性降至最低来提供数据集大小的上 限，随着数据集大小增加，该可能性将增加。为了将故障匹配的发生降至最低，数据集Y的大 小优选地保持最小。在计算出的平均R0P和已知的传感器间隔距离L的帮助下，可限定逼近 时间窗，在逼近时间窗中，期望后传感器S2经过前传感器S1的位置。
[0145] 假定平均R0P = RP并且传感器距离是L。用P指示R0P中的不确定性。可通过以下给 出用于搜索模式的时间窗Ta:
[0147]如果通过ω给出用于进行各测量的时间窗，则通过以下给出数据集Y的长度
[0149] 模式（即，数据集）的大小还可以是待优化的参数。模式优选地长得足以捕获足够 的变量，从而有助于找寻数据集Y中的特有对应匹配。然而，计算上，功率和时间的考虑导致 数据集X的大小的上限。模式大小可被选定为针对特定地层的常数。对于另一个地层，可使 用算法重复地优化数据集X的大小。
[0150] 例如，使用来自两个传感器S1和S2的响应的线性回归分析，可估计数据集X的最佳 大小。因为来自两个传感器S1和S2的测量值之间的差异可被表现为高斯噪声，所以针对数 据集X的最小代表性模式大小将优选地提供数据集Y和数据集X的测量值之间的线性关系。 即，这两个传感器在预定容差内测量特定地下位置的近似数量的被计数的伽玛射线。一旦 满足这个线性关系，在数据集X中添加更多数据点就将在模式中添加值。
[0151] 为了估计通过基于混合匹配和/或交叉相关的算法计算出的噪声水平、传感器漂 移和总计数的效果，运行各种模拟。针对每种情况，采集时间上的绝对误差。选择25个数据 点(数据集X)的模式大小和100个数据点（数据集Y)的设置大小。为了避免由于不匹配导致 的不切实际的高RP，针对R0P，设置150米/小时的上限。
[0152] 通过模拟得到的关键发现是：
[0153] 1.交叉相关保持不受偏移的影响并且产生相同的结果；
[0154] 2.随着总计数的增大，最大误差减小；
[0155] 3.总体上，随着噪声级的增大，误差增大；
[0156] 4.由于R0P本身的平均而导致的误差可高达大约22m;
[0157] 5.误差更快速地在钻孔的更多水平片段上增大。在水平片段上，伽玛射线读数的 变化更受限制，这样会导致认为错误匹配被用于深度计算。
[0158]以上，伽玛射线测量被描述为模式匹配的可能测量，但其他测量也会是合适的。例 如，还可出于增量深度测量的相同目的，考虑密度测量。伽玛射线和密度测量提供了良好的 可重复性。
[0159] 总体上，基于密度数据估计的深度可以与基于伽玛射线的结果一样好或者更好。 基于密度的测量的最大偏差可以是相对于地表测得深度的大约l〇m(基于，杆片段的数量， SP，钻柱的长度）。基于伽玛射线的深度测量的最大偏差可以是大约15m至20m。实际上，密度 测量需要有源中子源。
[0160] 在用记录相关性计算出的深度的两个误差源：
[0161] 1.错误识别的匹配；以及
[0162] 2.R0P 平均。
[0163] 由于因地下测量的相关，深度计算有固有限制以使测得深度的整体误差降至最 低，导致辅助系统优选地与提议的系统并行使用。辅助系统所需的测量频率取决于主系统 的误差和可接受误差的限制值。接下来，说明用于此辅助系统的两种可能方法。与提议的系 统并行地，添加一个或多个辅助系统可提高深度估计的整体质量。
[0164] 在钻具表数据的帮助下进行校正
[0165] 在估计深度中引入校正的最早解决方案之一是从地表起下行。钻柱的总长度以钻 杆片段数量的形式可得自地表。例如，钻杆的各片段可具有例如大约l〇m的长度。钻柱的总 长度是大约等于片段数量的10倍。钻孔的增量深度将大约等于钻柱的这个总长度，至少在 钻头接合钻孔的井下端时。钻杆的总长度可被发送并且存储到井下系统中。可在形成钻杆 的各片段时触发信号24。然后，信号4可通过泥浆柱发送到井下(参见图2)。
[0166] 在接收到信号24时，井下系统可识别钻柱中添加了另一个杆接头。逐渐地添加将 各管状片段的长度使得井下系统能够计算井中的钻柱的总长度，因此计算钻孔的相关增量 深度。在图8中说明用于此系统的算法。
[0167] 算法开始于步骤100
[0168] 在第二步骤102中，井下控制单元20检查是否添加了钻杆的额外片段。实际上，控 制单元检查它是否已经接收到指示添加钻杆片段的泥浆脉冲信号24(参见图2)。
[0169] 如果没有接收到信号24,则算法继续到步骤106,其中，计算井下深度。随后，在步 骤112中更新增量深度。
[0170]如果控制单元20确实接收到信号24,则在步骤104中，算法指示添加钻杆片段。
[0171] 在步骤108中，算法得到钻柱中包括的杆片段的最近数量。所述数量提供了钻柱总 长度的度量。
[0172] 在步骤110中，算法通过将一个钻杆片段的预定长度与最近存储的增量深度相加， 计算钻孔的增量深度。即，当前增量深度=最近深度+1个钻杆片段长度。
[0173] 在步骤112中，存储井下的当前增量深度，并且取代最近深度。
[0174] 尽管这个系统利用泥浆脉冲遥测技术来传递信号24,但由于不经常进行传递(取 决于接头长度和R0P)，整体带宽要求将受到限制。然而，基于地表的测量的误差也渗透到井 下测量系统，并且在不引入用于估计与基于地表的测量关联的误差的任何算术模型的情况 下，不可进行绝对校正。
[0175] 通过人造标记跟踪进行校正
[0176] 基于地理标记相关性的构思，增量深度的井下估计的另一个构思是基于在地下形 成标记。原理基本上是关于在时间tl使用源120(N)在地层2中形成标记122(图9A)，然后在 时间t2用后检测器124(D)检测标记122。
[0177] 当检测器124检测到标记122时，推导出钻柱已经前进了距离L，该距离L等于源122 和检测器124之间的间隔L。增量井下深度可被计算为标记-检测器间隔L的多倍。也可通过 测量形成标记和检测标记之间的时间差A t = (t2-tl)来计算平均R0P。
[0178] 源120可以是核源。基于诸如Californium-252的放射性核素或加速剂的系统可用 作记录的核源，即，能够测量诸如本体密度和孔隙率的地下属性。在用中子轰击地下元素 时，可激活这些元素，使核素不稳定，核素在去激活时产生具有特定特性的辐射。可在可切 换或脉动中子源的帮助下，激活钻孔的特定片段。中子源可以是商购的脉动中子发生器。脉 动中子装减? PND?_):谱仪是例如由国际有限公司市售的。
[0179] 地层中的被激活元素的辐射可包括任何核辐射、α、β或γ射线。实际上，辐射将是 γ射线辐射。本文中的γ射线指示电磁辐射，具有近似于X射线的波长，但通常具有较短波 长，进而具有较高频率。本文中的伽玛辐射源自伽玛衰减，即，将经常发射单伽玛光子的被 激活核的衰减。
[0180] 通过识别源自不稳定核素的特征发射，可识别这种发射的位置。这样提供了地下 标记-检测器系统使用中子发生器作为中子源120以形成被激活地带作为标记122并且使用 核辐射检测器124检测被激活地带发射的核辐射的构思。检测器124可以是基于伽玛射线探 测器的检测器。
[0181] 地层通常包括巨量的各种元素，这些元素例如构成岩石地层的化合物和相关流 体。然而，并非所有地下元素适于本发明的标记-检测器系统。元素的实用性取决于三个主 要因素：
[0182] 1.中子俘获剖面:这意味着给定温度下针对给定中子能量的任何核素被激活的概 率；
[0183] 2.地下的丰度:浓度越高，被激活原子的产量将越高并且源自去激活的特征辐射 越多。
[0184] 3.衰减常数:因俘获中子而形成的放射性核素的半衰期意味着任何给定时间的衰 减原子的数量。在极短时间内高度不稳定衰减的核素会不足以使检测器达到被激活位置并 且俘获衰减光子。另一方面，几乎稳定的核素将不衰减充足的数量，以提供可检测信号。
[0185] 在考虑到地下最常见元素的半衰期时，最有利元素中的一些是氯同位素 C1-37、硅 同位素 Si-30和天然钠 Na-23。氯原子和钠原子二者都在含盐地层水中作为离子存在。硅常 常出现在砂质和泥质地层中的娃酸盐矿物质中。即使氯和娃本质是丰富的，C1-37和Si-30 同位素也经常以较低浓度(分别地，24%和3%)出现并且具有较低的中子俘获剖面(分别 地，0.433和0.11靶）。会需要进行详细计算来评价被激活原子和所得光子的数量。
[0186] 对于任何被激活元素，在任何给定时间点，衰减的量取决于放射性核素的浓度。随 着时间增加，被激活核素的浓度呈指数减小，源自衰减的辐射也如此。因此，为了得到更强 的信号，优选地，检测器和源之间的间隔小。但是来自源的中子在所有可能方向上前进并且 不限于源正上方的地层。由于在标记-检测器间隔较小的情况下，激活区域较大，因此精确 定位对应于激活的源位置是非常困难的。因此，需要决定源和检测器之间的距离，使得被激 活核的浓度保持足够高，而位置保持不可区分。另外，检测器可与中子源屏蔽。
[0187] 为了研究构思的可行性和来自地下被激活核的所得信号，执行核模拟来模拟地下 的中子和光子传输。用MCNP执行这些模拟。这是够通过整体几何形状进行基于模拟时间的 核粒子(诸如，中子、光子和电子)的传输的基于蒙特卡洛的通用核转运代码。基于蒙特卡洛 的计算涉及用各粒子的概率分布在使用随机数字进行统计学取样的帮助下各粒子的行为 模拟。它由实际上遵循在来自粒子源的各粒子在其从生到灭的整个过程中的结束类别（吸 收、逸出等）。使用传输数据对概率分布进行随机取样，以确定其寿命中各阶段的结果。构思 验证需要两个问题的模拟：
[0188] 1.激活合适的同位素
[0189] 2.检测因去激活而生成的光子
[0190] 这两个问题都需要在MCNP代码中进行分开的模拟。在来自中子发生器的中子的第 一模拟转运中，在地层中执行中子的俘获。第一问题的输出给出因中子俘获而得到的被激 活同位素的位置。使用这个输出来形成针对第二问题的输入。因为光子的起始位置是被激 活核素的位置并且是已知的，为了完成第二问题的输入，使用在随机数字和核素的衰减常 数的帮助下对衰减时间的统计学取样将衰减的时间与各核素关联。在第二模拟中使用这个 间隔和时域光子源来模拟光子在地下的运转，以对可到达检测器的光子的数量进行计数。 这些光子将检测器中的信号的强度定量并且其数量应该大得足以被检测到并且与背景噪 声区分开，以便实施该构思。
[0191] 在第一模拟中，形成地层的简化模型。包括岩石块(地层）的模型包括核源处于其 中的钻孔。岩石块是5m X 5m X 5m。钻孔的直径是20cm。核源是5cm高。模型化岩石类型是具有 30 %孔隙率和30，OOOppm的盐分的沙石储层。
[0192] 假设盐分主要是由于有氯化钠。本文中，大概氯原子的大约25%是天然的C1-37同 位素。
[0193] 用水模型化填充直径20cm的钻孔，以使用基于水的钻井流体来模拟钻井。中子源 和检测器优选地被包括在BHA中。因此，源和检测器被认为是位于钻孔的轴上。中子源是 14. IMeV的、每秒大约108个中子的中子产量的D-T源。这是脉动中子发生器的工业标准规 格。源被建模为5cm长且直径为3cm的圆柱体(图10B) 〇
[0194] 在第二模拟中，检测器124被建模为具有5cm的高度Η并且具有6.5cm的内径ri和 8cm外径r。的井型设计的Nal闪烁晶体，以得到较大的俘获量（图10A)。检测器将在向前钻井 的同时随着BHA移动。
[0195] 为了确定关于信号强度的检测器位置的效果，将四个检测器124建模(如图10B中 所示）。一个检测器124(1)被建模在源120的位置。其他三个检测器被设置成与源120的间隔 距离是10cm(检测器124(2))、15cm(检测器124(3))和20cm(检测器124(4))。针对各检测器， 计算随时间的光子通量，以识别其中信号强度充分强的时间窗。
[0196] 针对108个中子执行中子运转模拟。在所谓的PTRAC文件中得到各被激活核的位 置。在后续处理PTRAC文件时，发现被激活原子Na-23、Cl-37和Si-30的数量分别是14050、 3222和5555。因为蒙特卡罗模拟使用的是概率分布的随机取样，所以结果的不确定性常常 随着工作用源离子的数量(nps)的增大而减小。当源离子的数量增大10倍时，检测器4中的 105秒（27.8小时）内的每个源离子的总通量的标准偏差减小至31 % (从0.0325减小至 0.0100)。
[0197] 当中子从源均匀地在所有方向上开始时，所关注的是检验沿着钻孔的轴被激活的 核的扩展。测试表明，被激活核的轴向扩散限于相对于源位置的大约20cm。最大数量的被激 活核将处于源的位置。在源位置的两侧，根据高斯分布，被激活核的数量将下降。测试表明， 沿着钻孔轴的被激活核的数量下降至小于相对于源位置的50%，在+/-5cm内。被激活核相 对于标记122的轴向中心的轴向扩散限于大约20cm。可通过使用准确检测器来提高检测准 确度。
[0198] 由于轴线扩散受限制，导致可以以相对良好的准确度来确定标记122的轴向位置。 准确度是例如大约+/_5至10cm。使用本发明的标记-检测器检测核标记能实现的准确度可 以是大约+/-5cm或更小。因此，本发明的标记系统可提高增强深度测量的准确度，例如当使 用本发明的核标记系统作为辅助系统时。
[0199] 检测器中的信号强度取决于
[0200] 1.检测器与被激活位置的接近度；以及
[0201 ] 2.在激活核标记之后过去的时间。
[0202]放射性核素遵循指数衰减。如果每个时间单元产生的光子的数量保持恒定，则检 测器中的信号应该在接近被激活位置122的同时由于被激活核的浓度增大而增大。但是，由 于指数衰减，导致核标记发射的光子的数量和关联信号在时间上呈指数减小。因此，检测器 124中的实际信号将是由于检测器接近源位置而导致信号线性增大（即，随着在钻井期间在 图10B中检测器从位置3移动到2到1到4,信号强度增大并且随后减小)和核发射信号本身呈 指数减小的结果。
[0203]光子转运的分析已经表明，在接近被激活地带的同时，通量增大例如从20cm至0cm 的距离的大约5倍。假设最适宜的R0P是大约8米/小时，可在大约90秒内钻出20cm的距离。在 这90秒期间，核发射的信号强度减小不到5%。对于大约5至15英尺/小时（1.5至5米/小时） 的范围内的R0P，可在大约8分钟至2.4分钟内钻出20cm的距离。在一些典型元素的半衰期 内，这是很好的，从而使得能够在所指示的准确速率下使用本系统。例如，氯C1-38具有大约 37分钟的半衰期，而钠 Na-24具有大约15小时的半衰期。根据提供核衰减信号的特定同位 素，核发射的信号前度减小减少(差不多）不足5%。当检测器接近被激活地带122时的信号 强度增大是相比于典型指数核衰减的主导因素。
[0204] 依据上文，在具有大约8米/小时或更多的机械转速的典型钻井操作期间，尽管有 核衰减，但以上(即，大约+/-5cm或更小）的准确性将是能实现的。
[0205] 针对入射伽玛射线生成的光子数量基于检测器中使用的闪烁晶体的类型而变化。 然而，对于每种晶体，将入射伽玛射线转换成可见光光子的效率取决于日射辐射的几何形 状和能量。
[0206] 横越检测器的大部分光子的能量是大约0.5MeV或更小。随着闪烁晶体的效率随着 光子能量减小而增大(例如，Nal，参见图10A)，这将导致检测器中的信号较高。与Si-31、C1-38和Na-24的特征能量相关的峰值分别是大约1.3MeV、2. IMeV和2.75MeV。
[0207] 峰值相对小，所以检测器124优选地相对敏感。例如，在实施例中，检测器可包括多 个检测器元件，它们一起提供标记122的检测信号。在另一个实施例中，检测器可相对大。实 际上，本文中相对大的检测器将意味着包括相对大检测表面或检测晶体的检测器。以下，更 详细地说明后者。
[0208] 井下深度估计的动机是应用井下深度估计进行自动化地下导航。准确深度估计可 有助于基于闭环进行自动化地下导航。井下可用的深度估计还可消除由于基于表面的R0P 和深度测量导致的误差。
[0209] 在后续部分中，提出针对自动化地下导航的方案。本文中，使用井下深度和测量数 据来遵循预定的钻孔轨迹结合可操纵的钻井系统。
[0210] 地下导航的主要目的是：
[0211] 1.在保持最佳井轨迹的同时，到达储层；
[0212 ] 2.在实现地质导向目的的同时，遵循储层;和/或
[0213] 3.避免在去往储层的路上的地理危害。
[0214] 图11示出使用自动化可操纵钻井到达靶储层的简化方案。
[0215] 用于自动化钻井的方案开始于第一步骤140。
[0216] 在步骤142中，形成合成记录。本文中的合成记录是基于预期地层的预期属性的预 先计算出的记录结果。所述预期属性和预期地层是使用在开始井下钻井之前指定的地理测 试的结果通过专业地理软件提供的。合成记录是在钻井开始之前预先计算出的，并且被存 储在联接到控制单元2 0的存储介质中。在自动化钻井中，所述存储介质（未示出）和控制单 元20被包括在BHA(参见图3)中。
[0217] 如步骤146所指示的，通常包括在BHA 12中的随钻记录(LWD)工具(未示出）提供了 随钻的地层记录结果。
[0218] 在后续步骤146中，该方案将LWD工具的记录结果与合成记录结果进行比较。
[0219] 在下一个步骤148中，该方案检查记录结果是否指示钻孔的井下端是否已经到达 预定地层。本文中，该方法可考虑特定阈值特性(诸如，深度、地层压力、密度、电阻率等）。
[0220] 如果步骤148的结果指示钻孔在相对于储层的预定阈值距离内，则该方案切换至 步骤150。本文中，该方案切换至地质导向方案。
[0221] 钻井工程师使用可用地震测量、记录数据和地理信息来形成井计划，以避免与现 有井碰撞并且实现最佳井眼轨迹到达储层。定向钻机使用这些井计划进行井布置。在自动 化钻井环境中，可使用井计划进入储层，而此后可调配地质导向算法对储层片段进行钻井 来满足井靶。
[0222] 井计划包括测得的深度、倾斜度和方位值。还用定向测量来测量倾斜度和方位值 二者。通过连续监测计划好的测得的井路径并且计算偏差，可实时采取校正动作来遵循计 划好的钻孔轨迹。可预先在源自偏移井的记录数据的帮助下形成合成日志响应以识别正靠 近储层并且相应地可创建条目。类似地，还可使用此合成日志响应来确定用于设置套管靴 的关键地理深度并且因此可在原始井计划中进行校正。
[0223] 可在地表将典型的可操纵钻井系统编程。该程序可包括用于实现所需累积和转动 速率的工具设置和工具面角度。程序提供了能够使三维钻孔轨迹复杂的灵活性和经度。对 于不同工具，可用工具设置的数量是不同的。在任何设置下可实现的狗腿严密性取决于诸 如岩石属性、孔的倾斜度、BHA设计、钻孔参数等许多参数。在这些参数之中，诸如各向异性 和地层强度的岩石属性具有相对高不确定性并且经常造成相对于预测的行为和轨迹有偏 差。
[0224] 钻井自动化的井下系统使得能够利用钻头性能的实时井下测量和井下钻井参数 来重复地校正钻头行为模型并且使用测量来调节工具设置。
[0225] 为了实现自动化地下导航，依据井下工具的角度构建能力将井下工具的性能建 模。可使用简化的算术模型来模仿通过复杂建模软件预测的工具性能。在简化模型中，可假 设能实现的狗腿严密性⑴)与钻头13上的侧切割力(F s)成比例。
[0226] 图12示出具有钻头13的钻柱4的井下端。钻柱的中线170相对于垂直172成角度α， 垂直172通常是对齐重力的笔直向下的线。狗腿严密性(β)可被表达为：
[0227] P = kFs (18)
[0228] 其中，k是比例常数[deg/Nm]，可引起诸如位设计、岩石属性和BHA设计的因素。
[0229] 就"推动钻头"设计的情况下对钻头13的侧向力Fs可包括由于孔倾斜导致的钻柱 重量W的和由于钻柱弯曲导致的钻柱压缩的分量。另外，侧向力Fs可包括操纵装置174施加 的力。操纵装置可包括接合井下壁的可回缩焊盘，以提供所述力。
[0230] 可使用稳定器的各种布置来提供BHA所需的构建、下降或保持趋势。这些步骤允许 BHA弯曲成期望方向，成为通过稳定器支承于井下壁的弹性梁。可操纵钻井系统可例如使用 液压操作焊盘来形成侧向力和/或使用流体喷射。
[0231 ]传统上，从地表起通过定向钻井控制侧向力。定向钻机在泥衆脉动遥测计的帮助 下与井下系统通信，如图2中所示。
[0232] 在简化布置中，可假设靠近钻头稳定器的人形成小摆测效应。这将增加下降速率 并且减小建造速率。基于这个假设，建造速率可被简化成
[0233] P = k(FP+ff sina) (19)
[0234] 本文中，钻孔以相对于垂直轴172的角度a入射，如图12中所示。
[0235] 为了验证模型并且用比例常数k、重量W和侧向力匕的真实值来调节它，使用可操 纵钻孔数据，如在钻出特定钻孔的同时针对特定可操纵钻井工具得到的。可操纵钻井工具 可具有多个(例如，10至15个)工具偏转的设置，从而得到不同的侧向力F P。
[0236] 为了模拟相同行为，可假设侧向力匕随着工具设置的增大线性增大，从而忽略其 他效应。使用k= ldeg/kNm、W= 1.4kN和FP = 5kN的初始值，可形成计算机程序预测的BHA行 为。
[0237] 在证实井下深度测量与可操纵钻井系统的控制结合时使用相同的模型来遵循给 定的钻孔轨迹。实际上，实际性能和模型预测之间的任何偏差可引起可基于实时型评价来 更新的比例常数k。模型的这种简化使自动化钻井系统更稳健并且操作更容易。
[0238] 检查并且保持与预定轨迹的依附。井眼计划（即，存储的预定轨迹)和实际井眼轨 迹的实时比较导致要调节轨迹，以尽可能符合计划好的轨迹。
[0239] 然而，可操纵钻井系统对地层的响应和钻井参数会偏离预测行为。在任何给定的 工具设置，建造或下降角度会不同并且导致与井计划有偏差。表面上的定向钻井改变可操 纵钻井工具的工具设置，以使基于经历的偏差最小。
[0240] 实际上，由于当前泥浆脉冲遥测技术系统中的带宽限制和关联时滞的限制，导致 实际轨迹和计划好的轨迹之间将仍然有偏差。在自动化系统中，计划好的轨迹被存储在其 中将可得井下深度、倾斜度和方位角的BHA。通过将实际轨迹的实时测量与计划好的轨迹相 比较，可使偏差最小并且可提高钻井性能。
[0241]参照图13中的方案，描述根据本发明的自动化钻井。
[0242] 该方案开始于步骤180。
[0243] 在步骤182中提供钻孔特性的测量值。在钻孔的钻井期间，测量仪器(诸如，LWD工 具)提供例如钻孔倾斜度、方位角和深度的测量值。这些测量值可设置在沿着钻孔轨迹的任 何地下位置。可使用如上所述方法中的任一种来测量本文中的深度。可在基于表面的深度 测量的辅助下检查和更新井下深度测量。
[0244] 还提供相同钻孔特性的参考值。在步骤184中，基于钻孔的模型，提供参考值。如上 所述，参考值是与钻孔的模型，基于地表测量、基于测试的钻孔的模型化等。
[0245] 在步骤186中，该方案检查是否已经达到靶深度TD。如果达到了，则在步骤188中， 该方案结束。
[0246] 如果还未达到TD，则该方案在步骤190中计算计划好的轨迹的测量值和参考值之 间的偏差。
[0247] 在后续步骤192中，该方案计算所需的建造速率。本文中的建造速率指示倾斜度的 增大（即，角度a的变化）。建造速率可通常表达为每个距离单元的角度(例如，1°/1〇〇英尺或 l°/30m)〇
[0248] 执行测试已将实时井下钻井自动化系统的性能与基于地表的自动化钻井系统进 行比较，以执行计划好的钻孔轨迹。基于地表的系统取决于用于通信的泥浆脉冲遥测技术 并且每添加一个钻杆片段在设置中可能的变化次数是有限的。基于井下的系统可进行工具 设置的更频繁变化。井下系统可实现实时闭环控制。
[0249]实际上，井下系统可例如当已经钻出小于2m时频繁进行设置变化。井下系统可例 如每0.1m改变设置。基于地表的系统可只在已经钻了至少2.5m之后改变。
[0250]测试表明，由于计划好的轨迹的高扭转，实时井下控制和基于地表的离散控制二 者偏离计划好的轨迹。检查东边的相对于沿着深度的计划好的轨迹的偏离（A E)、沿着北边 的相对于计划好的轨迹的偏离（A N)和沿着垂直的相对于计划好的轨迹的偏离（△ TVD)。针 对井下系统的北边方向上的误差大大小于基于地表的系统的误差(2.4%)。对于东边方向， 基于地表的系统和井下系统的误差属于同一大小级别(0.5%至0.9% )。
[0251] 测试还用针对自动化钻井的控制创建了井下深度测量的可积分性。
[0252] 基于用于估计井下深度的测试，得出以下结论。
[0253] 1.可使用地下测量的相关性来估计井下深度；
[0254] 2.可将伽玛射线和密度测量成果用于进行相关；
[0255] 3.交叉相关性和混合模式匹配算法针对低噪声级别执行地足够好；
[0256] 4.通过优化模式大小，可依据准确性和计算效率，提高算法的性能；
[0257] 5.基于交叉相关性的方法对于传感器的漂移是稳健的；
[0258] 6. R0P变化是地下测量的相关性的误差的关键来源。可通过例如使R0P变化最小 和/或使源120和检测器124之间的间隔L最小将这个误差降至最小；
[0259] 7.对于包括两个伽玛射线传感器S1和S2的深度测量系统（图4A)，增量深度的误差 可在钻孔的水平片段上增大，这通常表明随着钻孔在同一地层内延伸，伽玛射线辐射的变 化越来越小；
[0260] 8.误差随着与各数据点关联的计数数量的增大而减小；
[0261] 9.可在辅助系统(诸如，核标记-检测器系统)的帮助下减小整体误差；
[0262] 10.核标记的信号强度取决于地下组成。可通过预先选择一个或多个靶原子来进 行优化，特定地层中的靶原子可相对丰富。然后，可针对靶原子的伽玛衰减来优化检测器。
[0263] 11.与核标记-检测器系统关联的误差独立于R0P变化；
[0264] 12.井下深度与RSS控制的结合可实现井下自动化井计划执行；
[0265] 13.检测器124优选地是为了收集较高计数率而优化的伽玛射线传感器。
[0266] 如上所述，本发明推荐在地层中使用人造核标记的系统作为辅助深度测量系统。 这种系统对地理特征和属性的依赖性低，致使系统适于整个钻孔，包括其倾斜和水平的片 段。可基于特定应用来选择测量分辨率。本发明使用可用工具来提供稳健标记系统并且不 需要开发硬件。可根据检测器的精度和/或敏感性来选择标记系统的准确度。
[0267] 市面上可得合适的伽玛射线检测器。由于伽玛射线检测器基于其放射性而具有的 低成本、可靠性和识别不同地理层的能力而例如常被包括在LWD系统中。
[0268] 伽玛射线发射是例如在不稳定同位素的放射性衰减期间出现的。对放射性衰减的 宏观观察涉及对有限时长内成千上万不稳定核素的衰减发射的数量进行计数。因为存在大 量事件（即，放射性衰减发射)和给定不稳定核发生衰减(根据半衰期）的小概率，所以可通 过泊松或高斯概率分布将放射性计数速率进行建模。对于这些概率分布中的任一个，分布 标准偏差等于估计的均方根。如果伽玛射线的总计数是X:
\并且相对
即，如果总计数X增加，则相对误差减小。因此，本发明的工具 i 和数据采集方案优选地被设计成提供属于信号中的各信号点的高总计数。后者将提高信号 的可重复性。
[0269] 改进总计数X的概率中的一些包括晶体大小的增大、由于包围检测器的保护圈的 厚度减小而导致的衰减减小、和/或增大每次测量的时间窗。
[0270] 核源可包括基于氘-氘或氘-氚聚变的电动脉动中子源。中子源可在高达150°C和 20，000psi (大约1400巴）的环境中操作。核源是使用氖-氖和/或氖-氣反应的相对紧凑的装 置。这些反应可分别产生大约2.5MeV和14. IMeV的能量的中子。以下，示出对应的反应。
[0273] 在公式（20)的反应中，氦原子具有大约0.82MeV的能量，中子具有大约2.45MeV的 能量。在公式(21)的反应中，氦原子具有大约3.5MeV的能量，氦具有大约14. IMeV的能量。其 他聚变反应也是可能的。
[0274] 装置的基本构造由用于产生带正电离子、离子加速器(例如高达110kV)和装载氘、 氚或二者的混合物的金属氢化物靶组成。
[0275] 伽玛射线检测器124可包括用于井记录的基于山锁的伽玛射线检测器（参见图 14)。检测器124可使用在与伽玛射线光子相互作用时产生光的闪烁晶体210。晶体可设置有 保护盖212。光电阴极214可接合晶体210,以将伽玛射线转换成可见光子。在闪烁晶体中添 加少量杂质(被称为活性剂），以增强可见光子的发射。通过活性剂杂质发光的一个重要结 果是，本体闪烁晶体210对于闪烁光而言是透明的。
[0276] 在进入晶体210时，伽玛射线可通过三个不同过程来产生快速电子：1)光电效应； 2)康普顿效应;和3)电子偶的产生(光子能量>1.02MeV)。这些快速电子噪声闪烁并且产生 可见光光子216。可见光光子随后经过光电倍增管218,光电倍增管218按比例增大信号，使 得最终以电脉冲（即，一个或多个阳极220和一个或多个倍增器电极222之间的电势）的形式 检测到。脉冲高度将入射伽玛射线能量定量。
[0277] 电脉冲的计数将与晶体210反应的伽玛射线光子的量定量。由于相比于电子的结 合能量，伽玛射线能量通常非常高，因此当伽玛射线由于光电效应而损失其能量时，几乎所 有与其关联的能量被传递到电子并且这些电子在光谱中形成光峰。另一方面，与康普顿散 射相关的光谱几乎是平坦的平台，因为康普顿散射电子的能量分布是恒定的，其能量范围 是从〇至被称为康普顿边界的最大能量损失。
[0278] 可使用各种有机和无机晶体和液体作为闪烁剂，诸如，BG0(锗酸铋）、LaCl3(Ce) (掺杂铈的氯化镧）、聚萘二甲酸乙二酯等。掺杂铊的碘化钠晶体[Nal(Ti)]由于其较高的光 产量是检测器的晶体的好选择。
[0279] 为了提高检测效率，晶体210可具有环形设计（图10A)。外径r。可以是大约8cm，内 径^可以是大约6.5cm，并且高度Η可以是大约5cm或更大。为了增大准确性和/或检测器灵 敏度，根据工具柱4或BHA 12中的可用空间，检测器晶体的外径和/或高度可增大。
[0280] 检测器计数效率η可被定义为检测器nd中计数的光子的数量与源ns发射的光子的 数量之比。可根据以下因素中的一个或多个，优化检测器计数效率：
[0281] i)晶体几何形状；
[0282] ii)用覆盖件212屏蔽晶体；
[0283] iii)晶体的光学效率，其将产生的总光子之中的可进行管理以通过晶体传输的光 子定量；
[0284] iv)光阴极的量子效率。这是所发射的光电子的数量与入射的伽玛射线光子的数 量之比；
[0285] v)倍增器电极的倍增因子。这是所发射的辅助电子与主要入射电子216的数量比。 辅助电子的发射是统计学过程，所以倍增因子不是常数，而是对于不同事件，围绕平均值进 行变化；
[0286] vi)阳极收集的效率。这是在阳极接收到的电子216的数量与最近倍增器电极222 处产生的总电子的数量比。
[0287] 在存在于地层中的众多元件之中，在优选实施例中，本发明的标记系统可特别地 基于一些元素的丰度和/或半衰期来瞄准一些元素。本文中，提供三个同位素的核属性，可 根据本发明的核标记系统来瞄准这三个同位素。
[0288] 氯-37。氯经常被包括在被限定在地层微孔内的地层水中。氯被溶解在水中，成为 氯离子。它具有33靶热俘获剖面并且可在衰减时产生1.5MeV至8.6MeV之间的密集伽玛射 线。C1-37是在天然氯中的具有24.2%的丰度的氯的天然同位素。可通过中子俘获将C1-37 转换成C1-38或Cl-39(氯的同位素）。对于C1-37,中子俘获剖面面积可以是大约0.433靶。
[0289] 与C1-38的衰减相关的伽玛射线的特征能量是1.6MeV和2. lMeVXl-38的半衰期是 大约37.24分钟。用于检测因衰减产生的伽玛射线的时间窗口优选地被设置在0分钟至大约 37.2分钟(50 %耗尽的时间）内。以下，总结对应的反应。
[0292]钠连同氯经常还存在于被限定在地层微孔内的含盐地层水中。钠-23具有大约 0.530靶的核俘获片段。在进行中子俘获时，Na-23通常转换成Na-24，Na-24具有大约14.96 小时的半衰期。与Na-24的β衰减相关的特性伽玛射线具有大约1.369和2.754MeV的能量。对 应反应是
[0295] 硅是泥土中最丰富的元素之一。天然硅包含Si-28(大约92.23%)、Si-29(大约 4.67% )和Si-30(大约3.1% )。本文中，Si-30可以是中子激活所关注的。Si-30具有大约 0.11靶的种子俘获剖面并且可被转换成Si-31，Si-31具有大约157.3分钟的半衰期。源自衰 减的特征伽玛射线的能量是大约1.266MeV。
[0298]放射性同位素通常遵循指数衰减。通过衰减常数来管理核的衰减时间。在从激活 后的任何时间t，通过以下给出时间t的剩余激活核N相对于时间to的原始激活核No:
[0302] 其中，λ是衰减反应的激活常数(^1)，并且T1/2是放射性核素的半衰期。
[0303] 在优选实施例中，用于检测的时间窗口限于大约靶向核素的半衰期。通过距离L结 合机械钻速来确定时间窗口。后者取决于各种方面，但通常可落入大约5米/小时至10米/小 时的范围（大约15英尺/小时至35英尺/小时的范围）内。考虑以上提供的靶物质的半衰期时 间，可在(几乎)〇和最大值Lmax之间的范围内选择源120和检测器124之间的距离L:
[0305] 在以上的表中，提供Lmax的示例性值。本文中，基于预期的R0P来计算Lmax:Lmax = (被激活靶物质的半衰期）X (预期平均R0P)。可在同时考虑到所需准确度的情况下，在〇(非 常接近源)和Lmax之间选择L的值。
[0306]在实施例中，本发明的系统可包括多个检测器D1、D2等，各检测器相对于源位于对 应的轴向距离L1、L2等。优选地，L1、L2等中的每个等于或小于Lmax，如以上表中提供的。可 分别基于所述L1、L2等的倍数来计算增量深度。可比较使用各个检测器的信号计算出的深 度值并且求平均，以进一步提高准确度。
[0307]本发明不限于上述实施例，其中，在随附权利要求书的范围内能料想到各种修改 形式。可例如组合各个实施例的特征。
【主权项】
1. 一种在地下钻孔中进行井下测量的系统，所述系统包括： -中子源，用于发送中子脉冲，以激活包围所述钻孔的地下地层中的靶物质； -核辐射检测器，位于与所述中子源相距轴向距离L处，用于检测由被激活的所述靶物 质发射的核辐射;以及 -井下钻井过程监测单元，被构造成基于核辐射的检测和所述轴向距离L来计算钻孔中 的钻井组件的增量深度。2. 根据权利要求1所述的系统，其中，所述井下钻井过程监测单元布置在钻井组件的底 部钻具组合BHA中并且被构造成向自动化井下钻头导航系统提供关于BHA的方位角和/或倾 斜度、增量深度、关联的机械钻速ROP的实时数据。3. 根据权利要求2所述的系统，其中，所述自动化井下钻头导航系统布置在所述BHA中 并且所述井下钻井过程监测单元被构造成借助所述BHA中的信号传输组件将所述实时数据 传输到所述自动化钻头导航系统，而不需要将所述实时数据传输到位于地表的数据处理设 备。4. 根据权利要求3所述的系统，其中，所述自动化钻头导航系统包括用于基于存储在所 述系统中的关于以下的数据将所述BHA向着地下钻井靶导航的装置： -所述钻孔的井口和所述地下钻井靶的三维地理坐标；以及 -基于所述钻井过程监测单元收集的关于所述钻孔轨迹的整个长度内所述BHA的方位 角和倾斜度、增量深度的实时数据进行的所述钻孔轨迹的估计。5. 根据权利要求1所述的系统，其中，所述检测器是圆柱形的。6. 根据权利要求1所述的系统，其中，所述检测器具有至少5cm的高度和/或被构造成检 测靶物质发射的核辐射，所述靶物质包括氯-37、钠-23和/或硅-30。7. 根据权利要求1所述的系统，还包括被布置成与所述中子源相距预定轴向距离LdPL2 的至少两个检测器。8. 根据上述权利要求中的任一项或任何组合所述的系统，其中，各轴向距离。。和"等 于或小于最大值Lmax。9. 根据权利要求8所述的系统，其中，各轴向距离L、LjPL2是至少20cm并且大于所述靶 物质的被激活核沿着所述钻孔的轴的扩散和/或L max在3米和6米之间。10. -种在地下钻孔中进行井下测量的方法，所述方法包括以下步骤： -通过中子源发送中子脉冲，以激活包围所述钻孔的地下地层中的靶物质；以及 -致使位于与所述中子源相距轴向距离L处的核辐射检测器检测由被激活的所述靶物 质发射的核辐射;以及 -致使井下钻孔过程监测单元基于所述核辐射的检测和所述轴向距离L来计算所述钻 孔中的钻井组件的增量深度。11. 根据权利要求10所述的方法，其中，所述核辐射包括伽玛射线和/或所述检测器检 测靶物质发射的核辐射，所述靶物质包括氯-37、钠-23和/或硅-30。12. 根据权利要求10所述的方法，其中，所述井下钻孔过程监测单元布置在钻井组件的 底部钻具组合BHA中并且将关于所述BHA的方位角和/或倾斜度、增量深度、关联的机械钻速 ROP的实时信息提供到自动化井下钻头导航系统。13. 根据权利要求12所述的方法，其中，所述自动化钻头导航系统布置在所述BHA中并 且借助所述BHA中的信号传输组件将所述实时数据传输到所述自动化钻头导航系统，而不 需要将所述实时数据传输到位于地表的数据处理设备。14. 根据权利要求13所述的方法，其中，所述自动化钻头导航系统包括基于所述系统中 存储的关于以下的数据将所述BHA向着地下钻井靶导航的装置： -所述钻孔的井口和所述地下钻井靶的三维地理坐标；以及 -基于所述钻井过程监测单元收集的关于所述钻孔轨迹的整个长度内所述BHA的方位 角和倾斜度、增量深度的实时数据进行的所述钻孔轨迹的估计。15. 根据权利要求14所述的方法，其中，当所述BHA已经达到所述地下靶时，所述自动化 钻井过程终止，由此，从所述钻孔取出所述钻井组件，所述钻孔完成并且被构造为烃流体生 产井并且由所述井产生烃流体。
【文档编号】G01V5/10GK105874354SQ201480071991
【公开日】2016年8月17日
【申请日】2014年12月19日
【发明人】P·杰因
【申请人】国际壳牌研究有限公司