控制方法

控制方法
【专利摘要】提供了一种用于控制共振增强旋转钻机的方法，该共振增强旋转钻机包括钻进模块和用于控制钻进模块的一个或更多个钻进参数的控制系统，该方法包括以下步骤：(a)采用被钻材料的一个或更多个初始特性和/或一个或更多个初始钻进参数，以控制钻进模块；(b)测量一个或更多个当前钻进参数，以获得一个或更多个测量得的钻进参数；(c)采用一个或更多个测量得的钻进参数，作为控制系统中的输入，以便获得来自控制系统的输出，该输出包括被钻材料的一个或更多个计算得的特性；(d)采用被钻材料的一个或更多个计算得的特性和/或一个或更多个测量得的钻进参数，作为控制系统中的输入，以便获得来自控制系统的输出，该输出包括一个或更多个计算得的钻进参数；(e)将一个或更多个计算得的钻进参数可选择地应用于钻进模块；(f)可选择地重复步骤(b)、(c)、(d)以及(e)。
【专利说明】控制方法
[0001 ] 本发明设及冲击增强旋转钻进（percussion enhanced ;ro1:a;ry化illing)，并且 更具体地，设及共振增强钻进(RED:resonance enhanced化illing)。本发明的实施方式致 力于：用于控制共振增强旋转钻进W提高钻进性能的方法和设备，更具体地致力于使用模 糊逻辑控制的方法和设备。运里所描述的进一步实施方式致力于可W根据运些方法和设备 控制的共振增强旋转钻进设备。本发明的特定实施方式可应用于任意尺寸的钻机或要钻的 材料。特定的更具体的实施方式致力于钻过石油工业、气体工业W及采矿工业中的深孔钻 进应用中可能遭遇的岩层(具体是可变组成的岩层）。
[0002] 冲击旋转钻进本身已知。冲击旋转钻机拥有引起对旋转钻头的轴向载荷的机构。 该机构提供对被钻材料的冲击力，W便破碎帮助旋转钻头穿透材料的材料。
[0003] 共振增强旋转钻进是如下的特殊类型的冲击旋转钻进，其中，具有W高频振动的 振荡器，W便实现与被钻材料的共振。运导致有效使用旋转钻头处的机械能，由此，在与标 准冲击旋转钻进相比时提高钻进效率。
[0004] US3，990，522公开了一种冲击旋转钻机，该冲击旋转钻机使用用于钻螺栓孔的旋 转钻机中所安装的液压键。公开了可变行程和频率的冲击周期可W被应用并调整为被钻材 料的自然频率，W产生钻头的尖端处所施加的压力的放大。伺服阀维持冲击控制，并且转而 伺服阀由操作员借助经由电导体连接到伺服阀的电子控制模块来控制。操作员可W通过控 制增压流体到和从致动器的流动，从0至2500圈每分钟（即，0至42化）选择性地改变冲击频 率，并从0至1/8英寸（即，0至3.175mm)选择性地改变钻头的行程。描述了通过选择具有等于 被钻岩层的天然频率或共振频率的频率的冲击行程，岩层中由冲击力所存储的能量将导致 钻头的尖端处所施加的压力的放大，使得固体材料将巧塌并离开原位，并且允许钻率在3至 4英尺每分钟的范围内。
[0005] 存在用上述结构识别的并且在下面讨论的几个问题。
[0006] 使用US3,990,522的设备无法实现高频，该设备使用较低频率液压周期冲击器(描 述为振荡器，但并不在严格地技术意义上振荡）。因此，虽然US3,990,522讨论了共振的可能 性，但看来经由US3,990,522的冲击器可实现的低频不足W实现穿过许多硬材料的共振增 强钻进。
[0007] 不管上述的频率问题如何，在使用US3,990,522的布置的任意情况下无法容易地 实现并维持共振，特别是如果钻机穿过具有不同共振特性的不同材料时。运是因为US3， 990，522的布置中的冲击频率和行程的控制由操作员手动实现。由此可见，难W控制设备持 续调整冲击力的频率和行程，W随着钻机穿过不同类型的材料而维持共振。运对于如US3， 990,522所述的钻浅螺栓孔而言可能不是运种主要问题。操作员仅能对于内部要钻螺栓孔 的材料选择适当的频率和行程，然后操作钻机。然而，对于穿过许多不同层岩石恶化了问 题。位于深钻孔上方的操作员无法看见正在钻过什么类型的岩石，并且随着钻机从一个岩 石类型到另一个岩石类型穿过无法容易地实现并维持共振，特别是在岩石类型频繁变化的 区域中。
[000引上述问题中的一些已经被本发明人如WO 2007/141550所述地解决。本专利申请书 描述了一种共振增强旋转钻机，该共振增强旋转钻机包括自动化反馈和控制机构，该自动 化反馈和控制机构可W持续调整冲击力的频率和行程，W随着钻机穿过不同类型的岩石而 维持共振。钻机设置有：调整装置，该调整装置响应于钻机正在穿过的材料的情况;和井下 位置中的控制装置，该控制装置包括用于进行材料特性的井下测量的传感器，借此，设备可 在闭环实时控制下井下操作。
[0009] US2006/0157280提出了一种振荡器的井下闭环实时控制。描述了传感器和控制单 元可W在监测关键钻进效率参数(诸如前进率(ROP))的同时初始扫描频率的范围。然后，可 W控制振荡装置W提供最佳频率下振荡，直到进行下一频率扫描为止。频率扫描的模式可 W基于钻进操作的一个或更多个元素(诸如地层(formation)的变化、所测量ROP的变化、预 定时间段或来自地面的指令）。详细实施方式使用将扭转振荡应用于旋转钻头的振荡装置， 并且提及扭转共振。然而，还描述了应用于钻头的振荡的示例性方向包括跨所有自由度的 振荡，并且为了引发被钻材料的裂缝而不被使用。更确切的说，描述了钻头的旋转引起被钻 材料的初始碎裂，然后，为了确保旋转钻头保持与碎裂材料接触，应用瞬间振荡(momentary oscillation)。似乎没有提供如下的振荡器的任何公开或提议：为了如根据如WO 2007/ 141550所述的共振增强钻进所需的、引发旋转钻头穿过的材料的裂缝，该振荡器可W向钻 头输入足够高的轴向振荡载荷。
[0010] 现有技术没有提供关于如何监测轴向振荡的任何详情。US2006/0157280中和WO 2007/141550中总体公开了传感器，但未讨论运些传感器相对于部件(诸如隔振单元和振动 传递单元)的位置。
[0011] 本发明人的已公布国际申请书W02011/032874希望提供对现有技术中所描述的方 法和设备的进一步改进。该申请书公开了一种用于将共振增强旋转钻机中的振荡器的频率 控制在由Us(被钻材料的压缩强度)限定的范围内。该方法还设及将振荡器的动态力(Fd)控 制在由Us限定的范围内。Us使用传感器来测量。
[0012] WO 2007/141550、1152006/0157280^及胖02011/032874中所述的解决方案具体集 中于振荡器的振荡性质。运些仅是所设及的自由度的子集，并且运些方法不考虑Rm)的其他 重要方面。另外，材料特性的估计太简单，运依赖于使用传感器进行的压缩强度的非精确测 量，并且根本不解释其他材料特性。因此，期望对其中描述的方法和设备进行进一步改进。 本发明的实施方式的目的是为了进行如下改进:其在限制对设备的磨损W便增加设备的寿 命的同时，提高钻进效率。进一步目的是，特别是在钻过快速变化的岩石类型时更精确地控 制共振增强钻进。
[0013] 因此，本发明提供一种用于控制共振增强旋转钻机的方法，该共振增强旋转钻机 包括钻进模块和控制系统，该控制系统用于控制所述钻进模块的一个或更多个钻进参数， 该方法包括W下步骤：
[0014] (a)采用被钻材料的一个或更多个初始特性和/或一个或更多个初始钻进参数，W 控制所述钻进模块；
[0015] (b)测量一个或更多个当前钻进参数，W获得一个或更多个测量得的钻进参数；
[0016] (C)采用所述一个或更多个测量得的钻进参数，作为所述控制系统中的输入，W便 获得来自所述控制系统的输出，该输出包括所述被钻材料的一个或更多个计算得的特性；
[0017] (d)采用所述被钻材料的一个或更多个计算得的特性和/或所述一个或更多个测 量得的钻进参数，作为所述控制系统中的输入，W便获得来自所述控制系统的输出，该输出 包括一个或更多个计算得的钻进参数；
[0018] (e)将所述一个或更多个计算得的钻进参数可选择地应用于所述钻进模块；
[0019] (f)可选择地重复步骤(b)、（c)、（d似及（e)。
[0020] 在当前上下文中，共振增强旋转钻机不受特别限制，并且可W是领域中已知的任 意共振增强旋转钻机，只要它包括钻进模块和控制系统即可。钻进模块不受特别限制，并且 可W是用于共振增强钻进的领域中已知的任意钻进模块（诸如所公布的专利申请书 W02012/076401中所公开的钻进模块）。控制系统不受特别限制，只要它能够执行本发明的 方法即可。为了实现该目的，控制系统应当能够计算被钻材料的一个或更多个特性和一个 或更多个钻进参数。虽然特定具体方法确实提供明显优点，但用于进行该过程的方法不受 特别限制，并且将在下面详细描述方法。控制系统可W包括一个或更多个控制器。由此，在 一些实施方式中，可W采用单个控制器，W计算所有参数。然而，在其他实施方式中，可W有 两个或更多个控制器。在一些实施方式中，例如控制系统可W包括:用于计算被钻材料的特 性的控制器和；用于计算钻进参数的进一步控制器。一个或更多个或两个或更多个控制器 一起构成控制系统。
[0021] 发明人已经发现能够将钻进参数用于准确且精确地计算被钻材料的具体特性、然 后将运些特性连同测量得的钻进参数(实时钻进参数）用于计算可W应用于钻进模块的更 新后的钻进参数的过程。一个或更多个实时钻进参数(前进率、钻压(W0B:weight on bit) 等)在控制过程的各循环迭代中来测量(实时钻进参数随着被钻材料变化而变化，所W需要 保持实际钻进参数的跟踪的测量），并且如果必要(例如因为岩石类型变得更硬或更软），则 将新计算得的钻进参数应用于钻机。通常如果计算得的钻进参数中的一个或更多个不同于 所测量钻进参数中等效的一个或更多个(或在所选范围内充分不同），则可W进行将一个或 更多个计算得的钻进参数应用于钻进模块的步骤(e)。
[0022] 使用更新后的钻进参数重复控制周期将使钻进更有效。由此，虽然通常前进率 (ROP)或"钻进的速度"和/或钻压(WOB)为最期望优化的参数，但也可W使所期望的任意钻 进参数最大化(或最小化或优化）。
[0023] 之前，由于与普通钻进(钻机的振荡将额外的自由度引入到系统中）相比的更复杂 的RED性质（诸如与振荡关联的振荡频率和振幅W及动态力）使得任意预测或计算更加复 杂，所W控制方法在共振增强钻进(RED)过程中不是特别有效。已知方法尝试将振荡考虑在 内，但仍然需要改进。使用本方法，发明人能够将更多参数考虑在内，并且还能够更加有效 地将被钻材料的性质考虑在内。因此，发明人解决了许多现有问题，并且使得共振增强钻进 比已知系统中更快且更便宜，并且减少对材料的磨损。
[0024] 将仅参照附图经由示例描述本发明。
[0025] 图1示出了根据本发明的实施方式的钻进模块。
[0026] 图2的(a)和图2的（b)示出了对于被钻材料的各种振动质量和各种压缩强度将频 率例示为振幅的函数的图表。
[0027] 图3示出了给定固定动力供应的情况下对于各种振动质量将频率例示为振幅的函 数的图表。
[0028] 图4示出了本发明的控制系统的示意性表现。在该方案中，用于计算被钻材料的特 性的控制器被标记为"岩石识别机构"，用于计算钻进参数的控制器被标记为"控制器"，并 且如控制方法所采用的示出各种示例性材料特性和钻进参数。
[0029] 图5示出了用于计算或识别被钻材料的特性的控制器中用于信号识别的过程(例 如，岩石识别过程）。
[0030] 图6示出了人工神经网络(ANN)的神经元模型。
[0031] 图7示出了对于花岗岩和砂岩在特定钻进条件下生成的训练信号（位置和加速度 信号）。
[0032] 图8示出了对于ANN训练(位置和加速度信号)采用的各种测量。
[0033] 图9示出了 ANN训练过程的执行。
[0034] 图10示出了识别矩阵。
[0035] 图11示出了作为由阻尼器和弹黃支撑的无质量防护板的软冲击模型。
[0036] 图 12示出 了对于mi = lkg、h〇= lm、g = 9.8m/s2、ts = 9.5Ns/m、ks = 20 X l〇3N/m而计 算得的单个冲击的时间关系图表。质量mi的位移由上线来示出，并且质量的加速度由下线 来标记。
[0037] 图 13示出了具有对于mi = lkg、h〇 = lm、g = 9.8m/s2、ts = 9.5Ns/m、ks = 30 X l〇3N/m (左下点组）、ks = 32 X 103N/m(下一更高点组）、ks = 35 X 103N/m(下一更高点组）、ks = 40 X IO3NAi(下一更高点组）W及ks = 50X IO3NAi(下一更高点组)而计算得的各种刚度系数的冲 击图。
[003引图14示出了具有对于ml = lkg、g = 9.8m/s2、ts = 9.5化/m而计算得的各种刚度系数 和不同冲击速度的冲击图。
[0039] 图15示出了漂移振荡器的物理模型。
[0040] 图 16的（a)示出了对于a = 0.3、b = 0.16、C = 0.05、《=0.1、g = 0.02?及 cp=Jr/2 而 计算得的质量对时间的位移x(实线)和连接器底部对时间的位移V(虚线）。图16的（b)示出 了具有冲击面(直线）的相位图（曲线）。
[0041 ]图 17的(a)示出了对于a = 0. l、b = 0、C = 〇.〇5、O =0. l、g = 0.02W及:巧=巧。而计算 得的质量的位移(实线）、滑动器顶部的位移(虚线）W及滑动器底部的位移（点划线），图17 的(b)示出了对于a = 0. l、b = 0、C = 0.05、W =0. l、g = 0.02W及9=71/2而计算得的质量的加 速度。
[0042] 图 18示出 了对于a = 0. l、b = 0、C = 〇.〇5、O = 0. l、g = 0.02W及巧=征/2而计算得的 具有冲击面(垂直线）的质量加速度对相对位移的轨迹（曲线）。
[0043] 图 19示出 了对于a = 0. l、b = 0、C = 0.05、《=0.1、6 = 0.02、爭=瓦沒:、4 = 0.5(最右 线）、k=l(右中线及k = 5(左中线）而计算得的具有冲击面(最左线）的质量加速度对相 对位移的轨迹。
[0044] 图20示出了冲击振荡器的数学模型（参见Ing,J. ,Pavlovskaia,E. ,Wiercigroch, M.,Soumitro,B."Experimental study of impact oscillator with one-sided elastic constraint". Phi I. Trans.R.Soc. A,366(2008),679-705)。
[0045] 图21示出了模糊控制算法的示意图。
[0046] 图22示出了本模糊控制过程的流程图。
[0047] 图23示出了刚度斜率的模糊化的示意图。
[004引图24示出了刚度斜率的从属函数。
[0049] 图25示出了静态载荷b和激励幅度a的从属函数。
[0050] 图26示出了对于111=1、口: = 1、| = 0.05、《=0.1、邑=0.02、.(()=況/2、4 = 3、6 = 0.13、曰 = 0.2(下线）W及b = 0.1432、a = 0.3126(上线)而计算得的位移的时间关系图表。
[0051] 图27示出了对于111=1、口: = 1、| = 0.05、《=0.1、邑=0.02、啤=瓜'2''、4 = 3、6 = 0.13、曰 = 0.2(下线）W及b = 0.1432、a = 0.3126(上线)而计算得的平均前进的时间关系图表。
[0052] 图28示出了（a)0.75kN、（b)1.0kN、（c)1.2化N、（d)1.5kN、（e)1.75kNW及（f)2.0kN 静态载荷和2.化N动态力幅度时的位移的时间关系图表。（左顶线为常规钻进，而右底线为 R邸）。
[0053] 图29示出了（a)改进因子与(b)2.5kN动态幅度的前进W及静态力之间的关系。
[0054] 图30示出了（a)1.5kN、（b)2.5kN、（c)3.化NW及(d)5kN激励幅度和lkN静态载荷时 的位移的时间关系图表。（左上线为常规钻进，而右下线为RED)。
[0055] 图31示出了激励幅度和(a)穿透率W及(b)改进因子的变化。
[0056] 图32示出了（曰）150化、（6)180化、（。）190化、（(1)200化^及（6)210化激励频率和 1.75kN静态载荷时的位移的时间关系图表。（左上线为常规钻进，而右下线为RED)。
[0057] 图33示出了激励频率和(a)穿透率W及(b)改进因子的变化。
[005引 图;34不出了（a)30rpm、（b)60rpm、（。）90巧111^及（(1) 120rpm、200Hz激励频率和 1.75kN静态载荷时的位移的时间关系图表。（左上线为常规钻进，而右下线为RED)。
[0059] 图35示出了旋转速度和(a)穿透率W及(b)改进因子的变化。
[0060] 现在将更详细地描述本发明。
[0061] 如上所述，本发明提供了一种用于控制共振增强旋转钻机的方法，该共振增强旋 转钻机包括钻进模块和用于控制钻进模块的一个或更多个钻进参数的控制系统，该方法包 括W下步骤：
[0062] (a)采用被钻材料的一个或更多个初始特性，和/或一个或更多个初始钻进参数， W控制钻进模块；
[0063] (b)测量一个或更多个当前钻进参数，W获得一个或更多个测量得的钻进参数；
[0064] (C)采用一个或更多个测量得的钻进参数，作为控制系统中的输入，W便获得来自 控制系统的输出，该输出包括被钻材料的一个或更多个计算得的特性；
[0065] (d)采用被钻材料的一个或更多个计算得的特性和/或一个或更多个测量得的钻 进参数，作为控制系统中的输入，W便获得来自控制系统的输出，该输出包括一个或更多个 计算得的钻进参数；
[0066] (e)将一个或更多个计算得的钻进参数可选择地应用于钻进模块；
[0067] (門可选择地重复步骤化）、（(3)、（(1)^及(6)。
[0068] 在本发明的一些实施方式中，在步骤(d)中，除了测量得的钻进参数之外，还可W 采用来自控制过程的上次迭代的一个或更多个计算得的钻进参数，作为进一步输入，W确 定计算得的钻进参数。运可W例如在迭代中未测量所有钻进参数时发生。另外或另选地，在 计算得的材料特性的值不可用的情况下，反而可W使用来自上次迭代的计算得的值。
[0069] 在步骤(a)中，步骤"采用……W控制钻进模块"可W是指，将关于被钻材料的一个 或更多个初始特性的信息"应用"或"发送"到控制系统和/或钻进模块。因此，本发明还提供 一种用于控制共振增强旋转钻机的方法，该共振增强旋转钻机包括钻进模块和用于控制钻 进模块的一个或更多个钻进参数的控制系统，该方法包括W下步骤：
[0070] (a)将被钻材料的一个或更多个初始特性和/或一个或更多个初始钻进参数应用 于控制系统和/或钻进模块；
[0071] (b)测量一个或更多个当前钻进参数，W获得一个或更多个测量得的钻进参数；
[0072] (C)采用一个或更多个测量得的钻进参数，作为控制系统中的输入，W便获得来自 控制系统的输出，该输出包括被钻材料的一个或更多个计算得的特性；
[0073] (d)采用被钻材料的一个或更多个计算得的特性和/或一个或更多个测量得的钻 进参数，作为控制系统中的输入，W便获得来自控制系统的输出，该输出包括一个或更多个 计算得的钻进参数；
[0074] (e)将一个或更多个计算得的钻进参数可选择地应用于钻进模块；
[0075] (f)可选择地重复步骤化）、（(3)、（(1)^及(6)。
[0076] 在本发明中，钻进参数不受特别限制，并且可W是在控制钻进过程中有用的任意 参数。参数的数量也不受特别限制，只要采用至少一个。在一些实施方式中，参数的数量可 W是2个或更多、3个或更多、4个或更多、5个或更多、6个或更多、7个或更多、8个或更多、9个 或更多或者10个或更多。通常，钻进参数从W下参数选择：
[0077] (a)对被钻材料的轴向钻机力（也称为"钻压"（W0B:weight on bit)，或"静态 方'）;
[0078] (b)钻机(或钻头或钻进模块）的速率或速度，并且通常为穿过被钻材料的速率或 速度(也称为前进率(ROP);
[0079] (C)钻机(或钻头或钻进模块)的加速度，并且通常为穿过被钻材料的加速度；
[0080] (d)钻机的振荡频率；
[0081] (e)钻机的振幅；
[0082] (f)对被钻材料的振荡轴向钻机力(也称为"动态方'）；
[0083] (g)钻机的旋转速率或旋转速度；
[0084] 化)对被钻材料的、钻机的旋转力或扭矩；
[0085] (i)流体流率;W及
[0086] (j)钻头的相对位移。
[0087] 在该上下文中，流体流率包括穿过钻头、在钻头周围和/或在钻头上的流体量。通 常，流体穿过钻进区域，W便润滑、冷却和清除岩屑用途。在一些情况下，（诸如具有ROP的情 况），钻进参数无法直接应用于钻进模块，而是可W通过应用与间接参数有关的其他参数中 的一个或更多个来间接应用的参数。
[0088] 根据条件、钻进方法、材料(岩石)类型W及过程的目的，来使运些参数中的一个或 更多个可W最大化、最小化或优化。然而，通常对于最大化和优化，期望间接钻进参数(诸如 ROP)。例如，应当最小化的另一个运种参数可W为钻头磨损。
[0089] 在本发明中，被钻材料的特性不受特别限制，并且可W是控制钻进过程中有用的 任意特性。特性的数量也不受特别限制，只要采用至少一个即可。在一些实施方式中，特性 的数量可W是2个或更多、3个或更多、4个或更多、5个或更多、6个或更多、7个或更多、8个或 更多、9个或更多或10个或更多。通常，被钻材料的特性从W下特性选择：
[0090] (a)材料的压缩强度；
[0091] (b)材料的刚度或有效刚度；
[0092] (C)材料的屈服强度；
[0093] (d)材料的冲击强度；
[0094] (e)材料的疲劳强度；
[00M] (f)材料的拉伸强度；
[0096] (g)材料的剪切强度；
[0097] 化)材料的硬度；
[0098] (i)材料的密度；
[0099] (j)材料的杨氏模量;W及
[0100] 化)材料的泊松比。
[0101] 通常，步骤(a)中的被钻材料的一个或更多个初始特性从经验信息(优选地从数据 库)获得。进一步通常地，步骤(a)中的一个或更多个初始钻进参数从经验信息(优选地从数 据库)获得。数据库通常由在(各种尺寸、直径、重量等的）钻进模块对各种岩石类型和岩层 的测试期间获得的信息构成。
[0102] 步骤(C)中的被钻材料的一个或更多个计算得的特性可W使用计算或估计W任意 适当的方式来获得，但通常它们使用一个或更多个模型(诸如一个或更多个经验模型和/或 数学模型)来获得。模型不受特别限制，只要它们是用于计算讨论中的参数的有效模型（即， 它们能够进行所期望参数的合理估计）即可。通常，模型为经验模型（诸如从经验性数据而 不是从基于物理或数学的定理开发的模型）。优选地，一个或更多个模型包括神经网络。在 一些实施方式中，一个或更多个模型可W包括NARMAX(具有外因输入的非线性自回归移动 平均，Nonlinear AutoRegressive Moving Average with eXogenous inputs)模型和/或 基于冲击的模型和/或基于穿透的模型。若需要，所用的任意模型可W连同一个或更多个其 他模型一起使用。在一些情况下，几个不同模型一起或单独使用可W提供更准确的结果。
[0103] 在本发明中，步骤(d)中的一个或更多个计算得的钻进参数可W使用计算或估计 W任意适当的方式来获得，但通常它们使用模糊逻辑模型来获得。所采用的模糊逻辑模型 不受特别限制，并且可W采用领域中已知的任意模糊逻辑模型。通常，本发明中所用的模糊 逻辑模型包括：
[0104] (i)处理步骤(d)输入参数，W产生模糊输入参数；
[0105] (ii)处理模糊输入参数，W产生模糊结论；W及
[0106] (iii)处理模糊结论，W给出步骤(d)输出参数。
[0107] 在一些实施方式中，步骤（i)包括采用从属函数(membership化nction)，W将一 个或更多个输入参数转换成一个或更多个模糊输入参数。优选地，步骤(ii)包括推理机制 和规则库，其中，推理机制用于基于输入参数和/或步骤(i)的模糊输入参数，从规则库选择 要应用的一个或更多个规则，并且其中，一个或更多个规则应用于一个或更多个模糊输入 参数，W将一个或更多个模糊输入参数转换成一个或更多个模糊结论。通常但不唯一地，步 骤(iii)包括采用从属函数和/或重屯、法，W将一个或更多个模糊结论转换成步骤(d)输出 参数。
[0108]在本发明的一些优选实施方式中，被钻材料的特性的模型采用神经网络，并且用 于钻进参数的模型使用模糊逻辑。
[0109] 在步骤(b)中测量一个或更多个当前钻进参数(存在于应用所计算参数之前的过 程的当前迭代中的钻进参数）的方法不受特别限制。然而，优选地，运使用在钻进模块中所 并入的一个或更多个传感器来进行。所述一个或更多个传感器不受特别限制，只要它们能 够测量所设及的具体参数即可，但在优选实施方式中，传感器从测压元件（load cell)、满 流传感器、加速计、用于测量模块中电流和/或电压变化的传感器，W及运些传感器中的每 一个的一个或多个的任意组合选择。
[0110] 方法中还可W使用进一步的参数(诸如与钻进设备有关的、不会改变的固定参数 (重量、尺寸(钻头的直径、钻孔的直径、钻进模块的长度)等)和可W或可W不固定的外部参 数(诸如溫度、压力等））。通常，但不唯一地，该进一步参数从W下参数选择：
[011U (a)钻头的直径；
[0112] (b)钻进模块的质量；
[0113] (C)被钻材料的溫度；
[0114] (d)被钻材料的地静压力（geostatic pressure)或孔隙压力（pore pressure);
[0115] (e)钻进模块的长度；
[0116] (f)钻孔的宽度；W及
[0117] (g)钻孔的深度。
[0118] 本发明还提供一种钻进的方法，该方法包括W下步骤:操作共振增强钻机，W对材 料进行钻进，其中，共振增强钻机根据如上所定义的方法来控制。
[0119] 本发明还提供一种用于共振增强钻进的控制系统，该控制系统包括：
[0120] (a)用于确定被钻材料的一个或更多个特性的控制器，W及
[0121] (b)用于确定一个或更多个钻进参数的控制器；
[0122] 其中，控制系统被构造为执行如上所定义的方法。
[0123] 本发明更进一步提供一种用于共振增强钻进的设备，该设备包括：
[0124] (a)钻进模块;和
[0125] (b)如上所定义的控制系统。
[01%]共振增强钻进
[0127] 如上面所强调的，本发明设及共振增强钻进过程或设备的控制。本发明的方法可 W应用于常规R邸法(诸如W02007/141550和W02011/032874所公开的那些方法）。在运种已 知的方法中，通常基于将传感器用于直接确定被钻岩石的压缩强度的反馈机构，来控制振 荡器的动作的参数。在本方法中，如上所述，采用更加精细的反馈。
[0128] 通常，在已知Rm)法中，在操作期间，旋转旋转钻头，并由振荡器向钻头施加轴向定 向的动态载荷，W产生裂纹扩展区域，W帮助旋转钻头穿透材料。在已知方法中，通常如 W02011/032874所述来控制振荡器。在运些方法中，频率和动态力的范围基于下面提供的分 析。在本发明中，下面的细节代表作为更精细控制过程可W应用到的过程的指南。
[0129] 地层的压缩强度给出必需冲击力的下限。最小所需的动态力的幅度如下计算：
[0130]
[0131] Deff为作为根据与被钻材料接触的钻头的形态(fraction)而计算的钻头的直径D 的旋转钻义的巧满首巧-由此，有效直径Deff可W被定义为：
[0132]
[013引其中，Scnntact为对应于接触被钻材料的钻头的形态的换算系数。例如，估计出仅钻 头表面的5%与被钻材料接触，有效直径Deff可W被定义为：
[0134]
[0135] 上述计算提供振荡器的动态力的下限。使用比该下限更大的动态力在操作期间在 钻头前面产生裂纹扩展区域。然而，如果动态力太大，那么裂纹扩展区域将远离钻头延伸， 危害钻孔稳定性并降低钻孔质量。另外，如果旋转钻机上由振荡器给予的动态力太大，那么 可能产生加快的且灾难性的工具磨损和/或故障。
[0136] 参数
[0137] 如上面使清楚的，本控制过程的重要部分是，将测量得的钻进参数和模型化的材 料特性用作输入来进行钻进参数的建模。
[0138] 为了对各种钻进参数对R抓的影响进行充分建模，需要将钻进模块(特别是钻头） 和材料(例如，岩石）的相互作用考虑在内的物理模型。运种模型应当优选地不仅捕获材料 破碎过程，还模拟系统的钻进动力学。运可W如下所示例地完成，假如对材料的力学响应充 分建模。那么，运些响应可W在钻进参数计算控制器中，例如用于为模糊逻辑过程建立适当 的规则库和推理机制。
[0139] 与钻孔产生的任意其他方法一样，共振增强钻进(RED)的执行取决于所应用的输 入参数。例如，期望钻压(W0B)、旋转速度、激励的力幅W及频率对实现的常规钻进上的前进 率和百分比改进具有重大影响。对运些参数如何独立地影响R邸的结果的理解，对于确定最 优执行所需的组合是重要的。
[0140] 发明人针对砂岩执行基于实验确定的特性的模拟。钻头和剩余钻柱被认为是刚性 元件。
[0141] 与更早的实验研究相比，研究并观察对于W0B、旋转速度、激励幅度W及激励频率 的影响而获得的结果。
[0142] WOB的影响
[0143] 钻压(WOB)提供常规旋转钻进中的岩石破碎所需的力。然而，在R邸中，力为静态力 (或W0B)和动态力的组合。因此，需要理解WOB对过程的作用。为了研究WOB对R抓的执行的影 响，进行有限元模拟。保持材料特性、旋转速度、激励幅度W及激励频率恒定，而静态力可W 在最小值和最大值之间(诸如从0.75kN至2kN)变化。
[0144] 图28示出了针对各所用载荷的钻头的位移的时间关系图表。据观察，在模拟的常 规钻进阶段期间的时间关系图表的斜度(gradient)总是低于R抓期间经历的时间关系图表 的斜度。运定性地类似于实验室中的实验结果。
[0145] 在图29(a)中，描述了旭0的改进因子的变化。所获得的值对于所用的所有WOB还算 恒定。值得一说的是运不对应于实验。运是由于如下事实，实验中所用的真正岩石为非均质 的，而该模型中的岩石具有同质的特性。因此，对于改进因子期望获得还算恒定的值。图29 (b)是对于不同WOB的、WOB与前进之间的关系的图形表现。据观察，前进随着增大的WOB而增 加，直到1.75kN(其之后观察到下降）为止。
[0146] 激励幅度的影响
[0147] 当RED中钻头冲击岩石时可用的能量为谐波激励的幅度函数。然而，由于可用于 R抓致动器的功率，所W存在对可达到的最大值的限制。由此，有用的是建立其如何影响钻 进过程的整个执行。为了研究幅度的变化的作用，在保持其他WOB(静态力)和其他钻进参数 恒定的同时，W不同的力幅度进行=个模拟。
[0148] 图30呈现了常规钻进和R抓期间钻头的位移的时间关系图表。注意的是，在R抓阶 段期间获得的平均穿透率总是大于在模拟的常规钻进阶段期间实现的平均穿透率。
[0149] 在图31(a)中，呈现了前进率对激励幅度的图表。前进率随着幅度而增大，改进因 子也是如此(图31(b))。观察结果符合从实验获得的结果。
[0150] 激励频率的影响
[0151] 动态力的激励频率确定每单位时间内可用于岩石破碎的能量率。因此，理解频率 如何影响RED的效率将帮助选择钻进操作中要使用的最优值。
[0152] 频率对执行的影响也使用相同有限元模型来研究。保持其他模型参数恒定的同时 改变频率的值。从150Hz、180Hz、190Hz JOOHzW及210Hz进行模拟。在图32中的钻头的位移 的时间关系图表方面来描述结果。
[0153] 图33(a)示出了频率与穿透率之间的关系，而在图33(b)中用图表呈现了对于不同 速度获得的改进因子。对于所用范围，穿透率随着增大的频率而增大。还据观察，图表的斜 度在超过180Hz时变得明显更高。运提示如果进行更高频率的模拟，则可W发现最优执行所 需的最佳值。对于如图33(b)所示的改进因子获得类似的结果。在R抓实验中观察运些类型 的结果。
[0154] 旋转速度的影响
[0155] 钻柱的旋转速度有助于钻进过程的效率。已知在常规旋转钻进中，旋转速度在岩 石破碎和钻屑移除中起到重要作用。有用的是，研究旋转速度对R邸的执行的影响。
[0156] WSOrmp的梯级（St邱)对从30rpm至12化pm变化的旋转速度进行研究。分析的结果 在钻头的位移的时间关系图表方面用图表呈现在图34中。R抓阶段的平均穿透率总是超过 常规旋转钻进的平均穿透率。
[0157]根据图34和图35(a)，还注意的是，钻进组件对于旭0和常规钻进运两者随着增加 的旋转速度而更快地移动。图35(b)是对于各种速度而获得的改进因子的图形简况。观察到 改进因子对于所应用的所有旋转速度保持还算恒定。运些结果符合从实验获得的结果。
[015引参数影响的总结
[0159] 具体地，示出了WOB和动态力运两者有助于岩石破碎所需的力。由此可见，WOB和动 态力各自的量影响钻进过程的效率。WOB被示出具有用于要实现的最佳性能的最优值。增大 所产生的激励幅度增大R邸的穿透率和改进因子。
[0160] 数据示出了激励频率在可用于岩石破碎的每单位时间内的能量的量中起到了重 要作用。另外，旋转速度被示出为确定从钻头的面的材料移除率。穿透率被示出为随着增加 的旋转速度而增大。
[0161] 从所进行的有限元模拟获得的所有结果定性地类似于从R邸实验观察的结果。
[0162] 总之，为了对各种钻进参数对R抓的影响充分建模，开发了考虑钻进模块(例如，钻 头)和材料(例如，岩石）相互作用的物理(或经验性)模型。运些模型被示出为，适于确定岩 石和其他材料的最佳钻进参数，并且可W用于作为本发明的控制过程的部件的控制器中。 在用于计算钻进参数的模糊逻辑控制器中，模型可W用于建立适当的规则库和推理机制， 或者可W用于其他控制策略中。
[0163] 控制系统
[0164] 图4中示出了根据本发明的控制系统的示例。在该方案中，用于计算被钻材料的特 性的控制器被标记为"岩石识别机构"，用于计算钻进参数的控制器被标记为"控制器"，并 且如控制方法中所采用的示出了各种示例性材料特性和钻进参数。
[0165] 用于计算钻进参数的控制器
[0166] 由发明人进行的模拟研究展示了，给出最优穿透率的最佳激励幅度和WOB根据岩 石刚度而变化。本发明的模糊控制策略被优选为，对于不同的岩石(或其他材料)条件调节 RED钻进模块的最佳频率和幅度(或其他钻进参数）。
[0167] 图21呈现了所提出方法的示意图，该方法含有模糊化、推理机制、规则库W及去模 糊化。模糊化是运样的过程:将清晰（crisp)输入转换成语言术语，然后基于预定义的规则 库将它们置于推理机制中。稍后，将从推理机制获得的语言决策转换成用于漂移振荡器的 清晰控制参数。如图21所示，本控制系统具有至少一个输入(例如刚度斜率，如附图中的）和 至少一个输出（附图中为两个输出，在该示例中为激励幅度和静态力）。
[0168] 图22中给出了本控制过程的流程图。如从该图所示，该过程W初始化开始，并且一 旦模拟开始则记录针对一个推进时期(period of forcing)的数据。然后，模糊系统评估当 前控制参数(b和a)是否为可W在当前刚度条件下产生最优效率的最优控制参数。如果是， 则系统将继续使用当前控制参数;如果不是，则将采用由模糊控制系统计算出的清晰控制 参数(更新后的、或计算得的钻进参数）。
[0169] 模糊化
[0170] 模糊控制器具有至少一个输入信号(在该示例中诸如为刚度斜率）。推理机制将基 于模糊后的输入选择控制规则。图23中呈现了刚度斜率的模糊化。如图23中所看出的，大 于-0.5024的刚度斜率定义为"软"，-0.5024至-5.0237之间的刚度斜率定义为"中等"，并且 小于-5.0237的刚度斜率定义为"硬"。表4中给出了模糊化的简介，并且图24中给出了刚度 斜率的从属函数。
[0171] 表4-刚度斜率的模糊化
[0172]
[0173] 规则库
[0174] 因为对于该示例假定静态载荷固定，所W模糊控制器的目的可W被转换成，发现 (例如）产生最优平均前进的最佳激励幅度。如下给出用于该示例的基于马丹尼模型 (Mamdan i mo d e 1)的模糊控制规则：
[0175] ?如果刚度为软，则WOB为大，并且激励幅度为小；
[0176] ?如果刚度为中等，则WOB为中等，并且激励幅度为中等；
[0177] .如果刚度为硬，则WOB为小，并且激励幅度为大；
[0178] 可W按需应用用于参数之间的其他关系的其他规则。当存在几个输入时，可W对 于不同的输入参数应用不同的规则。类似地，当需要优化不同的参数时，可W应用进一步不 同的规则。推理机制可W适于确保正确的规则与所采用的输入和所优化的参数匹配。
[0179] 去模糊化
[0180] 图25示出了用于来自推理机制的模糊结论的去模糊化的静态载荷和激励幅度的 从属函数。特别优选的是，可W使用如下的重屯、法，借助给重屯、法，使用各隐含模糊集合的 总面积的中屯、来计算清晰输出。
[0181] 数值结果
[0182] 发明人通过使用W下参数来进行对所提出模糊控制器的模拟研究:m=l、Pf = l、写 =0.05、《=0.1、邑=0.02^及9=巧/2。图26中示出了使用4 = 3、6 = 0.13^及曰=0.2的、由下 面的一组线来指示的模拟结果。模糊控制器指示最佳参数为b = 0.1432和a = 0.3126,并且 图26中示出了使用最佳参数的、由上面的一组线指示的模拟结果。可W清晰的看出，显著增 大使用最佳参数的位移，并且可W在图27中看出使用最佳参数的改进后的平均前进。
[0183] 用于计算被钻材料的特性的控制器
[0184] 现在将更详细地描述用于计算被钻材料的特性的控制器。当前的控制系统和方法 可W用于许多类型的钻进中，因此，不限于大型钻进(诸如钻过岩石）。然而，通常，本发明很 好地适于大型钻进(诸如采矿工业和石油工业中的岩石钻进）。因此，被钻材料通常为一种 类型的岩石。然而，在W下内容中，对岩石的提及也可W当作适用于其他材料。
[0185] 当材料为岩石或任意其他材料时，可W在数学优化方面用公式来表达材料识别问 题:给定所观察的信号x(t)(从钻进测量获得的），有用的是发现模型函数y(t，p)(p为模型 的参数），使得误差：
[0186]
[0187] 为最小。为了解决该问题，数个算法(例如图5中所示的基于反分析的算法)可用。 然而，运种算法的成功应用需要存在如下的鲁棒数学模型：其W可W预测信号y(t，p)的运 种方式来描述要被研究的现象。
[0188] 之前控制系统的问题和本发明潜在的问题是，迄今在物理学方面没有可用于描述 经受钻进过程的岩层的具体特性的准确模型。因此，发明人确定应用不同的方法。由发明人 所探索的另选途径包括W下经验模型：
[0189] ?遗传算法-基于一代又一代的遗传库过渡，源于达尔文的进化论。
[0190] ?蚁群算法-基于昆虫的行为。
[0191] ?禁忌研究(Tabu research)-探索所期望临界值的本地邻域。
[0192] ?粒子群优化-基于组(诸如天空中飞翔的鸟或建筑物内部寻找出口的人)的社会 行为。
[0193] .人工神经网络-受模式识别问题欢迎，其中，该方法模拟中枢神经系统的学习过 程。
[0194] 所有运些方法为系统提供物理理论的有用另选方案。总之，本发明中优选的是采 用运种另选经验模型。运可W是或可W不是数学模型。
[0195] 用于识别材料(或岩石)特性的其他潜在有用的模型包括基于冲击的识别、基于穿 透的识别W及NARMAX模型。
[0196] 基于神经网络的识别
[0197] 在本发明中，虽然可W采用上述经验模型中的任意一个(或任意其他适当模型）， 但对于解决材料特性(岩石特性)识别问题而言，人工神经网络(ANN)是最优选的。虽然人工 神经网络对于其他参数也同样有用，但它们在材料的刚度是要采用的一个参数时特别有 用。Ar^N呈现不需要讨论中的现象的显式数学模型（或物理学模型）的显著优点。Ar^N的关键 点是，对于所选网络进行适当的学习过程。图6中示出了 Ar^N的神经元模型。
[0198] Ar^N为对生物神经系统的操作(例如，学习、分类W及模式识别)进行建模的一般数 学计算范式。中屯、思想是定义如下的合适特征向量(措施），其捕获要被研究过程的特性(在 运种情况下为钻进条件下的材料(岩石)特性）。一旦定义了一组特征向量，则Ar^N经受学习 过程，在该过程中，为了正确识别岩石特性，将调整其系数。在该阶段，重要的是，W训练 试样代表实际应用中所期望的钻进条件的运种方式，具有足够大范围的训练试样。
[0199] 在本发明中，存在可W用于MN识别过程的若干可用信号。具体地，为钻头的位置 和加速度。运允许定义对于位置和加速度信号运两者而采用的各种特征测量(诸如动力、峰 值、傅里叶系数等）。运可W由实验来说明。考虑W下钻进条件下产生的一组训练信号：
[0200] 表1用于训练的实验数据的参数 「0勺011
LUAJ^」 刈 于化冈若和做'若，凶丫甲不出J上巧巧巧杀件h广生的训珠信号U亞置和观化巧 信号）。
[0203] 在图8中，示出了对于Ar^N训练而采取的各种测量。该方法的一个关键方面是，采取 能够区分感兴趣的特性的测量，该感兴趣的特性在该示例中为岩石类型(花岗岩/砂岩）。如 可W在图中看出的，对于各个岩石类型，所选测量有效地累积，由此，给出它们识别什么类 型的岩石正被钻的能力的指示。
[0204] 凭借运些测量，可W进行选择监督式学习作为实施方案的学习算法(Matlab A順 工具箱所用于的）。在图9中可W观察该过程的执行。如可W看出的，算法收敛非常快，运意 味着潜在问题是适定的(well-posed)。因为问题的适定性与识别技术的鲁棒性密切相关， 所W运是要考虑的非常重要的方面。运鉴于必须仔细识别岩层的可能苛刻条件而是重要 的。在运点上，该方法示出了非常好的收敛特性，运是该技术为计算实际钻进条件的很大改 进的指示。
[0205] 一旦已经训练ANN，则为了验证它的识别能力，可W测试网络。运可W通过W所有 过程用同组的实验数据来进行的运种方式将试样数据划分为两组来进行：训练组和测试 组。然而，为了说明所提出的识别技术的鲁棒性和有效性，如下面表2中所详细说明地使用 单组实验数据。
[0206] 表2用于的1试的实验数据的参数 「rm7l
[0208] 如可W看出的，测试数据的参数不同于训练数据的参数，但当然它们合理地类似。 由此，使用该新组实验数据可W把位置和加速度当作主要信号，并获得之前为了测试所训 练网络而定义的测量。
[0209] 图10中可W看出结果。运是所谓的混淆矩阵，其示出了试样的绝对数及其对应于 正确和非正确识别信号的百分比。在运种情况下，可W看出，所有试样累积在混淆矩阵的对 角线中，运意味着对应于种类1(花岗岩)的试样已经被网络识别为属于种类1等。因此，为了 实施图5中所示的岩石识别控制器，可W应用并扩展该方法。
[0210] 基于冲击的识别
[0211] 也可W使用基于钻头的加速度的方法来实现材料特性的识别。图11中所示的简单 冲击模型，其中，质量mi与被建模为轻防护板的基础碰撞，该轻防护板由具有恢复系数ks的 轻弹黃和具有阻尼系数ts的轻粘性阻尼器支撑。模型的动力学由W下方程式来管理：
[0219] 在质量虹-防护板接触开始时，并且在时间Ti之后的质量虹-防护板接触结束时，振 荡器(IUh干平锥术杰.所W接触在等于振荡周期的一半的时间期间发生：
[0212]
[0213]
[0214]
[0215]
[0216] 则获得I
[0217]
[021 引
[0220]
[0221] 图12呈现了在与防护板第一次冲击时质量mi的等效响应。如可W从图看出的，对 于第一次冲击标记冲击持续时间Tl和冲击幅度出。。，并且运两个重要参数用于构造可W识 别岩石的特性的冲击图.
[0222] 如果保持控制参数恒定，对于不同类型的岩石，冲击持续时间和冲击幅度不同。对 于软型岩石，冲击持续时间更长，并且冲击幅度更小。对于硬型岩石，冲击持续时间更短，并 且冲击幅度更大。
[0223] 图13示出了使用相同控制参数的具有各种刚度系数的冲击图。当刚度增大(从右 底部至左顶部)时，冲击持续时间变得更短，并且冲击幅度变得更大。所W，通过在该图中监 测冲击的位置来识别岩石的刚度是明确的。
[0224] 为了更加理解冲击动作，图14示出了具有各种刚度系数和不同冲击速率的冲击 图。发现，具有相同刚度系数的情况下，冲击持续时间在冲击速率不同时保持不变。还可W 推断的是，通过应用更大的冲击速度而产生的岩石的刚度的变化（由最上圈标记的），比通 过应用更小的冲击速率而产生的变化（由最下圈标记的)更明显。因为该现象展示了在共振 条件下操作R邸的要求，所W它是有趣的。
[0225] 基于穿透的识别
[0。6]冲击识别与该方法之间的主要区别是：前者通过使用加速单个（acceleration single)的冲击持续时间和冲击幅度来考虑冲击动作的几何学，而后者考虑加速单个的冲 击幅度和钻头到岩石中的穿透。
[0227] 考虑如图15所示的具有漂移的冲击振荡器(从化vlovskaia E. ,Wiercigroch M., Grebogi C.:(2001),Modelling of an impact system with a drift,Phy Rev E,64, 056224改编），其中，质量m由含有幅度Pd、频率Q W及相移夺的谐波分量和静态分量Ps的外 力来驱动。假定在初始时刻t = 0时，间隙G是质量与无重量滑动器顶板之间的距离，该无重 量滑动器顶板经由具有刚度k的线性弹黃和具有阻尼系数C的黏性阻尼器连接到滑动器底 板。Xm、Xt、Xb分别表示质量、滑动器顶部W及滑动器底部的绝对位移。假定重力被包括在静 态力中，或者模型在水平面中操作。
[0228] 所考虑的系统被写到一组一阶微分方程：
[0229]
[0230]
[0231]
[0232] 其中，H(.)是亥维赛化eaviside)函数，其被描述为：
[0233] Pi = Pi(x,z) =H(x-z-g),
[0234] P2 = P2(z，z'，v) =H(2|z'+z-v)，
[0235] P3 = P3(z，z'，v)=H(2|z'+z-v-1)，
[0236] P4=P4(v' ) =H(v')，。
[0237] 假定未知线性弹黃的刚度k，则该方法的目的是，估计模拟钻进过程期间岩石的刚 度的漂移振荡器的刚度。
[0238] 图16的(a)中示出了所研究系统的位移的时间关系图表，实线表示质量的位移，并 且虚线表示滑动器底部的位移。如可W看出的，系统在周期-1响应中，其中每推进周期质量 具有与滑动器顶部的一个碰撞。在化vlovskaia E. ,Wiercigroch M. ,Grebogi C. : (2001), Modelling of an impact system with a 化ift,Phy Rev E,64,056224中，发现如图 16的 (b)所示使用质量速度y对相对位移X-V的相位图给出该冲击行为的深刻理解。然而，从相位 图，难W估计与滑动器顶部的冲击期间的质量的加速度。该冲击行为用于模拟钻头与岩石 之间的冲击。如果在冲击期间识别岩石的刚度，则为了实现最佳的穿透率，可W应用适当的 动态力和静态力。本发明提供了一种如下的方法，该方法用于通过使用质量加速度和相对 位移构建相位图而进行刚度识别。
[0239] 图17中的数字结果示出了，仅具有冲击的系统的时间关系图表，图17的(a)呈现了 质量、滑动器顶部W及滑动器底部的位移，并且图17的（b)呈现了冲击期间质量的加速度。 图18示出了质量的加速度对相对位移的轨迹。图17和图18中标记了S个临界点A、BW及C， 它们分别指示冲击的开始、相对位移的最大值W及冲击的结束。如可W看出的，从点A开始， 质量开始接触滑动器顶部，并和滑动器顶部一起向下移动，并且来自滑动器顶部作用在质 量的反作用力立刻增大。一旦质量由反作用力停止，则质量的加速度沿着向上方向实现最 大值，并且质量与滑动器底部之间的相对位移也最大。如果反作用力更大，则沿向上方向的 加速度的最大值更大，并且相对位移的最大值更短。因此，我们可W推断质量的加速度和相 对位移的最大值取决于线性弹黃的刚度系数。图19通过使用=个不同的刚度系数模拟系统 来展示了上述结论。看出的是，具有刚度系数k = 0.5的响应示出了最大相对位移和最小质 量加速度，而具有刚度系数k = 5的响应示出了最小相对位移和最大质量加速度。
[0240] 基于NARMAX的识别
[0241] 已经广泛研究了经受正弦输入的分段线性振荡器的行为（参见Wiercigroch M.， Wojewoda J.,Krivtsov A.M.:(2005),Dynamics of ultrasonic percussive drilling of hard rocks,J.Sound and V;Lbration,280,739-757)。然而，大多数研究集中于，表示在 经由分支图（bifurcation map)改变输入幅度或频率时观察的复杂动力学体制的特征(参 见Ing,J. ,Pavlovskaia,E. ,Wiercigroch,!. ,Soumitro,B. ('Experimental study of impact oscillator with one-sided elastic constraint''.Phil.Trans .R.Soc.A,366 (2008) ,679-705)。虽然该信息提供了关于工作在已知环境中时的钻机响应的有用桐察，但 是它无法用于识别与物理上有意义的刚度变化关联的新场景。发明人通过使用非线性系统 建模和分析来解决该问题。
[0242] 该方法包括使用来自图20所示的冲击振荡器系统的输入-输出数据建立多项式 NARMAX(具有外因输入的非线性自回归移动平均，Nonlinear AutoRegressive Moving Average with eXogenous inputs)模型（参见Leonaritis , I. J. ,Billings, S.A. ('Input-output P曰rametric models for nonlinear systems,Part I:deterministic nonlinear systems". International journal of Control,41(1985),303-328)。然后，使用基于 ALE (关联线性方程）的新型通用框架（参见Fei joo ,J.A.V. ,Worden,K.和Stanway , R . "Associated linear equations for Volterra operators''.Mechanical Systems and Signal processing, 19(2005) ,57-69)，在感兴趣的特定频率范围内，通过提取系统NOFRF (非线性输出频率响应函数，Nonlinear 0u1:put Frequen巧Response F^mctions)在频域 中分析所识别模型的响应（参见Lang,Z.Q. ,Billings,S.A/'Ene;rgy hansfer properties of non-linear systems in the frequency domain". International journal of Control ,78(2005) ,345-362)。已经进行模拟研究。结果展示了基于简单NOFRF能量的指数 可W明显区分刚度变化的不同场景，W便能够用于监测系统的操作条件。
[0243] 考虑冲击振荡器的无量纲模型：
[0244]
[0245] 所乂用的参数力
庚中，m二 1, g二0.00126, fn二9.38, C二 1.3, Xo二 0.001。如勺不同值用于不同的模拟中。测试W下值:0=1、5、6、8、l〇W及20。
[0246] 系统使用正弦输入来识别：
[0247]
[0248] 其中，A = O.005,并且f = 8。因为运些值对于所研究的e的范围产生可区别的简单 (没有分支)输出，所W选择运些值(更小的频率产生太相似的输出，更大的频率产生难W识 别的动力学体制）。
[0249] 系统凭借固定Ts = O.06集成步骤使用第四阶龙格-库塔(Runge-Kutta)求解程序 来求积分。各模拟运行，直到Tf = 294.62巧000点）为止。
[0250] 具有2个试样的最大延时和多达第S度的非线性的NARMAX模型使用基于ERR的正 交估计器来识别。用于移动平均数(MA)和非MA项的阔值为1(T 8。通过将各NARMAX模型分成对 应的关联的线性方程式(ALE)而对于正弦输入估计NOFRF:
[0251]
[0巧2] 第n阶NOFRF通过划分由11。的DFT输出的第n阶ALE的DFT来获得。运些估计对于作为 正弦输入的频率范围的1 <k< 1000是可靠的。
[0253 ] 允许不同情况之间的区分的测量为第n阶NOFRF的能量：
[0 巧 4]
[0巧5] 其中，Gn化)为频率分量k下的第n阶NOFRF。
[0256] 表3中呈现了示出了指数对于不同的刚度值如何变化的识别结果。
[0257] 表3-各种刚度值的识别结果 「0巧引 L U之。7」 t 口化佚W
[0260] 本发明之前所述的实施方式可应用于任意尺寸的钻机或被钻材料。特定更具体的 实施方式致力于钻过石油工业、气体工业W及采矿工业中的深孔钻进应用中可能遭遇的岩 层(具体是可变构成的岩层）。本发明有助于确定适于钻过运种岩层的各种参数的数值。
[0261] 在构造用于与本发明一起使用的钻进模块时，将所设想的钻进类型中设及的条件 考虑在内。例如，用于石油钻进的钻进模块很可能与用于建筑工业的钻进模块显著不同。当 在岩层中钻进时，岩层的压缩强度具有大变化(从砂岩的大约Us = 70MPa-直到花岗岩的Us = 230MPa)。在大型钻进应用(诸如石油工业）中，钻头直径范围从90mm到800mm(3i/2到32")。 只要钻头表面的大约5%与岩层接触，那么所需动态力的最低值可W被计算为大约20kN(使 用90mm钻头穿过砂岩）。类似地，所需动态力的最大值被计算为大约GOOOkN(使用800mm钻头 穿过花岗岩）。由此可见，对于钻过岩层，根据钻头的直径，动态力优选地被控制维持在范围 20kN到6000kN内。因为驱动具有6000kN动态力的振荡器将消耗大量的动力，所W有利的是， 将具有中到小直径钻头的本发明用于许多应用。例如，90mm至400mm的钻头直径产生20kN至 ISOOkN的操作范围。进一步缩小钻头直径范围给出了20kN至1000 kN的动态力的优选范围， 更有选地为40kN至500kN，仍然更优选地为50kN至300kN。
[0262] 对于需要的振动位移幅度的更低估计是:具有比来自由于岩层中的非均质性而产 生的随机小型尖端反弹的位移明显更大的振动。运样，该振幅有利的为至少1mm。因此，振荡 器的振幅可W维持在范围Imm至1 Omm内，更优选地为1mm至5mm。
[0263] 对于大型钻进设备，模块的振动质量可W处于IOkg至1000 kg级别。运种大型钻进 设备的可行频率范围不伸展到高于几百赫兹。由此可见，通过选择之前描述的范围内的钻 头直径、振动质量W及振幅的适当值，振荡器的频率(f)可W在提供对于不同岩石类型的范 围产生裂缝扩展区域的足够动态力的同时被控制维持在范围IOOHz至500化内，并且该频率 为实现共振效应的足够高的频率。
[0264] 图2的(a)和图2的(b)示出了对于具有150mm直径的钻头将频率例示为振幅的函数 的图表。图表(a)适合于振动质量Hi= 10kg，而图表(b)适合于振动质量m=30kg。下曲线对于 更弱岩层有效，而上曲线适合于具有高压缩强度的岩石。如可W从图表看出的，曲线上方区 域中IOOHz至500Hz的工作频率将提供足够高的频率，W使用范围Imm至10mm(0.1 cm至Icm) 内的振幅在所有岩石类型中产生裂缝扩展区域。
[0265] 图3示出了给定固定动力供应的情况下对于各种振动质量将频率例示为振幅的函 数的图表。对于可W由用于驱动钻头的旋转动作的泥浆马达或满轮在井下产生的30kW的动 力供应，而计算该图表。上曲线适合于IOkg的振动质量，而下曲线适合于50kg的振动质量。 如可W从图表看出的，100化至500化的频率范围可用于范围Imm至10mm(0.1 cm至Icm)内的 振幅。
[0266] 可行的是在地面上设置计算机，该计算机处理来自钻孔下面的传感器的信号，然 后，将控制信号发送回钻孔下面，W便控制钻头。然而，因为地面与钻孔的底部之间的信令 不是直通的，并且还可能非常缓慢，所W运在实践中对于深钻孔钻进而言将难W实现。另选 地，可W将反馈机构的感测、处理W及控制元件定位于钻孔的下面但在钻头组件的外部。然 而，在实践中，钻孔下面可能具有小空间，同样地，该机构可能经受杂乱物理条件。
[0267] 因此，提供反馈控制的最优结构是使反馈机构的所有感测、处理W及控制元件定 位于井下组件内（例如，在钻进模块内）。该结构最紧凑，提供对共振条件的变化的更快的反 馈和更快速的响应，并且还允许钻进模块用集成在钻进模块内部的必需反馈控制器来制 造，使得模块可W被改造成适配现有钻柱，而不需要替换整个钻进系统。
[0268] 对于大型钻进设备，振荡器有利地包括具有机械放大的压电致动器、磁致伸缩致 动器、气动致动器、或电驱动机械致动器。发现，运些致动器可W实现所期望的频率、动态 力、振幅W及动力消耗范围，W便与之前所述的方法一起使用。
[0269] 本发明的实施方式的用途包括:钻井(例如，油井钻进）；采矿(例如煤、钻石等）；表 面钻进(例如，道路工程等）；W及手持钻机(例如，家庭用途的DIY钻机、牙钻等）。
[0270] 本发明的实施方式的优点包括:提高的钻进速度;更好的钻孔稳定性和质量;导致 更长寿命的、对设备更小的应力；W及降低能量成本的更大效率。
【主权项】
1. 一种用于控制共振增强旋转钻机的方法，该共振增强旋转钻机包括钻进模块和控制 系统，该控制系统用于控制所述钻进模块的一个或更多个钻进参数，该方法包括以下步骤： (a) 将被钻材料的一个或更多个初始特性和/或一个或更多个初始钻进参数应用于所 述控制系统和/或所述钻进模块； (b) 测量一个或更多个当前钻进参数，以获得一个或更多个测量得的钻进参数； (c) 采用所述一个或更多个测量得的钻进参数，作为所述控制系统中的输入，以便获得 来自所述控制系统的输出，该输出包括所述被钻材料的一个或更多个计算得的特性； (d) 采用所述被钻材料的所述一个或更多个计算得的特性和/或所述一个或更多个测 量得的钻进参数，作为所述控制系统中的输入，以便获得来自所述控制系统的输出，该输出 包括一个或更多个计算得的钻进参数； (e) 将所述一个或更多个计算得的钻进参数可选择地应用于所述钻进模块； (f) 可选择地重复步骤(b)、（c)、（d)以及(e)。2. 根据权利要求1所述的方法，其中，在步骤(d)中，采用来自所述控制过程的上次迭代 的一个或更多个计算得的钻进参数作为进一步输入，以确定所述计算得的钻进参数。3. 根据权利要求1或2所述的方法，其中，所述钻进参数包括以下各项中的一个或更多 个： (a) 对所述被钻材料的轴向钻机力，也称为"钻压WOB"，或"静态力"； (b) 穿过所述被钻材料的钻头和/或钻进模块的速率或速度； (c) 穿过所述被钻材料的所述钻头和/或钻进模块的加速度； (d) 所述钻头和/或钻进模块的振荡频率； (e) 所述钻头和/或钻进模块的振幅； (f) 对所述被钻材料的振荡轴向力，也称为"动态力"； (g) 所述钻机的旋转速率或旋转速度； (h) 对所述被钻材料的、所述钻机的旋转力或扭矩； (i) 流体流率；以及 (j) 所述钻头的相对位移。4. 根据任意前述权利要求所述的方法，其中，所述被钻材料的所述特性包括以下各项 中的一个或更多个： (a) 所述材料的压缩强度； (b) 所述材料的刚度或有效刚度； (c) 所述材料的屈服强度； (d) 所述材料的冲击强度； (e) 所述材料的疲劳强度； (f) 所述材料的拉伸强度； (g) 所述材料的剪切强度； (h) 所述材料的硬度； (i) 所述材料的密度； (j) 所述材料的杨氏模量；以及 (k) 所述材料的泊松比。5. 根据任意前述权利要求所述的方法，其中，步骤(a)中的所述被钻材料的所述一个或 更多个初始特性从经验信息获得，优选地从数据库获得。6. 根据任意前述权利要求所述的方法，其中，步骤(a)中的所述一个或更多个初始钻进 参数从经验信息获得，优选地从数据库获得。7. 根据任意前述权利要求所述的方法，其中，步骤(c)中的所述被钻材料的所述一个或 更多个计算得的特性使用一个或更多个模型来获得，优选地使用一个或更多个经验模型 和/或一个或更多个数学模型来获得。8. 根据权利要求7所述的方法，其中，所述一个或更多个模型包括神经网络。9. 根据权利要求7或8所述的方法，其中，所述一个或更多个模型包括NARMAX模型和/或 基于冲击的模型和/或基于穿透的模型。10. 根据任意前述权利要求所述的方法，其中，步骤(d)中的所述一个或更多个计算得 的钻进参数使用模糊逻辑模型来获得。11. 根据权利要求10所述的方法，其中，所述模糊逻辑模型包括： (i)处理步骤(d)输入参数，以产生模糊输入参数； (i i)处理所述模糊输入参数，以产生模糊结论；以及 (iii)处理所述模糊结论，以给出步骤(d)输出参数。12. 根据权利要求11所述的方法，其中，步骤（i)包括如下步骤:采用从属函数，以将所 述一个或更多个输入参数转换成一个或更多个模糊输入参数。13. 根据权利要求11或12所述的方法，其中，步骤(ii)包括推理机制和规则库，其中，所 述推理机制被用于基于步骤（i)的所述输入参数和/或所述模糊输入参数从所述规则库选 择要应用的一个或更多个规则，并且其中，所述一个或更多个规则应用于一个或更多个模 糊输入参数，以将所述一个或更多个模糊输入参数转换成一个或更多个模糊结论。14. 根据权利要求11至13中任意一项所述的方法，其中，步骤(iii)包括如下步骤:采用 从属函数和/或重心法，以将所述一个或更多个模糊结论转换成步骤(d)输出参数。15. 根据任意前述权利要求所述的方法，其中，如果所述计算得的钻进参数中的一个或 更多个不同于所述测量得的钻进参数中等效的一个或更多个，则进行将所述一个或更多个 计算得的钻进参数应用于所述钻进模块的步骤(e)。16. 根据任意前述权利要求所述的方法，其中，使用所述钻进模块中所并入的一个或更 多个传感器，来进行测量一个或更多个当前钻进参数以获得一个或更多个测量得的钻进参 数的步骤(b)。17. 根据权利要求16所述的方法，其中，从测压元件、涡流传感器、加速计、用于测量电 流或电流变化的传感器、以及用于测量电压或电压变化的传感器，选择所述一个或更多个 传感器。18. 根据任意前述权利要求所述的方法，其中，所述方法包括以下步骤:在所述控制方 法的一个或更多个步骤中使用进一步参数。19. 根据权利要求18所述的方法，其中，所述进一步参数从以下各项选择： (a) 所述钻头的直径； (b) 所述钻进模块的质量； (c) 所述被钻材料的温度； (d) 所述被钻材料的地静压力或孔隙压力； (e) 所述钻进模块的长度； (f) 钻孔的宽度；以及 (g) 所述钻孔的深度。20. -种钻进的方法，该方法包括以下步骤： 操作共振增强钻机，以对材料进行钻进，其中，根据如权利要求1至19中任意一项所述 的方法来控制所述共振增强钻机。21. -种控制系统，该控制系统用于共振增强钻进，该控制系统包括： (a) 用于确定被钻材料的一个或更多个特性的控制器，以及 (b) 用于确定一个或更多个钻进参数的控制器； 其中，所述控制系统能够执行如权利要求1至19所述的方法。22. -种设备，该设备用于共振增强钻进，该设备包括： (a) 钻进模块；以及 (b) 如权利要求21所述的控制系统。
【文档编号】G06K9/66GK105829645SQ201480067604
【公开日】2016年8月3日
【申请日】2014年10月9日
【发明人】M·维尔茨洛奇, Y·刘
【申请人】Iti 苏格兰有限公司