CN103270514B - 用于生成所关注的地质体积的、受所关注的地质体积的基于过程的模型约束的地质统计模型的系统和方法 - Google Patents
用于生成所关注的地质体积的、受所关注的地质体积的基于过程的模型约束的地质统计模型的系统和方法 Download PDFInfo
- Publication number
- CN103270514B CN103270514B CN201280004208.XA CN201280004208A CN103270514B CN 103270514 B CN103270514 B CN 103270514B CN 201280004208 A CN201280004208 A CN 201280004208A CN 103270514 B CN103270514 B CN 103270514B
- Authority
- CN
- China
- Prior art keywords
- model
- kernel
- based methods
- close attention
- paid close
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 193
- 230000003750 conditioning effect Effects 0.000 claims abstract description 76
- 230000008569 process Effects 0.000 claims description 54
- 230000008878 coupling Effects 0.000 claims description 40
- 238000010168 coupling process Methods 0.000 claims description 40
- 238000005859 coupling reaction Methods 0.000 claims description 40
- 238000005259 measurement Methods 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 16
- 238000013179 statistical model Methods 0.000 claims description 13
- 230000007704 transition Effects 0.000 claims description 8
- 238000013139 quantization Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 125000006850 spacer group Chemical group 0.000 claims 1
- 241001269238 Data Species 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 description 17
- 238000003860 storage Methods 0.000 description 17
- 238000012545 processing Methods 0.000 description 11
- 239000011435 rock Substances 0.000 description 9
- 238000009826 distribution Methods 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 5
- 241001597008 Nomeidae Species 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 239000008187 granular material Substances 0.000 description 4
- 230000008021 deposition Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 241001484259 Lacuna Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- -1 flow velocity rate Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 210000002706 plastid Anatomy 0.000 description 1
- 230000001376 precipitating effect Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000008080 stochastic effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
- Feedback Control In General (AREA)
- General Factory Administration (AREA)
Abstract
生成一种所关注的地质体积的基于过程的模型。该基于过程的模型利用与所关注的地质体积相关联的条件化信息来条件化。从该基于过程的模型生成局部地表示该所关注的地质体积的参数的统计数据。这些统计数据被用于约束所关注的地质体积的一个或更多个地质统计模型。
Description
技术领域
本公开一般涉及生成受所关注的地质体积的基于过程的模型约束的地质统计（geostatistical）模型。
背景技术
用于生成所关注的地质体积的基于过程的模型的系统和方法是已知的。生成基于过程的模型，以获得对被形成以创建所关注的地质体积的过程的理解。基于过程的模型可以被条件化成条件化信息，举例来说，如地震信息、井数据信息、模拟信息、钻心数据、以及和/或其它信息，但条件化通常是近似和/或不实用的。
用于生成所关注的地质体积的地质统计储集层模型的系统和方法是已知的。地质统计模型可以被条件化成条件化信息，但通常缺乏任何过程信息。
发明内容
本公开的一个方面涉及一种系统，该系统被配置成生成所关注的地质体积的地质统计模型。在一些实施例中，该系统包括：一个或更多个处理器，该一个或更多个处理器被配置成执行条件化信息模块、过程模型模块、统计模块、以及统计模型模块。该条件化信息模块被配置成获取与所关注的地质体积相关联的条件化信息，其中，该条件化信息包括从在所关注的地质体积处或其附近所进行的测量导出的信息。该过程模型模块被配置成生成所关注的地质体积的基于过程的模型，其中，该基于过程的模型符合与所关注的地质体积相关联的条件化信息。该统计模块被配置成，根据该基于过程的模型生成忠实于该基于过程的模型中的过程信息的统计数据，其中，该统计数据中的单个统计数据在该基于过程的模型内局部地确定，以使给定统计数据的值作为该基于过程的模型内的位置的函数而改变。该统计模型模块被配置成，生成该所关注的地质体积的、受由该统计模块所生成的统计数据约束的地质统计模型。
本公开的另一方面涉及一种生成所关注的地质体积的地质统计模型的方法。在一些实施例中，该方法包括：获取与所关注的地质体积相关联的条件化信息，其中，该条件化信息包括从在所关注的地质体积处或其附近所进行的测量导出的信息；生成所关注的地质体积的基于过程的模型，其中，该基于过程的模型符合与所关注的地质体积相关联的条件化信息；根据该基于过程的模型生成忠实于该基于过程的模型中的过程信息的统计数据，其中，该统计数据中的单个统计数据在该基于过程的模型内局部地确定，以使给定统计数据的值作为该基于过程的模型内的位置的函数而改变；以及生成所关注的地质体积的、受由统计模块所生成的统计数据约束的地质统计模型。
本公开的又一方面涉及存储机器可读指令的非暂时电子存储介质，该机器可读指令使一个或更多个处理器执行导致生成所关注的地质体积的地质统计模型的操作。在一些实施例中，该操作包括：获取与所关注的地质体积相关联的条件化信息，其中，该条件化信息包括从在所关注的地质体积处或其附近所进行的测量导出的信息；生成所关注的地质体积的基于过程的模型，其中，该基于过程的模型符合与所关注的地质体积相关联的条件化信息；根据该基于过程的模型生成忠实于该基于过程的模型中的过程信息的统计数据，其中，该统计数据中的单个统计数据在该基于过程的模型内局部地确定，以使给定统计数据的值作为该基于过程的模型内的位置的函数而改变；以及生成所关注的地质体积的、受由该基于过程的模型所生成的统计数据约束的地质统计模型。
当参照附图考虑下面的描述和所附权利要求书时，在此公开的系统和/或方法的这些和其它目的、特征以及特性，和结构与组合部分的相关部件的操作与功能的方法以及制造的经济性将变得更清楚，其全部形成了本说明书的一部分，其中，相同标号指定各个图中的对应部分。然而，应当明白，附图仅仅是出于例示和描述的目的，而非旨在作为对本发明的限制的解说。如在本说明书和权利要求书中使用的，单数形式“一（a）”、“一（an）”，以及“该/所述（the）”包括多个指代物，除非上下文另外清楚地规定。
附图说明
图1例示了生成所关注的地质体积的地质统计储集层模型的方法。
图2例示了被配置成生成所关注的地质体积的地质统计储集层模型的系统。
具体实施方式
本技术可以按一系统和要通过一计算机执行的计算机方法的一般背景来描述和实现。这种计算机可执行指令可以包括程序、例程、对象、组件、数据结构、以及可以被用于执行特定任务和处理抽象数据类型的计算机软件技术。本技术的软件实现可以按针对多种计算平台和环境中的应用的不同语言来编码。应当清楚，本技术的范围和基本原理不限于任何特定计算机软件技术。
而且，本领域技术人员应当清楚，本技术可以利用硬件和软件构造中的任一个或组合来实践，包括但不限于，具有单一和/或多处理器计算机处理器系统的系统、手持式装置、可编程消费类电子设备、迷你计算机、大型计算机等。本技术还可以在其中通过经由一个或更多个数据通信网络链接的服务器或其它处理装置执行任务的分布式计算环境中实践。在分布式计算环境中，程序模块可以位于包括存储器存储装置的本地和远程计算机存储介质两者中。
而且，用于与计算机处理器一起使用的制造品（如CD、预记录盘或其它等同装置）可以包括计算机程序存储介质和记录在其上的、用于引导计算机处理器以易于实现和实践本技术的程序装置。这种装置和制造品也落入本技术的精神和范围内。
下面，参照附图，对本技术的实施例进行描述。本技术可以按许多方式来实现，例如包括系统（包括计算机处理系统）、方法（包括计算机实现方法）、装置、计算机可读介质、计算机程序产品、图形用户接口、门户网站、或者有形地固定在计算机可读存储器中的数据结构。下面，对本技术的几个实施例进行讨论。附图仅例示了本技术的典型实施例，并由此不应被视为对其范围和宽度的限制。
图1例示了生成所关注的地质体积（例如，包括地下体积）的地质统计模型的方法10。下面呈现的方法10的操作旨在例示。在一些实施例中，方法10可以利用未描述的一个或更多个附加操作和/或不利用所讨论操作中的一个或更多个来完成。另外，其中图1例示和下面描述的方法10的操作的次序不旨在进行限制。
在一些实施例中，方法10可以在一个或更多个处理装置（例如，数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或用于电子地处理信息的其它机构）中实现。该一个或更多个处理装置可以包括响应于电子地存储在电子存储介质上的指令来执行方法10的一些或全部操作的一个或更多个装置。该一个或更多个处理装置可以包括通过要具体设计用于执行方法10的一个或更多个操作的硬件、固件、以及/或软件配置的一个或更多个装置。
在操作12，获取与所关注的地质体积相关联的条件化信息。该条件化信息包括从在所关注的地质体积处或其附近所进行的测量导出的信息。该条件化信息可以指定所关注的地质体积的、针对该所关注的地质体积内的特定位置的性质和/或参数。这种信息例如可以包括：井数据、钻心数据、地震数据、模拟信息、专门（expert）信息、和/或其它信息。该性质和/或参数可以包括以下的一个或更多个：测量值、概括性统计数据（包括：平均值、比例、方差、相关性等）、单变量、多变量以及/或空间分布（纹理、协方差、走向分布等）、几何结构（厚度、斜率、波动等）和任何趋向、组合、和/或衍生物、以及/或其它性质或参数。该条件化信息例如可以从电子存储部获取，该电子存储部被配置成，非暂时地存储先前从用户输入、通过网络、以及/或从其它源获取的条件化信息。
在操作14，生成一种所关注的地质体积的基于过程的模型。如在此使用的，基于过程的模型可以指这样的所关注的地质体积的模型：其通过从初始储集层（或其它地质结构）开始，接着对初始储集层和/或先前建模的过程顶部的各个过程（例如，流、坡身不稳（slopefailure）、局部侵蚀和沉积、剥蚀和淤积、分解和沉淀、胶结、生长、压实和构造修改、以及/或其它过程类型）建模来生成。所建模的过程在形成期间生成并且应用于基于过程的模型，以试图表示和/或重建随着时间发生的原始过程，以创建所关注的地质体积。在操作14生成的基于过程的模型表示针对遍布所关注的地质体积的位置的岩石类型（包括：粒度大小、趋向、类型等）和储集层性质（包括：孔隙度、渗透性）及其相关联的分布、趋向以及几何结构。
单个过程的一个或更多个参数可以推测地和/或确定性地确定。例如，可以推测地和/或确定性地确定各个过程的初始盆形、裂源（origin）、路径、沉积物流量、类型以及位置（和相关定时）、压实和构造历史、局部基准面变化、以及/或其它参数。通过非限制例的方式，于2009年12月23日提交的、题名为“System And Method ForModeling A Geologic Volume of Interest”的美国专利申请（其全部内容通过引用而并入本申请中）提供了对基于过程的建模的描述。
该基于过程的模型中的过程被生成为针对一算法的输出。该算法包括多个输入，举例来说，如初始地形（包括形状和梯度）和相关联的初始盆地岩石性质（粒度大小、内聚性、侵蚀性、岩石类型等）中的一个或更多个。而且包括针对相关联的过程的输入，如流（例如、排出、粒度大小以及岩石类型等）、岩石生成（例如，碳酸盐生长、分解和沉淀等）、构造变化（例如，断层、速率、速率变化）、海面升降周期（例如，幅度和频率以及速率变化）、压实（例如，基于岩石类型和其它信息的速率）、成岩蚀变（例如，在沉积之后的岩石性质作为岩石类型和局部状态等的函数的变化）、以及其它运输机制（例如，涌流、风力、坡身不稳等）的参数化。这些可以基于物理定律、经验规则或者可以通过专门知识直接施加在该模型上（例如，代替利用压实法，可以随着时间步映射压实并且施加），并且可以针对时间和局部性而全部被认为是静态或可变的。可设想其它输入。
生成基于过程的模型可以包括：针对该模型一次生成并应用一个过程，和/或在形成所关注的实际地质体积期间生成并且应用与某时间增量相对应的一系列过程。用于生成单个过程和/或一系列过程的输入可以包括：对于该模型的多个过程/系列相同的一个或更多个输入，和/或在该模型的多个过程/系列之间调节的一个或更多个输入（例如，用于考虑在形成所关注的实际地质体积期间随着时间的状态变化）。
在一些实施例中，用于在操作14生成基于过程的模型的技术被配置成针对在操作12获取的条件化信息条件化该模型。这种条件化可以包括：根据条件化信息约束对一过程的一个或更多个参数的推测性确定，基于条件化信息调节针对被用于生成该模型的算法的一个或更多个输入，以及/或用于条件化该过程的其它机制。该条件化可以是近似的和/或准确的。准确是这样一种模型状态，即，条件化在数据位置上精确地遵守，并且对于基于过程的模型通常不可能或不实际。这意指在所测量参数被包括在条件化数据中的特定位置处，在基于过程的模型中该特定位置处的同一参数的值被约束成该条件化数据中的测量值。近似是这样一种状态，即，在该条件化中观察到的相同的一般特征在该基于过程的模型中如通过任何局部统计数据（包括部分19或任何其它部分中列出的那些）确定的同一位置中再现。
在操作16，获取匹配标准。该匹配标准建立在操作14生成的基于过程的模型与在操作12获取的条件化信息之间的必需对应。如果具有基于过程的模型的结构的真实地质体积被希望生成和在操作12获取的条件化信息相同的条件化信息，则认为该基于过程的模型与条件化信息“匹配”。通过非限制例的方式，该匹配标准可以建立所需的匹配程度、必须匹配的参数、应当超出或不超出的阈值、应当存在的状态、以及/或其它标准。必须匹配的参数可以包括以下的一个或更多个：局部地层相关样式、用于约束储集层网格和体积的隔室和边界、具有唯一统计数据的局部区域或总体（populations）以及特征化这种区域之间的过渡的统计数据、包括几何结构的本地可变空间连续性模型、半方差图参数和连续性的方向、针对可推测协同仿真（cosimulation）的局部可变相关系数、局部可变过渡概率、二变量和多变量分布、多点统计数据、以及和/或其它参数。
各种标准可以基于每个过程（或其它局部化区域）、每个条件化信息源、以及/或其它分段来建立。在一些实施例中，匹配程度通过量化基于过程的模型与条件化信息之间的匹配的匹配度量的阈值级来设置。该匹配标准可以通过自动化确定和/或基于用户输入来获取。匹配标准可以例如基于通过方法10生成的地质统计模型的当前和/或将来应用来确定和/或调节。例如，出于体积测定的目的，在不考虑异方差性的情况下可以应用简单几何形状匹配标准。在另一个实施例中，针对详细的井布局，匹配标准可以考虑在条件化数据中可获得的高分辨率异质性。
在操作18，量化基于过程的模型与条件化信息之间的匹配。该量化可以包括：生成指示基于过程的模型与条件化信息之间的匹配的紧密度的一个或更多个匹配度量。该匹配度量可以基于每个条件化信息源、每个条件化信息类型、以及/或其它特定信息来确定。该匹配度量在该基于过程的模型内局部地确定。这包括针对各个过程、该基于过程的模型内（例如，单一过程内和/或多个过程内）的其它局部化区域、基于过程的模型内的各个位置、以及/或该基于过程的模型内的其它方位来确定该匹配度量。
在操作20，确定该基于过程的模型是否充分匹配条件化信息。例如，在操作18确定的量化匹配可以与在操作16获取的匹配标准相比较，以确定该基于过程的模型是否充分匹配条件化信息，以满足该匹配标准。在一些实施例中，这包括比较在操作18确定的一个或更多个匹配度量与在操作16获取的一个或更多个匹配阈值。响应于确定该基于过程的模型不充分匹配条件化信息，方法10可以返回至操作14以生成另一基于过程的模型，在操作22调节该基于过程的模型，以及/或采取其它动作。
在一些实施例中，部分地或整个地基于由用户/专家人工检查该基于过程的模型来执行操作20。该检查可以补充对该基于过程的模型与条件化信息之间的匹配的自动化分析，和/或可以代替这种分析。
应当清楚，图1的例示和在此讨论的针对基于过程的模型执行的操作14、16、18、以及20总体上不是进行限制。在一些实施例中，操作14、16、18、以及/或20可以作为单个过程迭代地执行、或作为过程系列执行，被添加至该模型并且在下一个过程之前通过操作18和20验证、或者添加过程系列。
在操作24，根据该基于过程的模型生成忠实于该基于过程的模型中的过程信息的统计数据。所生成的统计数据忠实于该基于过程的模型中的过程信息，因为推测地生成的统计数据被约束成对应于和/或不抵触该基于过程的模型中的过程信息。通过非限制例的方式，从储集层模型和过程模型两者直接提取变量，如方位和/或接近度、厚度、储集层特性（如岩石类型、相图、孔隙度和渗透性、取向、取向或厚度变化、形状或几何结构、多点统计数据、连通性、流响应、流障和管道）、以及从储集层模型推断并且与该过程模型进行一致性比较的变量（如梯度、梯度变化、纹理、相距源的距离或过程过渡、过程变化、与条件化数据附近的侵蚀或沉积特征相关联的古（paleo）设置（包括水深、流速率、沉积物源序列、位置和类型、粒度大小分布、浓度等）），以及这些或其它的任何组合和/或衍生物。
该统计数据在该基于过程的模型内局部地确定，以使得给定统计数据的值作为基于过程的模型内的位置的函数而变化。在生成地质统计储集层模型中使用这些统计数据可以被视为将过程信息集成到储集层模型中。这些统计数据可以在操作24直接从基于过程的模型或者基于过程的模型的任何衍生物（包括解释、数值处理、或求导）提取。
用户可以确定该提取的特定性（例如，通过调节局部窗口的尺寸以平均化，和/或通过其它机制），以从该基于过程的模型提取被认为对于该模型的应用重要的、并且被认为信息丰富且局部准确的信息级别。在一些实施例中，在该基于过程的模型内的给定位置或区域处提取（或可用于提取）的局部统计数据的分辨率（或特定性）可以取决于针对给定位置或区域与条件化信息的匹配程度。例如，所提取的或可用于提取的统计数据的分辨率可以是在操作18针对该给定位置或区域确定的一个或更多个匹配度量的函数。
所提取的统计数据可以包括分别特征化为局部比例或针对分类或连续储集层性质的手段的趋势。其它局部统计数据可以包括但决不限于：局部地层相关样式、用于约束储集层网格和体积的隔室和边界、具有唯一统计数据的局部区域或总体以及特征化这种区域之间的过渡的统计数据、局部相分类、局部变量连续和分类分布以及所有关联的分布参数、包括几何结构的本地可变空间连续性模型、连续性的半方差图参数和方向、针对可推测协同仿真的局部可变相关系数、局部可变过渡概率、二变量和多变量分布、多点统计数据、图案、以及/或其它局部统计数据。另外，任何其它空间统计数据（例如，基于点（如最近邻域和Rippley-K函数等）、基于光栅（如数学形态算子和滤波器得分、缺项）、基于表面（如梯度、曲率、皱曲以及/或其它））或连通性或流代理（flow proxy）可以被计算以用于储集层模型的约束。以及这些统计数据在任何标度的任何组合和衍生物，或表示这些统计数据在标度或位置之间的变化的统计数据。这些统计数据的多个集合（和/或对应约束）可以根据多个基于过程的模型和/或通过统计数据推断的多个实现来构造，以提供相关联的约束下的不确定性模型。
操作24提取的统计数据中的一个或更多个可以保存来自该基于过程的模型的过程信息。这可以包括：局部梯度、古流（paleoflow）、水深测量（bathymetry）、水深、以及和/或其它处理信息。虽然这些可以不直接通知地质统计模型，但这些参数与地质统计特性之间的关系被建模成全部用于局部约束的进一步改进。
在操作26，一个或更多个地质统计储集层模型利用来自在操作24生成的基于过程的模型统计数据的局部约束进行构造。在一些实施例中，多个地质统计模型可以根据多个统计数据输入来构造，表示这些输入中的不确定性。可以仿真针对每一组统计约束的多个推测性实现，以表示随机效应的影响。对于某些统计数据来说，如网格、体积、全局比例、局部变量比例、半方差、几何结构、以及/或其它统计数据，与适当的建模算法集成是容易的，因为这些是标准输入。对于其它非标准统计数据（如过渡概率、其它空间统计数据和流代理、以及/或其它统计数据）来说，可以应用多种技术以改进地质统计模型与生成的统计数据的统计匹配。对于基于规则的其它约束来说，可以应用另一些基于规则的建模技术。
在操作26，如果生成地质统计模型提供了有关过程的新见解，则方法10可以返回至操作14，以在新的基于过程的模型中包括这种见解。确定是否应当执行这种迭代（其在图1中未示出）可以由用户进行，和/或自动地进行。
在操作28，可以验证在操作26生成的地质统计模型。集成多个信息源和/或复杂工作流可以引入不期望的模型特征、数值伪像或者增加严重错误的机会。这种验证可以包括：人工和/或自动检查该基于过程的模型与地质统计模型之间的匹配、条件化信息与地质统计模型之间的匹配，搜索条件化信息附近和/或远离条件化信息的区域中的伪像，以及和/或其它验证过程。
响应于确定地质统计模型包括使得产生对该地质统计模型的准确度的怀疑的一个或更多个伪像，进行方法10。如果伪像出现在多个模型上，则应当检查在操作12获取的条件化数据的矛盾和不准确，并且有可能修正或去除部分条件化数据，并且在操作14重新生成过程模型。如果伪像出现在有限数量的模型中，则可以检查操作14的对应过程模型的不准确和伪像，并且可以将其去除，否则，可以人工地或自动地修正或者废弃单个模型。在过程信息与条件化之间可以存在折衷。该工作流程考虑到生成多个模型（例如，通过遍历，例如，操作14、16、18、20、22、24、26、以及/或28），其按不同级别遵守过程和条件化，以相对于约束模型结果和对模型应用的影响来估定条件化的值和过程信息。这可以包括：从共同条件化信息生成多个基于过程的模型，从共同基于过程的模型生成多个地质统计模型，以及/或遍历所述操作。响应于确定该地质统计模型可靠，方法10结束，和/或进行至下一操作（图1中未示出）。
图2例示了被配置成生成所关注的地质体积的地质统计模型的系统40。在一些实施例中，系统40包括：一个或更多个处理器42，其被配置成执行以下模块中的一个或更多个：条件化信息模块44、过程模块46、匹配标准模块48、匹配量化模块50、匹配验证模块52、统计数据模块54、统计模型模块56、模型验证模块58、以及/或其它模块。
该条件化信息模块44被配置成获取与所关注的地质体积相关联的条件化信息。该条件化信息包括从在所关注的地质体积处或其附近所进行的测量导出的信息，该条件化信息可以自动地和/或基于用户输入来获取。在一些实施例中，条件化信息模块44被配置成执行与操作12相似或相同的操作（如图1所示和在此描述的）。
该过程模块46被配置成生成所关注的地质体积的基于过程的模型。该基于过程的模型符合与所关注的地质体积相关联的条件化信息。在一些实施例中，过程模块46被配置成执行与操作14相似或相同的操作（图1所示和在此描述的）。
该匹配标准模块48被配置成获取指定基于过程的模型与和所关注的地质体积相关联的条件化信息之间的可接受相关程度或匹配的匹配标准。该匹配标准可以自动地和/或基于用户输入（例如，通过用户接口60）来获取。在一些实施例中，匹配标准模块48被配置成执行与操作16相似或相同的操作（图1所示和在此描述的）。
该匹配量化模块50被配置成量化基于过程的模型与条件化信息之间的匹配。这可以包括确定量化匹配的一个或更多个匹配度量。该匹配度量根据基于过程的模型和条件化信息来确定。在一些实施例中，匹配量化模块50被配置成执行与操作18相似或相同的操作（图1所示和在此描述的）。
该匹配验证模块52被配置成验证基于过程的模型与条件化信息之间的匹配的紧密度是否满足匹配标准（例如，通过匹配标准模块48获取的匹配标准）。该验证可以包括：确定由匹配量化模块50确定的匹配度量是否满足该匹配标准，接收对基于过程的模型与条件化信息之间的匹配的人工估定的结果（例如，经由用户接口60），和/或其它活动。在一些实施例中，匹配验证模块52被配置成执行与操作20相似或相同的操作（图1所示和在此描述的）。
该统计模块54被配置成根据基于过程的模型生成针对所关注的地质体积的统计数据。所生成的统计数据忠实于该基于过程的模型中的过程信息。该统计数据在该基于过程的模型内局部地确定，以使得给定统计数据的值作为该基于过程的模型内的位置的函数而改变。在一些实施例中，统计模块54被配置成执行与操作24相似或相同的操作（图1所示和在此描述的）。
该统计模型模块56被配置成生成所关注的地质体积的一个或更多个地质统计储集层模型。该统计模型模块56被配置成使得该地质统计模型受由统计模块54生成的统计数据约束。在一些实施例中，统计模型模块56被配置成执行与操作26相似或相同的操作（图1所示和在此描述的）。
该模型验证模块58被配置成验证通过统计模型模块56生成的地质统计模型。这可以包括检查针对基于过程的模型的符合性、针对条件化信息的符合性、条件化信息附近和/或远离其的伪像、以及和/或其它验证。在一些实施例中，模型验证模块58被配置成执行与操作28相似或相同的操作（图1所示和在此描述的）。
处理器42被配置成在系统40中提供信息处理能力。同样地，处理器42可以包括数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、以及/或用于电子地处理信息的其它机构中的一个或更多个。尽管处理器40在图2中被示出为单一实体，但这仅仅是出于例示性目的。在某些实现中，处理器42可以包括多个处理单元。处理器42可以被配置成通过软件；硬件；固件；软件、硬件、以及/或固件的某一组合；以及/或用于配置处理器42上的处理能力的其它机制来执行模块44、46、48、50、52、54、56、以及/或58。
应当清楚，尽管模块44、46、48、50、52、54、56、以及/或58在图2中被例示为共同位于单一处理单元内，但在其中处理器42包括多个处理单元的实现中，模块44、46、48、50、52、54、56、以及/或58中的一个或更多个可以相对于其它模块远程定位。下面描述的、对由不同模块44、46、48、50、52、54、56、以及/或58提供的功能的描述出于例示性目的，而非旨在进行限制，因为模块44、46、48、50、52、54、56、以及/或58中的任一个可以提供比所述的更多或更少的功能。例如，模块44、46、48、50、52、54、56、以及/或58中的一个或更多个可以消除，并且其一些或全部功能可以由模块44、46、48、50、52、54、56、以及/或58中的其它模块提供。作为另一实施例，处理器38可以被配置成执行可以执行下面归因于模块44、46、48、50、52、54、56、以及/或58中的一个的一些或全部功能的一个或更多个附加模块。
系统40可以包括可通过处理器42存取的电子存储部60。在一些实施例中，电子存储部60包括电子地存储信息的电子存储介质。电子存储部60的电子存储介质可以包括与系统40集成地（即，基本上不可去除）设置的系统存储部和/或例如经由端口（例如，USB端口、固件端口等）或驱动器（例如，盘驱动器等）可去除地连接至系统40的可去除存储部中的一个或两个。电子存储部60可以包括一个或更多个光学可读存储介质（例如，光盘等）、磁可读存储介质（例如，磁带、磁硬盘驱动器、软盘驱动器等）、基于电荷的存储介质（例如、EEPROM、RAM等）、固态存储介质（例如，闪存驱动器等），以及/或其它电子可读存储介质。电子存储部60可以包括虚拟存储资源，如经由云和/或虚拟专用网络提供的存储资源。电子存储部60可以存储软件算法、由处理器42确定的信息、经由用户接口60接收的信息、以及/或使得系统40能够适当地起作用的其它信息。电子存储部60可以是系统40内的分离组件，或者电子存储部60可以与系统40的一个或更多个其它组件（例如，处理器42）集成地设置。
用户接口60被配置成，在系统40与用户之间设置用户可以通过其向系统40提供信息和接收来系统100的信息的接口。这使得统称为“信息”的数据、结果、以及/或指令和任何其它可传送项能够在用户与系统40之间传送。适于包含在用户接口60中的接口装置的示例包括：小键盘、按钮、开关、键盘、按钮（knob）、杠杆、显示屏、触摸屏、扬声器、麦克风、指示灯、声音报警器、以及打印机。
应当明白，其它通信技术（硬布线或无线）也可以被本发明设想为用户接口60。例如，本发明设想用户接口60可以与由电子存储部62提供的可去除存储接口集成。在这个示例中，可以将信息从可去除存储部（例如，智能卡、闪存驱动器、可去除盘等）加载到系统40中，使得用户能够定制系统40的实现。适于与系统40一起使用作为用户接口60的其它示例性输入装置和技术包括但不限于：RS-232端口、RF链路、IR链路、调制解调器（电话、线缆或其它）。简单地说，用于与系统40传送信息的任何技术被本发明设想为用户接口60。
尽管基于当前所考虑的是最有用且优选的实现，出于例示的目的，对本公开的系统和/或方法进行了详细描述，但要明白的是，这种细节仅出于该目的，并且本公开不限于所公开实现，而且，正相反，其旨在覆盖处于所附权利要求书的精神和范围内的修改例和等同布置。例如，要明白的是，本公开设想，在尽可能的情况下，可以将任何实现的一个或更多个特征与任何其它实现的一个或更多个特征相组合。
Claims (11)
1.一种被配置用于生成所关注的地质体积的地质统计模型的系统，该系统包括：
一个或更多个处理器，被配置成执行计算机程序模块，所述计算机程序模块包括：
条件化信息模块，被配置成获取与所关注的地质体积相关联的条件化信息，其中，该条件化信息包括从在所关注的地质体积处或其附近进行的测量导出的信息，
过程模型模块，被配置成生成所关注的地质体积的基于过程的模型，其中，该基于过程的模型符合与所关注的地质体积相关联的条件化信息；
统计模块，被配置成从该基于过程的模型生成忠实于该基于过程的模型中的过程信息的统计数据，其中，该统计数据中的单个统计数据在该基于过程的模型内局部地确定，以使得给定统计数据的值作为该基于过程的模型内的位置的函数而改变；
统计模型模块，被配置成生成该所关注的地质体积的、受该统计模块所生成的统计数据约束的地质统计模型；
匹配量化模块，被配置成确定匹配度量，所述匹配度量量化由过程模型模块生成的基于过程的模型与和所关注的地质体积相关联的条件化信息之间的匹配的紧密度；
匹配验证模块，被配置成基于该匹配度量验证该基于过程的模型与该条件化信息之间的匹配的紧密度满足匹配阈值；
其中，该匹配量化模块被配置成使得该匹配度量在该基于过程的模型内局部地确定，以使得该匹配度量作为该基于过程的模型中的位置的函数而改变。
2.根据权利要求1所述的系统，其中，该条件化信息模块被配置成，使得与所关注的地质体积相关联的条件化信息包括井数据、钻心数据、模拟信息、专门信息、或地震数据中的一个或更多个。
3.根据权利要求1所述的系统，其中，该匹配量化模块被配置成使得该匹配度量针对包括在该基于过程的模型中的各过程和/或各时间步而分别确定。
4.根据权利要求1所述的系统，其中，该过程模型模块被配置成使得，响应于该基于过程的模型与该条件化信息之间的匹配的紧密度不满足该匹配阈值，该过程模型模块调节该基于过程的模型以更好地匹配该条件化信息，或者生成新的基于过程的模型。
5.根据权利要求1所述的系统，其中，该统计模块被配置成使得统计数据包括以下的一个或更多个：空间连续性、纹理、连通性、几何结构、性质范围或特定类型的比例、类型或性质范围之间的过渡概率、或者几何结构的相对位置或堆叠。
6.根据权利要求1所述的系统，其中，所述计算机程序模块还包括模型验证模块，该模型验证模块被配置成验证该统计模型。
7.根据权利要求6所述的系统，其中，该模型验证模块被配置成使得验证该统计模型包括：分析该统计模型以识别数值伪像和/或非预期的模型特征。
8.一种生成所关注的地质体积的地质统计模型的方法，该方法包括：
获取与所关注的地质体积相关联的条件化信息，其中，该条件化信息包括：从在所关注的地质体积处或其附近进行的测量导出的信息，生成所关注的地质体积的基于过程的模型，其中，该基于过程的模型符合与所关注的地质体积相关联的条件化信息；
从该基于过程的模型生成忠实于该基于过程的模型中的过程信息的统计数据，其中，该统计数据中的单个统计数据在该基于过程的模型内局部地确定，以使得给定统计数据的值作为该基于过程的模型内的位置的函数而改变；
生成所关注的地质体积的、受从该基于过程的模型生成的统计数据约束的地质统计模型；
确定匹配度量，所述匹配度量量化基于过程的模型与和所关注的地质体积相关联的条件化信息之间的匹配的紧密度；以及
基于该匹配度量验证该基于过程的模型与该条件化信息之间的匹配的紧密度满足匹配阈值，
其中，该匹配度量被在该基于过程的模型内局部地确定，以使得该匹配度量作为该基于过程的模型中的位置的函数而改变。
9.根据权利要求8所述的方法，其中，所获取的与所关注的地质体积相关联的条件化信息包括以下的一个或更多个：井数据、钻心数据、模拟信息、专门信息、或地震数据。
10.根据权利要求8所述的方法，其中，该局部统计数据包括以下的一个或更多个：空间连续性、纹理、连通性、几何结构、性质范围或特定类型的比例、类型或性质范围之间的过渡概率、或者几何结构的相对位置或堆叠。
11.根据权利要求8所述的方法，还包括：验证该统计模型。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US13/170,462 | 2011-06-28 | ||
| US13/170,462 US8954303B2 (en) | 2011-06-28 | 2011-06-28 | System and method for generating a geostatistical model of a geological volume of interest that is constrained by a process-based model of the geological volume of interest |
| PCT/US2012/037779 WO2013002898A1 (en) | 2011-06-28 | 2012-05-14 | System and method for generating a geostatistical model of a geological volume of interest that is constrained by a process-based model of the geological volume of interest |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103270514A CN103270514A (zh) | 2013-08-28 |
| CN103270514B true CN103270514B (zh) | 2016-10-26 |
Family
ID=47391457
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201280004208.XA Active CN103270514B (zh) | 2011-06-28 | 2012-05-14 | 用于生成所关注的地质体积的、受所关注的地质体积的基于过程的模型约束的地质统计模型的系统和方法 |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US8954303B2 (zh) |
| EP (1) | EP2727035B1 (zh) |
| CN (1) | CN103270514B (zh) |
| AU (1) | AU2012276056B2 (zh) |
| BR (1) | BR112013010516A2 (zh) |
| CA (1) | CA2819053A1 (zh) |
| EA (1) | EA201490172A1 (zh) |
| WO (1) | WO2013002898A1 (zh) |
Families Citing this family (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2573166C2 (ru) | 2010-05-28 | 2016-01-20 | Эксонмобил Апстрим Рисерч Компани | Способ сейсмического анализа углеводородных систем |
| US20130132052A1 (en) * | 2011-11-18 | 2013-05-23 | Chevron U.S.A. Inc. | System and method for assessing heterogeneity of a geologic volume of interest with process-based models and dynamic heterogeneity |
| US9798027B2 (en) | 2011-11-29 | 2017-10-24 | Exxonmobil Upstream Research Company | Method for quantitative definition of direct hydrocarbon indicators |
| US9097821B2 (en) * | 2012-01-10 | 2015-08-04 | Chevron U.S.A. Inc. | Integrated workflow or method for petrophysical rock typing in carbonates |
| EP2815255B1 (en) | 2012-02-13 | 2017-03-01 | Exxonmobil Upstream Research Company | System and method for detection and classification of seismic terminations |
| AU2013266865B2 (en) | 2012-05-23 | 2015-05-21 | Exxonmobil Upstream Research Company | Method for analysis of relevance and interdependencies in geoscience data |
| US10395215B2 (en) * | 2012-10-19 | 2019-08-27 | International Business Machines Corporation | Interpretation of statistical results |
| US10422900B2 (en) | 2012-11-02 | 2019-09-24 | Exxonmobil Upstream Research Company | Analyzing seismic data |
| US9915742B2 (en) | 2012-12-20 | 2018-03-13 | Exxonmobil Upstream Research Company | Method and system for geophysical modeling of subsurface volumes based on label propagation |
| US9348047B2 (en) | 2012-12-20 | 2016-05-24 | General Electric Company | Modeling of parallel seismic textures |
| WO2014099204A1 (en) | 2012-12-20 | 2014-06-26 | Exxonmobil Upstream Research Company | Method and system for geophysical modeling of subsurface volumes based on computed vectors |
| WO2014099200A1 (en) | 2012-12-20 | 2014-06-26 | Exxonmobil Upstream Research Company | Vector based geophysical modeling of subsurface volumes |
| US9529115B2 (en) | 2012-12-20 | 2016-12-27 | Exxonmobil Upstream Research Company | Geophysical modeling of subsurface volumes based on horizon extraction |
| US9297918B2 (en) | 2012-12-28 | 2016-03-29 | General Electric Company | Seismic data analysis |
| WO2014150580A1 (en) | 2013-03-15 | 2014-09-25 | Exxonmobil Upstream Research Company | Method for geophysical modeling of subsurface volumes |
| WO2014149344A1 (en) | 2013-03-15 | 2014-09-25 | Exxonmobil Upstream Research Company | Method and system for geophysical modeling of subsurface volumes |
| US9952340B2 (en) | 2013-03-15 | 2018-04-24 | General Electric Company | Context based geo-seismic object identification |
| US9824135B2 (en) | 2013-06-06 | 2017-11-21 | Exxonmobil Upstream Research Company | Method for decomposing complex objects into simpler components |
| US9804282B2 (en) | 2014-02-17 | 2017-10-31 | General Electric Company | Computer-assisted fault interpretation of seismic data |
| WO2016118223A1 (en) | 2015-01-22 | 2016-07-28 | Exxonmobil Upstream Research Company | Adaptive structure-oriented operator |
| US10139507B2 (en) | 2015-04-24 | 2018-11-27 | Exxonmobil Upstream Research Company | Seismic stratigraphic surface classification |
| WO2018022649A1 (en) * | 2016-07-29 | 2018-02-01 | Exxonmobil Upstream Research Company | Method and system for generating a subsurface model |
| CN106844999A (zh) * | 2017-02-13 | 2017-06-13 | 深圳泛华工程集团有限公司 | 基坑范围内坚硬岩石的体积分析方法及系统 |
| US11604909B2 (en) * | 2019-05-28 | 2023-03-14 | Chevron U.S.A. Inc. | System and method for accelerated computation of subsurface representations |
| CN113933909B (zh) * | 2021-08-30 | 2024-03-19 | 核工业北京地质研究院 | 一种快速确定裂谷盆地边界的方法 |
| GB2620556A (en) * | 2022-07-05 | 2024-01-17 | Equinor Energy As | Internal imaging and model construction |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6480790B1 (en) | 1999-10-29 | 2002-11-12 | Exxonmobil Upstream Research Company | Process for constructing three-dimensional geologic models having adjustable geologic interfaces |
| AU2005262438B2 (en) | 2004-07-01 | 2009-12-24 | Exxonmobil Upstream Research Company | Hydrodynamics-based gridding geologic modeling (Hydro-Gridding) |
| WO2007070413A1 (en) * | 2005-12-12 | 2007-06-21 | Allergan, Inc. | Treatment of corneal ulcers with topical gatifloxacin |
| US8463457B2 (en) * | 2008-06-13 | 2013-06-11 | Schlumberger Technology Corporation | Feedback control using a simulator of a subterranean structure |
| US8606555B2 (en) | 2009-10-23 | 2013-12-10 | Chevron U.S.A. Inc. | System and method for modeling a geologic volume of interest |
| US8234073B2 (en) | 2009-10-23 | 2012-07-31 | Chevron U.S.A. Inc. | System and method for estimating geological architecture of a geologic volume |
-
2011
- 2011-06-28 US US13/170,462 patent/US8954303B2/en active Active
-
2012
- 2012-05-14 EA EA201490172A patent/EA201490172A1/ru unknown
- 2012-05-14 CA CA2819053A patent/CA2819053A1/en not_active Abandoned
- 2012-05-14 BR BR112013010516A patent/BR112013010516A2/pt not_active IP Right Cessation
- 2012-05-14 CN CN201280004208.XA patent/CN103270514B/zh active Active
- 2012-05-14 EP EP12804365.0A patent/EP2727035B1/en active Active
- 2012-05-14 AU AU2012276056A patent/AU2012276056B2/en active Active
- 2012-05-14 WO PCT/US2012/037779 patent/WO2013002898A1/en active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| EP2727035A1 (en) | 2014-05-07 |
| AU2012276056A1 (en) | 2013-04-11 |
| EP2727035B1 (en) | 2017-06-21 |
| CA2819053A1 (en) | 2013-02-01 |
| EA201490172A1 (ru) | 2014-04-30 |
| BR112013010516A2 (pt) | 2016-08-02 |
| AU2012276056B2 (en) | 2015-07-09 |
| CN103270514A (zh) | 2013-08-28 |
| EP2727035A4 (en) | 2015-10-07 |
| US8954303B2 (en) | 2015-02-10 |
| US20130006591A1 (en) | 2013-01-03 |
| WO2013002898A1 (en) | 2013-01-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103270514B (zh) | 用于生成所关注的地质体积的、受所关注的地质体积的基于过程的模型约束的地质统计模型的系统和方法 | |
| US10062207B2 (en) | Method and system for reconstructing a three-dimensional model of point clouds | |
| CA2879773C (en) | Multi-level reservoir history matching | |
| US8838428B2 (en) | Methods and systems to volumetrically conceptualize hydrocarbon plays | |
| Abbaszadeh Shahri et al. | A hybrid ensemble-based automated deep learning approach to generate 3D geo-models and uncertainty analysis | |
| Tzu-hao et al. | Reservoir uncertainty quantification using probabilistic history matching workflow | |
| EP2780747A1 (en) | System and method for assessing heterogeneity of a geologic volume of interest with process-based models and dynamic heterogeneity | |
| CN105706089A (zh) | 伪相生产模拟：在储层流动模拟中经由连续相似阶跃函数相对渗透率受控模型来访问准多相流生产的信号处理方法 | |
| CN107221028A (zh) | 一种基于地震解释数据的地质体闭合曲面三维重建方法 | |
| Torrado et al. | Opening new opportunities with fast reservoir-performance evaluation under uncertainty: Brugge field case study | |
| Colombera et al. | Models for guiding and ranking well-to-well correlations of channel bodies in fluvial reservoirs | |
| Yeh et al. | A practical workflow for probabilistic history matching and forecast uncertainty quantification: Application to a Deepwater West Africa reservoir | |
| Zhao et al. | An efficient model to estimate the soil profile and stratigraphic uncertainty quantification | |
| Zhang et al. | An initial guess for the Levenberg–Marquardt algorithm for conditioning a stochastic channel to pressure data | |
| Giraud et al. | Utilisation of probabilistic MT inversions to constrain magnetic data inversion: proof-of-concept and field application | |
| US20100017131A1 (en) | Method for evaluating measured electromagnetic data relating to a subsurface region | |
| He et al. | IOR pilot evaluation in a brown-field fractured reservoir using data analytics of reservoir simulation results | |
| Demyanov et al. | Reservoir modelling with feature selection: a kernel learning approach | |
| CN103299214B (zh) | 改进地球模型、速度模型、图像体积准确度的方法和系统 | |
| Giraud et al. | Utilisation of probabilistic magnetotelluric modelling to constrain magnetic data inversion: proof-of-concept and field application | |
| Yang | Holistic strategies for prediction uncertainty quantification of contaminant transport and reservoir production in field cases | |
| CN108876916A (zh) | 辫状河训练图像生成方法及设备 | |
| Yu et al. | Machine Learning Empowered Well Placement in a Large Unconventional Tight Gas Field in China | |
| Julio et al. | Spatial constraints for the stochastic modeling of fault networks in the presence of large structural uncertainties | |
| Aly et al. | Integration of Production Data Analysis in Ensemble History Matching |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |