用于井数据分析的系统和方法转让专利

申请号 : CN200880127730.0

文献号 : CN101960087A

文献日 :

基本信息:

PDF:

法律信息:

相似专利:

发明人 : 尤利安·J·波普琼-马克·福利尼

申请人 : 普拉德研究及开发股份有限公司

摘要 :

本发明公开了用于分析在地层测试期间遇到的井数据的技术的示例。测试的一些部分可以显示在测试期间可能发生的异常行为、故障、误差或事件的指示。在执行测试期间或在执行测试之后可以识别一个或多个置信度标记。可以对一个或多个置信度标记进行分析以确定在测试期间这种异常行为、故障、误差或事件是否已经发生。然后,这些置信度标记可以用于确定由所执行的测试得到的结果和/或所述测试结果潜在的数据和解释的置信度水平。

权利要求 :

1.一种用于确定由随钻测试工具获得的测量值的置信度的方法,所述随钻测试工具位于穿透地下地层的井眼中,所述方法包括以下步骤:建立所述随钻测试工具的压力传感器与所述地层之间的压力耦合;

利用所述随钻测试工具执行第一压力下降;

利用所述压力传感器测量指示压力的数据;

根据压力数据确定至少一个置信度标记;以及显示所述至少一个置信度标记。

2.根据权利要求1所述的方法,还包括第二压力下降。

3.根据权利要求2所述的方法,其中,所述第二压力下降的参数至少部分地基于所述至少一个置信度标记。

4.根据权利要求1-3中任一项所述的方法,其中,使用参数比较技术确定所述置信度标记。

5.据权利要求1-4中任一项所述的方法,还包括以下步骤:终止所述压力耦合。

6.根据权利要求1-5中任一项所述的方法,还包括以下步骤:识别井下条件。

7.根据权利要求6所述的方法,还包括以下步骤:显示识别的所述条件。

8.根据权利要求1-7中任一项所述的方法,其中,使用趋势分析技术确定所述至少一个置信度标记。

9.根据权利要求8所述的方法,其中,所述趋势分析技术包括在压力恢复期端点之外映射压力数据。

10.根据权利要求1-9中任一项所述的方法,使用噪点离散分析技术确定所述至少一个置信度标记。

11.根据权利要求8-10中任一项所述的方法,所述趋势分析技术是整体分析技术。

12.根据权利要求1-11中任一项所述的方法,还包括以下步骤:使用压力比较技术、参数预测技术、模型相关技术、量规比较技术、和增压技术中的一种确定第二置信度标记。

13.根据权利要求1-12中任一项所述的方法,其中,所述至少一个置信度标记表征压力恢复期结束时的压力曲线的斜率。

14.一种用于确定由位于井眼中的工具获得的测量值的置信度的方法,所述井眼穿透地下地层,所述方法包括以下步骤:选择多个井下条件;

使不同的数值与所述井下条件中的每一个相关联;

执行井下测量;

根据所述井下测量识别所述井下条件中的一个;

将与识别的所述井下条件相关联的数值发送到地面显示器;

在地面显示器处接收所述数值;以及显示指示被识别的所述井下条件的记号。

15.根据权利要求14所述的方法,其中,根据所述井下测量识别所述井下条件中的一个的步骤包括:确定多个置信度标记,并且将所述置信度标记中的每一个与预定阈值进行比较。

16.根据权利要求14或15所述的方法,还包括以下步骤:确定至少一个置信度标记,并且显示所述置信度标记。

17.根据权利要求16所述的方法,还包括以下步骤:识别与井下工具的操作相关联的多个事件;

选择用于由所述井下工具发送的数据点,且所述数据点被选择作为所述事件的函数和增长函数;

确定与所述事件和用于由所述井下工具发送的数据点相关联的值;

将确定的所述值发送到地面显示器;以及在测井图上显示被发送的所述数据。

说明书 :

用于井数据分析的系统和方法

背景技术

[0001] 地层测试器通常包括适于下入到井眼内并在井眼内被定位在邻近于期望数据的地下地层的深度处的细长工具。一旦定位在井眼内,这些工具被放置成与地层进行流体连通以从地层采集数据。通常,探头、连通管或其它装置可与井壁密封接合以建立这种流体连通。
[0002] 地层测试器通常用于尤其测量诸如井筒压力、地层压力和地层流动性的井下参数。地层测试器还可以用于从地层采集样品,使得可以确定地层中所含有的流体类型及其它流体特性。在地层测试期间确定的地层特性是确定井的商业价值和可以从井开采油气的方式的重要因素。此外,通过随钻测量(MWD)确定的地层特性可能在控制进一步的钻井操作中是非常有用的。
[0003] 可以参照图1A和图1B中所示的传统的电缆式地层测试器的结构更加容易地理解地层测试器的操作。如图1A中所示,电缆式地层测试器100从钻机2被下入到填充有在本行业中通常被称作为“泥浆”的流体的井眼3内。井眼被覆盖有在钻井操作期间沉积到井眼的壁上的泥饼4。井眼3穿过地层5。美国专利No.4,860,581和No.4,936,139中更详细地说明了具有多个相互连接的模块的传统的模块化电缆式地层测试器的操作。图2示出了在用于确定诸如地层压力的参数的传统的电缆式地层测试操作期间由地层测试器测量的压力迹线随时间的图示。
[0004] 如图1A和图1B中所示,地层测试器100通过钢丝电缆6被下入到井眼3内到达井眼内的期望位置。然后,可以通过打开平衡阀(未示出)使地层测试器中的流动管线119中的压力与井眼中的流体的流体静压相等。压力传感器或压力计120用于测量井眼中的流体的流体静压,所述流体静压在图2中沿线103被示意性地示出。然后,可以利用液压致动活塞将地层测试器100固定在适当的位置,从而抵靠井眼的井壁定位探头112以建立与地层的流体连通,并且闭合平衡阀以使工具的内部与井液隔离。图2中在105处示意性地示出了在探头与地层直接进行密封并建立流体连通的、被称作为“工具安置”点的点。然后,通过将活塞118缩回预测试室114内以在流动管线119中生成小于地层压力的压降而将来自地层5的流体吸入到地层测试器100内。在图2中沿线107示意性地示出了被称作为“压力下降”期的这种体积膨胀期。
[0005] 当活塞118停止缩回(在图2中的点111处被示出)时,来自地层的流体继续进入探头112,直到在给定充足的时间下流动管线119中的压力与地层5中的压力相同,如图2中的115处所示。被称作“压力恢复”期的这种周期在图2中沿线113被示出。如图2中所示,通常被称作为“井底压力”的在115处表示的最终恢复压力通常被假定为与地层压力的良好近似。
[0006] 由压力迹线生成的曲线的形状和相对应的数据可以用于确定不同的地层特性。例如,在压力下降(图2中的107)和压力恢复(图2的中的113)期间测量的压力可以用于确定地层流动性,所述地层流动性是地层渗透率与地层流体粘度的比值。当地层测试器探头(图1B中的112)与井壁分离时,流动管线119中的压力由于流动管线中的压力与井筒压力平衡而迅速增加,如图2中的线117所示。在地层测量循环已经完成之后,可以卸除地层测试器100并将所述地层测试器重新定位在不同的深度处,并且随意地重复地层测试循环。
[0007] 在用于电缆输送工具的这类测试操作期间,在井下采集的压力数据通常通过有线通信系统以电子的方式被通信给地面。在地面上,操作员监视流动管线119中的压力并且电缆测井系统实时记录压力数据。可以在井位计算机处实时或随后在数据处理中心分析在压力下降和压力恢复期期间记录的数据以确定诸如地层流体压力、泥浆过平衡压力(井筒压力与地层压力之间的差)、和地层的流动性的关键的地层参数。
[0008] 电缆式地层测试器允许高数据率通信用于通过使用有线遥测术实时监测和控制测试和工具。这类通信系统能够使现场工程师当发生测试测量时对测试测量的质量进行评价,并且如有必要采取立即行动以在试图进行另一个测量之前放弃测试过程和/或调节预测试参数。例如,通过观察预测试压力下降期间的数据,工程师可以进行选择以改变诸如压力下降速度和压力下降体积的初始预测试参数,以在试图进行另一个测试之前使所述初始预测试参数与地层特征进行更好地匹配。例如,在美国专利No.3,934,468;No.4,860,581;No.4,936,139;和No.5,969,241中说明了现有技术电缆式地层测试器和/或地层测试方法的示例。
[0009] 地层测试器还可以在钻井操作期间使用。美国专利No.6,230,557;No.5,803,186;No.7,114,562;和No.5,233,866中说明了适于在钻井操作期间从地下地层采集数据的示例性井下钻井工具。
[0010] 已经发展了用于执行专业的地层测试操作、或预测试的各种技术。例如,美国专利No.5,095,745和No.5,233,866说明了通过分析压力与线性压力下降偏离的点来确定地层参数。美国专利No.6,932,167;No.7,011,155;No.7,234,521;和No.7,178,392中提供了其它的示例。
[0011] 尽管对用于执行预测试的方法做了改进,但是仍然需要消除预测试过程中的延迟和误差,并提高由这种测试获得的参数的精度。因为地层测试操作在整个钻井操作中都使用,因此测试的持续时间和与工具的非实时通信是必须要考虑的重要约束。与用于这些操作的实时通信相关联的问题很大程度上是由于通常在钻井操作期间使用的遥测技术(例如,泥浆脉冲遥测术)的当前限制。用于大多数随钻测井(LWD)或随钻测量(MWD)工具的诸如上行链路和下行链路遥测数据速度的限制使得在井下工具与地面之间产生较慢的信息交换。例如,在工程师将指令发送到井下以根据发送的数据收回探头之后将预测试压力迹线发送到地面的简单过程可能产生往往对钻井操作产生不利影响的相当大的延迟。
[0012] 延迟还增加工具在井眼内卡钻的可能性。为了减少卡钻的可能性,通常建立基于主要岩层和钻井条件的钻井操作规范以指示在给定井眼中可以固定多长的钻柱。在此规范下,仅可以允许钻柱固定持续有限的时期以布置探头并执行压力测量。因此,由于与遥测技术带宽相关联的限制,实时发送在测试期间获得的所有数据可能是不可行的,因此适当的数据分析和/或控制可能是不可以的。
[0013] 其中实施两个阶段测试协定的地层压力随钻(FPWD)测量说明了需要实时地层测试数据通信。例如,FPWD预测试可以包括假定包括压力下降期和压力恢复期的、作为调查研究阶段实施的第一阶段,和假定也包括压力下降期和压力恢复期的、作为测量阶段被实施的第二阶段。来自调查阶段的数据可以用于构造/执行测量阶段。如果来自调查阶段的数据不是在井口处提供的,则可能不能相对于构造测量阶段、继续测试等进行适当分析和/或控制。类似地,如果来自测量阶段的数据不是在井口处提供,则可能不能相对于继续钻井操作、进一步测试等进行适当分析和/或控制。在16比特/样品的情况下具有15Hz采样率的五分钟时限预测试例如产生72,000比特每数据通道。然而,在实施泥浆脉冲遥测中,通信信道容量通常局限于0.5到12比特/秒之间。这种通信信道通常不足以实时传输FPWD预测试数据。
[0014] 在发展用于地层测试的方法中已经有了进步,但是仍然需要提高对在井下测试期间生成的数据的评价和/或通过测试数据质量控制改进测试程序。例如,需要对在测试过程中产生的影响测试结果的误差进行评价。此外,恶劣的井下条件可能影响设备的性能、井下参数的测量和/或可能影响所提供的整个数据的各种其它因素。由于错误的测试结果可能会进行不正确的判定。因此,期望提供用于检测数据中的潜在问题或误差的技术。此外,期望这种系统提供用于分析井下测量以确定结果的精度和/或结果的置信度的测量的(自动或手动)技术。

发明内容

[0015] 公开了用于分析在地层测试期间遇到的压力迹线的技术的示例。测试的一些部分可以显示在测试期间可能发生的异常行为、故障、误差或事件的指示。在执行测试期间或在执行测试之后可以识别一个或多个置信度标记(confidence token)。可以对一个或多个置信度标记进行分析以确定在测试期间这种异常行为、故障、误差或事件是否已经发生。然后,这些置信度标记可以用于确定由所执行的测试和/或所述测试的潜在的数据和解释得到的结果的置信度水平。
[0016] 因此,本公开提供了一种用于确定由随钻测试工具获得的测量值的置信度的方法。所述方法包括以下步骤:建立测试工具的压力传感器与地层之间的压力耦合;利用测试工具执行第一压力下降;利用压力传感器测量指示压力的数据;根据压力数据确定至少一个置信度标记;以及显示至少一个置信度标记。
[0017] 本公开还提供了一种用于确定由测试工具获得的测量值的置信度的方法。所述方法包括以下步骤:建立测试工具的压力传感器与地层之间的压力耦合;利用测试工具执行第一压力下降;利用压力传感器测量指示压力的数据;使用趋势分析或噪点离散分析技术根据压力数据确定至少一个置信度标记;以及显示至少一个置信度标记。
[0018] 本公开还提供了一种用于确定由井下工具获得的测量值的置信度的方法。所述方法包括以下步骤:选择多个井下条件;使不同的数值与井下条件中的每一个相关联;执行井下测量;根据井下测量识别井下条件中的一个;将与识别的井下条件相关联的整数发送到地面显示器;在地面显示器处接收所述整数;以及显示指示被识别的井下条件的记号。
[0019] 以上已经概述了本公开的相当广泛的一些特征和技术优点,以便可以更好地理解随后的详细说明。以下描述形成权利要求的主体的另外的特征和优点。本领域的技术人员应该认识的是可以容易地使用所公开的原理和具体实施例作为修改或设计用于实施相同目的的其它结构和/或方法。本领域的技术人员还应该认识的是这些等效结构和/或方法没有背离本公开的如所附权利要求所述的保护范围。但结合附图时,将从以下说明更好地理解本公开。然而,特别要理解的是提供每一幅图用于仅用于说明和图示,并且不旨在作为本公开的限制的限定。

附图说明

[0020] 参照附图进行以下说明,其中:图1A、1B和2示出了现有技术的多个方面,而其余的附图示出了本公开的多个方面。
[0021] 图1A显示设置在井眼中的电缆式地层测试器;
[0022] 图1B是图1A的测试器的横截面图;
[0023] 图2对于使用地层测试器执行典型的预测试程序显示了压力对时间图;
[0024] 图3是预测试方法的流程图;
[0025] 图4是地层测试器模块的示意图;
[0026] 图5对于图3的预测试显示压力对时间图;
[0027] 图6是图3的方法的调查阶段的流程图;
[0028] 图7是图5的图表的调查阶段部分的、显示压力下降终止的详细视图;
[0029] 图8是图5的图表的调查阶段部分的、显示压力恢复终止的确定的详细视图;
[0030] 图9是图3的方法的测量阶段的流程图;
[0031] 图10是并入泥浆压缩系数阶段的预测试的流程图;
[0032] 图11A对于图10的预测试显示压力对时间图;
[0033] 图11B显示与图11A相对应的体积变化率;
[0034] 图12是图10的方法的泥浆压缩系数阶段的流程图;
[0035] 图13是并入泥浆滤失预测试的流程图;
[0036] 图14A对于图13的预测试显示压力对时间图;
[0037] 图14B显示与图14A相对应的体积的变化率;
[0038] 图15显示被修改用于与泥浆滤失阶段一起使用的图12的修改的泥浆压缩系数阶段;
[0039] 图16A-16C是图13的方法的泥浆滤失阶段的程序框图,其中,图16A显示泥浆滤失阶段,图16B显示具有重复压缩循环的修改的泥浆滤失阶段,而图16C显示具有减压期的修改的泥浆滤失阶段;
[0040] 图17A对于执行包括修改的调查阶段的预测试显示压力对时间图;
[0041] 图17B显示与图17A相对应的体积的变化率;
[0042] 图18是图17A的被修改的调查阶段的流程图;
[0043] 图19A对于执行包括修改的调查阶段的预测试显示压力对时间图;
[0044] 图19B显示与图19A相对应的体积的变化率;
[0045] 图20是图19A的被修改的调查阶段的流程图;
[0046] 图21是可以用于当在不同的温度和/或压力下执行原始泥浆压缩系数时提供校正的泥浆压缩系数的流体压缩系数校正图;
[0047] 图22显示由地层测试器生成的压力对时间图;
[0048] 图23是用于提供数据压缩和通信的流程图;
[0049] 图24是用于构造图23的抽取/压缩数据步骤的流程图;
[0050] 图25和图26是用于数据抽取以便进行数据压缩的程序框图;
[0051] 图27是大致与图22的预测试相对应的、与用于通信的数据集相关联的曲线;
[0052] 图28是大致与图22的预测试测量阶段相对应的、与用于通信的数据集相关联的曲线;
[0053] 图29是用于量化技术的流程图;
[0054] 图30是提供非均匀量化的数据压缩扩展器的操作的图示;
[0055] 图31显示由地层测试器生成的、显示沿压力恢复的数据点的压力对时间图;
[0056] 图32显示由地层测试器生成的、显示沿压力恢复的压力区间的压力对时间图;
[0057] 图33A显示了用于确定曲线的选定点处的平滑值的过滤器的示例;
[0058] 图33B显示了用于确定选定点处的曲线斜率的平滑值的过滤器的示例;
[0059] 图34是确定预测试的置信度的方法的流程图;
[0060] 图35是使用压力比较技术确定置信度的方法的流程图;
[0061] 图36A显示说明了丧失密封的压力对时间图;
[0062] 图36B显示说明了阻塞流的压力对时间图;
[0063] 图37是使用参数比较技术确定置信度的方法的流程图;
[0064] 图38是使用参数预测技术确定置信度的方法的流程图;
[0065] 图39是使用曲线分析技术确定置信度的方法的流程图;
[0066] 图40A对于非扩展区间显示压力对时间图;
[0067] 图40B对于扩展区间显示压力对时间图;
[0068] 图40C显示说明了虚拟数据集的压力对时间图;
[0069] 图41是使用数据方差技术确定置信度的方法的流程图;
[0070] 图42是使用模型相关技术确定置信度的方法的流程图;
[0071] 图43显示说明了参数函数与数据的后期拟合的压力对时间图;
[0072] 图44A显示说明了在压力恢复期间的渗漏的压力对时间图;
[0073] 图44B显示在图44A的压力曲线下方的面积与时间图;
[0074] 图45是使用量规比较技术确定置信度的方法的流程图;
[0075] 图46是使用增压技术确定置信度的方法的流程图;
[0076] 图47是分析如图34中所述的置信度标记的方法的流程图;以及
[0077] 图48是显示例如利用图47的方法被识别的井下条件的方法的流程图。

具体实施方式

[0078] 将要理解的是以下公开提供许多不同的实施例、或示例,用于执行各种实施例的不同特征。部件和装置的具体示例被描述如下以简化本公开。然而,这些仅仅是示例性的并且不旨在进行限制。此外,本公开可以在各种示例中重复附图标记和/或字母。这种重复出于简单和清楚的目的并且本身不表示所述的各种实施例和/或结构之间的关系。此外,在以下说明中第一特征经过或在第二特征上的形成可以包括其中第一和第二特征被形成为直接接触的实施例,并且还可以包括其中另外的特征可以形成在第一与第二特征之间使得第一和第二特征没有直接接触的实施例。
[0079] 图3中示出了根据本公开的一个或多个方面的用于评价地层特性(例如,地层压力和流动性)的方法。所述方法包括调查阶段13和测量阶段14。可以利用本领域公知的任意地层测试器(例如,相对于图1A和1B所述的测试器,以及美国专利No.4,860,581;No.4,936,139;No.6,230,557;和/或No.7,114,562所述的设备)实施所述方法。
[0080] 图4中示出了可与这种地层测试器一起使用的探头模块101。模块101包括探头112a、包围探头的封隔器110a、和从探头延伸到模块内的流动管线119a。流动管线119a从探头112a延伸到探头隔离阀121a,并且具有压力计123a。第二流动管线103a从探头隔离阀121a延伸到取样管线隔离阀124a和平衡阀128a,并且具有压力计120a。预测试室114a中的可逆预测试活塞118a也从流动管线103a延伸。引出管线126a从平衡阀128a延伸并延伸出来到达井眼,并具有压力计130a。样品流动管线125a从取样管线隔离阀124a延伸并通过工具。在流动管线125a中取样的流体可以被捕获、冲洗或用于其它目的。
[0081] 探头隔离阀121a隔离流动管线119a中的流体与流动管线103a中的流体。取样管线隔离阀124a隔离流动管线103a中的流体与取样管线125a中的流体。平衡阀128a隔离井眼中的流体与工具中的流体。通过操纵阀以选择性地隔离流动管线中的流体,压力计120a和123a可以用于确定各种压力。例如,当探头与地层流体连通同时最小化连接到地层的工具体积时,通过闭合阀121a,可以通过压力计123a获取地层压力。
[0082] 在另一个示例中,在平衡阀128a打开的情况下,可以借助于预测试活塞118a将泥浆从井眼吸入到工具内。在闭合平衡阀128a、探头隔离阀121a和取样管线隔离阀124a时,可以在工具内将流体捕集在这些阀与预测试活塞118a之间。压力计130a可以用于在工具的整个操作期间连续监测井筒流体压力,并且与压力计120a和/或123a一起可以用于直接测量泥饼两端的压降并监测泥饼两端的井眼干扰的传输,用于随后在对测量的井底压力进行校正时使用。
[0083] 预测试活塞118a的一个功能是从地层抽取流体或将流体注入到地层内或者压缩或膨胀捕集在探头隔离阀121a、取样管线隔离阀124a和平衡阀128a之间的流体。预测试3 3
活塞118a优选地具有在低流量(例如,0.01cm/秒)和在高流量(例如,10cm/秒)下操
3
作的能力,并且具有能够在单个冲程中抽取较大体积(例如,100cm/秒)的能力。此外,
3
如果需要在不收回探头的情况下从地层采出超过100cm 的流体,可以重复利用预测试活塞
118a。优选地,可以连续监测并强制控制预测试活塞118a的位置,并且当预测试活塞118a静止时,可以“锁定”所述预测试活塞的位置。在一些实施例中,探头112a还可以包括过滤阀(未示出)和过滤活塞(未示出)。
[0084] 阀、预测试活塞和探头的各种操纵允许根据所述方法对工具进行操作。本领域的技术人员将认识到虽然这些说明限定了优选的探头模块,但是在不背离本公开的保护范围的情况下可以使用其它说明。虽然图4示出了探头型模块,但是将要认识的是在假定在一些修改的情况下可以使用探头工具或封隔器工具。以下说明假定使用探头工具。然而,本领域的技术人员将认识到类似的程序可以与封隔器工具一起使用。
[0085] 这里公开的技术也可与装有流动管线的其它装置一起使用。如这里所使用的术语“流动管线”应该表示用于在地层与预测试活塞之间建立流体连通和/或用于允许流体在所述地层与所述预测试活塞之间流动的导管、空腔、或其它通道。其它这种装置可以包括例如其中探头和预测试活塞是一体的装置。美国专利No.6,230,557和No.6,986,282中公开了这种装置的示例。
[0086] 如图5中所示,调查阶段13表示获得诸如地层压力和地层流动性的地层参数的初始估计值。然后,这些初始估计值可以用于设计测量阶段14。如果期望并允许,则可根据这些参数执行测量阶段以生成改进的地层参数的估计值。图5示出了当执行图3的方法时显示压力随时间的变化的相对应的压力迹线。将要认识的是虽然可以通过图4的设备执行图5的压力迹线,但是也可以通过诸如图1A和1B的测试器的其它井下工具执行所述压力迹线。
[0087] 图6中更加详细地显示了调查阶段13。调查阶段包括在工具安置之后开始压力下降并持续持续时间Ti到t3的步骤310,执行压力下降的步骤320、终止压力下降的步骤330、执行压力恢复的步骤340以及终止压力恢复的步骤350。为了开始根据步骤310的调查阶段,探头112a被放置成与地层流体连通,并被固定到适当的位置,并且工具的内部与井眼隔离。通过使活塞118a在预测试室114a中移动而执行压力下降的步骤320。为了结束压力下降的步骤330,停止活塞118a。在步骤340处,压力将开始在流动管线119a中恢复直到在步骤350处压力恢复结束。调查阶段持续时间TIP。也可以如先前相对于图1B和图2所述的执行调查阶段,且在调查阶段开始之前预先限定压力下降流量和压力下降终点。
[0088] 图7中更加详细地示出了调查阶段13的压力迹线。可以从对由调查阶段的压力迹线获得的数据进行分析来确定诸如地层压力和地层流动性的参数。例如,终点350表示地层压力的临时估计值。可选地,可以通过使用本领域的技术人员所公知的技术对在压力恢复340期间获得的压力变化趋势进行外推来更精确地估算地层压力,且允许与已经获得的压力相对应的外推压力使得可允许压力恢复无限持续。这种过程可能需要另外的处理以达到地层压力。
[0089] 还可以从由线340表示的压力恢复阶段确定地层流动性(K/μ)1。本领域的技术人员所公知的技术可以用于由在压力恢复340期间压力随时间的变化速度来估计地层流动性。这种过程可能需要另外的过程以达到地层流动性的估计值。
[0090] 可选地,在Goode等人的题目为“Multiple Probe Formation Testing and Vertical Reservoir Continuity”(SPE 22738,准备在1991年10月6日到9日在美国德克萨斯州达拉斯举办的1991年石油工程师协会年技术会议和展览会上介绍)出版物中所述的工作暗示了由阴影区所示并由附图标记325表示的面积(这里由A表示)可以用于预测地层流动性。此面积由从终点350(表示在终点处估算的地层压力P350)水平延伸的线321、压力下降线320和压力恢复线340界限。此面积可以通过使用以下公式确定并与地层流动性的估计值有关:
[0091]
[0092] 其中(K/μ)1是地层流动性的第一估计值(D/cP),其中K是地层渗透率(达尔西,由D表示),而μ是地层流体粘度(cP)(因为由地层测试器确定的数是地层渗透率与地层3
流体粘度的比值,即,流动性,因此不需要粘度的确切值);V1(cm)是在调查预测试期间从地层开采的体积,V1=V(t7+T1)-V(t7-T0)=V(t7)-V(t7-T0),其中V是预测试室的体积;rp是探头半径(cm);而εk是通常对于具有流动性大于1mD/cP的地层来说为非常小(小于百分之几)的误差项。
[0093] 解释有限尺寸井眼对探头的压力响应的影响的变量ΩS可以由在F.J.Kuchuk的 题 目 为“Multiprobe Wireline Formation Tester Pressure Behavior in Crossflow-Layered Reservoirs”(In Situ,(1996)20,1,1)的出版物中所述的以下公式确定:
[0094]
[0095] 其中rp和rw分别表示探头的半径和井的半径;ρ=rp/rw,η=Kr/Kz;而Kr和Kz分别表示径向渗透率和垂向渗透率。
[0096] 在说明公式1中所述的结果中,已经做了以下假设:地层渗透率是各向同性的,即,Kr=Kz=K,在测试期间的流态是“球形的”,和保持确保达尔西关系的有效性的条件。
[0097] 在图7中,可以对调查阶段的压力下降步骤320进行分析以确定随时间的压降,从而确定压力迹线的各种特征。由沿压力下降线320的点获得的最佳拟合线32被示出为从起始点310延伸。可以沿曲线320确定偏差点34,所述偏差点表示曲线320到达距离最佳拟合线32的最小偏差δ0的点。偏差点34可以用作“流动的开始”的估计值,即,在调查阶段压力下降期间流体从地层被输送到工具内的点。
[0098] 可以由诸如美国专利No.5,095,745和No.5,233,866中所公开的技术的公知技术确定偏差点34,所述专利公开了一种用于由与使用来自预测试的压力下降阶段的数据点产生的最佳拟合线从偏差的点估计地层压力的技术。可选地,可以通过对最近获得的点进行测试以观察是否当获得连续压力数据时所述最近获得的点保持在表示流动管线膨胀的线性趋势上来确定偏差点。如果不是,可以终止压力下降,并允许稳定压力。还可以通过对在320期间记录的压力相对于时间求导数来确定偏差点。当导数改变(假定变小)2%-5%时,取相对应的点以表示来自地层的流动的开始。如有必要,为了确认与膨胀线的偏差表示来自地层的流动,则可以执行进一步更小体积的预测试。
[0099] 其它技术可以用于确定偏差点34,例如,用于确定偏差点34的另一种技术基于泥浆压缩系数并且以下相对于图9-11说明所述技术。
[0100] 一旦确定了偏差点34,压力下降持续超过点34直到满足一些规定的终止判据。这种判据可以基于压力、体积和/或时间。一旦已经满足判据,终止压力下降,并且到达终点330。理想的是终点330在相对于与图7的偏差点相对应的偏差压力P34的给定压力范围ΔP内发生在给定压力P330处。可选地,可以期望的是在确定偏差点34之后在给定时间周期内终止压力下降。例如,如果偏差发生在时间t4处,终止可以被预先设置成为在时间t7前发生,且时间t4与时间t7之间过去的时间被表示为TD并且被限制到最大持续时间。用于终止预测试的另一种判据是在已经识别偏差点34之后限制从地层抽取的体积。这种体积可以由预测试室114a的体积的变化来确定(图4)。最大体积变化可以被指定为用于预测试的限制参数。
[0101] 限制判据、压力、时间和/或体积中的一个或多个可以单独或结合在一起使用以确定终点330。例如,在高渗透地层的情况下,不能满足诸如预定压降的期望的判据,预测试的持续时间还可能受限于一个或多个其它判据。
[0102] 在达到偏差点34之后,压力沿线320持续下降直到膨胀在点330处终止。然后,闭合探头隔离阀121a和/或停止预测试活塞118a,并且调查阶段压力恢复340开始。流动管线中的压力的恢复持续直到在点350处发生压力恢复的终止。
[0103] 压力恢复变得充分稳定的压力通常作为地层压力的估计值。监测恢复压力以提供用于由恢复压力的递增稳定估计地层压力的数据。具体地,获得的信息可以在设计测量阶段瞬变时使用,使得在压力恢复结束时可获得地层压力的直接测量。保留应该允许调查阶段压力恢复持续多久以获得地层压力的初始估计值的问题。
[0104] 从先前的论述清楚地了解到在压力已经恢复到可识别与流动管线减压的偏差的水平(即,由图7上的P34表示的压力)之前不会终止压力恢复。在一种方法中,一组时限可以用于压力恢复的持续时间T1。T1可以被设定为一些数字,例如为从地层流动的时间的2-3倍。可以考虑其它技术和判据。
[0105] 如图5和图7中所示,终点350示出了压力恢复的结束,调查阶段的结束和/或测量阶段的开始。一些判据可以用于确定终点350应该何时发生。用于确定终点350的一种可能的方法是允许测量的压力稳定。为了建立可以相对迅速地得到终点350处的地层压力的合理精确估计值的点,可以使用用于确定关于建立何时终止的判据的程序。
[0106] 如图8中所示,一种这样的程序包括建立在压力下降终止的点330处开始压力增量。例如,这种压力增量可以是压力计分辨率的较大的倍数、或压力计噪音的倍数。当获得压力恢复数据时,连续压力点将落入一个这样的区间。选择在每一个压力增量中的最高压力数据点并且在相对应的时间之间构造差以产生时间增量Δti(n)。压力恢复持续直到两个连续时间增量的比值大于或等于预定数,例如2。在满足这种判据时,在最后的区间内最后被记录的压力点是计算的终点350。这种分析可以以数学的方式表示如下:
[0107] 起始于t7,调查阶段的压力恢复开始,得到下标的顺序:
[0108] i(n)>i(n-1),n=2,3,.......,使得对于n≥2,i(1)=1,并且[0109]
[0110] 其中nP是具有等于或大于例如4通常为10或更大数值的数字,δP是压力测量仪的额定分辨率;而εP是压力表噪音的较小倍数(例如,2倍)-诸如在泥浆压缩系数实验期间在安置工具之前可以确定的量。
[0111] 本领域的技术人员将认识到在不背离本公开的保护范围的情况下可以根据期望的结果选择其它nP和εP值。如果除了基点之外在由右手侧的公式(3)定义的区间内不存在任何点,可以使用在区间外部的最接近点。
[0112] 定义Δti(n)≡ti(n)-ti(n-1),当满足以下条件时可以终止压力恢复:pi(n)≥p(t4)=P34(图7)并且
[0113]
[0114] 其中mP是大于或等于例如2的数字。然后,地层压力的第一估计值被定义为:
[0115] p(ti(max(n)))=p(t7+T1)=P350 (5)
[0116] 在近似项中,当在压力恢复期间的压力大于与偏差点34相对应的压力并且压力的增加速率下降至少2个因子时,终止根据当前判据的调查阶段预测试。地层压力的近似值作为在压力恢复期间测量的最高压力。
[0117] 公式(3)和(4)一起设定在调查阶段期间确定地层压力的精度:公式(3)定义误差的下限,而mP大致定义估计值与实际地层压力的接近程度。mP的值越大,地层压力的估计值将越接近实际值,并且调查阶段的持续时间越长。
[0118] 用于终止调查阶段压力恢复的另一种判据可以基于压力恢复曲线的平直度,例如,通过比较压力恢复点的范围的平均值与压力计噪音的较小倍数(例如,2或4)来确定所述平直度。将要认识的是这里所公开的判据中的任一个单独或组合在一起可以用于终止调查阶段压力恢复(图5中的340)、测量阶段压力恢复(图5中的380并且如下所述)、或一般地说任意压力恢复。
[0119] 如图7中所示,终点350示出了在完成压力恢复阶段340之后调查阶段13的结束。然而,具有其中需要或期望终止预测试的情况。例如,在诸如当探头被堵塞时、测试是干的或地层流动性非常低使得测试基本上是干的、泥浆压力与地层压力精确平衡、检测到假破坏(false breach)、测试到非常低的渗透率地层、检测到流动管线流体的压缩系数的变化或发生其它问题的过程中的问题可以在完成整个循环之前证实预测试终止的正确。
[0120] 一旦期望在调查阶段期间终止预测试,则可以停止预测试活塞或闭合探头隔离阀121(如果存在),使得流动管线119中的体积被减小到最小值。一旦已经检测到具有问题,可以终止调查阶段。如果期望,可以执行新的调查阶段。
[0121] 参照回到图5,在完成调查阶段13时,可以判定条件是否允许或实现测量阶段14的期望性能。这种判定可以手动执行。然而,优选地,进行自动判定并且基于所设定的判据。
[0122] 可以使用的一种判据是时间。可能需要确定是否具有充足的时间TMP来执行测量阶段。在图5中,具有充足的时间来执行调查阶段和测量阶段。换句话说,执行两个阶段的总时间Tt小于被分配给循环的时间。通常,当TIP小于总时间Tt的一半时,具有充足的时间来执行测量阶段。
[0123] 可以用于确定是否继续进行测量阶段的另一个判据是体积V。也可能需要或期望确定测量阶段的体积是否至少与在调查阶段期间从地层开采的体积一样大。如果没有满足一个或多个条件,可以不必执行测量阶段。其它判据还可以确定是否应该执行测量阶段。可选地,尽管无法满足任意判据,但是在其余的规定时间继续调查阶段直到结束,使得调查阶段由于默认变成调查阶段和测量阶段。
[0124] 将要认识的是虽然图5示出了单个调查阶段13之后为单个测量阶段14,但是根据本公开的一个或多个方面可以执行不同数量的调查阶段和测量阶段。在极端情况下,因为在调查阶段压力恢复期间的压力增加可能非常缓慢使得分配给测试的整个时间被这种调查阶段耗费掉,因此调查阶段估计值可以是可获得的唯一估计值。这通常是地层具有非常低的渗透率的情况。在其它情况下,例如,在其中地层压力的恢复相对迅速的适当渗透地层到高渗透地层的情况下,在不和规定时间约束相冲突的情况下可以执行多个预测试。
[0125] 仍然参照图5,一旦确定了要执行测量阶段14,则调查阶段13的参数用于设计测量阶段。由调查阶段获得参数(即,地层压力和流动性)用于规定测量阶段预测试的操作参数。具体地,期望的是使用调查阶段参数求出测量阶段预测试的体积及其持续时间,并因此求出相对应的流量。优选地,以此方式确定测量阶段操作参数以优化在测量阶段预测试期间使用的体积,从而在给定的范围内生成地层压力的估计值。更具体地,期望的是开采正好够的体积,优选地,开采比在调查阶段期间从地层开采的体积大的体积,使得在测量阶段结束时,压力恢复到实际地层压力pf的期望范围δ内。优选地,选择在测量阶段期间开采的体积,使得也可以满足时间约束。
[0126] 令H表示地层对如先前所述由探头工具产生的流量的单位阶跃的压力响应。在测量阶段结束时测量压力在实际地层压力范围δ内的条件可以被表示为:
[0127]
[0128] 其中 是分配给调查阶段和测量阶段的总时间减去流动管线膨胀所花费的时间,即,在图5中 (在执行测试之前被指定-秒);T0是在调查阶段期间的地层流动的近似持续时间(在采集期间确定的-秒);T1是在调查阶段期间压力恢复的持续时间(在采集期间确定的-秒);T2是在测量阶段期间压力下降的持续时间(在采集期间确定的-秒);T3是在测量阶段期间压力恢复的持续时间(在采集期间确定的-秒);q1和q2分别表示调查阶段和测量阶段的恒定流量(在采集之前指
3
定并且在采集期间被确定-cm/秒);δ是在测量阶段期间确定地层压力的精度(被指定-大气压),即,pf-p(Tt)≤δ,其中pf是实际地层压力;φ是地层孔隙度,Ct是地层总压缩系数(在采集之前由地层类型和孔隙度的知识通过标准相互关系指定-1/大气压);
其中n=t,0、1、2表示无因次时间,而 表示时间常数;以及,
r*是由 定义的有效探头半径,其中
κ是具有模数 的第一类完全椭圆积分。如果地层是各向同性的,则r*=2rp/(πΩS)。等效地,测量阶段可以通过指定第二预测试流量与第一预测试流量的比值和测量阶段预测试的持续时间T2并因此指定其体积来限制。
[0129] 为了完整地详细说明测量阶段,理想的是根据附加条件进一步限制测量阶段。因为在完成调查阶段之后已知测量阶段的持续时间,即, 因此这种条件可以基于指定测量阶段的压力下降部分的持续时间相对于完成整个测量阶段可获得的总时间的比值。例如,希望允许使压力测量阶段的压力恢复的两倍(或大于两倍)的时间用于压力下降,则T3=nTT2,或, 其中nT≥2。公式(6)则可以求解测量
值与调查阶段预测试流量率的比值,并因此可以求解测量阶段的体积V2=q2T2。
[0130] 完成测量阶段预测试参数的说明的另一个条件将限制在测量阶段压力下降期间的压降。在与公式(6)中使用的相同符号和相同控制假设的情况下,此条件可以被写成为:
[0131]
[0132] 其中Δpmax(单位为大气压)是在测量阶段可允许的压力下降压降的最大值。
[0133] 用于确定测量阶段预测试参数的公式(6)和(7)的应用最好地说明了具体简单但重要的情况。为了进行说明,如前所述,假设调查阶段预测试和测量阶段预测试在精确控制的流量下实施。此外,假设工具库(tool storage)对压力响应的影响可以忽略,在压力下降和压力恢复中的流型是球状的,地层渗透率是各向同性的,并且确保满足达尔西关系的有效性。
[0134] 在以上假设下公式(6)具有以下形式:
[0135]
[0136]
[0137] 其中erfc是补余误差函数。因为误差函数的自变量通常较小,因此使用通常的平方根近似值时几乎不损失精度。在重新排列一些项之后,公式(8)可以被示出为具有以下形式:
[0138]
[0139]
[0140] 其中λ≡T2+T3,一旦已经完成调查阶段预测试,则测量阶段的持续时间是已知的。
[0141] 一旦在左手侧的括号中的表达式被进一步近似以获得用于测量阶段预测试的期望体积的公式,则此关系的实用性如下非常清楚:
[0142]
[0143] 在由公式(6)对公式(8)做相同的假设的情况下,公式(7)可以被写成为:
[0144]
[0145]
[0146] 其中,在对补余误差函数和重排项应用平方根近似值之后,可以被表示如下:
[0147]
[0148]
[0149] 合并公式(9)和(12),得到:
[0150]
[0151]
[0152] 因为在最后两个括弧/括号表达式中的项每一个都非常接近于一,因此公式(13)可以被近似为:
[0153]
[0154] 其中公式(14)给出了用于确定测量阶段压力下降的持续时间的公式,并因此与以上用于测量阶段预测试体积的结果结合,给出测量阶段预测试流量的值的公式。为了由公式(14)获得T2的实际估计值,应该保持以下条件:
[0155]
[0156] 公式(15)表示近似于最终压力的目标邻域应该大于由调查阶段留下的残余瞬态值的条件。
[0157] 总之,由公式(10)和(14)得到的V2和T2的估计值可以在使用公式(8)和(11)的更加复杂的参数估计方案中用作起始值。虽然公式(8)和(11)已经用于说明计算测量阶段参数的过程中的步骤,但是将要认识的是诸如工具库、地层复杂性等其它影响可以容易地并入判断过程中。如果已知地层模型,则在参数估计过程中使用更加通用的地层模型公式(6)和(7)。
[0158] 用于确定测量阶段预测试的上述方法假设在可以估算最佳预测试体积和持续时间之前将指定一些参数。这些参数包括:地层压力测量的精度δ;可允许的最大压力下降(Δpmax);地层孔隙度φ-其通常可从裸眼测井图获得;和总压缩系数Ct-其可以由又取决于岩性和孔隙度的已知的相互关系获得。
[0159] 在确定测量阶段预测试参数的情况下,应该在分配给整个测试的时间范围内获得地层压力和地层流动性的改进的估计值。
[0160] 在点350处,调查阶段结束而测量阶段可以开始。由调查阶段确定的参数可用于计算确定执行测量阶段14的参数所需的流量、预测试持续时间和/或体积。以下使用由在调查阶段估算的原始地层参数确定的改进的一组参数执行测量阶段14。
[0161] 如图9中所示,测量阶段14包括执行第二压力下降的步骤360、终止第二压力下降的步骤370、执行第二压力恢复的步骤380和终止压力恢复的步骤390。如先前根据图6的调查阶段13所述执行这些步骤。优选地,根据调查阶段的结果预先确定测量阶段的参数(例如,流量、时间和/或体积)。
[0162] 重新参照图5,测量阶段14优选地在调查阶段的终点处350开始,并且持续由测量阶段规定的持续时间TMP直到在点390处终止。优选地,执行调查阶段和测量阶段的总时间落入规定数量的时间内。一旦完成测量阶段,可以估算地层压力,并且工具被收回用于另外的测试、井下操作或从井眼移除。
[0163] 以下参照图10,示出了用于评价地层特性的并入泥浆压缩系数阶段11的上述方法的可选实施例。此实施例包括泥浆压缩系数阶段11、调查阶段13和测量阶段14。泥浆压缩系数的估计可以用于改进调查阶段过程,从而由调查阶段13和测量阶段14产生更好的参数估计值。图11A示出了与图10的方法相对应的压力迹线,而图11B显示预测试室体积的变化率的相关图示。
[0164] 在此实施例中,图4的地层测试器可以用于执行图10的方法。根据此实施例,隔离阀121a和124a可以与平衡阀128a结合使用以捕集流动管线103a中的大量液体。此外,隔离阀121a可以用于减少工具库体积效应以有助于进行迅速压力恢复。另外,平衡阀128a允许容易地冲洗流动管线以排出诸如气体的不期望的流体并且有助于利用井内流体重新填充流动管线段119a和103a。
[0165] 泥浆压缩系数测量可以例如通过以下步骤执行:首先借助于预测试活塞118a通过平衡阀128a从井眼将大量泥浆吸入到工具内,通过闭合平衡阀128a和隔离阀121a和124a隔离流动管线中的大量泥浆;通过借助于预测试活塞118a调节预测试室114a的体积压缩和/或膨胀被捕集泥浆的体积,以及借助于压力计120a同时记录被捕集流体的压力和体积。
[0166] 例如,借助于图4中未示出的适当的线性电位计或通过其它沿用已久的技术通过测量预测试活塞的位移可以非常精确地测量预测试室的体积。可以精确控制预测试活塞的速度以给出对预测试活塞流量qp的期望控制的装置在图4中也没有示出。用于实现这些精确流量的技术是公知的,例如,通过将活塞连接到正确形式的丝杠的使用,可以容易地实现本方法所需的齿轮箱和计算机控制的马达的速度。
[0167] 图11A和图12更加详细地示出了泥浆压缩系数阶段11。在安置工具之前并因此在实施调查和测量阶段之前执行泥浆压缩系数阶段11。具体地,工具没有必要安置抵靠在井眼上也没有必要固定在井眼中来实施泥浆压缩系数测试,从而减小工具由于固定钻柱而被卡的风险。然而,优选的是在靠近测试点的点处对井内流体进行取样。
[0168] 图12中更加详细地示出了用于执行压缩系数阶段11的步骤。这些步骤也与沿图11A的压力迹线的点相对应。如图12中所示,泥浆压缩系数测试的步骤包括:开始泥浆压缩系数测试的步骤510;从井眼将泥浆吸入到工具内的步骤511;隔离流动管线中的泥浆的体积的步骤512;压缩泥浆体积的步骤520;以及终止压缩的步骤530。接下来,在步骤540处,开始泥浆体积的膨胀,在步骤550处,泥浆体积膨胀一段时间直到在步骤560处终止。在步骤561处,打开流动管线与井眼的连通,并且在步骤570处,使流动管线中的压力与井筒压力相等,直到在步骤575处终止。再循环的预测试活塞以下可以在步骤580处开始。在步骤581处,从流动管线将泥浆排出到井眼中,并且在步骤582处,重复利用预测试活塞。当期望执行调查阶段时,则在步骤610处安置工具,并且在步骤620处终止打开流动管线与井眼的连同。
[0169] 泥浆压缩系数与流动管线流体(通常为整个钻井泥浆)的压缩系数有关。泥浆压缩系数的知识可以用于更好地确定直线32的斜率(如先前相对于图7所述),直线32的斜率又提高了对表示来自地层的流动的偏差点34的确定。泥浆压缩系数的值的知识因此产生更加有效的调查阶段13并且提供进一步改进由调查阶段13获得的估计值并最终提高由测量阶段14获得估计值的另外的手段。
[0170] 可以通过分析图11A的压力迹线和相应生成的压力和体积数据确定泥浆压缩系数Cm。具体地,泥浆压缩系数可以从以下公式确定:
[0171] 或 等 价 地3
[0172] 其中Cm是泥浆压缩系数(1/psi),V是被捕集的泥浆的总体积(cm),p是测量的流动管线压力(psi),是测量的流动管线压力的时间变化率(psi/秒),而qp表示预测试3
活塞流量(cm/秒)。
[0173] 为了获得泥浆压缩系数的精确估计值,理想的是在泥浆压缩系数测量期间采集多个数据点以限定压力-体积趋势的每一个路程。在使用公式(16)确定泥浆压缩系数时,已经进行通常的假设,具体地,压缩系数恒定,并且在测量中使用的增加的预测试体积与捕集在流动管线中的泥浆的总体积V相比较小。
[0174] 以下说明在获得更加精确的偏差点34a时测量泥浆压缩系数的实用性。通过使调查阶段13的压力下降数据的初始部分的已知斜率的线32a与数据拟合开始所述方法。线32a的斜率由先前所确定的泥浆压缩系数、流动管线体积、和预测试活塞压力下降速度固定。因为在固定和精确的控制速度下操作压力下降,并且流动管线流体的压缩系数是已经由上述实验确定的已知常数,因此以下给出说明具有已知斜率a的这种线的公式:
[0175]
[0176]
[0177] 其中V(0)是在膨胀开始时流动管线体积,Cm是泥浆压缩系数,qp是活塞减压流量,+p 是在膨胀过程开始时的视压力。假设V(0)由于预测试室的膨胀而远远大于体积的增量。
[0178] 因为斜率a现在是已知的,需要被指定以完全限定公式(17)的唯一参数是截距+ +p,即,b。总之,p 是未知量,然而,当属于流动管线膨胀的线性趋势的数据点与具有斜率a的线拟合时,所述数据点应该都产生类似的截距。因此当流动管线膨胀的线性趋势被识别+
时,截距p 的值将形成。
[0179] 可识别落到具有限定斜率a的线上的数据点到给定精度内的延伸。此线表示实际泥浆膨胀压降趋势。本领域的技术人员将认识到在对数据点与线进行拟合时,没有必要使所有点精确地落在线上。相反,足够的是,在基于工具特征和操作参数所选择的精度限制范围内使数据点与线拟合。利用这种方法,可以避免与早期的数据点相关联的不规则趋势,即,在预测试活塞压力下降开始周围的那些点。最后,在限定与线显著偏离(或超过精度限制)的直线的点之后的第一点34a是其中发生压力下降趋势的偏离的点。偏差点34a通常发生在比通过对直线进行外推而预测到的更高的压力下。这点表明泥饼的破坏。
[0180] 不同过程可用于识别属于流动管线膨胀线的数据点。然而,任意过程的细节取决于希望如何确定流动管线膨胀线、如何选择最大区间,以及如何选择精度的测量等。
[0181] 如下给出两种可能的方法以说明细节。在进行此之前,定义以下项:
[0182]
[0183] 以及 (19)
[0184]
[0185] 其中,通常,N(k)<k表示选自获得的k个数据点(tk,pk)的数据点的数量。基于上下文,N(k)可以等于k。公式(18)和(19)分别表示具有固定斜率a的最小二乘线和通过N(k)数据点的具有固定斜率a的最小绝对偏差的线,而公式(20)表示关于固定斜率线的数据的方差。
[0186] 用于限定跨度最长时间区间的具有斜率a的线的一种技术是当获得单独的数据点时拟合所述单独的数据点与固定斜率a的线。这种拟合产生一组截距{bk},其中单个{bk}由以下公式计算:bk=pk+atk。如果bk的连续值变得逐渐接近并且最终落入窄带内,与这些符号相对应的数据点用于拟合最终的线。
[0187] 具体地,所述技术包括以下步骤:
[0188] (i)由给定的一组截距{bk}确定中位数
[0189] (ii)得到属于集合 的符号,其中nb是诸如2或3的数字,并且其中εb的可能选择由以下公式定义:
[0190]
[0191] 其中最后一个表达式由时间测量值是精确的假设产生。对εb的其它非自然选择可以例如为εb=Sp,k;
[0192] (iii)拟合固定斜率a的线与具有属于 的符号的数据点;以及
[0193] (iv)得到产生 的第一点(tk,pk),其中 或 取决于用于拟合线的方法,而nS是诸如2或3的数字。由图11A中的34a表示的此点往往表示泥饼的破坏和来自地层流动的开始。
[0194] 可选的方法基于当拟合线遇到实际流动管线膨胀数据时数据关于恒定斜率的线的方差序列(sequence of variance)应该最终变得或多或少恒定。因此,根据本公开的一个或多个方面的方法可以实施如下:
[0195] (i)固定斜率a的线首先与累加直到时间tk的数据拟合。对于每一组数据,由确定线,其中,通过公式(18)来计算;
[0196] (ii)在N(k)=k的情况下使用公式(20)构造方差 的序列;
[0197] (iii)发现连续符号属于以下集合:
[0198]
[0199] (iv)固定斜率a的线与具有在 中的符号的数据拟合。令N(k)是集合中的符号的数量。
[0200] (v)确定与一系列固定斜率的线中的具有在上述集合中的标记的最后一条线分离的点作为满足 的第一点,其中nS是诸如2或3的数字;
[0201] (vi)定义
[0202] (Vii)获得 的点的子集从而
[0203] (viii)通过具有 中的符号的点拟合具有斜率a的线;以及
[0204] (ix)定义泥饼的破坏作为第一点(tk,pk),其中 当在先前的选择中,由图11A上的附图标记34a再次表示的此点表示泥饼的破坏和来自地层的流动的开始。
[0205] 一旦确定了最佳拟合线32a和偏差点34a,可以相对于图7如上所述确定终点330a、压力恢复370a和压力恢复的终止350a。测量阶段14则可以由在图11A的调查阶段
13生成的改进的参数来确定。
[0206] 以下参照图13,示出了用于评价地层特性的并入泥浆滤失阶段12的方法的可选实施例。在此实施例中,所述方法包括泥浆压缩系数阶段11a、泥浆滤失阶段12、调查阶段13和测量阶段14。图14A中示出了相对应的压力迹线,而图14B中示出了相对应的预测试体积的变化率的图示。与相对于图10的方法的相同的工具也可以与图13的方法结合使用。
[0207] 图14A和图14B更加详细地示出了泥浆滤失阶段12。在安置工具之后并且在执行调查阶段13和测量阶段14之前执行泥浆滤失阶段12。在泥浆滤失阶段12之前执行修改的泥浆压缩系数阶段11a。
[0208] 图15中更加详细地示出了修改的压缩系数测试11a。修改的压缩系数测试11a包括与图12的压缩系数测试11同样的步骤510-580。在步骤580之后,重复泥浆压缩系数测试的步骤511和512,即,在步骤511a处,从井眼将泥浆吸入到工具内,并且在步骤512a处隔离流动管线与井眼。以下在步骤610处可以安置工具,并且在步骤620处,在设定的周期终止时可以隔离流动管线,以便为泥浆滤失阶段、调查阶段和测量阶段做准备。
[0209] 图16A中更加详细地示出了泥浆滤失阶段12。在步骤710处开始泥浆滤失阶段,在步骤711中压缩流动管线中的大量泥浆,直到在点720处终止,并且在步骤730处流动管线压力下降。在初始压缩之后,在步骤751处打开井眼内的流动管线的连通,在步骤752处,工具和井眼内的压力平衡,并且在步骤753处,流动管线与井眼隔离。
[0210] 任选地,如图16B中所示,可以执行修改的泥浆滤失阶段12b。在修改的泥浆滤失阶段12b中,在打开流动管线的连通之前在步骤751处执行第二压缩,所述第二压缩包括以下步骤:在步骤731处开始对流动管线中的泥浆重新压缩,在步骤740处将流动管线中大量泥浆压缩到更高压力,在步骤741处,终止再压缩。然后在步骤750处允许流动管线压力下降。可以如相对于图16A所述执行步骤751-753。图14A的压力迹线显示图16B的泥浆滤失阶段12b。
[0211] 在图16C中所示的另一种选择12c中,可以在第一压缩步骤711中的流动管线压力下降的步骤730之后执行减压循环,所述减压循环包括以下步骤:在步骤760处开始对流动管线中的泥浆进行减压,在步骤770处减压到在井筒压力以下的适当压力,以及在步骤780处终止减压。然后,在步骤750中允许流动管线压力下降。步骤751-753则可以如先前相对于图16A所述而被重复。图14A的压力迹线显示图16C的泥浆滤失阶段12c。
[0212] 如图14A中的压力迹线所示,图16A的泥浆滤失方法12可以在图16B的泥浆滤失阶段12b或图16C的泥浆滤失阶段12c的情况下执行。任选地,图16A-C中所示的一个或多个技术可以在泥浆滤失阶段期间执行。
[0213] 泥浆滤失涉及通过沉积在井壁上的泥饼的基液泥浆的滤失和在现有井眼条件下滤失的体积流量的确定。假设泥饼特性在测试期间保持不变,通过泥饼的滤失流量由一下简单的表达式给出:
[0214]
[0215] 其中Vt是被捕集的泥浆的总体积(cm3),而qf表示泥浆滤失流量(cm3/秒);Cm表示泥浆压缩系数(1/psi)(其中Cm在修改的泥浆压缩系数测试11a期间确定,或被输入);表示当在图14的步骤730和750期间测量时压力下降速度(psi/秒)。公式(22)中的体积Vt表示在如图4中所示的阀121a、124a和128a之间所含有的流动管线的体积。
[0216] 对于不足以密封井壁的泥饼来说,泥浆滤失的流量可以是在调查阶段的流动管线减压期间的预测试活塞流量的显著分数,并且如果不考虑,则可以使被检测为来自地层的流动开始的点的点(图7的34)产生误差。在流动管线减压阶段期间使用以检测来自地层的流动开始的点(即,图7的偏差点34)的固定斜率线的斜率a在在这种情况下使用以下公式确定:
[0217]
[0218]
[0219] 其中V(0))是在膨胀开始时的流动管线体积,Cm是泥浆压缩系数,qp是活塞减压+流量,qf是从流动管线通过泥饼滤失进入地层内的滤失流量,而p 是如先前所述在确定偏差点34的过程期间确定的膨胀过程开始时的视压力。
[0220] 一旦已经确定了泥饼滤失流量qf和泥浆压缩系数Cm,则可以继续在其中通过泥饼的泥浆滤失是显著的情况下由调查阶段13估计地层压力。
[0221] 在本公开的保护范围内的实施例可以以自动方式实施。此外,所述实施例可以适用于井下钻井工具,和由诸如钻柱、钢丝电缆、连接的油管、或挠性管的任意类型的工作柱被输送到井下的电缆式地层测试器。有利地,本公开的方法可以允许井下钻井工具以时间最有效的方式执行时间约束地层测试,使得可以最小化或避免与停止钻具相关联的潜在问题。
[0222] 以下参照图17A、图17B和图18说明执行调查阶段测量的另一个实施例。在安置地层测试器的步骤805之前,优选地如上所述确定泥浆压缩系数(未示出)。在确定泥浆压缩系数之后并且在安置地层测试器之前,在步骤801处,由工具测量的压力是井筒流体、泥浆静水压。在步骤805中安置工具之后,在步骤810中,启动如图4中所示的预测试活塞118a以精确和固定的流量抽取流体,从而在步骤814中,在期望的时间815中获得指定压降。优选地,如果超平衡几乎已知,则期望的压降(Δp)与所述深度处的预期的超平衡同量级但是小于所述预期的超平衡。超平衡是泥浆流体静压与地层压力之间的压力差。可选地,期望的压降(Δp)可以是大于“流动初始压力”的最大预期值(例如,200psi)的某个值(例如,300psi)。实际地层压力是否在此范围内对于本公开的实施例是不重要的。因此,以下说明假定地层压力不在所述范围内。
[0223] 根据本公开的实施例,用于实现此限定压降(Δp)的活塞压力下降流量可以由以下公式估计:
[0224]
[0225] 其中Cm是被假定为与井内流体相同的流动管线流体的压缩系数;Vt是在图4中所示的阀121a、124a和128a之间在流动管线103中的被捕集流体的体积;Δp是期望的压降,而tpi是预测试压力下降的持续时间。
[0226] 参照图17A、图17B和图18,根据本公开的实施例的执行调查阶段13b的方法包括开始压力下降的步骤810和执行控制压力下降的步骤814。优选地,活塞压力下降速度被精确控制,使得可以很好地控制压降和压力变化的速度。然而,不是必须要在低流量下实施预测试(活塞压力下降)。当已经达到指定增加的压降(Δp)时,在步骤816处,预测试活塞停止,并且终止压力下降。然后在步骤817处允许压力平衡,并在步骤818处持续周期所述周期 可以比在步骤817处的压力下降周期tpi长,例如, 在压力已经平衡之后,使点820处的稳定压力与在点810处的压力下降开始时的压力进行比较。此时,如图18中被显示为步骤819所示,进行判定是否重复循环。进行判定的判据是:是否平衡压力(例如,在点820处)与压力下降(例如,在点810处)开始时的压力差别大致与预计的压降(Δp)一致的量。如果是这样的话,则重复此流动管线膨胀期。
[0227] 为了重复流动管线膨胀期,例如,如上所述重新启动预测试活塞期和重复压力下降期,即,在820处开始预测试,在824处在大致相同的速度下使压力精确地下降相同的量(Δp),对于先前的循环持续时间826,在825处压力下降终止,并在830处稳定。再次,比较820和830处的压力以判定是否重复循环。如图17A中所示,这些压力明显不同,并且基本上与由流动管线中的流体的膨胀而产生的预计的压降(Δp)一致。因此,重复所述循环,即,830-834-835-840。重复“流动管线膨胀”循环直到连续稳定的压力差基本上小于施加/规定的压降(Δp),例如如图17A中的840和850所示。
[0228] 在连续稳定的压力差基本上小于施加/规定的压降(Δp)之后,可以再一次重复“流动管线膨胀”期,在图17A中被显示为850-854-855-860。如果在850和860处的稳定压力基本一致,例如在小倍数的压力计可重复性(gauge repeatablity)内,两个值中更大的一个值作为地层压力的第一估计值。本领域的普通技术人员将认识到图17A、17B和图18中所示的过程仅仅是为了说明。本公开的实施例不限于执行多少个流动管线膨胀期。此外,在连续稳定的压力差基本上小于施加/规定的压降(Δp)之后,任选地,重复一次或多次循环。
[0229] 从流动管线流体膨胀到来自地层的流动发生转变的点在图17A中被识别为800。如果850和860处的压力在规定的稳定时间结束时一致,有利地,可以允许860处的压力继续恢复并且使用先前部分所述的过程(见对图8的说明)以终止压力恢复,从而获得更好的地层压力的第一估计值。先前的部分说明了进行判定以继续调查阶段或执行测量阶段
864-868-869从而获得地层压力870的最终估计值的过程。在870处完成测量阶段之后,探头与井壁脱离,并且在874处压力在时间周期895中恢复到井筒压力并且在881处达到稳定。
[0230] 一旦在图17A和18中所示的调查阶段13b中获得地层压力的第一估计值和地层流动性,则如此获得的参数可以用于建立将在规定的测试时间内产生更加精确的地层参数的测量阶段14预测试参数。在先前的部分中已经说明了用于使用在调查阶段13b获得的参数设计测量阶段14预测试参数的过程。
[0231] 在图17A、图17B和图18中所示的实施例中,指定在流动管线膨胀阶段期间的压降(Δp)的大小。在可选的实施例中,如图19和图20中所示,指定在流动管线膨胀阶段期间的体积增量(ΔV)的大小。在此实施例中,流体的被固定和精确调节的体积的(ΔV)在控制流量下在每一个步骤中被提取,以产生可以由以下公式估算的压降:
[0232]
[0233] 在此实施例中使用的过程类似于图17A、17B和图18中所示的实施例所述的过程。在安置地层测试器之前,优选地确定泥浆压缩系数(没有示出)。在确定泥浆压缩系数之后并且在安置地层测试器之前,由工具测量的压力是井筒压力或泥浆流体静压201。
[0234] 参照图19A、19B和图20,在205处安置工具之后,启动图4中所示的预测试活塞118a。根据本公开的一个实施例,用于执行调查阶段13c的方法包括开始压力下降的步骤
210,在精确并且固定的流量下抽取流体直到预测试室114a的体积增加指定量ΔV的步骤
214。预测试室的体积的增量变化可以例如为0.2-1立方厘米的量级。本领域的普通技术人员将认识到指定的体积增量(ΔV)不局限于这些示例性体积,并且应该根据被捕集流体的总体积来选择。流动管线流体的最终膨胀产生流动管线中的压降。
[0235] 当已经达到预测试室体积的指定增量时,在215处停止预测试活塞118a并且终止压力下降。然后在217处允许流动管线中的压力平衡持续时间toi218,所述时间比216处的压力下降期tqi长,例如,toi=2tqi。在压力已经稳定之后(在图19A中的点220处被示出),在步骤219处,进行判定是否重复“流动管线膨胀”循环(图20中所示)。用于进行判定的判据类似于图17A和图18中所示的实施例所述的判据。即,如果在稳定或平衡之后的压力(例如,在点220处)明显不同于压力下降开始时的(例如,在点210处)的压力,并且压差基本上与由在流动管线中的流体的膨胀产生的预计的压降一致,则重复“流动管线膨胀”循环。
[0236] 为了重复“流动管线膨胀”循环,例如,在步骤220处重新启动预测试活塞,在步骤224处使流动管线精确地膨胀相同的体积ΔV,以及在步骤230处允许稳定压力。再次,如果220和230处的压力明显不同并且基本上与由流动管线中的流体的膨胀产生的预期的压降一致,则重复循环,例如230-234-235-240。重复“流动管线膨胀”循环直到连续稳定的压力(例如,如图19A中所示的230和240处的压力)差基本上小于由于流体在流动管线中的膨胀而引起的预计的压降。
[0237] 在连续稳定的压力差基本上小于预计的的压降(Δp)之后,可以再一次重复“流动管线膨胀”期,在图19A中被显示为240-244-245-250。如果在240和250处的稳定压力基本一致,取两个值中更大的一个值用于表示地层压力的第一估计值。本领域的普通技术人员将认识到图19A、19B和图20中所示的过程仅仅是为了说明。本公开的实施例不限于执行多少个流动管线膨胀期。此外,在连续稳定的压力差基本上小于预计的压降(Δp)之后,任选地,重复一次或多次循环。
[0238] 从流动管线流体膨胀到来自地层的流动发生转变的点在图19A中被识别为300。如果240和250处的压力在规定的稳定时间结束时在选择的范围(例如,小倍数的压力计可重复性)内一致,有利地,可以允许250处的压力继续恢复并且使用先前部分所述的过程(见图8)以终止压力恢复,从而获得更好的地层压力的第一估计值。在先前的部分中说明了进行判定以继续调查阶段或执行测量阶段250-258-259-260从而获得地层压力260的最终估计值的过程。在260处完成测量阶段之后,探头从井壁脱离,并且在264处压力在时间周期295内恢复到井筒压力并且在271处达到稳定。
[0239] 一旦在图19A和图20中所示的调查阶段13c中获得地层压力的第一估计值和地层流动性,则如此获得的参数可以用于建立将在规定的测试时间内产生更加精确的地层参数的测量阶段14预测试参数。在先前的部分中已经说明了用于使用在调查阶段13c获得的参数设计测量阶段14预测试参数的过程。
[0240] 在先前的部分中,概括了用于确定泥浆压缩系数的方法。泥浆压缩系数取决于其组分和流体的温度和压力。因此,泥浆压缩系数通常随深度变化。因此,理想的是在靠近其中将要执行测试的地方的位置处现场测量泥浆压缩系数。如果工具结构不允许如上所述确定泥浆压缩系数,现场泥浆压缩系数可以通过以下所述的迭代方法来估算。
[0241] 在根据本公开的一个或多个方面的方法中,地层测试器可以例如靠近套管鞋安置在套管中以建立与套管的流体密封。借助于图4中所示的预测试活塞118a执行被捕集在测试器流动管线中的井液的压缩和减压。以上参照图11A和11B说明了用于执行泥浆压缩系数测试的过程。一旦预测试活塞流量qp、压力变化的速度 和被捕集体积V是已知的,则可以由 估算泥浆压缩系数。
[0242] 在此具体实施例中,在执行压缩系数测量的实际垂直深度(因此温度和压力)可能明显不同于将要测量的地层压力的深度。因为钻井液的压缩系数受温度和压力的影响,因此必须对如此测量的压缩系数进行校正以估算钻井泥浆在将要执行测试的深度处的压缩系数。
[0243] 在根据本公开的一个或多个方面的方法中,在测量开始之前,即在图17A中所示的点801处,使用传统的压力和温度传感器获得井筒压力和温度信息。基于已知的钻井泥浆特性和现场温度和压力测量值,如图21中所示的图表可以被构造成用于实施温度和压力校正。可选地,本领域公知的分析方法可以用于计算校正因子,当应用到原始压缩系数测量值时,所述校正因子将提供在将要被测量的地层的深度处的现场流动管线流体压缩系数。例如,参见E.Kartstad和B.S.Aadnoy的“Density Behavior of Drilling of Drilling Fluids During High Pressure High Temperature Drilling Operations “IADC/SPE47806,1998。
[0244] 在本公开的一个或多个方面的另一个方法中,测量在地面得到(例如,泥浆槽)的样品在期望的井下温度和压力条件的范围内的压缩系数。井下条件下的现场泥浆压缩系数的估计值则可以根据本领域公知的方法由泥浆密度和泥浆压力以及泥浆温度之间的已知关系来估算。例如,参见图21和E.Kartstad和B.S.Aadnoy的“Density Behavior of Drilling Fluids During High Pressure High Temperature drilling operation”,IADC/SPE 47806,1998。
[0245] 图21对于油基泥浆和水基泥浆示出了流体压缩系数(Cm)与流体压力(p)之间的典型关系。实线10对于典型的油基泥浆示出了泥浆压缩系数随井筒压力的变化。虚线11对于典型的水基泥浆示出了泥浆压缩系数的相对应的变化。在地面上的油基泥浆的压缩系数由附图标记7表示。套管鞋处的油基泥浆的压缩系数由附图标记8表示。在套管鞋下方给定测量深度处的油基泥浆的压缩系数由附图标记9表示。压缩系数校正值ΔC表示套管鞋处的油基泥浆的压缩系数8与在测量深度处的油基泥浆的压缩系数9之间的差。在套管鞋8处得到的压缩系数测量值可以通过压缩系数校正值ΔC来调节以确定测量深度处的压缩系数9。如虚线11所示,水基泥浆的压缩系数的变化和相对应的压缩系数校正值可明显地小于油基泥浆的由实线10所示的校正值。
[0246] 如上所述,可以在本公开的实施例中使用直接在现场测量或由其它测量值外推的井下条件下的泥浆压缩系数,以提高由例如图11A中如图所示的调查阶段和/或测量值阶段获得的地层特性的估计值的精度。
[0247] 图22示出了在预测试操作期间获得的压力(P)与时间(t)图2200。此压力迹线类似于以上在图5中所述的预测试,但是具有更多的细节。以下参照图22说明用于预测试的通用过程,尽管要注意这种说明是示例性的,并且在不背离本公开的保护范围的情况下可以使用其它过程。
[0248] 在预测试开始之前,诸如探头(图4中的112a)的流体连通装置位于缩回位置,使得工具的内部受到井筒压力或流体静压Ph1,在2201处被示出。为了执行预测试,将流体连通装置压靠在井壁上以形成密封并建立与地层的流体连通。当探头接合井壁时,‘安置’流体连通装置,并且流动管线中的压力增加。这种压力增加是由于当将探头压入到井壁上的泥饼内时流动管线中的流体的压缩而产生的。此‘安置’动作具有安置压力(Pset)并且在图22中在2203处被显示。如所述,尽管可以不必总是这种情况,但是安置压力(Pset)可以比在2201处的井筒压力(Ph1)高。安置压力(Pset)相对于流体静压(Ph1)的相对位置对以下所述的适用性是不重要的。
[0249] 在图22中,点2204表示调查阶段的压力下降阶段2205的开始。这被称作膨胀压力(Pex),因为所述膨胀压力是在膨胀阶段刚刚开始之前测量的压力。点2204可以在井筒压力(Ph1)以上,或者点2204可以在安置工具之后降回到井筒压力(Ph1)或甚至降回到井筒压力(Ph1)以下。
[0250] 在压力下降阶段,位于工具内部并连接到流动管线(例如,图4中的119a)的测试活塞(例如,图4中的118a)移动,使得流动管线的体积增加。在这种情况下,以稳定并且已知的速度发生增加,但是如果期望则可以变化。当体积增加并且执行压力下降时,流动管线中的压力下降。此‘压力下降阶段’2205从2204延伸以在下降压力2209处终止压力下降。
[0251] 在第一压力下降期间的某点处,期望的是在工具的探头内被隔离的井壁上的泥饼(图1的附图标记4)将破裂,这将能够使流体从地层流入到探头流动管线内。当泥饼破裂时,并且如果地层具有充分的流动性,则流动管线中的压力可以产生细微压力恢复,在2206处被示出。通常,这发生在比稳定的2240处的井底压力(Psf)低的压力下,所述井底压力在泥饼破裂时通常对于操作者是未知的。因此,在2206处的泥饼破裂的压力(PMC)提供在*2240处稳定的井底压力(Psf)以及最终地层压力(Pf或P)存在的范围的初始指示。
[0252] 一旦泥饼破裂,如附图标记2206所示,压力下降沿2207持续直到流动管线中的压力在2209处达到下降压力(Pd1)。要注意的是除了在2206处泥饼破坏之外,大多数压力下降阶段(即,2205、2207)非常接近以上相对于图7所述的压力的线性下降。在2208处靠近压力下降阶段终止时,压力变化趋势变成非线性的。这是因为流体从地层流入到工具内,并且来自地层的流体的流量开始与由于活塞的运动施加的变化的体积流量相匹配。
[0253] 在被称作为‘压力下降阶段’2205的压力下降期间的最低压力被称作“下降压力”(Pd1)2209。具有多个方法用于确定何时压力下降停止。以上相对于图7说明了用于确定压力下降的终止的技术的一些示例。
[0254] 可以用于选择下降压力(Pd1)2209的一种技术基于如果泥饼破裂被检测则检测泥饼破裂的压力(PMC)2206。例如,如果检测到泥饼破裂,下降压力(Pd1)2209可以设定为在泥饼压力(PMC)2206以下的给定或预先选定值。
[0255] 在其它情况下,根本不会具体选择下降压力(Pd1)2209。相反,例如,根据在泥饼破裂2206之后探头流动管线的有效体积的变化来终止压力下降阶段。例如,压力下降阶段可以通过在泥饼破裂2206之后使活塞移动以空出选定体积来限定。在其中没有检测到泥饼破裂2206的情况下,可以基于通过移动活塞而被空出的流体的总体积来终止压力下降阶段。因此,可以指定固定流量和总体积。压力下降阶段将继续使活塞在固定流量下移动直到达到规定的总体积。此时,停止活塞,并且下降压力(Pd1)2209将取决于地层输送流体的能力和为预测试所选择的操作参数。
[0256] 一旦下降压力(Pd1)达到2209,工具中的活塞停止移动,并且工具中的压力传感器监测由于地层流体流动到工具而产生的压力恢复。这种压力恢复或压力恢复阶段2210从下降压力2209延伸直到达到最终的压力恢复2216。在压力恢复阶段2210期间,压力以渐近的方式朝向虚线2240处稳定的井底压力(Psf)恢复。要注意的是第一压力恢复阶段2210结束时的最终恢复压力(Pb1)2216被示出为小于稳定的井底压力(Psf)2240,但是所述最终恢复压力可以更大。压力恢复阶段2210可以在压力完全稳定之前(例如,当仅给预测试规定了短时期时)被终止。
[0257] 如图22中所示,执行两个连续预测试。如刚刚所述的称为‘调查阶段’的第一预测试在图22中的压力曲线从2204横跨到2216。调查阶段可以类似于例如以上相对于图2所述的预测试。可以在如上所述的第一预测试或调查阶段之后执行第二预测试或‘测量阶段’。另外的预测试可以如所期望地执行。
[0258] 第二预测试或‘测量阶段’在图22中从2216延伸到2231。如上所述,可以限制所述阶段的持续时间,并且可以基于判据设定这些阶段的终点。通常,与测量阶段(一个或多个)相比,调查阶段的持续时间较短,并且通常用于提供地层参数的估计值和/或设计用于实施测量阶段的判据。测量阶段可以具体地适于基于由调查阶段获得的结果实现预测试目标。通常,所述测量阶段的持续时间比调查阶段长,并且可以提供更加精确的结果。
[0259] 如以上相对于图7所述,预测试可以用于生成地层压力的估计值(Pf)和地层中的流体的“流动性”的估计值。流动性说明了地层流体是如何易于流入地层中。这可在评价从井开采油气的经济能力时是有用的。流动性被定义为地层的渗透率除以地层中的流体的粘度。因此,流动性M被定义为M=K/μ,,其中K是地层渗透率,而μ是地层流体的粘度。
[0260] 如先前相对于图7所述,地层的流动性的估计值可以由在水平延伸通过最终恢复压力的线以下和在压力下降和压力恢复曲线上方的面积(在图7中由325所示)确定。例如,在图22中,在通过第一恢复压力(Pb1)2216的水平线2242以下和在压力下降2205和压力恢复2210曲线的至少一部分以上的面积2251是流动性的指示器。例如,可以使用以上公式(1)估算流动性(K/μ)1,并且其中V1是线2242与线2205的交点(例如,在图22中的2261)与压力恢复阶段的端点(例如,图22中的2216)之间的预测试室体积的变化,而A是在曲线以下的面积(例如,图22中的面积2251)。
[0261] 图22还显示了在点2216与2231之间测量阶段的压力曲线。除了测量阶段可以而不是必须地具有更大压降并且通常具有更长的用于压力恢复阶段的时间之外,图22中所示的测量阶段类似于调查阶段(2204-2216)。可以基于如这里先前所述的调查阶段的结果设计用于测量阶段的判据。
[0262] 第二压力下降在图22中的点2216处开始,并且持续直到压力达到第二下降压力(Pd2)2219。类似于第一压力下降的最后部分2208,第二压力下降2217的最后部分2218显示出非线性。如同第一压力下降阶段2205一样,第二压力下降2217可以通过本领域公知的任意方法被终止。例如,第二压力下降阶段2217可以在膨胀预先选定的体积之后停止。此外,一旦已经达到预先选定的第二下降压力(Pd2),则可以终止第二压力下降2217。可以基于已经已知的关于井和地层的信息、来自由先前的预测试获得的信息、或来自在试验井中执行的测试的信息、或使用上述任意判据选择第二下降压力(Pd2)2219。
[0263] 可选地,可以基于如上所述的在调查阶段期间获得的信息终止第二压力下降2217。例如,可以基于在调查阶段2204-2216期间获得的压力数据选择为第二压力下降
2217所选择的体积流量和总体积。在另一个示例中,可以基于对在调查阶段2204-2216中获得的压力数据的分析具体地选择第二下降压力2219。用于终止第一和第二压力下降阶段的方法不旨在限制本公开。
[0264] 可以通过移动活塞以膨胀工具中的流动管线内的体积来产生第二压力下降2217。虽然可以使用另一个活塞,但是优选地用于测量阶段的活塞是与用于调查阶段相同的活塞。另外,可以使用如本领域所公知的用于降低压力的其它方法。用于执行压力下降的方法不旨在限制本公开。
[0265] 在压力下降阶段2217在点2219处终止之后,可以停止活塞,并且允许流动管线中的压力增加。这是第二压力恢复阶段2220。优选地,当执行多个预测试时,第二压力恢复阶段2220的持续时间比第一压力恢复阶段2210长。第二压力恢复阶段2220中的压力恢复直到在2231处的第二恢复压力(Pb2)。此第二恢复压力可以用作稳定的井底压力(Psf)2240的第二指示器。
[0266] 如同调查阶段一样,测量阶段的图表上位于第二恢复压力(Pb2)2231以下并在第二压力下降阶段2217和第二压力恢复阶段以上的面积2252可以用作地层中的流体的流动性的指示器。面积2252的值和点2216与点2231之间的预测试室体积的变化一起可以用于估算流动性。例如,以上公式(1)可以用于估算地层中的流体的流动性。可选地,本领域公知的任意其它方法可以用于确定流动性。
[0267] 在测量阶段之后(即,在第二压力恢复阶段2220在2231处终止之后),预测试活塞通常被部分伸出,平衡阀打开,并且流体连通装置从井壁缩回。然后再次使流动管线受到井筒压力。流动管线中的压力上升(在2232处)到井筒压力(Ph2)2233。
[0268] 在大多数情况下,在预测试开始时测量的井筒压力(在2201处的Ph1)类似于在预测试结束时测量的井筒压力(在2233处的Ph2)或与在预测试结束时测量的井筒压力(在2233处的Ph2)相同。要注意的是基于多种情况可能具有差异。例如,温度的变化可能影响压力测量。另外,如果在钻井时执行预测试,井眼中的液体动压力在泥浆泵运行时执行预测试的情况下可能会波动。其它因素可能影响井筒压力测量值(Ph1,Ph2)。
[0269] 要注意的是当在钻井操作期间执行预测试时,即使泥浆流可能产生噪音并且使井筒压力产生波动,但是期望保持泥浆泵运行。泥浆泵提供泥浆通过钻柱的流动,这允许使用泥浆脉冲遥测术。通过在执行预测试时使泥浆泵运行,可以发生与地面的至少某种水平的通信。
[0270] 在根据本公开的多个方面的操作中,使用数据压缩技术以利用将被通信的数据(例如,上述预测试数据等)填充预定通信信道容量(例如,在上述泥浆脉冲遥测术通道中可用于数据传输的带宽)。使用这种数据压缩技术,即使在数据通信信道例如由于低速数据传输而使带宽严重受限和/或带宽由于其它/另外的数据的传输而被耗费掉,也可以实时或近似实时地提供测试数据(例如,在钻井测试时由地层压力获得压力对时间数据)的稳固井口通信。例如,使用本公开的数据压缩技术,相对于图22在以上所述的足以精确地表示图中所示的图表的预测试的数据可以被实时或近似实时地通信给地面。
[0271] 可以使用稳固的数据通信以有助于在不需要移除地层测试工具并因此移除钻柱的情况下分析和/或控制钻井操作,和/或根据由预测试的结果获得信息等允许连续进行钻井操作和/或迅速修改钻井操作。然而,本公开不仅仅局限于上述预测试压力数据或以上刚刚所述的预测试数据的通信。例如,如果期望,本方法尤其可以用于通信预测试压力导数数据、预测试马达转速和体积、在取样操作期间液压泵的体积、来自流体分光仪的光密度、取样流的流体密度和/或粘度、以及与诸如收回和安置管线的压力(retract and setline pressure)的工具操作有关的信息或与工具的内部状态有关的信息。在地层测试工具不适于自动利用调查阶段数据构造测量阶段测试时,可以使用本公开的数据压缩技术以将足以精确表示图22中所示的图的调查阶段的数据实时或近似实时地通信到地面。可以在地面上分析此数据,用于在钻井操作停止以便实施预测试所规定的有限时间内构造预测试的测量阶段。在泥浆循环泵关闭的情况下实施压力测试、此时在执行测试期间在工具与地面之间不进行通信的情况下,使用这里所述的方法可以具有很大的优点。通过使用这些方法,在这种泵关闭测试期间由工具获得的数据的精确表示可以在允许关于工具的操作和井的状态做即时的判定之后被有效地传输到地面。虽然上述示例参考井口数据传输,但是应该认识的是可以相对于井下或其它数据通信应用本公开的原理。
[0272] 直接参考图23,示出了根据本公开的多个方面的提供数据压缩和通信的操作的高级流程图。如图23的流程图所示,在步骤3702处采集数据,所述数据例如可以包括上述预测试数据。例如,地层测试工具可以如上所述执行一个或多个测量以提供期望的数据。
[0273] 之后,在步骤3704处,优选地,使用以下更加全面所述的技术对采集的数据(例如,表示相对于测试过程所感兴趣的一部分的数据)的全部或选定部分进行抽取/压缩。要注意的是“抽取”在这里使用其最广泛的意思,包括减小信号离散序列或数据流中的样品的数量,并且不旨在局限于整个的十分之一(“十分之一”)。
[0274] 在步骤3704处提供数据抽取/压缩时,在用于通过数据通信信道通信的测试数据内优选地识别事件数据点。数据抽取器优选地利用这些事件数据点以识别采集数据内的另外的数据点,例如设置在用于优选地通过数据通信信道进行通信的事件数据点之间的曲线上的具体的数据点。优选地,选择另外的数据点以使事件数据点、另外的数据点、和相对于采集数据的通信所使用的开销数据以几乎尽可能地填充通信信道中的所有可用带宽。根据本公开的多个方面被填充的通信信道中的带宽可以是通信信道的整个带宽或没有以其它方式被利用、预定、或不能用于上述数据通信的信道带宽的一部分。
[0275] 在步骤3706处,对被抽取/压缩的数据进行编码,用于在通信信道内进行传输。对数据进行编码可以包括将位分组或量化和分配给数据,处理数据以提供误差检测和/或校正,将数据封入适当的输送容器中等。此外,如在步骤3706处提供的对数据进行编码的步骤可以包括将被抽取/压缩的数据追加到将通过通信信道被通信的其它数据,或使被抽取/压缩数据与所述其它数据交错(interleave)。
[0276] 在步骤3708处使用通信信道传输被编码数据。这种传输可以包括对载波的调制或用于将数据放置在介质上用于进行传输的其它公知的技术。在优选的结构中,编码数据被调制为脉冲,用于通过泥浆脉冲遥测通信信道进行传输。
[0277] 在步骤3710处,通过与通信信道通信的系统接收编码数据。例如,在其中地层测试工具已经执行测试(已经从所述测试采集了数据)的情况下,诸如耦合到通信信道的井口接收器的地面系统可以接收所述数据。在步骤3710处的接收可以包括对载波信道的解调或用于从传输介质提取数据的其它公知的技术。在优选的结构中,通过泥浆脉冲遥测通信信道的脉冲解调接收到的数据。
[0278] 在步骤3712处对接收到的数据进行解码。对数据进行解码可以包括拆包或解量化和重构数据,处理所述数据以检测和/或校正误差,展开或解封包来自输送容器内的数据等。而且,如步骤3712处提供的对数据进行解码的步骤可以包括分离期望的数据与已经通过通信信道被通信的其它数据。另外或可选地,在步骤3712处对数据进行解码的步骤包括诸如以下相对于图30所述的使用特殊函数将一个或多个反函数应用到被压缩数据。此外,对数据进行解码的步骤可以包括应用如下所述的由数据抽取器应用的“增长”函数的反函数。这样一个反函数的应用可以利用与应用到通过通信信道被通信的数据的函数有关的信息,或可以例如通过采用用于由发送数据的系统确定将被应用的函数的相同的算法而被独立确定。
[0279] 在步骤3714处,分析和/或使用解码或重构的数据。通常将解码数据添加到测井图。测井图可以具有位于钻机(例如,图1A中的钻机2)上的屏幕上的显示器的形式。测井图也可以具有存储在本领域公知的任意公知的存储装置中的印刷文献或电子记录的形式。例如,在其中地层测试工具已经执行了测试(已经从所述测试采集了数据)的情况下,诸如计算机或终端机的地面系统可以处理所述数据以将关于继续钻井操作、执行另外的测试、完成测试等的信息提供给井工程师或其它操作者。可选地或另外,随后可以存储和使用信息,例如用于确定油藏模型、评价储层的收益性、选择开采设备或用于其它应用。
[0280] 已经如图23中所示根据本公开的多个方面大致说明了提供数据压缩和通信的操作,以下参照图24-26,其中提供了关于优选的数据压缩技术的进一步细节。具体地,图24的流程图提供了关于图23的抽取/压缩数据步骤3704的优选的结构的细节。类似地,图25和图26的程序框图提供了关于图24的用于通信的抽取数据步骤3816的各种结构的细节。
[0281] 为了使读者更好地理解本公开,这里将参照地层预测试数据说明由图24-26的程序框图表示的操作,从而提供更加确实的示例性结构。然而,应该认识的是本公开不局限于与参考这里所述的示例性数据一起使用。
[0282] 以下参照图24,其中流程图在步骤3802处开始,在步骤3802处,选择或推导事件数据点用于进行通信。参照图22中所示的预测试数据,可以看出在所述预测试数据中示出了具体的事件。例如,预测试的测量阶段的数据包括与在预测试期间经历或与所述预测试相关联的具体事件相关联的数据点。具体地,数据点2216表示压力下降开始事件,数据点2219表示达到下降压力事件,而数据点2231表示恢复压力近似事件。可能与井预测试相关的其它事件包括以下事件:在测试之前井筒压力的识别2201、工具安置事件2203、预测试调查阶段的开始/流动管线膨胀的开始2204、泥饼破坏检测2206、调查压力下降的终止2209、调查恢复压力近似/预测试调查阶段的结束/预测试测量阶段的开始2216、测量阶段压力下降的终止2219、测量阶段恢复压力近似/预测试测量阶段的结束/最终地层压力的达到2231、以及在实施测试之后的井筒压力2233。这些事件可以不必存在于所有测试中,例如,可以不必存在于丧失密封或干测试等。这些及其它事件可以容易地被检测到(例如,开始或终止具体的测试操作,例如,接合预测试活塞、启动马达、获得马达的具体转速、接合工具、卸除工具等)或者在数据流中被相对容易地识别(例如,反向或快速变化的趋势、在一个或多个测量参数中的峰值或峰谷、达到一个或多个测量参数的稳定状态、达到超时等)。
[0283] 各种事件数据点可以被认为是相对于所执行的预测试尤其感兴趣的数据,或者另外可以表示在数据流中尤其感兴趣的数据。例如,上述事件数据点可以限定用于压缩和/或实时通信的数据的值或部分的区间。因此,图24的步骤3802优选地操作以选择或推导这些事件、事件点或数据点中的一个或多个,用于包含在表示全部数据流(例如,图22中所示的曲线的整个预测试或预测试测量阶段部分)的压缩通信中。
[0284] 在3804处,确定与上述事件点相关联的一个或多个值。例如,在事件数据点表示具体时间下的压力的情况下,可以确定每一个选定事件数据点的压力值和时间值用于进行传输。在另一个示例中,对获得的数据在取样时间或在取样时间之外进行外推以精确地确定趋势变化下的值或渐近值。在另一个示例中,通过使采集数据或采集数据中的走向平滑来确定选定数据点处的值,例如,如以下相对于图31-33详细所述。
[0285] 在图24中所述的操作中,在步骤3806处,每一个选定事件数据点的值被量化用于通过通信信道进行通信。例如,在传输之前,每一个事件数据点的值可以被量化用于进行编码。如果期望,可以提供非均匀量化的数据。例如,有利的是使用一个量化精度对位于一个区间内的数据点的值进行量化,而使用另一个量化精度对位于另一区间或多个区间中的数据点的值进行量化。可以使用压缩扩展器以基于期望的精度水平改变用于量化每一个事件数据点的值的量化精度。以下论述关于根据本公开的多个方面可以使用的压缩扩展器结构的细节。
[0286] 分配给由于对数据进行量化而产生的被抽取的数据点值的位的数量可以基于期望的精度。例如,在数据点表示压力和时间信息的情况下,由上述量化提供的位的数量可以根据以下法则计算:
[0287]
[0288]
[0289] 其中 表示比x大的最小整数,tacc和Pacc分别是期望的时间精度和压力精度,nbitstime和nbitspress分别是分配给被抽取的时间和压力的位的数量,而tmax和Pmax分别是最大压力值和最大时间值。
[0290] 在步骤3808处,相对于预定的信道容量(例如,在用于对预测试数据进行通信的通信信道中可获得的带宽)对被识别用于进行通信的数据(例如,事件点数据和与其传输相关联的任意开销数据)进行分析,以确定另外的预测试数据是否可以在通信信道内被通信。例如,泥浆脉冲遥测通信信道可以提供大约0.5位/秒到大约12位/秒,这取决于各种因素。例如可通过经验评价确定使用泥浆脉冲遥测可相对于任意具体井可获得的最大位速度。类似地,可确定完成数据通信的周期。例如,钻井操作可以被中断持续最大时间(例如,15分钟),如果期望在钻井操作恢复之前完成预测试操作和所有相关联的通信,预测试操作(从所述预测试操作捕获将被通信的数据)可能需要10分钟,从而留下大约5分钟用于数据通信(对于此示例忽略数据通信可以在预测试操作期间完成)。可选地,如有必要,数据传输可以与钻井的恢复同时发生。假设在此示例中泥浆脉冲遥测通信信道支持1位每秒并且假设此时没有其它数据正在通过通道被通信,300位的带宽能力可用于对预测试数据进行通信(假设5分钟传输时间)。在步骤3808处的操作优选地比较来自选定事件数据点的量化值的位和与所述选定事件数据点相关联的任何开销位(overhead bit)(例如,分组标题(packet header)、误差检测/校正位等)的数量与可用带宽容量,以确定容量是否保持用于另外的数据的通信。
[0291] 优选地,在步骤3810中确定与选定事件数据点的通信相关联的数据、和当前选定用于通信的任意其它数据的数量是否小于在用于这种通信的通信信道中可获得的容量。如果在通信信道中可获得另外的容量(或者如果具有超过足以允许对另外的数据进行通信的最小阈值数量的另外的容量),处理根据所示的流程图进行到步骤3816,在步骤3816中,优选地选择另外的预测试数据用于进行通信。以下图39和图40讨论了相对于可以用于选择这种另外的数据的各种数据抽取技术的细节。
[0292] 然而,如果在通信信道中不能获得另外的容量(或者如果没有足以允许对另外的数据进行通信的容量),处理根据所示的流程图进行到步骤3811,在步骤3811处,调节量化精度。例如,可以改变值的分辨率,以获得分配给数据点的较小位数和/或可以减小数据点的数量直到达到充分的带宽。
[0293] 在步骤3812处,对被选择位用于通信的数据(例如,选定事件数据点和选定另外的数据点)进行编码。相对于步骤3812的操作优选地与以上相对于图23所述的步骤3706相对应。在所述的结构的步骤3814处,传输编码数据。相对于步骤3814的操作与以上相对于图23所述的步骤3708相对应。
[0294] 以下参照图25,示出了根据本公开的多个方面的为将被通信的数据压缩提供数据抽取的操作的流程图。应该认识的是图25中所示的流程图的步骤可以作为用于在图24中所示的抽取数据进行通信步骤3816的一部分而执行。
[0295] 在执行图25中所示的流程图中使用的数据抽取器操作以通过使用用于选择进行通信的数据的一个或多个变量优化被选择为用于进行通信的具体数据和/或被通信的数据的数量。图25中所示的用于选择进行通信的数据的变量是压力变化(ΔP)和时间变化(ΔT),这与其中预测试数据点是压力对时间的示例相一致。然而,根据本公开的多个方面,可以使用其它变量来选择用于进行通信的数据作为压缩数据集。
[0296] 在步骤3902处,选择ΔP和ΔT的值。可以通过多种技术中的任一项选择这些变量的值。例如,初始可以为这些变量选择与数据的最高分辨率相关联的步长值(例如,与在采集测试数据中使用的采样率相对应),因为这种选择将提供最大信息。可选地,可以初始为这些变量选择考虑可能导致足以填充通信信道的容量的数据点选择的步长值。可以初始为这些变量选择考虑可能导致小于需要填充通信信道的容量的数据点的选择的步长值,使得迭代过程可以用于增加选定数据点的数量,以基本上填充通信信道的容量。换句话说,迭代过程可以包括选择、识别、和确定数据点以集中对备选数据点进行选择。可以使用历史信息、模拟、统计分析等选择这种步长值。压力步长值特别有利的初始选择为选择压力信道噪点的整数倍数,例如四倍或更多倍数,且由通过信号处理中所公知的方法压缩的压力迹线直接确定压力噪点。
[0297] 在优化ΔP和/或ΔT中,可以通过离散优化算法确定压力和/或时间步长,所述离散优化算法自动调节压力和/或时间步长大小以实现表示将被通信的预测试压力-时间迹线的位数的指定目标。
[0298] 在步骤3904处,选择在根据本公开的多个方面将被压缩的数据流中的数据点。在图25的结构中,被分析用于进行选择的数据点是位于来自参考点(这里是事件数据点)的步长值(这里是ΔP或ΔT)的曲线上的数据点。优选地,被分析用于进行选择的数据点是在两个选定事件数据点(例如,在图24的步骤3802处选择的事件数据点)之间的曲线上的数据点。因此,根据这种操作,可以在压缩数据流中容易地表示由两个选定事件数据点界限的一段数据集曲线。可以根据上述压缩多段曲线,从而提供预测试数据的分段压缩。以下通过参考图27可更加容易地认识上述原理。
[0299] 图27显示与用于根据本公开的多个方面进行通信的数据集相关联的、基本上与图22的预测试相对应的曲线。数据点4102-4138显示为初始作为被用于进行通信的压缩数据集。即,如果与通过通信信道通信这些数据点中的每一个相关联的位数小于或等于通信信道的有效容量(例如,这可以在图24的步骤3810处进行确定),数据点4102-4138将被选择为用于进行通信。数据点4102-4138包括例如可以已经在图24的步骤3802处选择的事件数据点4102、4112、4114、4124、和4136。在对曲线的一部分进行压缩用于进行通信的情况下,可识别用于界限曲线的所述部分的事件数据点,并且在与上述变量相关联的步骤中选择沿所述曲线的数据点用于进行通信。
[0300] 因此,在将要压缩由事件数据点4114和4124界限的曲线的所述部分的情况下,可以识别事件数据点4114并分析数据流,以选择具有ΔP或ΔT比事件数据点4114的对应值大或小的值的下一个数据点。在所示的示例中,数据点4116的压力值与事件数据点4114的压力值之间为ΔP(尽管时间变化保持小于ΔT)。使用选定数据点4116作为参考再次重复此,因此选择压力值与数据点4116的压力值之间为ΔP的数据点4118(再次,尽管时间变化保持小于ΔT)。数据点4122显示选择具有不同于前述选定数据点的时间值的时间值ΔT的数据点的示例(虽然压力变化保持小于ΔP)。应该认识的是可以根据上述容易地抽取整个数据集或所述整个数据集的多个部分。
[0301] 一旦已经选择了数据点,在步骤3906处确定与上述选定数据点相关联的一个或多个值,并且在步骤3908处所述值被量化用于通过通信信道进行通信。使用与相对于选定事件数据点使用的相同技术或通过使用另一种技术可以完成对所述值的量化(步骤3806)。
[0302] 因为本公开的上述结构的操作最大化了可在通信信道中获得的带宽内通信的数据的数量,因此使用上述变量对另外的数据点进行选择优选地是一种迭代过程。因此,在使用上述抽取技术选择另外的数据点之后所述的示例返回到图24的步骤3810,用于确定与选定事件数据点的通信相关联的数据、和选定的另外的数据点的数量是否小于在用于这种通信的通信信道可获得的容量。如果选定数据点的通信没有填充通信信道有效容量,优选地,通过对一个或多个上述变量进行调节(例如,减少步长ΔP和/或ΔT)重复抽取步骤,以增加另外的选定数据点的数量。类似地,如果选定数据点的通信将超过通信信道有效容量,优选地,通过对一个或多个上述变量进行调节(例如,增加步长ΔP和/或ΔT)重复抽取步骤,以减少另外的选定数据点的数量。
[0303] 用于进行调节的上述变量中的具体的一个变量的选择和所提供的调节量可以基于多个条件中的一个。例如,在这里所述的示例中,其中,压力和时间步长用于选择另外的数据点,可以期望的是在时间相关变量已经被选择为取样数据的最大或最小“超时”的函数的情况下调节压力相关变量。然而,可以根据本公开的原理以任意数量的方式调节任意或所有这种变量。此外,可以选择不同的变量用于在不同的时间进行调节(例如,连续迭代),和/或根据本公开的原理的调节不同量。
[0304] 以下参照图26,示出了根据采用本公开的原理的方法的另一方面的为将被通信的数据压缩提供数据抽取的操作的流程图。应该认识的是图26中所述的流程图的步骤可以作为用于在图24中所示的抽取数据进行通信步骤3816的一部分被实施。此外,应该认识的是参照图26所述的数据压缩技术可以用作参照图25所述的数据压缩技术的替换技术或与所述数据压缩技术结合使用。例如,图25的数据压缩技术可以用于曲线的一段,而图26的压缩技术可以用于曲线的另一段。可以使用最适于具体数据特征的压缩技术。
[0305] 在执行图26中所示的流程图中使用的数据抽取器操作以通过使用适当的“增长”函数选择用于通信的具体数据点和/或被通信的数据点的量来优化被选择用于通信的具体数据和/或被通信的数据的数量。通过这种数据抽取器实施的函数可以利用例如线性、对数、指数、球形、或几何级数,或任意其它适当的类时间函数,例如,时间或生成的体积。例如,在由数据点表示的曲线显示在曲线的开始时的迅速变化值并且所述值的变化率随后在曲线中下降的情况下,可以期望的是执行数据点选择技术以沿曲线展开选定数据点,从而避免在曲线中捕获不相称的大百分比的后期数据点和在曲线中捕获低百分比的早期数据点,其中正在发生大多数变化。数据抽取器对增长函数的应用可以通过选择在曲线的弧上更加均匀地分布的数据点来用于优化对用于通信的具体数据的选择。通过参照图28可更加容易地认识上述原理。
[0306] 图28显示了与用于根据本公开进行通信的数据集相关联的大致与图22的预测试测量阶段相对应的曲线。虽然在图28中示出了数据点4202-4230,但是此示例显示了相对于曲线的一部分对数据点进行选择。因此,数据点4212-4228初始被显示为被选择为用于进行通信的压缩数据集。即,如果与通过通信信道通信这些数据点中的每一个相关联的位数小于或等于通信信道的有效容量(或曲线的这部分的通信可获得的通信信道的有效容量)(例如,这可以在图24的步骤3810处进行确定),数据点4212-4228将被选择为用于进行通信。数据点4212-4228包括在图24的步骤3802处被选择用于进行通信的并因此可以不必重复通信的事件数据点4212和4228、和使用所述结构的抽取器选择的另外的数据点4214-4226。因此,在压缩通过事件数据点4212和4228界限的曲线的所述部分的情况下,优选地实施用于提供沿这些界限事件数据点之间的曲线部分的选定数据点的相对均匀分布的增长函数。
[0307] 在根据图26的示例提供数据抽取中,在步骤4002处确定将要被选择的多个另外的数据点。例如,在将要抽取曲线的具体部分的情况下,可以在步骤4002处确定界限曲线的所述部分的事件数据点之间的多个另外的数据点。可以通过从通信信道可用带宽减去选定事件数据点和相关联的通信开销以及与另外的数据点的通信相关联的任意开销来确定将要被选择的另外的数据点的数量。
[0308] 在步骤4004处,确定两个事件之间的期望区间。可选地,并且如这里所述,期望的区间可以由图28中所示的诸如t0和tn的两个时间来界限。此区间可以从一个事件跨度到另一个事件,或者可以跨度两个事件之间的任意部分。例如,在其中数据点包括压力和时间信息的上述示例中,可以选择时间间隔,当应用选定的增长函数时,所述时间间隔将有助于选择在步骤4002中确定的另外的数据点的数量。可以使用如图28中所示的时间步长Δt1确定区间的开始。例如,可以使用诸如1秒的预先选定时间延迟确定区间的开始。可以使用tn-t0的百分数确定区间的结束。可以以类似的方式确定区间的开始。
[0309] 应该认识的是本公开的操作不限于任意具体的参数或用于相对于使用增长函数选择另外的数据点的区间。然而,优选地相对于数据流的一部分实施使用增长函数的抽取,其中,数据点值单调增加或减少,以提供更加均匀分布的选定另外的数据点。
[0310] 在步骤4006处确定增长函数因子,所述增长函数因子将导致选择在步骤4002处确定的数据点的数量。在已经确定增长函数因子的情况下,所述示例的步骤4006进一步提供时间级数,从而识别与将被选择的另外的数据点相关联的时间。
[0311] 在步骤4008处,确定与在步骤4006中提供的时间级数相对应的数据点的压力值。应该认识的是通过应用这种增长函数,除了数据的抽取之外的数据压缩可以通过对部分数据集进行通信来实现。在上述示例中,在数据点表示压力对时间的情况下,可以使用上述几何级数以再现数据点的相关时间方面,从而仅允许对数据点的压力分量进行通信。
[0312] 因此,在所述示例的步骤4010处,量化与选定另外的数据点相关联的压力值和在确定时间级数中使用的增长函数因子,以便进行传输。如所期望的,在步骤4010处可以量化另外的或可选的信息。例如,在可以相对于数据抽取实施各种不同增长函数的情况下,可以量化指示所实施的具体的增长函数的信息。类似地,在选定数据点之间的期望区间,相对于增长函数所使用的具体数据点参数,等,对于通信的接收端来说是未知的的情况下,可以量化相对于这些参数的信息,以便进行通信。可以使用与相对于选定事件数据点相同的技术或使用另一种技术完成所述值的量化(步骤3806)。
[0313] 因为本公开的上述结构的操作最大化了可在通信信道中获得的带宽内通信的数据的数量,因此使用上述变量对另外的数据点进行选择优选地是一种迭代过程。例如,可以根据图24和图26所述的步骤对曲线的多个部分进行抽取。因此,在使用上述抽取技术选择另外的数据点之后所述的示例返回到图24的步骤3810,用于确定与选定事件数据点的通信相关联的数据、和选定的另外的数据点的数量是否小于在用于这种通信的通信信道可获得的容量。如果选定数据点的通信没有填充通信信道有效容量,优选地,相对于曲线的这部分或曲线的另一部分重复抽取步骤,以增加另外的选定数据点的数量。类似地,如果选定数据点的通信将超过通信信道有效容量,优选地,相对于曲线的这部分或曲线的另一部分重复抽取步骤,以减少另外的选定数据点的数量。
[0314] 在已经说明了根据本公开的原理如图24-26中所示的提供数据压缩和通信的操作的情况下,参照图29,其中示出了相对于如可以使用根据本公开的数据压缩扩展器实施的量化技术的细节。可以使用图29的流程图的步骤以例如在步骤3806(图24)、步骤3908(图25)、和步骤4010(图26)中的任一个中提供数据的量化。
[0315] 在步骤4302处开始根据图29的流程图对数据进行量化,在步骤4302处,确定在数据集的动态范围内的两个或更多个区间。在步骤4304处,例如相对于以上公式45选择期望的量化精度。在步骤4306处量化区间末端。在步骤4308处根据两个或更多个区间确定变换,并且在步骤4310处将所述变换应用到数据集中的至少一个点。在步骤4312处量化变换数据集的结果。
[0316] 继续参照诸如图22中所示的具有压力和时间值的地层预测试数据的示例,假定压缩扩展器用于量化包括井筒压力值Ph1和Ph2、探头安置压力值Pset、膨胀压力值Pex、恢复压力Pb1,Pb2、下降压力Pd1和Pd2以及泥饼破裂的压力值PMC的示例性数据集。如以下所使用的,术语“井筒压力”可以表示静水压和/或液体动压力中的任一个。此数据集中的值从等于Pd2的最小值Pmin延伸到等于Pset的最大值Pmax。
[0317] 由预测试提供的重要结果是井底稳定压力Psf的近似值。此压力的量化精度Pacc优选地被选择位用于量化此压力值、或者量化被识别为达到最终地层压力的至少一个事件(例如,在图27中的数据点4136或图28中的数据点4230)。更具体地,量化精度可以被设定为1磅/平方英寸(psi),用于在1磅/平方英寸分辨率下生成估算的井底稳定压力的测井图。
[0318] 在数据集中的数值分布分散在一个或多个区间的情况下,优选地使用根据图43的流程图操作的数据压缩扩展器以将位中的一小部分分配给在上述稀疏的一个或多个区间中的数据集的数值。继续参照其中地层预测试数据具有与所述地层预测试数据相关联的压力和时间值的示例,如图22中所示,可以认识的是除了调查阶段的压力下降部分中的初始流动管线膨胀之外,压力分布图基本上在Pmin(最小压力值)与Psf(例如与最终恢复压力Pb2近似)之间变化。可以有利的是缩小其中压力分布是稀疏的区间[PsfPmax],从而例如有效地表示在所述区间[Pmin Psf]中的压力。通过参考图30更加容易地认识上述原理。如图30中所示,可以通过多线性变换映射压力,所述多线性变换利用具有斜率小于1的线性函数压缩区间[Psf Pmax]。因此,在映射之后的区间[Pmin Psf]占据总区间的百分数比区间[Psf Pmax]大。在使用均匀量化器进行量化之后,整体效果与将大部分位分配给落入[Pmin Psf]中的值的非均匀量化器相同。
[0319] 更具体地,示例性变换基于涵盖示例性数据集的动态范围的两个区间[Pmin Psf]和[Psf Pmax]。这些区间的末端是优选地使用以上相对于公式45所述的多个位以精度Pacc被量化的Pmin(Pd2)、Psf(Pb2)、和Pmax(Pset)。示例性数据集中的其它数值首先通过图30的变换被映射。这种变换将区间[Pmin Psf]映射到区间[Vmin V]中,并将区间[Psf Pmax]映射到区间[V Vmax],其中例如Vmin等于0,V例如等于Psf-Pmin,而Vmax是通常小于Pmax-Pmin的确定值。
[0320] 将这种变换应用到示例性数据集的除了先前已经被量化的Pd2,Pb2和Pset之外的元素。优选地,以精度Pacc量化所变换的结果。要注意的是用于这些被变换数值的位数由以下公式给出:
[0321]
[0322] 本领域的技术人员将认识到如果期望则区间的数量可以大于二。此外,应该认识的是可选地可以使用除多线性变换之外的变换。例如,可以使用具有可变斜率的单个单调函数替换多线性函数,或者除了多线性函数之外,可以使用所述单个单调函数。也可以由一系列的数据点值对这种单调函数进行参数化。具体地,如果将被传输的变量具有涵盖多个数量级的值(例如,渗透率),可以将量化应用到变量的表示,而不是应用到变量本身的数值。在渗透率的情况下,一旦已经确定了希望涵盖的范围,则量化可以应用到数值的对数表示的指数。在这种情况下,指数的量化的精度是关键的。
[0323] 在已经说明了例如可以使用根据本公开的原理如图29-30中所示的数据压缩扩展器实施的量化技术的情况下,直接参考图31和图32,其中,图31和图32示出了相对于例如可以根据本公开的原理实施的数值确定技术的进一步细节。可以使用相对于图31-32所述的技术或以下所述的其它平滑技术以提供诸如在步骤3804(图24)、步骤3906(图25)、和步骤4008(图26)中的任一个内的选定数据点处的数值。
[0324] 图31显示了预测试的压力(P)对时间(t)图的压力恢复曲线2900的一个示例。这种压力恢复曲线2900可以与详细说明的图22的压力恢复2210或2220相同。压力恢复曲线表示在假想预测试中由井下工具随时间记录的所有压力数据点。由于压力传感器的操作、井底温度、和流体从地层流出的方式的变化,所述数据显示关于总体趋势的变化。然而,当观察时,数据表现为形成稍微平滑的压力恢复曲线。
[0325] 在一些情况下,有利的是沿压力恢复曲线的发展计算压力恢复曲线的选定点处的平滑压力值和压力导数或斜率。可以使用用于选择具体点的任意方法。在图31中,在压力恢复阶段中在时间为零处的第一数据点被选择为第一选定数据点2901。基于诸如压力步长值、时间步长值、时间增长函数等的各种判据选择其余的数据点。在此示例中,使用几何时间级数选择点2902-2907。可选地,所有采集的数据点可以用于进行分析。
[0326] 一旦选择了数据点,可以关于选定点确定压力的平滑值和导数(即,压力恢复曲线的斜率)。可能有用的是选择关于选定数据点的范围并且拟合曲线与在所述范围内的所有数据点。可以使用拟合曲线估算在选定数据点处的曲线的平滑值和导数。
[0327] 图32显示了压力恢复曲线2900的部分3000。选择数据点3001,将关于所述数据点估算平滑值和斜率值。数据点3001具有时间t0和压力P0。对压力区间(δ)选择关于数据点3001的范围。可以任意或通过多种不同的方法实现区间(δ)的选择。优选地,区间(δ)被选择为多个信号噪点。在其它情况下,区间(δ)可以被选择为多个压力传感器分辨率。通过在这些方法中的任一个中选择区间(δ),可以确保区间上的多个点之间的压差表示实际压力变化,而不是数据的统计变化。
[0328] 压力范围的上限和下限分别与压力PL和PH相对应,其中PL=P0-δ而PH=P0+δ。在图32中,压力PL和PH大约分别与压力恢复点3003和压力恢复点3005相对应。
[0329] 一旦定义了压力范围,通过区间拟合曲线。在一个示例中,平滑函数与所述范围中的数据相拟合。“平滑函数”是与数据进行拟合以产生接近所述范围中的数据的平滑曲线的任意函数。可以使用接近所述数据的任意函数。在一个示例中,平滑函数的数学表达式是诸如公式31中所示的时间的二次函数:
[0330] p(t)=a(t-t0)2+b(t-t0)+c (31)
[0331] 其中t0是选定数据点的时间,而a、b、和c是将被拟合的常数。用于拟合二次公式式的一种方法是如本领域所公知的稳固的最小二乘法。拟合所述公式以及特殊方式的公式的方法不旨在限制本公开。图32中的线3010表示已经与所述范围中的数据进行拟合的二次公式的曲线。
[0332] 在其中t=t0的点处,公式31中的压力将是常数c。此外,在得到公式31的解析导数的情况下,可以看出公式31在点t0处的导数是常数b。因此,通过拟合诸如公式31的二次公式与所述范围中的数据,压力恢复曲线在t0处的压力和斜率的“平滑”值可以分别被估算为常数c和b。因此,t0处的压力可以被估算为第三常数(即,公式31中的c),而t0处的压力导数可以被估算为第二常数(即,公式31中的b)。可以对用于压力恢复的数据集中的每一个选定数据点执行如图所示用于图32中的选定点3001的这种方法。例如,这种方法可以用于确定图30中的压力恢复曲线在点2902-2906处的“最典型的”压力值和斜率。这种方法还不局限于沿压力恢复曲线选定数据点,并且可以应用到在预测试曲线上的其它地方选择的其它数据点。
[0333] 可能重要的是知道预测试阶段中预测试曲线的端点数据点处的压力值和/或斜率。在一些情况下,选定数据点可以是曲线的最后记录的数据点(即,图30中的2907)。在其它情况下,选定数据点可以靠近压力变化趋势迅速变化的事件(即,图30中的2901)。
[0334] 应该认识的是输送到地面操作者用于并入测井图内的数值不局限于平滑值和斜率。例如,可以通过诸如曲率的曲线拟合确定其它数据,并且传输所述其它数据。此外,可以仅传输平滑值或斜率中的一个。另外或者可选地,如以下进一步所述,可将通过过滤技术确定的数值应用到关于数据点所选择的区间。
[0335] 再次参照图32,可以有利的是使用过滤技术关于选定点确定压力恢复曲线的平滑值和斜率。一旦关于选定数据点限定了压力范围或曲线部分,可以基于存在于所述范围中的数据点的数量选择过滤器。因此,可确定数据点3003(与选定区间的下限PL相关联)与选定数据点3001之间的数据点的数量NL。还可以确定选定数据点3001与数据点3005(与选定区间的更高界限PH相关联)之间的数据点的数量NH。可以基于NL、NH或NL和NH选择过滤器长度L。例如,可以根据NL和NH中的最小值Nmin(例如,由L=2Nmin+1给出)选择零相位有限脉冲响应(FIR)过滤器长度。
[0336] 过滤系数通常取决于选定的过滤器长度L。一些过滤器可以更加有效地过滤短区间上的数据,而另一些过滤器可以更加有效地过滤长区间上的数据,并因此选择所述过滤器。可以通过利用低通零相位FIR过滤器(例如,标准化锥形窗口或核加权过滤器)的褶积(convolution)获得选定数据点3001处的压力的“最典型的”压力值。更具体地,可以使用Welch窗口、Epanechenikov核或Savitsky Golay过滤器。图33中示出了可用于获得选定点处的压力的平滑值的过滤器的示例性示例。要注意的是过滤器可以包括正数值和负数值(未示出)。
[0337] 一旦选择了过滤器,使用同样在本领域中所公知的过滤法(即,褶积)关于选定数据点对记录的曲线进行过滤。在时间t0的被过滤的曲线的值则可以被传输。
[0338] 可选地或另外,可以通过过滤技术获得选定数据点3001处的压力导数或曲线斜率。例如,可以使用选定长度L的导数过滤器。导数过滤器通常具有基本上与所关心的信号的频带中的脉动jω成比例的频率响应H(ω)。例如,可以通过对低通滤波器求微分来得到导数过滤器。图33B中示出了FIR反对称导数过滤器的示例性示例。
[0339] 虽然仅在图33A-33B中示出了FIR过滤器,但是本领域的技术人员将要认识的是可以使用其它类型的过滤器。例如,无限脉冲响应(IIR)过滤器可以用于确定曲线平滑值、曲线斜率值或曲线的其它特征。此外,可以使用正向和反向过滤。另外,过滤可以用于对两个采集的时间之间进行数据内插。除了过滤之外,可以使用诸如异点检测和剔除的其它噪点剔除技术。
[0340] 图34-45示出了用于分析可能在地层测试期间遇到的压力迹线的技术的另外的示例。测试的某些部分可以显示已经在测试期间发生的异常特性、缺陷、误差或事件的指示。在测试的执行期间或所述测试的执行之后可以识别一个或多个置信度标记。可以分析这些置信度标记中的一个或多个以确定这种异常特性、缺陷、误差或事件是否在测试期间已经发生。然后,这些置信度标记可以用于确定由所执行的测试和/或所述测试的潜在的数据和解释获得的结果的置信度水平。
[0341] 通常,置信度标记用于识别在实际预测试期间测量的压力响应与理想条件或原型预测试中的相对应的期望响应之间的相似性。如这里所使用的,置信度标记可以例如用于检测这种相似性的程度。还可以获得关于测试条件或其它井下特征的另外的信息。
[0342] 图34是示出了用于确定归因于地层测试的置信度水平的方法2300的流程图。这种方法包括执行如上所述的至少一个预测试的步骤2302。可以在观察或不观察全部预测试的情况下执行评价。在一些情况下,可以执行一个或多个预测试。在其它情况下,在不执行另外的预测试的情况下可以终止测试,以便可以尤其在测试持续时间内达到总测试目标。
[0343] 在预测试期间,可以在步骤2304处确定一个或多个置信度标记。如以下更加全面地所述,具有多个不同类型的置信度标记和用于确定置信度标记的技术。基于所确定的置信度标记,在步骤2306处发现通常导致不可能恢复的灾难性事件已经发生。例如,标记可以表示工具已经失效、或者不显著地,工具与井壁之间的液体密封已经丧失。如果是这样的话,在步骤2308处,可以进行判定以在接收所述信息之后尽可能快地终止测试。可以终止测试,并且随后可以不再执行其它的测试,可以终止测试和重新开始所述测试,或者可以允许所述测试继续。
[0344] 可以在一个或多个预测试期间识别一个或多个置信度标记,并然后在步骤2310处分析所述一个或多个置信度标记。可以分析在单个预测试期间获得的多个标记。可选地,可以分析一个或多个预测试两端的一个或多个标记。在步骤2312处,此分析可以用于确定一个或多个预测试的置信度的总体测试水平。
[0345] 如果期望,可以在预测试2314期间或所述预测试之后调节预测试和/或井眼操作。例如,可以期望的是调节工具的操作或在井下重新构造工具以获得更好的测量值并且继续测试过程。在另一个示例中,第一预测试置信度评价的结果可以用于改变第二测试的参数。在一些情况下,可以期望的是根据由一个或多个置信度标记和/或预测试置信度获得的信息优化测试过程。在步骤2316处,可以利用调节的参数执行另外的预测试。
[0346] 如果期望,一个或多个置信度标记可以可选地在步骤2318处被编码,用于在地层测试器与地面之间的通信信道中进行传输。在步骤2320处使用通信信道传输编码标记。对标记进行编码的步骤可以包括附加诸如相对于图23-26中所述的压缩数据的其它数据,或与所述其它数据交错。在优选的结构中,通过泥浆脉冲遥测通信信道传输编码标记。在步骤2322处,通过与通信信道通信的系统接收编码标记。例如,耦合到通信信道的井口接收器可以接收传输信号。在步骤2324处对接收到的标记进行解码。如在步骤2324处提供的对数据进行解码的步骤可以包括分离期望的标记与已经通过通信信道被通信的其它数据。在步骤2326处,显示和/或使用解码标记。例如,在地层测试已经对已采集数据执行测试的情况下,诸如计算机或终端机的地面系统可以显示所述数据以将相对于执行另外的测试、完成测试等的信息提供给井工程师或其它操作者。通常将解码标记添加到测井图。地面操作者或其它操作者可以使用所述标记以调节如以上相对于步骤2314所述的测试操作,这因此可以通过地层测试器在井下完成或通过操作者或自动地面系统在井口完成,或者通过本领域所公知的其它方法完成。
[0347] 图35-46说明了用于确定如以上的方法2300的步骤2304中所述的一个或多个置信度标记的各种技术。每一个置信度标记提供关于预测试的一个方面的信息。可以使用各种技术确定这些置信度标记。
[0348] 图35和图36A-36B说明了用于根据压力比较技术确定置信度标记的方法2400。在此示例中,在预测试期间的不同时间处测量的压力的相对比较和/或次序可以用于检查预测试是否正在如所期望地执行。
[0349] 图35显示用于根据压力比较技术确定置信度标记的方法2400。置信度标记可以用于检查何时从最高变化到最低,预测试顺序中的每一个被识别事件处的压力应该具有具体的次序。在步骤2402处识别至少一个预测试中的至少两个点。则可以在步骤2404处确定每一个的相对应的压力。可以将所述点确定为预测试中的具体事件,例如,相对于图22所述的事件。
[0350] 例如,在图22中,安置压力(Pset)2203通常是最高压力。因此,预期正好在探头已经压缩泥饼层之后和在压力下降循环之前流动管线中的压力将是预测试期间测量的最高压力。预测试的公知特性可以用于识别诸如下降压力(2209)及其它的被识别事件。
[0351] 根据此信息,可以建立在预测试顺序中的每一个有效事件处测量的压力的次序。原型预测试以数学的方式被表示如下:
[0352] Pset>Pex>(Ph1≈Ph2)>(Pb1≈Pb2)>PMC>max(Pd1,Pd2) (26)[0353] 其中Pset是安置压力(例如,在图22的事件2203处测量的压力水平),Pex是膨胀压力(例如,在图22的事件2204处测量的压力水平),Ph1和Ph2是在测试之前和之后的井筒压力(例如,分别在图22的事件2201和2233处测量的压力水平),PMC是泥饼与井壁分离时的流动管线压力(例如,在图22的事件2206处测量的压力水平),Pb1和Pb2是最终的井底恢复压力(例如,在图22的事件2216和2231处测量的压力水平),而Pd1和Pd2是在压力下降结束时的压力(例如,在图22的事件2209和2219处的压力水平)。
[0354] 然后,在步骤2408处,对识别的压力进行比较以确定所述识别的压力是否以预计的次序发生。基于采集的预测试数据点是如何很好地与标准预测试的预期次序相对应,可以对置信度标记分配给定值。例如,可以根据如公式26中列出的次序的不满足或满足设定置信度标记。可选地,如以下进一步所述,可以根据预测试中被识别的点的测量的压力值设定置信度标记。
[0355] 可以进一步改进这些关系中的一些。例如,是否在安置工具时已经建立了密封的指示器可以用公式表示为:Pset-Ph1>D1,其中D1是具体的工具、地层、和泥浆类型的压力特性,并且可以具有预定水平的值。可以根据这种改进的压力比较技术设定改进的置信度标记。
[0356] 公式26的压力比较的改进的另一个示例可以基于以下关系:Pex-Ph1<(Pset-Ph1)/m,其中m是通常大于或等于2的预定数字。如果满足这种关系,“渗漏”泥饼可能被怀疑并因此可以设定另一个置信度标记。在这种情况下,可以利用如下限定的增压技术进一步检查恢复压力。
[0357] 在公式26的压力比较的改进的又一个示例中,井筒压力或流体静压(Ph1或Ph2)与恢复压力(Pb1或Pb2)的值的比较可以产生是否以过平衡或非过平衡的方式钻井的指示。基于对这种比较的无效或有效的又一个置信度标记可以基于(Ph1,Ph2)>(Pb1,Pb2),公式26中的不等式提供可以用于确定具体预测试是否有效的判据。
[0358] 在一些情况下,公式26中所表示的次序可以不被满足,但预测试仍然满足。例如,在欠平衡井中,在井筒压力或井眼中的钻井液的流体静压(Ph1,Ph2)通常低于地层压力(Pf)的情况下,井筒压力或流体静压(Ph1,Ph2)和恢复压力(Pb1,Pb2)的这些值将相反。此外,钻井操作可能使恢复压力(例如,Pb1,Pb2)高于井筒压力(Ph1,Ph2),从而指示潜在的危险操作条件。另外,如果泥浆泵正在一个点而不是在另一个点处运行时,在预测试的开始和结束时测量的井筒压力或流体静压(Ph1和Ph2)可以不同。因此,公式26中的压力比较提供可能的故障的指示。在一些情况下,另外的数据和/或分析(例如,改进的压力比较技术)可以提供充分的信息以得出在预测试中是否已经发生故障的结论。
[0359] 图36A-36B显示可能由在预测试期间遇到的问题产生的压力(P)对时间(t)的压力迹线的示例。这些问题可以使公式26中的条件不被满足。例如,图36A显示密封丧失的预测试的压力曲线2501。在压力下降期2502之后,压力恢复期2503开始。在压力恢复期2503期间,压力已经迅速恢复2504到井筒压力或流体静压2505。这表明建立在泥饼上的密封可能已经丧失,从而使井眼的压力渗漏到流动管线中。这里,在接近应该为恢复压力处测量的压力将基本上与在测试之前测量的井筒压力相同,并且使公式26不被满足。
[0360] 图36B显示预测试的另一个压力迹线2511。在压力下降期2512之后,压力保持在下降压力2514。这表明来自地层的流动没有进入工具内。这可能是由于流动管线阻塞或地层缺少可移动流体。再次,恢复压力将太低,并且基本上与下降压力相同,并且将使公式26不被满足。
[0361] 可以将置信度标记分配给图36A和图36B中所执行的预测试。图36A和图36B的迹线还可以指示突然故障。在这种情况下,可以期望的是如之前以上至少相对于图7所述的在测试进行以完成之前终止所述测试。在图34的任选步骤2308中,终止测试。当期望时可以重置工具并且可以重新执行测试。
[0362] 图37描述了用于根据参数比较技术确定置信度标记的方法2600。在此示例中,可以对不同预测试中的相同的测量参数(例如,噪点)进行相对比较以确认预测试是否正在预计的范围内执行。
[0363] 在此方法中,在步骤2602处识别来自第一预测试的至少一个参数。然后在步骤2604处识别来自至少一个另外的预测试的至少一个参数。然后在步骤2606处比较来自不同预测试的相对应的参数。然后在步骤2608处确定相对应的参数是否在预定范围内重复。
例如,限定噪点区、或其它传感器性能特征,并且比较来自不同预测试的相对应的参数以确认所述参数是否在限定的性能范围内重复。
[0364] 当执行多于一个的预测试时,预测试之间的比较可以提供关于将与预测试结果相关联的置信度水平。例如,如果2216处(图22)的第一恢复压力(Pb1)与2231处(图22)的第二恢复压力(Pb2)非常一致,所述一致性可以指示有效的测试。在这种参数比较技术中,可以根据第一和第二恢复压力在可接受范围内的条件设定置信度标记,例如:
[0365] |Pb1-Pb2|≤mmax(δ,η) (27)
[0366] 其中m是乘数,而max(δ,η)表示工具量规重复性(一个或多个)(δ)和与测量值相关联的噪声(η)中的最大值,所述测量值可以由在工具的操作期间获得的其它数据来确定。因为测量噪声通常大于传感器的固有噪声,因此通常需要通过本领域所公知的方法测量“不工作”的噪点。
[0367] 乘数m可以被设定为用于具体测试的适当数字。例如,m可以在其中泥浆泵正在运行并且噪声较高的情况下被设定为大约大于或等于2的数字。如果噪声非常高,m可以被设定为3或4。在其中泥浆泵关闭并且几乎没有噪声的情况下,m可以被尽可能低地设定为1。本领域的技术人员将认识到可以基于具体的测试情形修改乘数。此外,如果执行多于两个的压力恢复期,公式27可以被修改成包括除了第一和第二预测试之外的恢复压力。例如,如果执行三个压力恢复,公式27可以包括第一和第三压力恢复或第二和第三压力恢复。公式27中使用的具体压力不旨在限制本公开。
[0368] 可以在两个不同的预测试之间进行比较的另一个参数比较是压降响应的比较。第一预测试的压降响应是恢复压力(Pb1)与下降压力(Pd1)之间的差与压力下降速度(q1)的比值。因此,在这种第二参数比较技术中,可以根据两个压降响应之间的比较设定置信度标记,如以下被表示为:
[0369]
[0370] 其中e1和e2表示可以根据具体的测试情形选择的可接受的方差。要注意的是公式28中的中间项的第二半部分是第二压降响应的倒数。理想地,两个压降响应将近似相同,一个压降响应和另一个压降响应的倒数的乘积将近似等于一。通过使一个压降响应乘以另一个压降响应的倒数,可以将方差应用到所述乘积以评价预测试结果的置信度。
[0371] 可以在预测试之间进行的又一个比较是流动性之间的比较。如图22的描述中所述,可以使用公式1估算在第一压力下降顺序期间的流动性(K/μ)1和在第二压力下降顺序期间的流动性(K/μ)2。因此,在这种技术中,可以根据不满足或满足以下公式29中所述的条件来设定另一个置信度标记:
[0372]
[0373] 再次,e3和e4表示可以根据期望的结果和具体的测试情况选择的可接受的方差。要注意的是表示流动性由其估算的压力下降-压力恢复顺序(由所述压力下降-压力恢复顺序估计流动性)的数字的参考数字用作整个流动性项的下标。因为在这种流动性估算中没有区分渗透率或粘度,因此下标数字没有单独用于表示这些参数。
[0374] 例如,如果在第一和第二预测试期间几乎获得静流,则公式28和公式29中计算的比值可能非常类似。在这种情况下,比值可以接近一。然后可以选择置信度标记以指示高置信度水平。相反,在比值不接近一的情况下可以选择较低的置信度标记。在比值不接近一的情况下可以选择较低的置信度标记的情况下,可以合并多个置信度标记以选择参数的哪一个值可以最好地表示实际值。将容易认识的是例如由公式28和公式29表示的标记可以成对地应用到包括多于两个预测试的测试中。
[0375] 图38示出了用于确定与参数预测技术有关的置信度标记的方法2700。这种技术用于确定预测试是否如预计的那样执行。在一些情况下,可以期望的是使用由一个或多个测试获得的先验知识对预测试的参数的估计值进行预测。在图38中确定的估算的和/或计算的参数还可以用于确定各种井下条件。
[0376] 方法2700包括在步骤2702处对预测试的参数的估计值进行预测。可以选择诸如流动性、压力变化等的任意参数,然后在步骤2704处由在预测试期间采集的数据确定这种参数的计算值。然后在步骤2706处比较估算的参数和计算的参数。然后,在步骤2708处评价被比较的参数之间的差。可以根据所述评价赋予置信度标记。
[0377] 在一个示例中,参数预测技术可以用于确定可能影响测试结果的气体或其它可压缩流体在流动管线中的存在。参数预测技术的这种示例还可以被称为流动管线膨胀技术。如果工具中的流动管线(例如,图4中的119a)有诸如气体的可高度压缩的流体,则压力测量可能受到不利的影响。因
[0378] 有利的是确定这种流体是否存在于流动管线中。流动管线中的气体的一个指示器例如可由在流动管线膨胀期间的压力下降曲线的预测或估算斜率与压力下降曲线的实际斜率的比值来获得。例如可以根据流动管线中的泥浆的压缩系数的先验知识确定预测的斜率。当可以与图4中示意性地显示的工具一起使用时,泥浆压缩系数在预测试之前可以由数据库或泥浆特性与压力和温度的相关性已知,或者可以通过钻井液的分离测试来确定。
[0379] 这种比较以下可以使用公式30以数学的方式被示出。例如如果气体存在于流动管线中,压力由于流动管线流体压缩系数的预测的变化率与计算的压力的变化率的比值可以被表示如下:
[0380] (如果气体存在于流动管线中)
[0381] (30)
[0382] 其中 表示在压力下降期期间测量的压力的变化率, 表示估算的压力的变化率,Vt是在预测试期间流动管线的体积的估计值,例如流动管线的初始总体积加上在预测试中使用的体积的一半,Cm是钻井液的压缩系数,而qn是流动管线的变化的体积流量(例如,由于使连接到流动管线的诸如图4中的活塞118a的活塞移动而产生)。
[0383] 当左侧数量接近一时,流动管线中的流体的接近预计的钻井液的压缩系数。在这种情况下,如果有的话,在流动管线只有一点点气体。然而,如果在流动管线中具有大量气体,则测量的斜率将会比预测的斜率小很多。在这种情况下,公式30中的比值将明显地小于一。因此,在这种流动管线膨胀技术中,可以根据对公式30的不满足或满足来设定置信度标记,或者可选地,可以将置信度标记设定为预测的压力下降斜率与测量的压力下降斜率的比值。
[0384] 当在管线内检测到气体时,可以降低预测试结果的预计的置信度。在一些情况下,在已经从流动管线清除气体之后可以执行第二预测试。在其它情况下,执行另一个预测试可能是不切实际或不可能的。在这些情况下,操作者可以降低在气体存在于流动管线中的情况下所执行的预测试的结果的置信度或者重新评价所述结果。例如,如果对于不同垂直深度处的一系列测试中的一个预测试怀疑气体存在于流动管线中并且在所述深度处的地层压力的值表现出增加,操作者可以基于来自其它深度处的预测试的数据评价地层,而不是基于来自在流动管线中检测到气体的位置的数据评价所述地层。
[0385] 应该理解的是用于计算在流动管线膨胀期间测量的压力变化率(pressure rate)的方法不旨在限制本公开。可以由压力曲线斜率、压降等确定测量压力变化率。公知的技术包括曲线拟合、线性回归、代数计算等。此外,所述技术不局限于公式30的表达式。例如,可以由与图30的表达式在数学上等价的表达式确定置信度标记。
[0386] 图39示出了用于使用趋势分析技术确定置信度标记的方法2800。所述方法包括在步骤2802处沿预测试压力曲线的一部分选择一个或多个数据点。例如,可以使用靠近压力恢复曲线端部的数据点用于压力恢复的分析。优选地,在步骤2804处选择关于数据点的区间。在一些情况下,所述区间关于选定数据点定位。在其它情况下,所述区间通过将来自选定数据点的一侧的数据映射到数据点的另一侧来产生。以下进一步描述对区间的选择。在步骤2806处,例如,可以处理区间中的数据点,以分析噪点水平和/或从选定区间的数据剔除噪点。应该注意的是在不背离本公开的情况下可以使用各种处理技术,在步骤2808处,可确定压力曲线的诸如斜率、曲率等的一个或多个趋势特征,然后可以分析所述一个或多个趋势特征以确定置信度标记。
[0387] 利用这种技术,可分析预测试的一部分的诸如压力恢复的特征以确定在预测试中的一个或多个数据点处的压力变化趋势是否如表现出如所预计的特性。在一个示例中,关于测试的一部分的最后一个点的诸如斜率、和/或压力变化率(增加)的特征可以用于指示稳定性。在另一个示例中,可以分析关于测试的这部分分布的数据点的特征。
[0388] 例如,接近压力恢复结束的压力曲线可以在一些情况下相对水平或充分平坦,和/或压力变化率可以较小或接近零。这可以表示压力已经稳定并且达到地层压力,并且最终压力是地层压力的良好估计值。在其它情况下,压力变化率可能变大(增加或减小),这可以指示还没有达到地层压力。因此,可以根据接近压力恢复结束时的压力变化趋势将置信度标记赋予预测试。在这种局部应用到压力恢复期结束的示例性趋势分析技术中,置信度标记被设定为在接近压力恢复期结束的压力恢复曲线的斜率。可选地,将会根据接近压力恢复期结束时的压力恢复曲线的斜率与阈值之间的比较设定置信度标记。这种信息可以用于终止或继续测试,例如,直到达到稳定性。这种信息还可以用于确定预测试还没有达到稳定,并因此使得质量降低。
[0389] 因此,重要的是已知压力恢复曲线的在压力恢复期(即,图31中的2907)中的最后被记录的数据点处的斜率。在一些情况下,选定数据点可以在曲线的端部处。在这种情况下,可以期望的是实际上扩展数据以生成区间。用于扩展数据的一种方法是将数据扩展为奇函数。
[0390] 图40A、40B和40C显示了采用方法2800的压力恢复曲线的部分3100。压力恢复期的最后被记录的数据点3101包括在压力恢复部分3100中。数据点3101与时间t0处的压力P0相对应。在这种情况下,没有在选定数据点3101以上扩展的数据点。选定数据点3101以下的区间3107中的数据可以通过关于P0的压力范围(δ)来限定。压力的上限和下限在压力PL和P0之间,其中PL=P0-δ。存在大致与PL相对应的数据点3102。
[0391] 压力恢复部分3100实际上可以延伸超过数据点3101以产生关于选定数据点3101但同时适当地考虑数据中的噪点的区间。如图40B中所示,产生转折点3106。优选地,转折点可以被定义为时间t0处的压力的“平滑”值。例如,可以使用迭代最小二乘法拟合方法关于数据点3101形成第一切线3104。然后,通过关于转折点3106做对称通过选定数据点3101对数据点3102进行映射。映射的数据点3103定义压力PH。然后可以在如图40C中所示的数据点3101与数据点3103之间的区间产生数据的虚拟集。表示在范围的下端处的数据点3108的映射的虚拟数据点3105现在在关于选定数据点3101的范围的上端。
[0392] 例如,使用以上相对于图32或33A-33B所述的一个或多个平滑近似法,可以在最后的数据点3101处估算压力恢复的斜率。
[0393] 通过例如沿预测试的压力恢复部分分析多个数据点处的局部趋势可将如上所述的“局部”趋势分析技术自然地扩展到“整体”趋势分析技术,这种方法可以与如以上至少相对于图7所述的观察压力的次序一样简单。此外,这种方法可以包括沿压力恢复在选定点处观察相对于时间的次序压力导数或诸如先前至少相对于公式4所述的当量差。在这种情况下,可以通过如相对于图32、和33A-33B并且更加具体地相对于图40A-40C所述的局部趋势分析为压力恢复部分的端点获得所述每一个相对于时间的次序压力导数或所述当量差。
[0394] 例如,重新参照图31,可以有利的是沿压力恢复曲线的发展在选定点处计算压力恢复曲线的压力曲线趋势。可以使用用于选择具体点的任意方法。在图31中,在压力恢复期中时间为零处的第一数据点被选择为第一选定数据点2901。根据诸如压力、压力增量、时间、几何时间级数等各种判据选择剩余的数据点。在此示例中,使用几何级数选择点2902-2907。然后可以将局部趋势分析技术应用到选定数据点2901-2907,从而为每一个数据点提供一系列平滑的压力值、和平滑的斜率值。
[0395] 如果平滑压力显示单调增加趋势同时相对应点处的压力导数为正数并且在压力恢复结束时以非常小的值单调降低,将有良好的置信度,即,最终恢复压力(例如,Pb1)是稳定井底压力(Psf)的良好表示。整体趋势分析技术的一个示例可以以数学的方式被描述如下:
[0396] dp/dt(tk)>0 (49)
[0397] 其中tk是以上所述的选定时间;而dp/dt是例如利用相对于图32、33A-33B、40A-40C以及公式31所述的数据扩展和平滑近似法计算的相对于时间的压力导数。
[0398] 然而,如果导数显示为正数,并且几乎是恒定值,则可以怀疑具有渗漏。渗漏可能较小,使得通过压力迹线的肉眼检查而不容易被检测到。在这种情况下,可以将很少或几乎不将置信度赋予最终恢复压力的值。可以类似地诊断和评价表示异常动作的其它情况(例如,压力升高到最大值并然后以恒定负斜率降低的情况)。
[0399] 因此,在这些整体趋势技术中,可以根据在沿预测试的一部分的选定点处的一组局部趋势确定置信度标记。可选地,可以根据在预测试的一部分期间的压力的增加趋势和/或在预测试的一部分期间压力导数的降低趋势设定另一个置信度标记。如果传输相对较少的导数数据,则可以在地面上进行置信度标记的评价,或者在其中可获得更多数据用于进行分析的情况下,可以在工具的井下处理器中自动进行置信度标记的评价。
[0400] 图41示出了用于使用离散分析技术确定置信度标记的方法3200。所述方法包括沿预测试的一部分选择数据点的步骤3202。在步骤3204处选择关于数据点的区间。可以任意选择在所述方法中选择的压力范围或区间。另外,压力范围可以被选择为在选定数据点以上和以下的区间,其中,所述区间基于压力测量的噪点或压力计分辨率。压力区间可以被选择为噪点和压力计分辨率中的最大值的倍数;例如乘数可以是四。在不背离本公开的保护范围的情况下,可以使用用于选择关于选定数据点的压力范围的其它方法。然后在步骤3206处确定区间中的参考曲线。在步骤3208处关于参考曲线确定所述区间上的数据点的方差。在一些情况下,比较方差与阈值或噪点水平以确定所述方差是否满足一些判据,例如,以下由公式33定义的判据。可以根据如何更好地满足判据来赋予置信度标记。
[0401] 在方法3200的示例性实施例中,参考曲线是在区间内的中间压力水平的平直水平线。然后,利用例如以下公式32计算关于中间压力水平的区间上的方差。此方差是曲线关于选定数据点的平坦性的表征。可以分析区间中的关于中点的数据的方差以确定置信度标记。
[0402] 公式32显示了用于计算方差 的一种方法:
[0403]
[0404] 其中pk是在区间中的第k个点处的压力(即,关于t0位于中心的时间间隔),p(t0)是中间压力水平,而N是区间中的点的数量,优选地为奇数。公式33显示用于比较方差与阈值的一种方法:
[0405]
[0406] 其中m是乘数,而max(δ,η)表示工具分辨率的倍数(δ)和与测量值相关联的噪点(η)中的最大值。乘数m可以被设定为用于具体测试的适当的数字。在一个示例中,m被设定为4。本领域的技术人员将认识到可以选择m以适于具体的应用。因此,可以根据对公式33的不满足或满足来设定置信度标记。
[0407] 具体地,在接近压力恢复曲线的端部(例如,图7的350、或图22中的2210、2220)选择的区间的方差可以用作置信度标记。相对较低的方差表明压力几乎为常数,并且恢复压力(例如,Pb1,Pb2)近似于在井底压力(Psf)处被稳定。如以上相对于图7所述,这种置信度标记可以用于例如终止压力恢复期。
[0408] 在方法3200的另一示例性实施例中,可通过拟合诸如二次多项式函数的多项式函数与选定区间中的数据点来获得参考曲线。为了进行此,可以使用相对于图32或图40A-40C所述的技术。在方法3200的又一个示例性实施例中,可通过在选定区间上过滤测量曲线来获得参考曲线。为了进行此,可以使用相对于图33A所述的过滤器。本领域的技术人员将认识到多种方法可以用于平滑或去噪曲线的一部分并且这些方法可以用于确定参考曲线。然后可以计算关于参考曲线的方差以确定置信度标记。
[0409] 这里同样地,通过例如沿预测试的压力恢复部分分析多个数据点处的“局部”离散可以将如以上所述的“局部”离散分析法自然地被扩展成“整体”趋势分析技术。这种方法可以与如沿预测试的一部分观察在公式32中限定的方差的发展一样简单。例如,期望方差沿预测试的压力恢复期单调降低。当发生此时,置信度标记可以被设定为指示压力恢复如所期望地恢复并且预测试结果的置信度可以较高。
[0410] 图42提供了一种用于根据模型相关技术确定置信度标记的方法3300。所述方法包括在步骤3302处选择参数化系统响应函数。还可以在步骤3304处选择一个或多个参数化异常函数。例如,在诸如渗漏的异常被怀疑的情况下,可以选择期望的参数化异常函数。在3306处还选择成本函数。如下所述的成本函数定义参数化函数与数据之间的误差。
[0411] 在步骤3308处,可以优化参数化函数和/或异常函数的参数中的一个或多个以减小成本函数。在步骤3310处,可以比较优化的成本函数与预定值。在步骤3312处,还可以比较优化参数与预定值。这些比较中的一个或两个都可以用于确定一个或多个置信度标记。
[0412] 如果可以得到例如表示预测试压力恢复的参数化函数使得所述参数化函数近似表示实际压力恢复的特性,可以根据置信度标记解释如此得到的模型函数的参数。
[0413] 公式34显示用于模拟例如压力恢复的参数化函数的一个示例:
[0414] p(t)=P(t;Λ,Г)=F(t;Λ)+A(t;Г) (34)
[0415] 其中F(t;A)表示系统(例如,地层和工具)压力响应;Λ是表示系统的响应的一系列参数;A(t;Г)表示诸如渗漏、压力漂移等异常情况的模型;而Г是描述异常的一系列参数。例如,F(t;Λ)可以是表征地层中的球形流的综合影响的函数、本领域公知的函数、和工具库。可选地,用于系统压力响应的简单而不是必须精确的函数可以被写成为:
[0416]
[0417] 还可以选择参数化异常函数。在一个示例中,用于诸如逐渐渗漏的异常的模型可以被写成为:
[0418] A(t;Г)=γ(t-tγ)H(t-tγ) (36)
[0419] 其中H是在其自变量是负数的情况下值为零而在其自变量是非负数的情况下值为一的海维塞德阶跃函数。
[0420] 可以确定用于参数化函数和异常函数的参数。对于公式35和公式36来说,参数Λ,Г的列表被定义如下:
[0421] Λ={psf,po,tβ,β}
[0422] Г={γ,tγ}
[0423] 其中psf是估算的稳定的井底压力;p0是在压力恢复期开始时的压力;tβ是压力恢复期开始时的时间;γ是渗漏项的斜率(以下更详细地进行说明);tγ是估算开始以考虑渗漏的时间;而β是与地层和工具参数有关的压力恢复时间常数。可以由以下公式确定压力恢复时间β:
[0424]
[0425] 其中Ωs是考虑井眼的曲率对压力响应的影响的形状因子(参见公式2);rp是探头的半径;Vt是工具的总有效体积和流动管线体积加上预测试体积的一半;(K/μ)是地层中的流体的流动性;而Cm是占据工具流动管线的流体的压缩系数。
[0426] 在一些情况下,可能有用的是包括压力恢复项和渗漏项。在这种情况下,当井下工具中的探头在压力恢复期间的某一点处的没有与地层完全密封时,参数化函数可以与恢复压力更加近似地匹配。在这种情况下,井眼中的压力(Ph)使泥浆渗漏到流动管线中。这可以由除了正在被测量的井底压力(Psf)之外的源人为地增加流动管线中的压力。在其中探头与井壁进行有效密封的这些情况下,可以将渗漏参数(γ)简化为零。
[0427] 可以类似地识别和说明其它异常动作。例如,因为在测试之前或在所述测试期间停止循环,因此可观察到识别其中在压力恢复期间压力下降的动态滤失情况。在这种情况下,公式36中的渗漏参数γ是负数。
[0428] 一旦选择了参数化函数(例如,公式34-37),可以比较由参数化函数生成的作为时间函数的压力曲线与测量的压力数据。可以调节参数化函数中的参数,使得函数的曲线与压力数据更加近似地匹配。优选地,优化所述参数,使得参数化函数与数据尽可能近似地相匹配。
[0429] 参数优化算法的一个示例是当记录数据时最小化参数化函数的值与实际数据之间的误差。可以如公式38说明用于获得响应参数的优化过程:
[0430]
[0431] 其中Qk(Λ,Г)是进一步如下所述的成本函数,而N是所记录的数据的数量。
[0432] 优化可以包括改变可行或预测范围中的参数中的一个以确定哪一个参数值将产生最小误差,可以对所有参数重复此过程以进一步降低误差。在一些情况下,优化可以包括同时改变所有参数,并且由先前优化的值可以重复优化直到所有参数在指定范围内。优选地,使用诸如Levenberg-Marquardt程序的标准技术执行优化。还可以由本领域公知的其它优化技术估计方法确定模型函数的参数。
[0433] 公式39中示出了可以用于优化参数的成本函数的一个示例:
[0434]
[0435] 其中公式39中的成本函数的示例是数据(pk表示第k个压力数据点,而tk表示相对于测试部分的开始的与所述压力数据点相关联的时间)和参数化函数的函数(P(tk;Λ,Г),表示参数化函数在记录第k个压力数据点时的值)。公式39中的成本函数的示例还包括加权项w(tk)。差项(测量压力与参数化函数预测值之间的差)乘以加权项。在公式39中,可以选择加权项以将更多的权重给予数据的一些部分,例如选择w(tk)=(1+tk)强调了例如接近压力恢复结束时的参数化函数的方差。加权项允许在数据的早期部分的一些不拟合,但是却强调接近压力恢复结束时的参数化函数的拟合。通过在标记k上添加项(Ok(Λ,Г))以涵盖预测试的诸如压力恢复期的期望部分来计算成本函数的最终值。
[0436] 例如,图43显示了图表3400。在预测试期间采集的数据被显示为曲线3401。对参数化函数的最佳拟合被显示为曲线3403。曲线3403没有与数据中的早期部分的数据进行很好地拟合,但是曲线3403与后面部分数据非常近似地拟合。通过使用加权项,靠近数据结束的拟合比靠近数据开始的拟合更明显。
[0437] 在这些模型相关技术中,可以根据由如步骤3310所示的成本函数的最小值来设定置信度标记。例如,当成本函数的优化值较小时,测量数据点和在步骤3302处并可能在步骤3304选择的参数化函数近似匹配。这可以指示预测试的正在被研究的部分的确根据预计的函数表现。在这些情况下,置信度标记可以被设定为告知预测试的所述部分的形状已经以一定的置信度被辨别出。另外,成本函数的优化值中的较大值可以指示没有辨别出预测试的所述部分的形状,因此,置信度标记还可以被设定为不同的值。
[0438] 可选地,可以根据如步骤3312所示的最好地说明实际压力恢复的优化参数的值设定其它置信度标记。例如,当在压力恢复期间检测到渗漏时,压力恢复的最终压力可能是有问题的。渗漏的指示将是公式36中的参数的优化值(γ)的振幅。因此,可以根据参数的优化值(γ)设定置信度标记。
[0439] 此外,这种方法可以用于确定地层压力的改进的值。在一些情况下,与在压力恢复结束时记录的压力相比,参数化函数的优化能够更加精确地预测稳定的井底压力。例如,公式(35)的优化值psf可以是比恢复压力Pb1或Pb2更加精确的地层压力值。
[0440] 在方法3300的又一个实施例中,对如图44A和44B所示的由测量压力的数据点构造的在曲线以下的区域进行分析。这种可选的方法3300对于分析相对于对于很长的持续时间来说是负数或正数的系统响应具有偏差的异常是有利的。利用此方法,在系统响应周围波动的系统的扰动可以降低对分析的影响。
[0441] 图44A显示由地层测试器例如在压力恢复期期间获得的压力曲线4510。在此图中,假设在时间(tγ)时在封隔器两端正在发生渗漏。在没有渗漏的情况下,压力曲线4510已经朝向如由虚线4511所示的井底压力Psf稳定。在方法3300的这种实施例中,对曲线4510以下的区域进行分析。更具体地,选择多个压力恢复持续时间T。对于多个持续时间T中的每一个来说,计算曲线4510以下的横跨压力恢复的开始(被假定为t=0)与持续时间T之间的面积Ab(T)。例如,图44A显示对于持续时间T的具体值来说曲线4510以下的面积,所述持续时间大于被认为已经发生渗漏的时间(tγ)。此面积是面积A1 4530和面积A2 4540的总和。面积4530是在时间(tγ)时没有渗漏的情况下计算的面积。面积4540是由于在时间(tγ)时存在渗漏而计算的另外的面积。
[0442] 然后可以通过对多个持续时间T中的每一个将计算的面积Ab(T)绘制成为持续时间T的函数来构造曲线。图44B中示出了这种曲线的示例。曲线4520是在存在渗漏的情况下的曲线,而虚曲线4522是在没有渗漏(或至少当渗漏可以忽略或不能检测到渗漏)时产生的曲线。如图44B中所示,当在压力恢复期中不能检测到渗漏时,曲线将以渐近线的方式接近具有斜率psf和截距 的直线4521。此外,如图44B中所示,当可检测的渗漏发生在压力恢复期时,曲线4520与渐近直线4521分离。因此,可以通过分析通过将计算的面积Ab(T)绘制成为持续时间T的函数获得的曲线来检测渗漏。
[0443] 在步骤3302处,可以例如根据公式45选择参数化系统响应函数:
[0444]
[0445] 其中psf是估算的稳定的井底压力;p0是在压力恢复期开始时的压力;T是相对于压力恢复的开始所参考的压力恢复持续时间;而β是压力恢复时间常数。
[0446] 在步骤3304处,可以例如根据公式46选择参数化异常函数:
[0447]
[0448] 其中γ是渗漏项的斜率;而Tγ是开始进行估算以考虑渗漏的时间。
[0449] 在步骤3306处,选择成本函数O,例如:
[0450] O(psf,p0,β,Tγ,γ)=∑T(Ab(T)-A1(T,psf,p0,β)-A2(T,Tγ,γ))2 (47)[0451] 其中psf,p0和β是表示系统的响应的参数;γ和Tγ是描述异常的参数;Ab(T)是由测试期间测量的压力值计算的面积;而A1和A2是分别在步骤3302和3304处选择的函数。
[0452] 在步骤3308处,优化参数psf,p0,β,γ和Tγ的值以减少成本函数。可以使用任意的最优化算法。在一些情况下,利用参数的优化值计算的函数将与图44B的曲线4520近似匹配。
[0453] 图45中的方法3500描述了另一种用于使用量规比较技术确定预测试的置信度的技术。此方法包括在步骤3501处使用至少两个量规执行预测试。在步骤3502处沿由预测试生成的曲线选择区间。在步骤3504处确定每一个量规的噪点。可以沿所述区间确定噪点。可以使用来自一个或两个量规的噪点。在步骤3506处确定量规之间的方差。在步骤3508处比较方差与量规的噪点。可以根据结果赋予置信度标记。任选地,如果具有明显的总方差,所述方法可以包括识别引起诸如温度变化的方差的井下现象。
[0454] 在一个示例中,用于执行预测试的典型的地层测试工具可以包括应变压力计和石英压力计(例如,图4中的123a、120a可以包括应变压力计和石英压力计)。这两种类型的压力计具有不同的操作原理,因此所述压力计可以对相同的情况具有不同的响应。具体地,与石英压力计相比,应变压力计往往对压力变化可作出迅速反应,但应变压力计可能具有较差的绝对精度,并且可能通常噪声较高。此外,应变压力计可能受温度变化的影响较小。相反,石英压力计可以比应变压力计更加精确,但是石英压力计可以更加易于受到温度变化的影响;此外,与应变压力计相比,石英压力计可能对压力的变化反应更慢。如本领域所公知的其它类型的压力计可以代替石英和/或应变压力计而与此方法一起使用。
[0455] 可以在步骤3506处计算应变压力计与石英压力计之间的方差。例如,两个压力计之间的差CS,Q可以被定义
[0456]
[0457] 其中pk(Q)是通过石英压力计测量的第k个压力数据点,pk(S)是由应变压力计测量的相对应的压力数据点,n是可以由操作者选择的指数(例如,n可以是2),w(tk)是通常被选择以将更多的加权给予时间较晚的数据的加权函数,而N是例如在压力恢复结束时在步骤3502处所选择的区间中的数据点的数量。任选地,在步骤3506之前可以对多组数据点中的一组应用补偿。可以将补偿应用到来自第一压力计或第二压力计的数据,以更好地匹配(align)或覆盖压力计的测量的响应。补偿可以是压力计之间的压差的测量值。补偿可以可选地是压力计之间的时滞的测量值。因此,如果第一压力计数据被补偿,则以公式40中所示的同样的方式但是使用第一压力计的补偿压力数据代替实际第一压力计数据计算总的补偿方差。可以采用诸如互相关或本领域公知的其它方法的任意方法以自动确定最佳补偿。对哪一个数据进行补偿并且如何对所述数据进行补偿不旨在限制本公开。
[0458] 识别来自不同压力计的不同响应可以有助于确定压力计是否在井下失效。此外,如果压力计的响应在具体的区间上类似,则这将增加预测试的最终结果的可信程度。因此,应变压力计与石英压力计之间的方差可以用作预测试结果的置信度的指示器。如果CS,Q的值在选定值以下(例如,在步骤3504处计算的局部噪点或“最差”的压力计(通常为应变压力计)的压力计分辨率的小倍数以下),则可以认为预测试结果与压力计无关,因此可以认为所述预测试结果更加可靠。因此,在这种情况下,更高的置信度可以放入预测试结果中。在其它情况下,可以确定的是由不同的压力计进行的压力测量之间的差异具有原因。如果可以确定所述原因,可以将更高的置信度放入预测试结果中。
[0459] 图46描述了使用增压技术确定置信度标记的另一种方法。所述方法包括使用例如公式1计算流动性的步骤3602。如果流动性在某一水平以下,可以期望的是检查(check for)增压作用。可选地,可以不管流动性的水平而执行测试。然后在步骤3604处沿预测试选择点的子集。在步骤3606处确定子集中的每一个点处的球面导数。此外,在步骤3608处确定子集的几何平均数。可以使用几何平均数确定恢复压力是否被增压。例如,然后在步骤3610处可以比较几何平均数与预定的界限。
[0460] 这种方法用于评价稳定的井底压力是否是地层压力的良好表示。具有多种原因使得井底压力与地层压力不同,例如,泥浆滤液连续通过不良泥饼而渗漏到地层内的作用,这被公知为增压作用。这种现象大多数通常与“低”地层流动性相关联,其中低的定义取决于钻井实践、泥浆类型及其特征和执行预测试的条件,例如泥浆在测试期间是否正在被循环,并且如果泥浆在测试期间正在循环,则泥浆在多大流量下循环。被认为要增压的测量值可以认为具有比没有考虑增压的测量值具有更低的质量。
[0461] 在一个示例中,确定恢复压力是否被增压。优选地,首先使用任意预测试循环,例如通过使用相对于图7和/或公式1所述的技术计算流动性,可以如下比较这种流动性与阈值:
[0462]
[0463] 其中MS是通常小于1mD/cP到10mD/cP的流动性的界限,增压预期在所述界限以上。
[0464] 例如,如相对于图26所述沿压力恢复期选择数据点。优选地,在压力恢复开始之后在至少两个时间常数β下选择N个数据点(参见公式35和公式37),如以数学的方式表示如下:
[0465] 对于k=1,....,N来说,tk≥Mtβ其中Mt≥2
[0466] 可以例如利用相对于图32、图33A-33B、图40A-40C以及公式31所述的数据扩展和平滑近似法计算“普通的”压力导数。优选地,首先通常通过确保满足公式49来执行压力恢复的整体趋势分析以确保压力曲线表现为如所期望的。
[0467] 此外,计算如以下在公式41中进一步限定的在这些点处的压力的球面时间导数dp/dfS。对于持续时间tp的单个压力下降期通过以下公式给出球面导数:
[0468] dp/dfS(t)=2t3/2dp/dt(t){(1-τp)3/2/(1-(1-τp)3/2)} (41)[0469] 其中,τp≡tp/t,而dp/dt是如上所述确定的“普通的”压力导数。
[0470] 然后可以计算累加的球面导数的几何平均数。根据所述结果赋予置信度标记。在一些情况下,可以比较几何平均数与阈值。具体地,可以根据对公式系统42的不满足或满足设定置信度标记。如果满足以下公式,则稳定的井底压力(如由最终恢复压力Pb1或Pb2表示)被增压:
[0471]
[0472] 其中DS是对球面导数的几何平均数的界限,通常为100磅/平方英寸。
[0473] 图47更详细地示出了用于分析如在图34中的步骤2310中所述的置信度标记的方法4800。所述方法包括在步骤4810中确定置信度标记。所述方法可以利用已经在图34中的步骤2304中确定的置信度标记。所述方法进一步包括在步骤4820处比较多个置信度标记与阈值;在步骤4830处由所述比较确定多个指示值,并且在步骤4840处从所述多个指示值识别至少一个井下事件。
[0474] 在步骤4820处,比较先前确定的置信度标记与诸如噪点水平或特征值的阈值。这些阈值可以由测试条件的先验知识确定,例如,由钻井泥浆组分、由在同一或其它油藏中的先前的测试数据库等确定所述阈值。可选地,可以通过模拟诸如地层流动性的限制确定这些阈值,其中在所述地层流动性的限制处,对于具体的测试条件可以预计的增压作用。还可以通过诸如压力计噪点的实验确定这些阈值。最终,可以由诸如具体事件处的噪点测量值或压力水平的预测试数据计算阈值。已经相对于压力比较技术(例如,参见公式26及其改进公式)、参数比较技术(例如,参见公式27、28或29)、参数预测技术(例如,参见公式30)、趋势分析技术(例如,参见公式49)、离散分析技术(例如,参见公式33)、增压技术(例如,参见公式42和公式48)以及本公开所述的其它技术论述了各种比较。可选地,在此步骤中可以比较本领域公知的其它置信度标记与阈值。
[0475] 在步骤4830处,根据比较确定指示值。在步骤4830的一个示例中,指示值可以是基于比较的有效性的布尔数字。一般地,可以使用至少一个置信度标记的任意布尔评价函数。在另一个示例中,使用本领域公知的模糊逻辑原理推导指示值。所述值则可以是0与1之间的数,且0例如指示置信度标记很好地在阈值以下,而1指示置信度标记很好地在阈值以上,而0与1之间的值指示置信度标记稍微接近阈值。
[0476] 在步骤4840处,识别井下条件。井下条件可以是从操作者观点来看所关心的任何信息。在一个示例中,井下条件可与钻井操作有关。所述井下条件包括诸如“井是欠平衡的”和“井是超平衡的”的条件。在另一个示例中,井下条件可与工具状态有关。在这种情况下,所述井下条件包括诸如“流动管线被间歇地堵塞”、“探头没有到达井壁”等的条件。在又一个示例中,井下条件可与地层和井眼有关。所述井下条件包括诸如“地层是不可渗透的”、“泥饼正在渗漏”、“井底压力被增压”、“在流动管线中检测到气体”等的条件。在又一个示例中,井下条件可与预测试循环有关。所述井下条件包括诸如“在测试结束之前已经终止调查阶段”、“根据体积判据已经终止调查阶段”、“由调查阶段计算以设计测量阶段的测试参数在范围之外”等的条件。在又一个示例中,井下条件可与压力测量有关。在这种情况下,所述井下条件包括诸如“下降压力足以测量井底流压”、“测量阶段压力恢复期达到稳定”、“测量值是有噪点的”等的条件。以下说明这些条件中的一些及其它井下条件。
[0477] 将要认识的是各种置信度标记及其相关联的指示值如果单独考虑可以被含糊地解释为多于一个的井下条件。例如,与压力恢复事件结束相对应的压力水平,即Pb1或Pb2,其如公式26所要求的不小于静水压或井筒压力,且可以尤其被解释为在预测试循环期间丧失密封、被解释为探头没有充分伸出以到达井壁、或者被解释为以欠平衡的方式钻井。类似地,与压力恢复事件的结束相对应且几乎等于与压力下降的相应的结束相对应的压力水平(即,Pd1或Pd2)的压力水平可以尤其被解释为干测试(不渗透地层)或不足的压力下降。为了在更多必然性的水平下识别井下条件,有利的是分析多个指示值。
[0478] 更具体地,每一个井下事件可以与具有在步骤4830处确定的作为输入的指示值的真值表相关联。在一些情况下,仅一个指示值可以足以识别井下条件。在其它情况下,可能需要多个指示值以识别井下条件。参照以上示例,压力恢复期的整体趋势分析可以在以欠平衡的方式钻井的条件与丧失或不存在密封的条件下进行区分。通常,如果检测到压力恢复期并且所述压力恢复期具有适当级数的压力水平和/或压力曲线斜率,则可以确定丧失密封的井下条件。由此,可以使用真值表识别“井是欠平衡的”井下事件,所述真值表具有与压力恢复事件结束时的压力水平与井筒压力事件处的压力水平的比较相关联的第一指示值、和与应用到压力恢复期的整体趋势分析相关联的第二指示值。
[0479] 本领域的技术人员将认识到还可以在步骤4840中使用与不同井下条件相关联的其它真值表。例如,可以使用模型相关技术和压力比较技术识别渗漏泥饼。此外,可以使用参数预测技术以及识别缓慢压力恢复期的整体趋势分析技术识别流动管线中的气体。还将认识的是真值表的使用仅仅是用于执行步骤4840的示例性技术,并且可以代替地使用其它技术。具体地,可以使用模糊逻辑。
[0480] 图48显示了用于显示所识别的井下条件的方法4900。所述方法包括优选地在步骤4910处选择互斥的多个井下条件,以及在步骤4920处使不同的整数与多个井下条件中的每一个相关联。可以在将井下测试工具定位在井眼之前执行此步骤。还可以将这些步骤重复多次以包括可以独立于先前的集合所发生的井下条件的集合。所述方法4900进一步包括在步骤4930处执行井下测量,在步骤4940处识别多个井下条件中的一个,在步骤4950处发送以显示与所述条件相关联的整数,在步骤4960处接收发送的整数,以及在步骤4970处显示例如表示与接收到的整数相关联的井下条件的指示的句子。
[0481] 在一个示例性实施例中,在步骤4910处选择的井下条件包括:压力测试是正常的,井是超平衡的;压力测试是正常的,井是欠平衡的;压力试验是正常的,超平衡是不确定的;压力试验是干测试;在压力试验期间还没有实现密封(探头在井眼中被冲掉);在压力试验期间丧失密封;和/或压力试验是不可辨认的。这些井下条件优选地是互斥的。
[0482] 作为步骤4920的示例性实施例,以上所列的第一条件可以与整数0相关联,第二条件与整数1相关联等。因此,当以上所列的井下条件中的一个被识别为真时,可以对所述井下条件中的一个通过0与6之间的整数进行编码,所述整数可以被转换成适于3位或更多位的二进制字。
[0483] 在步骤4930处,可以执行任意测量。具体地,可以使用利用如上所述的测试工具和方法的压力测量。测量的类型不限于本公开。在步骤4940处,多个井下条件中的一个可以被识别为真。为了进行此,可以使用诸如方法4800的方法。可以代替地使用其它方法。
[0484] 在步骤4950处,发送与在步骤4940处识别的条件相关联的整数。例如,当识别到条件“压力试验正常,井是超平衡的”时,发送数字0,当识别到条件“压力试验正常,井是欠平衡的”时,发送数字1。将要认识的是如果多个条件是互斥的,则仅一个条件被识别为真,因此仅发送一个整数。因此,当遥测带宽受限时,这种编码方法是有利的。事实上,可以通过分析在井下采集并以压缩形式发送的大量数据来检测对地面操作者重要的信息。地面显示可以是能够接收数据并显示所述数据的例如在屏幕或印刷测井图上的任意系统。
[0485] 然后接收发送的整数。所述整数被解码并且指示与所述整数相关联的井下条件的句子被显示。参照以上示例,如果接收到整数0,句子“压力试验正常,井是超平衡的”可以显示给地面操作者。可以代替地显示类似意思的其它句子。
[0486] 如先前所述,还可以在步骤4910处选择第二组井下条件(例如,不可由第一组条件预测的)。与预测试的压力恢复期有关的第二组井下条件可以包括:沿压力恢复曲线降低压力方差和减小正斜率;沿压力恢复曲线降低压力方差并减小负斜率;沿压力恢复曲线忽略压力方差并且忽略斜率;沿压力恢复曲线几乎使压力方差和正斜率恒定;沿压力恢复曲线几乎使压力方差和负斜率恒定;沿压力恢复曲线增加压力方差和正斜率;和/或不可辨别形状(不是先前中的任一个)。在步骤4920处,与此第二组条件相关联的整数可以在0与6之间。在步骤4940处,可以例如使用如这里所述的整体趋势分析技术和整体离散分析技术识别与此第二组相关联的条件中的一个。
[0487] 在步骤4950处,与第二组被识别的条件相关联的整数可以在第二3(或以上)位二进制字中被编码。在一些情况下,可以有利的是链接与第一组条件相对应的二进制字和与第二组条件相对应的二进制字。在步骤4960处,地面解码器可以解除两个接收到的字的链接。在一些情况下,可以在步骤4970处显示与每一个字相对应的一个句子。在其它情况下,可以显示或多或少的句子。要注意的是还可以在不同组中任意重新组合所述条件。还要注意的是其它多组井下条件可以添加到上述示例。
[0488] 这里已经参照所述的具有压力和时间值的地层预测试数据说明了结构。然而,应该认识的是本公开的原理不限于具体的数据、数据源、数据通过其被发送的介质。此外,数据不必是压力数据。例如,数据可以包括来自压力传感器中的一个的温度、来自应变压力计的电压的温度。虽然温度和电压本身不是压力数据,但是所述温度和电压与压力测量值有关,并因此也可以应用到所述数据。
[0489] 此外,本公开不局限于具体的步骤、步骤的顺序、或以上所述的示例的结构。因此,可以添加或删去另外和/或可选的步骤。这里提供的一个或多个方法可以单独使用或组合使用。例如,可以期望使用一个或多个置信度标记方法为一个或多个预测试生成整体置信度标记。然后,可以使用置信度标记的结果调节预测试操作。在一些情况下,第一预测试的置信度标记可以用于有助于设计一个或多个随后的预测试。还可以使用其它的预测试设计判据。
[0490] 还应该认识的是本公开的原理不局限于具体的手动、可视或自动实施。此外,如果期望自动实施,这种实施可以由井下工具硬件、井口钻机硬件、客户办公室硬件、或所述井下工具硬件、井口钻机硬件、客户办公室硬件的任意组合来支持。
[0491] 应该认识的是通过使用本公开的原理,可以实时或近似实时地压缩和传输数据。例如,在数据包括地层预测试数据的情况下,可以在完成预测试之前,例如,在适当数量的事件数据点(例如一个或多个事件数据点)和另外的数据点(例如,在事件数据点之前或之后的一系列数据点)被捕获之后,执行压缩和传输。所述方法可以涉及由先前执行、和/或当前被测试的预测试获得的数据。
[0492] 以上概括了多个实施例的特征,使得本领域的技术人员可以更好地理解本公开的多个方面。本领域的技术人员应该认识到他们可以容易地使用本公开作为用于设计或修改用于实施相同目的和/或实现这里引入的实施例的相同优点的其它过程和结构的基础。本领域的技术人员还应该认识到这种等价结构不会背离本公开的保护范围,并且本领域的技术人员可以在不背离本公开的保护范围的情况下可以做各种变化、替换和改变。