CN102834739B - 估计地质体的地质结构的系统和方法 - Google Patents

Abstract

通过生成和/或选择作为感兴趣地质体内的位置的函数来表示感兴趣地质体的特性的感兴趣地质体的一个或多个数值模拟模型来估计感兴趣地质体的地质结构。感兴趣地质体的地质结构的估计可以在储层勘探和/或开发中实现。

Description

估计地质体的地质结构的系统和方法

技术领域

本发明涉及通过生成和/或分析数值模拟模型的组来估计感兴趣地质体的地质结构。

背景技术

在地质勘探领域中，储层勘探和开发需要结构不确定性模型。储层结构参数包括地质体的储层性质分布、走向、堆叠模式、和几何形状。这些参数往往对于储层对开采方案的响应具有关键性控制。通常，局部服从数据(例如，测井记录、岩心、地震数据等)不足以可靠地推断结构参数和它们的相关不确定性。可以应用宽泛的、初步的不确定性模型，但它们是高度主观的，忽略了重要的专业知识，并且在定量判定方面的使用是有限的。

发明内容

本发明的一个方面涉及一种利用计算系统来估计感兴趣地质体的地质结构的计算机实现方法，该计算系统包括被配置成执行一个或多个计算机程序模块的一个或多个处理器。在一个实施例中，该方法包含：在该计算系统可访问的一个或多个电子存储介质上获得感兴趣地质体的多个数值模拟，其中，感兴趣地质体的数值模拟代表作为感兴趣地质体内的位置的函数的感兴趣地质体的多种特性；识别作为感兴趣地质体内的位置的函数在数值模拟内表示的感兴趣地质体的多种特性之间的互相关性；以及指定各个数值模拟和感兴趣地质体之间的对应性的概率，其中，一个给定数值模拟和感兴趣地质体之间的对应性的概率是根据下述来指定的：作为感兴趣地质体内的位置的函数而在给定数值模拟中表示的多种特性的分布和确定的多种特性之间的互相关性。

本发明的另一个方面涉及一种利用计算系统来估计感兴趣地质体的地质结构的计算机实现方法，该计算系统包括被配置成执行一个或多个计算机程序模块的一个或多个电子处理器。在一个实施例中，该方法包含：在该计算系统可访问的一个或多个电子存储介质上获得感兴趣地质体的多个数值模拟，其中，感兴趣地质体的数值模拟代表作为感兴趣地质体内的位置的函数的感兴趣地质体的多种特性；在该计算系统可访问的一个或多个电子存储介质上获得作为感兴趣地质体内的位置的函数在数值模拟内表示的感兴趣地质体的多种特性之间的互相关性；获得局部服从数据，该局部服从数据包括与作为感兴趣地质体内的位置的函数在数值模拟中表示的多种特性中的一种或多种有关的测量值，该测量值是在感兴趣地质体处获得的；以及确定各个数值模拟和感兴趣地质体之间的对应性的概率，其中，一个给定数值模拟和感兴趣地质体之间的对应性的概率是根据下述来确定的：作为感兴趣地质体内的位置的函数在给定数值模拟中表示的多种特性的分布、局部服从数据、以及确定的多种特性之间的互相关性。

本发明的另一个方面涉及一种利用计算系统来估计感兴趣地质体的地质结构的计算机实现方法，该计算系统包括被配置成执行一个或多个计算机程序模块的一个或多个处理器。在一个实施例中，该方法包含：在该计算系统可访问的一个或多个电存储介质上获得感兴趣地质体的多个数值模拟，其中，感兴趣地质体的数值模拟代表作为感兴趣地质体内的位置的函数的感兴趣地质体的多种特性；在该计算系统可访问的一个或多个电存储介质上获得作为感兴趣地质体内的位置的函数在数值模拟内表示的感兴趣地质体的多种特性之间的互相关性；以及确定一个或多个测量参数以便从感兴趣地质体中采集局部服从数据，其中，局部服从数据包括与作为感兴趣地质体内的位置的函数在数值模拟中表示的多种特性中的一种或多种有关的测量值，该测量值是在感兴趣地质体处获得的，以及其中，所述一个或多个测量参数基于数值模拟和在数值模拟内表示的感兴趣地质体的多种特性之间的互相关性。

通过参考相同标号表示各个图形中的相应部件、形成本说明书一部分的附图对如下描述和所附权利要求书加以考察，可以更加清楚地了解本发明的这些和其它目的、特征、和特性，以及操作方法、结构的相关元件的功能、部件的组合和制造的经费。但是，显而易见，这些附图只是为了例示和描述的目的，并未打算作为限制本发明的定义。正如用在说明书和权利要求书中的那样，“一个”、“一种”、和“该”等的单数形式也包括复数指示物，除非上下文另有明确指出。

附图说明

图1例示了依照本发明的一个或多个实施例，被配置成估计感兴趣地质体的地质结构的系统；

图2例示了按照本发明的一个或多个实施例，量化感兴趣地质体的各种参数之间的相互作用的响应面；

图3例示了依照本发明的一个或多个实施例，量化感兴趣地质体的各种参数之间的相互作用的查找表；

图4例示了按照本发明的一个或多个实施例，示出只作为感兴趣地质体的各种参数的函数的对应性的概率的图；

图5例示了依照本发明的一个或多个实施例，局部服从数据如何可被用于区分通过不同数值模拟描述的潜在结构；

图6例示了依照本发明的一个或多个实施例的测量参数的确定；

图7例示了按照本发明的一个或多个实施例，估计感兴趣地质体的地质结构的方法；

图8例示了按照本发明的一个或多个实施例，估计感兴趣地质体的地质结构的方法；以及

图9例示了按照本发明的一个或多个实施例，估计感兴趣地质体的地质结构的方法。

具体实施方式

图1例示了被配置成估计感兴趣地质体的地质结构的系统10。通过生成和/或选择作为感兴趣地质体内的位置的函数表示感兴趣地质体的特性的感兴趣地质体的一个或多个数值模拟模型，来估计感兴趣地质体的地质结构。感兴趣地质体的地质结构的估计可在储层勘探和/或开发时实现。在一个实施例中，系统10包含电子存储器12、用户接口14、一个或多个信息资源16、一个或多个处理器18、和/或其它部件。

在一个实施例中，电子存储器12包含电子地存储信息的电子存储介质。电子存储器12的电子存储介质可以包括与系统10集成配备(即，基本上不可拆卸)的系统存储器和/或经由，例如，端口(例如，USB端口、火线端口等)或驱动器(例如，盘驱动器等)可拆卸地与系统10连接的可拆卸存储器之一或两者。电子存储器12可以包括光可读存储介质(例如，光盘等)、磁可读存储介质(例如，磁带、磁硬盘驱动器、软盘驱动器等)、基于电荷的存储介质(例如，EEPROM、RAM等)、固态存储介质(例如，闪速驱动器等)和/或其它电子可读存储介质中的一种或多种。电子存储器12可以存储软件算法、处理器18所确定的信息、经由用户接口14接收的信息、从信息资源16获得的信息和/或使系统10能够适当起作用的其它信息。电子存储器12可以是系统10内的分立部件，或电子存储器12可在单个设备(或一组设备)中与系统10的一个或多个其它部件(例如，处理器18)集成配备。

用户接口14被配置成提供用户可以向系统10提供信息和从系统10接收信息的系统10与一个或多个用户之间的接口。这使得可在用户与一个或多个电子存储器12、信息资源16、和/或处理器18中的一个或多个之间传送被统称为“信息”的数据、结果和/或指令以及任何其它可传送项目。适合包括在用户接口14中的接口设备的例子包括键板、按钮、开关、键盘、旋钮、手柄、显示屏、触摸屏、扬声器、麦克风、指示灯、声音报警器、和打印机。

应当明白，本发明也可以将其它通信技术，无论是硬连线的还是无线的，设想成用户接口14。例如，本发明设想可以将用户接口14与电子存储器12所提供的可拆卸存储接口集成。在本例中，可将信息从使用户能够定制系统10的实现的可拆卸存储器(例如，智能卡、闪速驱动器、移动盘等)装入系统10中。适合作为用户接口14与系统10一起使用的其它示范性输入设备和技术包括但不限于RS-232端口、RF链路、IR链路、调制解调器(电话、电缆等)。在一个实施例中，可将用户接口14配备在与执行本文分配给系统10的一些或所有功能的服务器可操作地通信的计算平台上。简而言之，本发明将与系统10传送信息的任何技术都设想成用户接口14。

信息资源16包括与感兴趣地质体和/或估计感兴趣地质体的地质结构的过程有关的一个或多个信息源。举非限制性的例子来说，信息资源16之一可以包括与感兴趣地质体的特性的分布有关的一组先前确定的规则。如下面进一步讨论的，这些规则可以包括一种或多种具体地质特性与一个或多个环境参数之间的一种或多种关系、多种地质特性之间的互相关性、对一种或多种地质特性的约束、和/或与感兴趣地质体的特性的分布有关的其它规则。所述规则可以包括对于所有(或几乎所有)模拟地质体通用的规则、和/或各种类型的储层类别、沉积设置、地质区、和/或其它多群或多组地质体专用的规则。所述规则可以包括由一个或多个用户输入和/或修改(例如，经由用户接口14)的规则、和/或自动确定(例如，如下所讨论，由处理器10或一些其它处理器)的规则。

作为信息资源16的另一个非限制性例子，信息资源16可以包括包括一个或多个地质体的局部服从数据的数据集。如本文所使用，“局部服从数据”指的是在地质体处获得的地质体的一种或多种特性的测量值。例如，“局部服从数据”可以包括从位于在地质体处或附近钻出的一口或多口井内的设备获得的测量值、在地质体处或附近的表面上采集的地震数据、和/或地质体的一种或多种特性的其它测量值。

处理器18被配置成在系统10中提供信息处理能力。这样，处理器18可以包括数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电子地处理信息的其它机构中的一种或多种。尽管处理器18在图1中被显示成单个实体，但这只是为了例示的目的。在一些实现中，处理器18可以包括多个处理单元。这些处理单元可以物理地处在同一个设备内，或处理器18可以代表协同运行的多个设备的处理功能。

如图1所示，处理器18可被配置成执行一个或多个计算机程序模块。所述一个或多个计算机程序模块可以包括模型模块20、规则模块22、服从模块24、概率模块26、选择模块28、测量模块30和/或其它模块中的一个或多个。处理器18可被配置成通过软件；硬件；固件；软件、硬件和/或固件的某种组合；和/或用于配置处理器18上的处理能力的其它机构来执行模块20，22，24，26，28和/或30。

应当懂得，尽管模块20，22，24，26，28和30在图1中被例示成共处在单个处理单元内，但在处理器18包括多个处理单元的实现中，模块20，22，24，26，28和/或30中的一个或多个可以处在远离其它模块的地方。对下面所述的不同模块20，22，24，26，28和/或30所提供的功能的描述是为了例示的目的，而不是为了限制的目的，因为模块20，22，24，26，28和/或30的任何一个都可以提供比所述多或少的功能。例如，可以省略模块20，22，24，26，28和/或30中的一个或多个，并且可由模块20，22，24，26，28和/或30中的其它几个来提供其一些或所有功能。作为另一个例子，处理器18可被配置成执行在下面分配给模块20，22，24，26，28和/或30之一的一些或所有功能的一个或多个附加模块。

模型模块20被配置成获得感兴趣地质体的多个数值模拟。模型模块20可将获得的多个数值模拟存储到电子存储器12中。每个获得的数值模拟代表感兴趣地质体的可能地质结构。通过数值模拟表示的作为位置的函数的一种或多种地质特性可以包括多孔性、渗透性、岩相或岩相比例、结构元类型、几何形状、粒度、沉积时的能量或速度、沉积速率、距震源的距离、相对于流轴线和边缘的位置、相对于流底部和顶部的位置、流动方向、流浓度、全粒度分布、流量、盆地外部影响海面升降/水深、沉淀物供给、构造/压缩参数和/或其它地质特性中的一种或多种。

通过数值模拟表示这些地质特性中的一种或多种定义了作为表示的地质特性的函数的相应感兴趣地质体的一个或多个地质特征。例如，数值模拟可以定义水道弯度、其它结构元的弯度、冲裂类型、冲裂频率或速率、淤积速率或频率、水道(或其它结构元)的分布、水道取向的空间分布、在数值模拟中描述的特性的分布、补偿指数、空隙度、连通性/曲折性、对任何特征序列的时序分析、堤防高度和宽度、张开频率、横向冲积速率、后退/推进速率、和/或其它地质特征中的一种或多种。

在一个实施例中，模型模块20被配置成从信息资源16之一获得感兴趣地质体的多个数值模拟中的至少一些。在一个实施例中，模型模块20被配置成实际生成数值模拟中的至少一些。数值模拟可以像在美国专利申请第12/140,901号中描述的那样由模型模块20生成，在此通过引用将其全文并入本公开中。在一些情况下，数值模拟的生成可能考虑局部服从数据。例如，在此通过引用将其全文并入本公开中的美国专利申请(申请人同时提交的美国专利申请参考文献第T-7823号)公开了部分根据局部服从数据来生成感兴趣地质体的多个数值模拟的系统和技术。像在上面通过引用并入的申请中描述的那些那样的技术可以使模型模块20能够获得在地质结构的估计中要作进一步处理的相对较大的一组数值模拟。举非限制性的例子来说，数值模拟的数量可以是几十个到几千个。

规则模块22被配置成获得由数值模拟表示的感兴趣地质体的多种特性之间的互相关性。这些规则提供可被用于约束结构不确定性和/或便于预测地质结构的地质特性之间的相互作用的量化。举非限制性的例子来说，由规则模块22获得的规则可以量化下列中的一个或多个之间的相互作用：淤积速率和净储层体的浓度(例如，较低的淤积速率趋向于导致较高的净储层体的浓度)、冲裂速率和连通性(例如，较高的冲裂速率趋向于导致较低的连通性)、横向阶跃和水道元内的潜在水道轴的保持(例如，在深水水道中，横向阶跃趋向于削弱水道元内的潜在水道轴的保持)、和/或其它相互作用。规则模块22可将获得的规则存储到电子存储器12中以便将来访问和/或处理。

由规则模块22获得的规则可包括一般性规则、敏感性、响应表面、查找表、多元回归模块、和/或量化地质性质之间的相互作用的其它规则中的一种或多种。举例来说，图2示出了量化结构元内(例如，水道内)的净储层体积、淤积速率、和冲裂频率之间的相互作用的响应表面。作为另一个例子，图3示出了量化净储层体积与无组织的水道设置的地质特性和/或特征之间的相互作用的查找表。

回头参照图1，在一个实施例中，规则模块22被配置成获得预定的感兴趣地质体的一条或多条规则。这样的规则可以从，例如，信息资源16之一中获得。预定规则可专用于一种类型的沉积设置和/或与感兴趣地质体相对应的储层类型，或可以更加通用。预定规则可由另一个系统根据以前对代表感兴趣地质体的局部服从数据和/或数值模拟的分析生成。规则模块22可以使一个或多个用户能够在实现之前修改或配置预定规则(例如，经由用户接口14)。

在一个实施例中，规则模块22被配置成根据对模型模块20所获得的感兴趣地质体的数值模拟的分析来生成一条或多条规则。所述规则可以通过在全体数值模拟上观察由数值模拟描述的特性之间的累积关系来生成。例如，在获得的数值模拟内通常可发现一种相对较高水平的特性与另一种相对较低水平的特性相结合。这种关系可由规则模块22以规则的形式来量化。要懂得的是，这种简单例子无意成为限制性的，由通过分析感兴趣地质体的数值模拟而创建的规则量化的这种特性所定义的两种或更多种特性和/或地质特征之间的更复杂关系也落在本公开的范围内。

上面通过分析为感兴趣地质体获得的数值模拟而描述的规则的生成(无论这样的分析实际上由规则模块22进行还是该规则随后由规则模块22从信息资源16之一获得)可以在估计感兴趣地质体的地质结构时提供各种增强。例如，该规则可以量化专用于感兴趣地质体和/或地质特性之间的相互作用的传统理解看起来关系不近或含糊的地质特性之间的互相关性。

在一个实施例中，规则模块22向用户呈现获得的规则(例如，经由用户接口14)。这使用户能够在实现之前审查所述规则，以便更详细地检查在用户看来似乎是由模型模块20获得的数值模拟中的统计异常的结果的规则。一旦用户审查了明显异常的规则和/或规则的基础，规则模块22可使用户能够舍弃该规则，以便在进一步的处理中将不使用该规则，或修改该规则。

服从模块24被配置成获得感兴趣地质体的局部服从数据。服从模块24可从信息资源16之一获得局部服从数据。服从模块24可将获得的局部服从数据存储到电子存储器12中以便将来访问和/或处理。

概率模块26被配置成确定模型模块20所获得的各个数值模拟和感兴趣地质体之间的对应性的概率。给定数值模拟的对应性的概率将给定数值模拟精确地描绘感兴趣地质体的地质结构的可能性表示成概率。对应性的概率由概率模块26根据规则模块22所获得的规则来确定。要懂得的是，概率模块26对呈现通过规则模块22所获得的规则量化的特性之间的互相关性的代表地质结构的数值模拟指定比未呈现通过规则量化的互相关性的代表地质结构的数值模拟高的概率。

举例来说，图4例示了只作为淤积速率和冲裂频率的函数的对应性的概率的图。图4进一步示出了感兴趣地质体的3种示范性结构(如由感兴趣地质体的3个分立数值模拟所表示的)，并且提供与基于该图的这些结构相关联的对应性的概率。在例示在图4中的实施例中，该图是与由如图1所示和如上所述的规则模块22获得的规则类似的规则的直观表示。要懂得的是，在一些情况下，对应性的概率的确定可以根据多条规则作出。在这样的情况下，可以将为例示的结构表示在图4中的概率与按照量化由正被分析的数值模拟表示的特征之间的互相关性的其它规则导出的其它概率结合。

回头参照图1，在一个实施例中，除了实现由规则模块22获得的规则之外，概率模块26在确定对应性的概率时也实现局部服从数据。在这个实施例中，将表示在数值模拟中的特征与局部服从数据相比较。表示的特征与局部服从数据具有较强关联性的数值模拟将趋向于具有比表示的特征与局部服从数据具有较弱关联性的数值模拟更高的对应性概率。

为了以这种方式使用局部服从数据，概率模块26可以实现根据与局部服从数据的关联性来区分由数值模拟表示的潜在结构的选择统计。例如，图5示出了例示局部服从数据如何可被用于区分由不同数值模拟描述的潜在结构的图。具体地说，在显示在图5中的图中，结构参数包括淤积速率和水道有序度，并且使用每口井的平均组元数来区分结构。要懂得的是，关于水道的参数和/或特征的特征描述无意成为限制性的。将相同的原理应用于感兴趣地质体内的其它类型结构元(例如，瓣面、沙洲等)的参数和/或特征都在本公开的范围之内。

回头参照图1，在一个实施例中，概率模块26通过空间自举以从数值模拟中重新取样局部服从数据，然后将来自测井记录的局部服从数据与在数值模拟中进行的重新取样的相比较，作为确定与测井记录相关联的局部服从数据和数值模拟之一之间的关联性的基础。在一个实施例中，局部服从数据包括地震数据。在这个实施例中，局部服从数据与数值模拟之间的关联性的确定基于局部服从数据的完全3D正向地震变换与由数值模拟定义的地质结构的比较。

选择模块28被配置成选择一个或多个数值模拟用于进一步处理。选择模块28根据概率模块26为数值模拟确定的对应性概率来选择一个或多个数值模拟用于进一步处理。在一个实施例中，选择模块28选择具有最高对应性概率的预定数量的数值模拟。所述预定数量可由一个或多个用户配置(例如，经由用户接口14)。在一个实施例中，选择模块28使用对应性概率来加权随机选择，随机地选择预定数量的数值模拟用于进一步处理。在一个实施例中，选择模块28将对应性概率与预定阈值相比较，并选择对应性概率大于预定阈值的数值模拟用于进一步处理。预定阈值可由一个或多个用户配置(例如，经由用户接口14)。

测量模块30被配置成确定用于从感兴趣地质体中采集局部服从数据的一个或多个测量参数。测量模块30所确定的一个或多个测量参数被设计成得到将增强系统10估计感兴趣地质体的地质结构的精度的局部服从数据。例如，一个或多个测量参数可被设计成得到在区分特定一组数值模拟时具有增强利用率的局部服从数据。一个或多个测量参数可包括为了获得服从数据而钻探的井的一个或多个参数、与用于收集局部服从数据的仪器(例如，地震源和/或检测器)的位置和/或取向有关的一个或多个参数、和/或与在感兴趣地质体处或附近采集局部服从数据有关的其它参数。

测量模块30被配置成根据要区分的特定一组数值模拟、规则模块22所获得的规则、和/或以前获得的局部服从数据，来确定一个或多个测量参数。通过分析要区分的该组数值模拟和规则模块22所获得的规则，测量模块30能够识别在确认或否认感兴趣地质体和各个数值模拟之间的对应性时具有增强影响力的局部服从数据。然后，测量模块30确定趋于得到为增强数量的数值模拟确认或否认对应性的局部服从数据的采集局部服从数据的测量参数。

在一个实施例中，测量模块30所确定的采集局部服从数据的测量参数将“最大化”局部服从数据在区分数值模拟时的影响。然而，要懂得的是，在一些实施例中，这种影响的“最大化”在计算上太昂贵(例如，从处理、时间、和/或信息存储的观点来看)，和/或可以不同地定义“最大化”(例如，可以为一个或多个数值模拟提供最大确定性，影响最大数量的数值模拟的对应性概率等)。这样，在无引导地采集局部服从数据期间为了区分数值模拟而增强局部服从数据的影响力未达到实际“最大化”的实施例仍然在本公开的范围之内。

举例来说，图6例示了通过数值模拟表示的几种潜在地质结构。在提供用于按照也显示在图6中的图进一步区分结构的这些地质结构中，示出了通过感兴趣地质体的五个建议井。这五个建议井的参数可由与如上所述和如图1所示的测量模块30相似或相同的测量模块来确定。在一个实施例中，测量模块对参数的确定可以由一个或多个用户引导。例如，在例示在图6中的例子中，一个或多个用户可以设置为采集数据钻探的井的最大数量(例如，5口)，然后，测量模块可以确定便于区分潜在地质结构的5口井的测量参数。

返回到图1，地质特性模块31被配置成确定感兴趣地质体的一种或多种感兴趣地质特性。地质特性模块31被配置成根据模型模块20所获得的和/或选择模块28所选择的数值模拟来确定一种或多种地质特性。感兴趣地质特性一般包括与位于感兴趣地质体中的矿物燃料储层有关的地质特性。例如，地质特性可以包括净毛比、连通性、结构几何形状、采收率和/或其它特性中的一种或多种。

这些地质特性可以不同于服从模块24和概率模块26用来确定数值模拟与感兴趣地质体之间的关联性的概率的地质特性。例如，一维地质特性可被服从模块24和概率模块26实现，以根据一维井数据来确定关联性的概率。这样的地质特性可以包括，例如，混合比、组元和组元间的频率分布、组元和组元间的厚度分布、一维空隙度、半方差参数、转移概率和/或其它特性中的一种或多种。但是，一旦加权和/或选择了(例如，由选择模块28)数值模拟，地质特性模块31使用感兴趣地质体的数值模拟来空间自举，以确定三维地质特性。

在一个实施例中，地质特性模块31通过分析多个数值模拟来确定感兴趣地质特性。这些数值模拟可以包括模型模块20所获得的所有数值模拟，或选择模块28所选择的一组数值模拟。举非限制性的例子来说，地质特性模块31可以扫描感兴趣地质体的每个数值模拟，并为每个数值模拟记录(i)服从模块24、概率模块26、和/或选择模块28用来加权和/或选择不同数值模拟的地质特性的值(例如，混合比、组元和组元间的频率分布、组元和组元间的厚度分布、一维空隙度等)，以及(ii)通过空间自举确定的感兴趣地质特征的值(例如，净毛比、连通性、结构几何形状、采收率等)。

一旦为感兴趣地质体的每个模型进行了这种扫描，地质特性模块31可以绘制(或数学模拟，或要不然关联)感兴趣地质特性的确定值与服从模块24、概率模块26、和/或选择模块28所使用的一种或多种直接测量的地质特性之间的关系。一旦为感兴趣地质体的每种模型完成了这种操作，地质特性模块31就可以实现这种绘制(或要不然关联)，以便绘制(或要不然表达)感兴趣特性的条件概率与可在局部数据中观察到的一种或多种地质特性(例如，来自井数据的一维特性)之间的关系。这样，地质特性模块31能够根据局部服从数据和模型模块20所获得的数值模拟，用公式表示结构不确定性模型。

图7例示了估计感兴趣地质体的地质结构的方法32。下面展示的方法32的操作旨在例示。在一些实施例中，方法32可以通过未描述的一个或多个附加操作，和/或没有所讨论的一个或多个操作而完成。另外，在图7中例示和下面描述方法32的操作的次序无意成为限制性的。

在一些实施例中，方法32可以在包含一个或多个处理器(例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电处理信息的其它机构)的计算系统中实现。一个或多个处理设备可以包括响应于电子地存储在电子存储介质上的指令来执行方法32的一些或所有操作的一个或多个设备。一个或多个处理器可以包括通过硬件、固件和/或软件配置成专门被设计用于执行方法32的一个或多个操作的一个或多个设备。

在操作34中，获得感兴趣地质体的数值模拟。获得数值模拟可以包括生成一个或多个数值模拟和/或获得一个或多个先前存储的模拟。在一个实施例中，操作34由与模型模块20(如图1所示和如上所述)相同或相似的模型模块来执行。

在操作36中，获得在数值模拟中描述的特性之间的互相关性。获得互相关性可以包括分析数值模拟以确定互相关性和/或获得先前存储的互相关性。互相关性可以包括一个或多个用户生成、修改和/或配置的互相关性。在一个实施例中，互相关性可以通过一条或多条规则被量化。在这个实施例中，获得互相关性可以包括获得量化互相关性的规则。在一个实施例中，操作36由与规则模块22(如图1所示和如上所述)相同或相似的规则模块来执行。

在操作38中，确定对应性的概率和/或将其指定给各个数值模拟。指定给给定数值模拟的对应性概率表达感兴趣地质体的实际地质结构对应于由给定数值模拟描述的地质结构的概率。通过将在操作36中获得的互相关性应用于数值模拟来确定数值模拟的对应性概率。在一个实施例中，操作38由与概率模块26(如图1所示和如上所述)相同或相似的概率模块来执行。

在操作40中，接受第一组数值模拟用于进一步分析和/或处理，以及拒绝第二组数值模拟用于进一步分析和/或处理。在操作40中，根据在操作38中确定的对应性概率将数值模拟分成第一组和第二组。在一个实施例中，操作40由与选择模块28(如图1所示和如上所述)相同或相似的选择模块来执行。

图8例示了估计感兴趣地质体的地质结构的方法42。下面展示的方法42的操作旨在例示。在一些实施例中，方法42可以通过未描述的一个或多个附加操作，和/或没有所讨论的一个或多个操作而完成。另外，在图8中例示和下面描述方法42的操作的次序无意成为限制性的。

在一些实施例中，方法42可以在包含一个或多个处理器(例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电处理信息的其它机构)的计算系统中实现。一个或多个处理器可以包括响应于电子地存储在电子存储介质上的指令来执行方法42的一些或所有操作的一个或多个设备。一个或多个处理器可以包括通过硬件、固件和/或软件配置成专门被设计用于执行方法42的一个或多个操作的一个或多个设备。

在操作44中，获得与感兴趣地质体的特性有关的局部服从数据。可以将获得的局部服从数据存储到电子存储器中。在一个实施例中，操作44由与服从模块24(如图1所示和如上所述)相同或相似的服从模块来执行。

在操作46中，获得感兴趣地质体的数值模拟。获得数值模拟可以包括生成一个或多个数值模拟和/或获得一个或多个先前存储的模拟。在获得的数值模拟中表示的至少一些特性与在操作44中获得的局部服从数据有关。在获得数值模拟包含生成至少一些数值模拟的实施例中，生成的数值模拟可以至少部分基于在操作44中获得的局部服从数据。在一个实施例中，操作46由与模型模块20(如图1所示和如上所述)相同或相似的模型模块来执行。

在操作48中，获得在数值模拟中描述的特性之间的互相关性。获得互相关性可以包括分析数值模拟以确定互相关性和/或获得先前存储的互相关性。互相关性可以包括一个或多个用户生成、修改和/或配置的互相关性。在一些实施例中，互相关性可以通过一条或多条规则被量化。在这样的实施例中，获得互相关性可以包括获得量化互相关性的规则。在一个实施例中，操作48由与规则模块22(如图1所示和如上所述)相同或相似的规则模块来执行。

在操作50中，确定对应性的概率和/或将其指定给各个数值模拟。指定给给定数值模拟的对应性概率表达感兴趣地质体的实际地质结构对应于通过给定数值模拟所描述的地质结构的概率。通过比较给定数值模拟所描述的地质结构与局部服从数据，并且将在操作48中确定的互相关性应用于给定数值模拟和/或与给定数值模拟有关的局部服从数据，来确定给定数值模拟的对应性概率。在一个实施例中，操作50由与概率模块26(如图1所示和如上所述)相同或相似的概率模块来执行。

在操作52中，接受第一组数值模拟用于进一步分析和/或处理，以及拒绝第二组数值模拟用于进一步分析和/或处理。在操作52中，根据在操作50中确定的对应性概率将数值模拟分成第一组和第二组。在一个实施例中，操作52由与选择模块28(如图1所示和如上所述)相同或相似的选择模块来执行。

图9例示了估计感兴趣地质体的地质结构的方法54。下面展示的方法54的操作旨在例示。在一些实施例中，方法54可以通过未描述的一个或多个附加操作，和/或没有所讨论的一个或多个操作而完成。另外，在图9中例示和下面描述方法54的操作的次序无意成为限制性的。

在一些实施例中，方法54可以在包含一个或多个处理器(例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或电处理信息的其它机构)的计算系统中实现。一个或多个处理器可以包括响应于电子地存储在电子存储介质上的指令来执行方法54的一些或所有操作的一个或多个设备。一个或多个处理设备可以包括通过硬件、固件和/或软件配置成专门被设计用于执行方法54的一个或多个操作的一个或多个设备。

在操作56中，获得代表感兴趣地质体的潜在地质结构的数值模拟、数值模拟的关联性的概率、和在数值模拟中表示的特性之间的互相关性。在一个实施例中，在操作56中还获得感兴趣地质体的局部服从数据。在一个实施例中，操作56包括方法32(如图6所示和如上所述)。在一个实施例中，操作56包括方法42(如图7所示和如上所述)。在这些实施例中，在操作56中获得的数值模拟是在操作40或52中接受的用于附加处理的第一组数值模拟。

在操作58中，确定用于从感兴趣地质体中采集局部服从数据的测量参数。测量参数被确定成导致便于区分在操作56中获得的数值模拟的局部服从数据的采集。测量参数的确定基于通过各个数值模拟描述的地质结构、通过数值模拟表示的特性之间的互相关性、和/或先前获得的局部服从数据。在一个实施例中，操作58由与测量模块30(如图1所示和如上所述)相同或相似的测量模块来执行。

在操作60中，获得在符合在操作58中确定的测量参数的一次或多次测量期间采集的感兴趣地质体的局部服从数据。这种局部服从数据是由一个或多个地球性质测量系统在感兴趣地质体处或附近采集的。局部服从数据可以在网络上经由可拆卸电子存储器，和/或通过将信息电子地传送给实现方法54的计算机系统的其它机构来获得。在一个实施例中，操作60由与处理器18(如图1所示和如上所述)相同或相似的处理器来执行，该处理器执行导致获得局部服从数据以便用在处理器上的附加处理中的一个或多个计算机程序模块。

在操作62中，提炼数值模拟的关联性概率和/或根据在操作60中获得的局部服从数据(和/或以前获得的在数值模拟中表示的特性之间的互相关性)拒绝至少一些数值模拟用于附加处理。在一个实施例中，操作62由与概率模块26和/或选择模块28(如图1所示和如上所述)相同或相似的概率模块和/或选择模块来执行。

在操作64中，如果在操作62中对关联性概率的提炼和拒绝数值模拟用于附加处理导致选择了最终一组数值模拟，则结束该方法54，并且可以从该最终一组被概率加权的数值模拟中估计感兴趣地质体的一种或多种特性。但是，如果操作64未导致最终一组数值模拟(例如，局部服从数据未提供足以将数值模拟减少到最终一组的明了性)，则方法54返回到操作58以便能够进一步减少数值模拟。

尽管为了例示的目的，根据当前认为最实用的优选实施例对本发明作了详细描述，但应该明白，这样的细节仅仅为了那个目的，本发明不局限于所公开的实施例，而是相反，旨在涵盖在所附权利要求书的精神和范围之内的所有修改和等效安排。例如，要明白的是，本发明设想，可以尽可能地将任何实施例的一个或多个特征与任何其它实施例的一个或多个特征组合。

Claims (6)

1.一种利用计算系统来估计感兴趣地质体的地质结构的计算机实现的方法，该计算系统包括被配置成执行一个或多个计算机程序模块的一个或多个电子处理器，该方法包含：

在所述计算系统可访问的一个或多个电子存储介质上获得感兴趣地质体的多个数值模拟，其中，感兴趣地质体的数值模拟代表作为感兴趣地质体内的位置的函数的感兴趣地质体的多种特性；

在所述计算系统的一个或多个电子处理器上执行被配置成识别作为感兴趣地质体内的位置的函数在数值模拟内表示的感兴趣地质体的多种特性之间的互相关性的一个或多个计算机程序模块；以及

在所述计算系统的一个或多个电子处理器上执行被配置成指定各个数值模拟与感兴趣地质体之间的对应性的概率的一个或多个计算机程序模块，其中，根据作为感兴趣地质体内的位置的函数在一个给定数值模拟中表示的多种特性的分布和确定的多种特性之间的互相关性，来指定所述给定数值模拟与感兴趣地质体之间的对应性的概率。

2.如权利要求1所述的方法，其中，识别作为感兴趣地质体内的位置的函数在数值模拟内表示的感兴趣地质体的多种特性之间的互相关性包含：识别对一种或多种特性的分布的约束，所述一种或多种特性在能够在将来模拟感兴趣地质体时或为了提炼感兴趣地质体的一种或多种现有模型而实现的数值模拟中表示。

3.如权利要求1所述的方法，其中，作为感兴趣地质体内的位置的函数在数值模拟内表示的多种特性包含下列中的一种或多种：多孔性、渗透性、岩相、岩相比例、结构元类型、几何形状、粒度、沉积时的能量或速度、沉积速率、距震源的距离、相对于流轴线和/或边缘的位置、相对于流底部和/或顶部的位置、流动方向、流浓度、全粒度分布、流量、盆地外部影响海面升降/水深、沉淀物供给、或构造/压缩参数。

4.如权利要求1所述的方法，其中，识别作为感兴趣地质体内的位置的函数在数值模拟内表示的感兴趣地质体的多种特性之间的互相关性包含：识别数值模拟的地质特征之间的互相关性，其中，数值模拟的地质特征是作为地质体内的位置的函数在数值模拟内表示的特性的函数。

5.如权利要求4所述的方法，其中，数值模拟的地质特征包含下列中的一种或多种：水道弯度、结构元弯度、冲裂类型、冲裂频率或速率、淤积速率或频率、水道的分布、在数值模拟中表示的一种或多种特性在各个水道中的分布、在数值模拟中表示的一种或多种特性在任何结构元外部的分布、补偿指数、空隙度、连通性/曲折性、对任何特征序列的时序分析、堤防高度和宽度、张开频率、横向加积速率、或后退/推进速率。

6.如权利要求1所述的方法，其中，在所述计算系统可访问的一个或多个电子存储介质上获得感兴趣地质体的多个数值模拟包含：在所述计算系统的一个或多个电子处理器上执行被配置成生成多个数值模拟的一个或多个计算机程序模块。