CN111971452A

CN111971452A - 使用互连储层区域图的渗透率预测

Info

Publication number: CN111971452A
Application number: CN201980025823.0A
Authority: CN
Inventors: 巴巴托普·嘉得
Original assignee: Saudi Arabian Oil Co
Current assignee: Saudi Arabian Oil Co
Priority date: 2018-02-21
Filing date: 2019-02-19
Publication date: 2020-11-20
Also published as: US11085291B2; WO2019164790A1; US20190257196A1; SA520420021B1; EP3755873A1

Abstract

用于执行包括以下操作的方法、系统及计算机可读介质：获得用于标识储层的多个井中的每个井的时间推移平均压力的数据；确定所述数据的多种均质趋势；基于所述均质趋势将所述数据聚类成多个分组，每个分组与对应的趋势相对应；以及生成包括多个互连储层区域(CRR)的互连储层区域图，每个CRR与所述多个分组中的一个分组相对应，每个CRR包括所述多个井的与该CRR的对应的分组相对应的子集，其中，每个CRR与均质压力响应相关联。

Description

使用互连储层区域图的渗透率预测

优先权要求

本申请要求于2018年2月21日递交的美国专利申请No.62/633,164的优先权，其全部内容由此通过引用并入。

技术领域

本公开涉及储层之中的渗透率预测，更具体地，生成用于促进这种预测的互连储层区域图。

背景技术

测地模型(geo-model)是储层的三维(3-D)原型。这种模型中定义的性质中的一种性质是储层渗透率。实际上，储层渗透率通常是在井穿透(well penetration)处或沿着井穿透竖直测量的。然而，取决于储层成熟度的状态，储层之中的井间隔可以在500m至1km(约0.31英里至0.62英里)的数量级。作为这种井间隔的结果，可能难以基于井穿透自身来测量井穿透之间的储层渗透率(也被称为“井间储层渗透率”)。为了克服这种缺陷，当前的测地建模方法使用在井穿透处获得的(例如，岩心或测井)数据、以及统计、地球物理或沉积分析技术来确定井间储层渗透率。

发明内容

本公开对使用时间推移平均关井压力数据来确定互连储层区域(CRR)进行了讨论。更具体地，本公开使用井数据和时间推移平均关井压力数据来确定CRR。诸如测井和岩心数据之类的井数据，提供非均质性的竖直分辨率，而时间推移关井压力数据指示了井间的性质的横向连续性的区域。CRR图可以用作用于创建不同的渗透率函数的基础，所述不同的渗透率函数用于已进行测井但未进行岩心采样(cored)的井的1维(1-D)渗透率预测。另外，CRR图可以用于3维井间渗透率预测或分配。

本说明书描述的主题的创新性方面可以体现在方法中，所述方法包括以下动作：获得对储层的多个井中的每个井的时间推移平均压力进行标识的数据；确定所述数据中的多种均质趋势；基于所述均质趋势将所述数据聚类成多个分组，其中，每个分组与对应的趋势相对应；以及生成包括多个互连储层区域(CRR)的互连储层区域图，每个CRR与所述多个分组中的一个分组相对应，其中，每个CRR包括所述多个井的与该CRR的对应的分组相对应的子集，并且其中，每个CRR与均质压力响应相关联。

这些方面的其他实施例包括：对应的被配置为执行所述方法的动作的系统、装置、以及在计算机存储设备上进行编码的计算机程序。

这些和其他实施例均可以可选地包括以下特征中的一个或多个。例如，针对所述多个CRR中的每个CRR，基于来自该CRR的岩心和测井数据，创建1维渗透率函数。此外，针对所述多个CRR中的每个CRR，在该CRR内分配渗透率包含仅使用该CRR内的井。确定所述数据的多种均质趋势还包含：确定所述多个井中的与平均压力的均质性相关联的两个或更多个井。生成所述互连储层区域图包含：创建多个沉积多边形，以将与每个分组相对应的井和与其他剩余的分组相对应的井划分开。所述每个井的时间推移平均压力是时间推移平均关井压力。基于所述互连储层区域图确定所述储层的方差图的范围。每个CRR的每个井液力连通到相应CRR中的每个其他的井。每个CRR没有大尺度渗透率阻挡层。

可以实现本说明书中描述的主题的特定实施方式，以通过提供已从井穿透观察到的关于储层性质的横向连续性的信息来改进储层测地模型。改进的储层测地模型可以包括，作为比当前实际获得的储层的表示更准确的表示的1维和3维渗透率分布。

在本文的附图和说明书中阐述了本公开中所描述的主题的一个或多个实施方式的细节。通过说明书、附图和权利要求书，所述主题的其他特征、方面和优点将变得明显。

附图说明

图1示出根据一些实施方式的储层的示意图。

图2示出根据一些实施方式的储层渗透率图的三种实现。

图3示出根据一些实施方式的用于生成互连储层区域图的系统。

图4示出根据一些实施方式的显示储层中的每个井的时间推移平均压力的图。

图5A-图5D示出根据一些实施方式的与储层中的每个井的时间推移平均压力的数据的均质趋势相对应的图。

图6示出根据一些实施方式的互连储层区域图的实现。

图7示出根据一些实施方式的互连储层区域图的另一实现。

图8示出根据一些实施方式的1维渗透率曲线图。

图9示出根据一些实施方式的储层的3维渗透率分布图。

图10A示出根据一些实施方式的储层的两种测地模型。

图10B和图10C各自示出根据一些实施方式的使用图10A的测地模型中的一种的仿真的结果。

图11示出根据一些实施方式的用于确定储层之中的3维渗透率分布的方法的流程图。

图12示出根据一些实施方式的用于实现本文描述的技术的示例计算环境。

具体实施方式

以下详细描述对用于改进储层中的渗透率预测的系统及方法进行了描述，并且被呈现为使得本领域技术人员能够做出和使用在一个或多个具体实施方式的上下文中所公开的主题。可以做出所公开的实施方式的各种修改、改变和置换并且对本领域普通技术人员将显而易见，并且在不背离本公开的范围的情况下，所定义的一般原理可以应用于其他实施方式和应用。因此，本公开不旨在被限制在所描述或所示出的实施方式，而是要符合与所公开的原理和特征一致的最宽的范围。

根据上述讨论，当前的测地建模技术使用可在井穿透处获得的数据、以及统计、地球物理或沉积分析技术来预测储层中的井间渗透率。虽然最终的预测在井穿透附近可能是准确的，但是在井间储层范围中的预测的准确性很差，这在某种程度上是由于所述预测基于统计方法。此外，在许多实例中，统计方法可以针对相同区域生成若干不同的等概率的井间渗透率预测。如果所述预测是等概率的，那么可能难以确定最准确的预测。为了说明这些缺陷，考虑使用当前的测地建模技术在图1的储层中做出的预测。

图1示出根据一些实施方式的包括井104的储层100。当前的测地建模技术使用来自井104的井数据，以及统计、地球物理或沉积分析来确定储层100中的渗透率。在井104处或附近的区域中，当前的技术使用井数据来确定井渗透率。在井间区域中，当前的技术使用统计、地球物理或沉积分析来确定井间渗透率。

图2示出使用当前的测地建模技术做出的储层100中的渗透率预测的三种实现200、202、204。每种实现表示使用来自井104的相同的井数据做出的渗透率预测。在这些实现中，相邻的绿色网格单元指示储层100中的渗透率的预测区域。如图2所示，在井穿透点处的渗透率预测是相同的，指示井104附近的渗透率预测是准确的。然而，在井间区域中，渗透率预测是不同的。如果所述预测是等概率的，则难以确定哪个预测是最准确的。因此，需要确定储层渗透率的(特别是用于储层的井间区域的)更准确的方法。

本文公开了一种方法，包括了图3所示的系统300，以确定储层之中的渗透率分布。在高层处，系统300使用来自井(例如，图1的储层100的井104)的时间推移压力数据的比较，来确定渗透率分布。具体地，系统300对时间推移压力数据进行比较，以生成具有类似的时间推移平均压力趋势的井分组的图。然后将已生成的图用于确定储层100中的渗透率分布。这种测定法的根本基础在于，在没有问题的储层中，相邻的井之间的时间推移压力趋势上的差别是储层渗透率上的不连续性的指示，而时间推移压力趋势上的相似性是储层连续性和连通性的指示。

如图3所示，系统300包括：计算设备302，其可以通过一个或多个网络与一个或多个其他计算设备(未示出)和数据储存器(未示出)通信。图3显示了被标记为(A)至(D)的多个步骤，其示出了数据流，并且可以以所示出的顺序或以另一顺序执行。

在步骤A中，计算设备302获得数据320，数据320标识多个井104中的每个井的时间推移平均压力。在一些实施例中，计算设备302从数据储存器(未示出)接收数据320，或当数据320被采集时而(例如，从多个井104)实时接收数据320。计算设备302可以通过一个或多个网络或其他通信方式接收数据320。

参考图4，图400显示多个井104中的每个井的时间推移平均压力。具体地，图400显示每个井的时间推移平均压力的非均质云404。在一些实施例中，每个井的时间推移平均压力是井的时间推移平均供油面积(drainage area)关井压力。在一些实施例中，数据320和图400指示(在大尺度上的)储层非均质性；然而，在所述非均质性范围内，可能存在显示均质性的井104的分组，本文中进行进一步描述。

在步骤B中，计算设备302确定数据320的多个均质趋势。具体地，计算设备320确定多个井104中的两个或更多个井与数据320的平均供油面积压力的均质性相关联。例如，图5A-图5D示出与数据320的均质趋势504a、504b、504c、504c(被统称为均质趋势504)分别相对应的图502a、502b、502c、502d(被统称为图502)。即，每个图502的每个均质趋势504与对应的平均压力数据的均质趋势相关联。

在步骤C中，计算设备302基于均质趋势504将数据320聚类成多个分组。具体地，计算设备302基于数据320的均质趋势504将数据320(包括非均质的井数据)分组成分组。即，计算设备302基于均质趋势504确定多个井104的与平均供油面积压力的相同的均质性相关联的子集。每个分组与对应的均质趋势504相对应。在一些实施例中，每个分组标识井的针对特定均质趋势504的子集；即，井的与特定均质趋势504的供油面积压力相关联的子集。

在步骤D中，计算设备302生成互连储层区域(CRR)图540，其包括多个互连储层区域(CRR)，每个CRR与多个分组中的一个分组相对应。

图6示出CRR图340(被示出为包括多个CRR 602的CRR图600)的示例。图6还包括与CRR 602并列的图502，指示均质趋势504和CRR 602之间的关系，本文中进行进一步描述。每个CRR 602与多个分组中的一个分组相对应，并且包括井104的与该分组相对应的子集。例如，CRR 602a与均质趋势504a相对应，并且包括与均质趋势504a相关联的井104的第一子集——即，与特定的均质趋势504a的供油面积压力相关联的井104。类似地，CRR 602b与均质趋势504b相对应，并且包括与均质趋势504b相关联的井104的第二子集——即，与特定的均质趋势504b的供油面积压力相关联的井104；CRR 602c与均质趋势504c相对应，并且包括与均质趋势504c相关联的井104的第三子集——即，与特定的均质趋势504c的供油面积压力相关联的井104；并且CRR 602d与均质趋势504d相对应，并且包括与均质趋势504d相关联的井104的第四子集——即，与特定的均质趋势504d的供油面积压力相关联的井104。在一些实施例中，井104的第一子集、第二子集、第三子集、第四子集中的每一个子集彼此不同。在一些实施例中，CCR 602中的每一个与均质压力响应相关联。

在一些实施例中，生成CRR图340和CRR图600包括：创建多个沉积多边形，以将与井104的每个分组相对应的井和与其他剩余的分组相对应的井划分开。即，通过对井104进行如此划分，CRR图340和CRR图600指示每个分组的每个井104液力连通到对应的CRR 602中的每个其他的井104。在一些实施例中，每个CRR 602没有大尺度渗透率阻挡层。因此，CRR图340和CRR图600可以用于确定井之间是否存在连续性、障碍或阻挡层。

在实施例中，CRR图可以用于确定储层中的渗透率分布。

在实施方式中，CRR图可以用于改进储层中的已进行测井但未进行岩心采样的井的1维渗透率的建模。作为背景技术，在新的储层中进行钻井的早期阶段，对大多数井穿透即进行测井也进行岩心采样，以获得用于储层建模和描述的足够的数据。然而，通常对随后钻取的井进行测井但不进行岩心采样。对于缺少岩心数据的这些井，可能难以仅基于测井记录准确地确定渗透率。

实际上，为了预测已进行测井但未进行岩心采样的井的渗透率，渗透率函数被创建用于预测渗透率。当前的技术通过将已进行测井和岩心采样的井用作训练集合来创建渗透率函数。具体地，训练数据集合用于建立岩心渗透率和测井响应之间的关系，然后将这种关系用于预测井的渗透率。将具有岩心和测井数据二者的至少一个井留在训练数据集合之外，使其可以用于评估在所谓的盲测试的过程中所得到的渗透率函数的预测能力。

然而，当前的这些技术使用来自储层的所有可用的岩心和测井数据来创建渗透率函数，而不区分储层中的各个非均质区域。结果，将储层质量好的和储层质量差的岩心和测井数据一起用于获得单个平均渗透率函数。这样的渗透率函数将导致低质量储层面积的乐观的预测和高质量储层面积的悲观的预测。

在实施例中，CRR图可以用作用于驱动1维渗透率函数建模的基础。在这个实施例中，不是使用储层中的所有已进行岩心采样和测井的井来创建储层中的所有已进行测井但未进行岩心采样的井的渗透率函数(这是当前的技术中的做法)，而是仅使用每个CRR中的数据集合为不同的CRR创建不同的渗透率函数。即，将指定CRR中的已进行岩心采样和测井的井用于创建该CRR中的已进行岩心采样但未进行测井的井的渗透率函数。当仅将每个CRR之中的数据用于创建该CRR的单独的渗透率函数时，结果将免受其他CRR的偏置，因此，更好地反映CRR的性质。

图7示出根据一些实施方式的储层的CRR图700。对CRR图700进行颜色编码，以示出储层中的不同的CRR。如图7所示，储层具有三个CRR 702a、702b、702c。图7也通过在CRR图700上添加的圆点示出了储层中的已进行测井和岩心采样的井。例如，CRR 702a包括井704a，CRR 702b包括井704b，并且CRR 702c包括井704c。在实施例中，使用每个CRR中的已进行测井和岩心采样的井，创建该CRR的渗透率函数。即，将已进行测井和岩心采样的井704a用于创建CRR 702a的渗透率函数，将已进行测井和岩心采样的井704b用于创建CRR 702b的渗透率函数，并且将已进行测井和岩心采样的井704c用于创建CRR 702c的渗透率函数。

图8示出根据一些实施方式的渗透率曲线图800。将渗透率曲线图800用作当前实际使用的渗透率函数和使用CRR图得到的渗透率函数的盲测试比较。为了比较不同的渗透率函数，图800描绘了通过CRR中的井的相同的分组的每个渗透率函数预测的渗透率。在图8中，圆点802指示盲测试井的岩心数据，线806是通过当前实际使用的渗透率函数预测的渗透率，而线804是通过CRR驱动的渗透率函数预测的渗透率。如图8所示，线804比线806更靠近圆点802。因此，CRR驱动的渗透率函数比当前实际使用的渗透率函数更准确。

在另一实施方式中，CRR图可以用作用于驱动储层中的3维渗透率分布的基础。在这种实施方式中，使用来自每个CRR中的井的井数据来确定该CRR的单独的3维渗透率分布。更具体地，在指定的CRR区域中，仅使用来自该CRR的1维井渗透率记录来确定3维渗透率分布。例如，在图6的示例中，对于每个CRR 602，仅使用CRR 602之中的井及其1维渗透率记录来确定CRR 602之中的渗透率分布。作为使用每个CRR中的井分别确定该CRR的渗透率分布的结果，来自相邻的低渗透率CRR 602的低渗透率数据没有统计分布到附近的高渗透率CRR602中。要注意，在包括已进行测井但未进行岩心采样的井的储层中，可以首先将CRR图用于确定那些井的1维渗透率记录，以便可以将1维渗透率记录用于确定3维渗透率分布。

附加地和/或备选地，当确定CRR的3维渗透率分布时，选择3维渗透率分布方差图以确保在该CRR之中没有分布低渗透率阻挡层。在实施方式中，选择方差图使该CRR之中没有出现低渗透率阻挡层。即，通过使用CRR图，并且结合储层中井的岩心和测井数据，选择方差图以避免在任意指定的CRR之中出现任何不想要的非均质性。在示例中，将方差图范围计算为，每个CRR之中的最大井间距离的一半。在这个示例中，指定CRR之中的渗透率完全相关，因此是一致的。

图9示出根据一些实施方式的储层中3维渗透率分布的示例图902。在实施方式中，储层中3维渗透率分布是使用储层的CRR图确定的。

图10A示出根据一些实施方式的储层的两种测地模型1002、1004。测地模型1002、1004各自描绘了储层的渗透率分布。两个测地模型是使用井数据的相同的集合生成的。然而，测地模型1002是使用当前实际使用的方差图确定的方法生成的，而测地模型1004是使用本文描述的CRR驱动的方法生成的。在这个示例储层中，当前实际使用的方差图确定的方法确定了小的方差图范围，因此导致严重的渗透率非均质性。另一方面，CRR驱动的方法确定了更大的方差图范围，因此导致储层渗透率的更好的连续性。此外，在这个示例中，从这个储层中的井观察到的时间推移压力数据指示：井间连通充分，表明渗透率均匀。

图10B示出根据一些实施方式的利用测地模型1002执行的储层仿真的结果。如图10B所示，所述结果示出了储层的井之间的连通很差，其不是由观察到的压力数据所暗示的情况。

图10C示出根据一些实施方式的利用测地模型1004执行的储层仿真的结果。如图10C所示，所述结果示出了井之间的连通极好，其与观察到的压力数据一致。

图11示出用于说明储层之中3维渗透率分布的方法1100的流程图。为了清楚地呈现，下文的描述在图2-图10的上下文中描述方法1100。例如，如图所示，在计算系统300上或在计算系统300处执行方法1100的特定步骤。然而，在适当时，方法1100可以例如通过任何其他合适的系统、环境、软件和硬件，或者系统、环境、软件和硬件的组合来执行。方法1100的操作可以包括一个或多个可选步骤，只包括正被执行的一个或多个步骤。

在步骤1102处，通过计算设备302获得，对储层102的每个井104的时间推移平均供油面积压力进行标识的数据320。在步骤1104处，通过计算设备302确定数据320的均质趋势504。在步骤1106处，基于均质趋势504将数据320聚类成分组。在一些实施例中，每个分组与相应的均质趋势504相对应。在步骤1108处，通过计算设备302生成互连储层区域图600。在一些实施例中，图600包括互连储层区域(CRR)602，每个CRR与多个趋势分组中的一个分组相对应。在一些实施例中，每个CRR 602包括井104的与CRR 602的相应分组相对应的子集。在一些实施例中，每个CRR 602与均质压力响应相关联。

图12示出结合此处描述的技术使用的通用计算机设备1200和通用移动计算机设备1250的示例。计算设备1200旨在表示各种形式的数字计算机，例如膝上型计算机、台式计算机、工作站、个人数字助理、服务器、刀片服务器、大型机和其他适当的计算机。计算设备1250旨在表示各种形式的移动设备，例如个人数字助理、蜂窝电话、智能电话以及其他类似的计算设备。在此处示出了组件，其连接及关系、以及其功能意在仅是示例性的，并且不意在限制本文档中描述和要求保护的本发明的实施方式。

计算设备1200包括处理器1202、存储器1204、存储设备1206、连接到存储器1204和高速扩展接口1210的高速接口1208、以及连接到低速总线1214和存储设备1206的低速接口1212。组件1202、1204、1206、1208、1210、以及1212中的每一个组件使用各种总线互连，并且被安装在公用的母板上或在适当时采用其他方式。处理器1202处理用于在计算设备1200之中执行的指令，包括被存储在存储器1204中或存储设备1206上的用于在外部输入/输出设备(例如，耦合到高速接口1208的显示器1216)上显示GUI的图形信息的指令。在其他实施方式中，在适当时，与多个存储器和多种类型的存储器一起使用多个处理器、多个总线、或二者。而且，多个计算设备1200与提供所需要的操作的部分的(例如，作为服务器组、刀片服务器组、或多处理器系统的)每个设备进行连接。

存储器1204存储计算设备1200内的信息。在一个实施方式中，存储器1204是一个或多个易失性存储器单元。在另一个实施方式中，存储器1204是一个或多个非易失性存储器单元。存储器1204也可以是另一形式的计算机可读介质，例如，磁盘或光盘。

存储设备1206能够为计算设备1200提供大容量存储。在一个实施方式中，存储设备1206包含计算机可读介质，例如软盘设备、硬盘设备、光盘设备、或磁带设备、闪存或其他类似的固态存储器设备、或设备(包括存储区域网络或其他配置中的设备)的阵列。计算机程序产品有形地体现在信息载体中。计算机程序产品也可以包含指令，当所述指令执行时，执行例如本文描述的那些的一种或多种方法。信息载体是计算机或机器可读介质，例如存储器1204、存储设备1206或处理器1202上的存储器。

高速控制器1208为计算设备1200管理带宽密集型操作，而低速控制器1212管理较低的带宽密集型操作。这样的功能分配仅是示例性的。在一个实施方式中，高速控制器1208耦合到存储器1204、显示器1216(例如，通过图形处理器或加速器)、以及容纳各种扩展卡(未示出)的高速扩展接口1210。在实施方式中，低速控制器1212耦合到存储设备1206和低速扩展端口1214。包括各种通信端口(例如，USB、蓝牙、以太网、无线以太网)的低速扩展端口耦合到一个或多个输入/输出设备，例如，键盘、指示设备、扫描仪、或例如通过网络适配器的联网设备(例如，交换机或路由器)。

计算设备1200以多种不同的形式实现，如图所示。例如，其被实现为标准服务器1220、或在这样的服务器组中实现多次。其也可以被实现为机架服务器系统1224的一部分。另外，其可以被实现在诸如膝上型计算机1222之类的个人计算机中。备选地，来自计算设备1200的组件与移动设备(例如，设备1250)中的其他组件(未示出)进行组合。这样的设备中的每一个设备包含计算设备1200、1250中的一个或多个，并且整个系统由相互通信的多个计算设备1200、1250组成。

计算设备1250包括处理器1252、存储器1264、输入/输出设备(例如，显示器1254)、通信接口1266和收发机1268、以及其他组件。设备1250也可以配备有用于提供附加存储的存储设备，例如微型驱动器或其他设备。组件1250、1252、1264、1254、1266、以及1268中的每一个组件使用各种总线互连，并且所述组件中的若干组件被安装在公用的母板上或在适当时采用其他方式。

处理器1252执行计算设备1250之中的指令(包括被存储在存储器1264中的指令)。处理器被实现为一组芯片，其中包括多个独立的模拟和数字处理器。例如，处理器提供设备1250的其他组件的协调，例如，用户接口的控制、通过设备1250运行应用、通过设备1250的无线通信。

处理器1252通过控制接口1258和耦合到显示器1254的显示器接口1256与用户通信。例如，显示器1254是TFT LCD(薄膜晶体管液晶显示器)或OLED(有机发光二极管)显示器、或其他适当的显示器技术。显示器接口1256包括用于驱动显示器1254以向用户呈现图形和其他信息的适当的电路。控制接口1258接收来自用户的命令，并且对它们进行装换，以提交给处理器1252。另外，外部接口1262与处理器1252进行通信，以便实现设备1250与其他设备的近区通信。例如，外部接口1262在一些实施方式中提供有线通信，或在其他实施方式中提供无线通信，并且可以使用多个接口。

存储器1264存储计算设备1250内的信息。存储器1264被实现为以下一种或多种：一个或多个计算机可读介质、一个或多个易失性存储器单元、或一个或多个非易失性存储器单元。也提供扩展存储器1274，并且通过扩展接口1272连接到设备1250，扩展接口1272可以包括例如SIMM(单列直插式存储器模块)卡接口。这样的扩展存储器1274为设备1250提供额外的存储空间，或者也为设备1250存储应用程序或其他信息。具体地，扩展存储器1274包括用于执行或补充本文描述的过程的指令，并且也包括安全信息。因此，例如，扩展存储器1274作为设备1250的安全模块提供，并且使用允许安全使用设备1250的指令进行编程。另外，经由SIMM卡提供安全应用程序、以及附加信息，例如，以无法被盗用的方式将标识信息置于SIMM卡上。

如本文所讨论的，存储器包括例如闪存、NVRAM存储器或二者。在一个实施方式中，计算机程序产品有形地体现在信息载体中。计算机程序产品包含指令，当所述指令执行时，执行例如本文描述的那些的一种或多种方法。信息载体是计算机或机器可读介质，例如，存储器1264、扩展存储器1274、处理器1252上的存储器、或例如通过收发机1268或外部接口1262接收到的传播信号。

设备1250通过必要时包括数字信号处理电路的通信接口1266进行无线通信。通信接口1266提供诸如GSM语音呼叫、SMS、EMS、或MMS消息收发、CDMA、TDMA、PDC、WCDMA、CDMA2000、或GPRS等之类的各种模式或协议下的通信。例如，这样的通信通过射频收发机1268进行。另外，例如，使用蓝牙、WiFi、或其他此类收发机(未示出)进行短距离通信。另外，GPS(全球定位系统)接收机模块1250向设备1250提供附加的导航和与位置有关的无线数据，在适当时，它们被设备1250上运行的应用程序使用。

设备1250也使用音频编解码器1260来进行可听通信，音频编解码器1260接收来自用户的口语信息，并且将其转换成可用的数字信息。音频编解码器1260同样诸如通过扬声器例如在设备1250的手持设备中生成用户可听见的声音。这样的声音包括来自语音电话呼叫的声音，包括经记录的声音(例如，语音消息、音乐文件等)，并且还包括由在设备1250上工作的应用生成的声音。

如图所示，计算设备1250以多种不同的形式实现。例如，其被实现为蜂窝电话1280。其也被实现为智能电话1282、个人数字助理、或其他类似的移动设备的一部分。

在数字电子电路、集成电路、专门设计的ASIC(专用集成电路)、计算机硬件、固件、软件和/或其组合中实现此处描述的系统和技术的各种实施方式。这些不同的实施方式包括采用在可编程系统上可执行的和可解释的一个或多个计算机程序的实施方式，所述可编程系统包括至少一个可编程处理器，专用的或通用的所述可编程处理器被耦合以从储存器系统、至少一个输入设备、以及至少一个输出设备接收数据和指令，以及向其发送数据和指令。

这些计算机程序(也被称为程序、软件、软件应用程序或代码)包括用于可编程处理器的机器指令，并且以高级过程和面向对象的编程语言和汇编/机器语言来实现。如本文所使用的，术语“机器可读介质”、“计算机可读介质”是指所使用的用于向可编程处理器提供机器指令和数据的任何计算机程序产品、装置和/或设备(例如，磁盘、光盘、存储器、可编程逻辑器件(PLD))，包括接收机器指令作为机器可读信号的机器可读介质。术语“机器可读信号”是指用于向可编程处理器提供机器指令和数据的任何信号。

为了提供与用户的交互，可以在具有用于向用户显示信息的显示设备(例如，CRT(阴极射线管)或LCD(液晶显示器))、以及用户可以通过其向计算机提供输入的键盘和指向设备(例如，鼠标或轨迹球)的计算机上实现所述系统和技术。也将其他类型的设备用于提供与用户的交互；例如，提供给用户的反馈是任意形式的传感反馈(例如，视觉反馈、听觉反馈或触觉反馈)；并且以任意形式(包括声音、语音或触觉输入)来接收来自用户的输入。

本文描述的系统和技术被实现在计算系统中，该计算系统包括后端组件(例如，数据服务器)、或包括中间件组件(例如，应用服务器)、或者包括前端组件(例如，具有用户通过其与本文描述的系统和技术的实施方式进行交互的图形用户界面或者web浏览器的客户端计算机)、或者此类后端组件、中间件组件或前端组件的任意组合。系统的组件通过数字数据通信的任何形式或介质(例如，通信网络)来互连。通信网络的示例包括局域网(“LAN”)、广域网(“WAN”)和互联网。

计算系统可以包括客户端和服务器。客户端和服务器一般相互远离并且通常通过通信网络交互。客户端和服务器的关系通过在相应计算机上运行并且相互具有客户端-服务器关系的计算机程序来产生。

虽然本公开包括一些细节，然而这些细节不应被考虑为对要求保护的本公开的范围的限制，而是作为本公开的示例实施方式的特征的描述。在单个实施方式中，还以组合方式来提供本公开中在独立实施方式的上下文中描述的特定特征。反之，在单个实施方式的上下文中描述的各种特征也在多个实施方式中分开地或以任何合适的子组合来实现。此外，尽管在本文将特征描述为以特定组合起作用并且甚至最初要求如此保护，但是在一些情况下，将来自所要求保护的组合的一个或多个特征从组合中删除，并且所要求保护的组合针对子组合或子组合的变化。

类似地，虽然在附图中以特定顺序描绘了操作，但这不应被理解为要求这些操作以示出的特定顺序或以连续的顺序执行，或者需要执行所有示出的操作来实现希望的结果。在某些环境下，多任务和并行处理是有利的。此外，在本文所描述的实施方式中的各系统组件的分离不应被理解为在所有实施方式中都要求这样的分离，并且应理解，所描述的程序组件和系统可以被一起集成在单个软件产品中或被封装成多个软件产品。

因此，已经描述了本公开的特定实施方式。其他实施方式在下文的权利要求的范围之内。例如，在权利要求书中记载的动作可以以不同顺序来执行，并且仍然实现希望的结果。已经描述了多个实施方式。然而，应理解，在不脱离本公开的精神和范围的情况下做出各种修改。例如，使用本文示出的对步骤进行重新定序、添加或删除的流程的各种形式。因此，其他实施方式在所附权利要求的范围内。

Claims

1.一种方法，包括：

获得用于标识储层的多个井中的每个井的时间推移平均压力数据的数据；

确定所述数据中的多种均质趋势；

基于所述均质趋势将所述数据聚类成多个分组，每个分组与相应趋势相对应；以及

生成包括多个互连储层区域CRR的互连储层区域图，每个CRR与所述多个分组中的一个分组相对应，其中，每个CRR包括所述多个井的与所述CRR的相应分组相对应的子集，并且其中，每个CRR与均质压力响应相关联。

2.根据权利要求1所述的方法，还包括：

针对所述多个CRR中的每个CRR，基于来自该CRR的岩心和测井数据创建1维渗透率函数。

3.根据权利要求1所述的方法，还包括：

针对所述多个CRR中的每个CRR，仅使用该CRR内的井在所述CRR内分配渗透率。

4.根据权利要求1所述的方法，其中，确定所述数据中的所述多种均质趋势还包括：确定所述多个井中的与时间推移平均压力的均质性相关联的两个或更多个井。

5.根据权利要求1所述的方法，其中，生成所述互连储层区域图包括：创建多个沉积多边形，以将与每个分组相对应的井和与其他剩余的分组相对应的井划分开。

6.根据权利要求1所述的方法，其中，每个井的时间推移平均压力是时间推移平均关井压力。

7.根据权利要求1所述的方法，还包括：基于所述互连储层区域图确定所述储层的方差图的范围。

8.根据权利要求1所述的方法，其中，每个CRR的每个井液力连通到相应CRR中的每个其他的井。

9.根据权利要求1所述的方法，其中，每个CRR没有大尺度渗透率阻挡层。

10.一种系统，包括：

一个或多个处理器；以及

非瞬时计算机可读存储介质，耦合到所述一个或多个处理器并存储用于由所述一个或多个处理器执行的编程指令，所述编程指令指示所述一个或多个处理器执行操作，所述操作包括：

确定所述数据中的多种均质趋势；

11.根据权利要求10所述的系统，所述操作还指示所述一个或多个处理器：

针对所述多个CRR中的每个CRR:(i)基于来自该CRR的岩心和测井数据，创建1维渗透率函数，以及(ii)仅使用该CRR内的井在该CRR内分配3维渗透率。

12.根据权利要求10所述的系统，其中，确定所述数据中的所述多种均质趋势还包括：确定所述多个井中的与时间推移平均压力的均质性相关联的两个或更多个井。

13.根据权利要求10所述的系统，其中，生成所述互连储层区域图包括：创建多个沉积多边形，以将与每个分组相对应的井和与其他剩余的分组相对应的井划分开。

14.根据权利要求10所述的系统，其中，每个井的时间推移平均压力是时间推移平均关井压力。

15.根据权利要求10所述的系统，所述操作还包括：基于所述互连储层区域图来定义所述储层的方差图的范围。

16.根据权利要求10所述的系统，其中，每个CRR的每个井液力连通到相应CRR中的每个其他的井。

17.根据权利要求10所述的系统，其中，每个CRR没有大尺度渗透率阻挡层。

18.一种存储指令的非瞬时计算机可读介质，所述指令用于使一个或多个处理器执行操作，所述操作包括：

确定所述数据中的多种均质趋势；

19.根据权利要求18所述的非瞬时计算机可读介质，所述操作还包括：

20.根据权利要求18所述的非瞬时计算机可读介质，其中，生成所述互连储层区域图包括：创建多个沉积多边形，以将与每个分组相对应的井和与其他剩余的分组相对应的井划分开。