CN105934689A - 地质特征拆分 - Google Patents
地质特征拆分 Download PDFInfo
- Publication number
- CN105934689A CN105934689A CN201480074060.6A CN201480074060A CN105934689A CN 105934689 A CN105934689 A CN 105934689A CN 201480074060 A CN201480074060 A CN 201480074060A CN 105934689 A CN105934689 A CN 105934689A
- Authority
- CN
- China
- Prior art keywords
- fracture
- convex closure
- instruction
- model
- split
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 claims abstract description 254
- 230000008569 process Effects 0.000 claims description 59
- 238000005194 fractionation Methods 0.000 claims description 33
- 238000004422 calculation algorithm Methods 0.000 claims description 29
- 230000015572 biosynthetic process Effects 0.000 claims description 20
- 238000005755 formation reaction Methods 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 3
- 206010017076 Fracture Diseases 0.000 description 176
- 208000010392 Bone Fractures Diseases 0.000 description 169
- 238000003325 tomography Methods 0.000 description 27
- 238000004458 analytical method Methods 0.000 description 22
- 239000010410 layer Substances 0.000 description 19
- 239000012530 fluid Substances 0.000 description 16
- 238000003860 storage Methods 0.000 description 16
- 241000220317 Rosa Species 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 13
- 238000004088 simulation Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 10
- 150000003839 salts Chemical class 0.000 description 9
- 230000011218 segmentation Effects 0.000 description 9
- 208000006670 Multiple fractures Diseases 0.000 description 8
- 239000003795 chemical substances by application Substances 0.000 description 7
- 238000007726 management method Methods 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 230000009471 action Effects 0.000 description 6
- 238000011049 filling Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000000151 deposition Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005553 drilling Methods 0.000 description 5
- 230000035699 permeability Effects 0.000 description 5
- 241000193935 Araneus diadematus Species 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 230000008021 deposition Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000000429 assembly Methods 0.000 description 2
- 230000000712 assembly Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 238000001311 chemical methods and process Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 235000001674 Agaricus brunnescens Nutrition 0.000 description 1
- 241000239290 Araneae Species 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 206010010149 Complicated fracture Diseases 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 241000257303 Hymenoptera Species 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 238000010796 Steam-assisted gravity drainage Methods 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 238000010923 batch production Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 239000012792 core layer Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 238000010422 painting Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000013396 workstream Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V99/00—Subject matter not provided for in other groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
- G01V2210/663—Modeling production-induced effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics And Detection Of Objects (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种方法可以包括:接收定义三维地下结构的信息;将三维地下结构拆分成各部分;针对各部分生成凸包;以及至少部分地基于凸包来生成离散裂隙网络。
Description
相关申请的交叉引用
本申请要求享有于2013年11月25日提交的第61/908,469号美国临时专利申请、于2013年11月25日提交的第61/908,480号美国临时专利申请以及于2014年11月24日提交的序列号为14/551,667的美国专利申请的权益和优先权,通过引用将其全部内容并入本文。
背景技术
可以处理和解释数据,例如以更好地理解地表下岩石的成分、流体内容、范围和几何形状。
发明内容
根据一些实施例,执行一种方法,其包括:接收定义三维地下结构的信息;将三维地下结构拆分成部分;生成所述部分的凸包;以及至少部分地基于凸包来生成离散裂隙网络。
在一些实施例中,方法的一方面包括:至少部分地基于包括离散裂隙网络的模型来模拟与地下地层相关联的现象。
在一些实施例中,方法的一方面包括:通过应用格雷厄姆扫描算法来生成凸包。
在一些实施例中,方法的一方面包括:通过实现角拆分参数进行拆分。
在一些实施例中,方法的一方面包括:拆分生成线段,其中,例如,每个线段包括相应的端点。
在一些实施例中,方法的一方面包括:拆分曲线,其中,曲线通过地下结构与平面的交叉来确定;以及例如沿着坐标方向重复针对平面的多个位置的拆分。
在一些实施例中,方法的一方面包括:用平面交叉三维地下结构来进行拆分以生成形成曲线的交叉点,然后将曲线拆分成部分,其中,例如,拆分可以包括用线段来表示每个部分,并且其中,例如,每个线段与相应的区域相关联。
在一些实施例中,方法的一方面包括:通过确定最佳拟合平面来生成至少一个凸包。
在一些实施例中,方法的一方面包括:输出离散裂隙网络作为至少部分地由凸包的集合表示的模型。
在一些实施例中,方法的一方面包括:将至少一部分地下结构镶嵌到多边形中;计算每个多边形的法线;以及基于多边形的法线来确定一个或多个倾角值。
根据一些实施例,提供一种计算系统,其包括至少一个处理器、至少一个存储器以及存储在至少一个存储器中的一个或多个程序,其中所述一个或多个程序被配置为由一个或多个处理器来执行,所述一个或多个程序包括用于以下的指令:接收定义三维地下结构的信息;将三维地下结构拆分成部分;生成所述部分的凸包;以及至少部分地基于凸包来生成离散裂隙网络。
在一些实施例中,计算系统的一方面包括至少部分地基于包括离散裂隙网络的模型来模拟与地下地层相关联的现象的指令。
在一些实施例中,计算系统的一方面包括至少部分地通过应用格雷厄姆扫描算法、针对地下结构的至少一部分生成凸包的指令。
根据一些实施例,提供一种计算机可读存储介质,该介质具有一个或多个程序的集合,所述一个或多个程序的集合包括指令,所述指令在被计算系统执行时使计算系统:接收定义三维地下结构的信息;将三维地下结构拆分成部分;生成所述部分的凸包;以及至少部分地基于凸包来生成离散裂隙网络。
在一些实施例中,计算机可读存储介质的一方面包括至少部分地基于包含离散裂隙网络的模型来模拟与地下地层相关联的现象的指令。
在一些实施例中,计算机可读存储介质的一方面包括通过应用格雷厄姆扫描算法、针对地下结构的至少一部分生成凸包的指令。
提供该概述以引入将在下面的具体描述中进一步说明的构思的选集。该概述不打算标识所要求保护的主题的关键或主要特征,也不打算用于帮助限制所要求保护的主题的范围。
附图说明
通过以下结合附图进行的描述,将更容易地理解所描述的实现方式的特征和优点。
图1图示包括用于对地质环境进行建模的各种组件的示例性系统;
图2图示地层的示例、关于倾角的规范的示例、数据获取的示例以及系统的示例;
图3图示方法的示例;
图4图示方法的示例;
图5图示方法的示例;
图6图示数据和经处理的数据的示例;
图7图示方法的示例;
图8图示格雷厄姆(Graham)扫描技术的示例;
图9图示方法的示例;
图10示出多边形集合的示例;
图11图示图表的示例;
图12图示方法的示例;
图13图示方法的示例;
图14图示方法的示例;
图15图示方法的示例;
图16图示方法的示例;
图17图示结构的示例;
图18图示系统的示例列、模块的示例以及裂隙网络的示例;
图19图示环境的示例;以及
图20图示系统和连网系统的示例性组件。
具体实施方式
以下描述包括用于实现所描述的实现方式的目前所想到的最佳模式。该描述不应以限制性的意义来领会,而是仅用于描述实现方式的一般原理的目的。所描述的实现方式的范围应当参照所提出的权利要求来确定。
作为示例,地下环境中的结构可以通过数据的获取和所获取的数据的处理来更好地理解。作为示例,数据可以是地震数据、影像数据和/或一个或多个其他类型的数据中的一种或多种。
作为示例,可以分析数据以揭示断裂、断层或其他结构。作为示例,可以分析数据以揭示包括交叉的多个断裂的裂隙网络。作为示例,裂隙网络的断裂可以在岩石例如响应于与板块构造活动相关联的力而被受压或者拉伸时形成。例如可以使用人工技术,例如通过水压断裂,来形成断裂。作为示例,裂隙网络可以包括自然断裂、人为断裂或者自然断裂和人为断断裂。作为示例,裂隙网络可以促成地层存储(例如经由多孔性等)和/或流体流(例如经由渗透性或可传性)。
作为示例,裂隙网络可以存在于所谓的非常规地层。作为示例,一个或多个因素可以确定地层是否可以被分类为非常规。例如,多孔性、渗透性、流体圈闭机理或其他特性中的一个或多个可以是对地层进行分类的分类方案的因素。作为示例,页岩气地层可以被分类为非常规地层。页岩储集层(reservoir)的质量可以取决于例如厚度和范围、有机质含量、热成熟度、深度和压力、流体饱和度、渗透性和/或一个或多个其他因素中的一个或多个。
如所提及的那样,裂隙网络可以包括自然断裂、人为断裂或者自然断裂和人为断裂。例如,包括页岩的地层中的裂隙网络可以包括自然断裂、人为断裂或者自然断裂和人为断裂。作为示例,关于自然断裂和/或存在的人为断裂的先验知识可以有助于一个或多个操作(例如钻井、注入、提取、断裂等)的计划、执行等。
地层建模将有助于理解例如关于地层可以如何表现(例如静态和/或动态表现)。作为示例,建模可以包括生成离散裂隙网络(DFN),其可以通过多孔性、渗透性、可传性等在空间上进行表征。作为示例,DFN的一个或多个特性可以随着时间而改变(例如,考虑时间模型)。作为示例,DFN可以帮助理解地层中的存储和/或流体流的方面。例如,动态模型可以包括DFN,其中流体的移动可以至少部分地经由DFN进行。作为示例,DFN可以包括断裂、断层或者断裂和断层。
作为示例,模型可以是包括定义模型元件的顶点或节点的多维模型,例如,其可以是多边形、多面体等。作为示例,包括多边形的模型可以服从一个或多个规则,诸如例如:(a)多边形是凸起的;(b)多边形与其他多边形连接;以及(c)多边形近似于结构并且具有地质学意义。例如,关于断裂,可能期望具有包括近似于断裂的多边形的模型,其可以是裂隙网络的一部分。
作为示例,在地层包括相对平坦的断裂的情况下,包括多边形的模型可能能够包括以将平坦断裂表示为例如在DFN中可以是单独断裂的表面的方式放置的多边形。然而,在断裂可能不同于这样的相对平坦的断裂的情况下,作为示例,对断裂进行建模的工作可能违反一个或多个规则。例如,可能存在断裂或多个断裂以不正确的位置、不正确的朝向等被建模的风险。作为示例,单个最佳拟合平面可能不足以充分地表示断裂(例如,取决于断裂的三维性质)。
作为示例,方法可以包括应用拆分算法以将断裂表面的基于数据的表示拆分成各部分。在这样的示例中,然后可以例如通过平面来近似至少一些部分。例如,可以将断裂表面拆分成部分,其中每个部分通过平面来近似。作为示例,邻近的平面可以表示断裂的邻近的部分,并且可以例如参照地点或多个地点以不同的角度来布置。
作为示例,方法可以包括将断裂和/或断层转换成凸包的集合。凸包可以根据点的集合以多个维度来定义,使得凸包是包含点的凸起集合的交集。
作为示例,可以采用所谓的格雷厄姆扫描来确定点的有限平坦集合的凸包。例如,给定平面中的点的有限集合,方法可以包括确定点的有限集合的凸包。格雷厄姆扫描可以确定点的集合中的哪些点可以是凸包的极值点,其可以定义凸包(例如,平面的边界)。作为示例,方法可以包括一个或多个过程,诸如:(i)找到平面中在点的有限集合的凸包的内部的点P;(ii)以极坐标来表示点的有限集合中的每个点,其中,原点在点P处,并且近似为空的极坐标θ在从点P开始的任意固定的射线的方向上;(iii)根据增加的θk来排序点的有限集合中的元素ρkexpiθk;(iv)基于一个或多个标准来删除点以得到点的简约集;以及(v)处理点的简约集以得到点的原始有限集合的子集,其中子集由凸包的所谓的极值点构成。
作为示例,给定点的集合,格雷厄姆扫描算法可以计算凸包,例如通过:找到极值点,其用作凸包上的支点,其可以被选择为具有最大坐标值(例如在坐标尺度上)的点;以围绕支点的增加角度的次序对点进行排序以得到星状多边形(例如,或者其一部分),其中支点可以“看到”星状多边形的分段和点(例如凸性条件);以及通过围绕星状多边形(例如,或者其一部分)行进、并且在向左转时添加边而在向右转时折返(例如,根据坐标系统的定义反之亦然)来建立凸包。
作为示例,断裂或断层可以由曲线、线、点等来指示。作为示例,断裂或断层可以根据角度标准来拆分。例如,可以实现角度以参照一个或多个标准来分析断裂或断层。在这样的示例中,在存在满足和/或超出角度标准的值的改变的情况下,例如,可以将断裂或断层的一部队发拆分成两个子部分。作为示例,处理可以导致多个子部分,其中每个子部分可以邻近于至少一个其他子部分。作为示例,一个或多个补片可以存在于例如两个子部分之间。
作为示例,子部分可以通过点的有限集合来定义。在这样的示例中,可以应用算法以生成凸包。作为示例,可以处理子部分以找到每个子部分的最佳拟合平面,其中可以应用算法以使用最佳拟合平面来生成凸包。
作为示例,例如作为工作流的一部分,可以使用凸包的集合来至少部分地形成DFN。作为示例,处理可以实现为从原始(raw)处理的地震数据的数据生成一个或多个属性立方体。作为示例,可以将处理实施为(例如,可选地,以自动化方式和/或半自动化方式)从一个或多个属性立方体来提取特征(例如断裂、断层等)。作为示例,工作流可以包括使用经处理的数据,可选地使用进一步地处理,来解释一个或多个结构。作为示例,工作流可以包括执行DFN变换以将断裂变换成DFN模型。作为示例,方法可以包括使用DFN模型来执行一个或多个模拟。作为示例,方法可以包括至少部分地基于模拟结果来执行一个或多个动作。
图1示出包括各种管理组件110以管理地质环境150(例如包括沉积盆地、储集层151、一个或多个断裂153等的环境)的各个方面的系统100。例如,管理组件110可以允许直接或间接地管理关于地质环境150的感测、钻井、注入、提取等。进而,关于地质环境150的另外的信息可以变得可作为反馈160可用(例如,可选地,作为对管理组件110中的一个或多个的输入)。
在图1的示例中,管理组件110包括地震数据组件112、附加信息组件114(例如井/测井数据)、处理组件116、模拟组件120、属性组件130、分析/可视化组件142以及工作流组件144。在操作时,针对每个组件112和114提供的地震数据和其他信息可以输入给模拟组件120。
在示例性实施例中,模拟组件120可以依赖于实体122。实体122可以包括诸如井、表面、储集层等地球实体或地质对象。在系统100中,实体122可以包括为了模拟的目的而被重构的实际物理实体的虚拟表示。实体122可以基于经由感测、观察等获取的数据(例如地震数据112和其他信息114)的实体。实体可以通过一个或多个属性来表征(例如,地球模型的几何柱状网格实体可以通过多孔性属性来表征)。这样的属性可以表示一个或多个测量(例如,所获取的数据)、计算等。
在示例性实施例中,模拟组件120可以依赖于诸如面向对象的框架这样的软件框架。在这样的框架中,实体可以包括基于预先定义的类的实体以便于建模和模拟。面向对象的框架的商业可用的示例是.NETTM框架(华盛顿,雷蒙德),其提供一组可扩展的对象类。在.NETTM框架中,对象类封装可重用的代码和相关联的数据结构的模块。对象类可以被用于实例化对象实例,以便在程序、脚本等中使用。例如,钻孔类可以定义用于基于井数据来表示钻孔的对象。
在图1的示例中,模拟组件120可以处理信息以符合由属性组件130(其可以包括属性库)指定的一个或多个属性。这样的处理可以在给模拟组件120的输入之前进行(例如,考虑处理组件116)。作为示例,模拟组件120可以基于由属性组件130指定的一个或多个属性对输入信息执行操作。在示例性实施例中,模拟组件120可以构造地质环境150的一个或多个模型,其可以取决于模拟地质环境150的表现(例如,响应于一个或多个活动,无论是自然的还是人为的)。在图1的示例中,分析/可视化组件142可以允许与模型或者基于模型的结果进行交互。作为示例,来自模拟组件120的输出可以按照工作流组件144所指示的那样输入给一个或多个其他工作流。
作为示例,模拟组件120可以包括诸如ECLIPSETM储集层模拟器(Schlumberger有限公司,休斯顿,德克萨斯州)、INTERSECTTM储集层模拟器(Schlumberger有限公司,休斯顿,德克萨斯州)等模拟器的一个或多个特征。作为示例,可以根据一个或多个增强的恢复技术来模拟储集层或多个储集层(例如,考虑诸如SAGD等热处理)。
在示例性实施例中,管理组件110可以包括诸如地震模拟软件框架(Schlumberger有限公司,休斯顿,德克萨斯州)这样的商业可用的模拟框架的特征。框架提供允许探测和开发操作的最优化的组件。框架包括地震模拟软件组件,其可以通过提高资产队伍生产力来输出例如供增加储集层性能使用的信息。通过使用这样的框架,各种专家(例如地球物理学、地质学家、储集层工程师)可以开发协作工作流,并且将操作集成到流水线处理。这样的框架可以被视为应用,并且可以被视为数据驱动的应用(例如,其中为了模拟地质环境的目的而输入数据)。
在示例性环境中,管理组件110的各个方面可以包括根据框架环境的规格来操作的附件或插件。例如,作为框架环境(Schlumberger有限公司,休斯顿,德克萨斯)出售的商业可用的框架环境允许将附件(或插件)集成到框架工作流中。框架环境利用工具(Microsoft公司,雷蒙德,华盛顿)并且提供用于高效开发的稳定、用户友好的接口。在示例性实施例中,各种组件可以被实现为符合并且根据框架环境的规格(例如,根据应用编程接口(API)规格等)操作的附件(或插件)。
图1还示出框架170的示例,框架170包括模型模拟层180以及框架服务层190、框架核心层195和模块层175。框架170可以包括商业可用的框架,其中模型模拟层180是容纳框架应用的商业可用的以模型为中心的软件包。在示例性实施例中,软件可以被视为数据驱动应用。软件可以包括用于模型建立和可视化的框架。这样的模型可以包括一个或多个网格(grid)。
模型模拟层180可以提供作为数据源184的域对象182,用于呈现186以及提供各种用户接口188。呈现186可以提供图形环境,其中,应用可以显示它们的数据,用户接口188可以提供应用用户接口组件的公共外观和感觉。
在图1的示例中,域对象182可以包括实体对象、特性对象,并且可选地,可以包括其他对象。实体对象可以被用于几何地表示井、表面、储集层等,同时特性对象可以被用于提供特性值以及数据版本和显示参数。例如,实体对象可以表示井,其中特性对象提供测井信息以及版本信息和显示信息(例如以显示井作为模型的一部分)。
在图1的示例中,可以将数据存储在一个或多个数据源(例如数据仓库、一般地物理数据存储设备)中,其可以处于相同或不同的物理场所处,并且可经由一个或多个网络来访问。模型模拟层180可以被配置为对项目进行建模。因此,可以存储具体的项目,其中所存储的项目信息可以包括输入、模型、结果和案例。因此,当完成建模会话时,用户可以存储项目。在稍后的时间,项目可以使用模型模拟层180来访问和恢复,这可以重新创建相关域对象的实例。
在图1的示例中,地质环境150可以包含多层(例如,层理),其包括储集层151并且可能被断层153横断的。作为示例,地质环境150可以配备各种传感器、检测器、致动器中的任何一种。例如,装备152可以包括通信电路以与一个或多个网络155接收和传送信息。这样的信息可以包括与井下装备154(其可以是用以获取信息的装备)相关联的信息,以帮助资源恢复等。其他装备156可以位于远离井地点之处,并且可以包括感测、检测、发射或其他电路。这样的装备可以包括储存器和通信电路以存储和传达数据、指令等。作为示例,可以提供一个或多个卫星用于通信、数据获取等目的。例如,图1示出与可以被配置用于通信的网络155通信的卫星,注意,卫星可以附加地或者替代地包括用于影像(例如,空间的、光谱的、时间的、辐射测量的,等等)的电路。
图1还将地质环境150示出为可选地包括与井相关联的装置157和158,井包括可以与一个或多个断裂159交叉的基本上水平的部分。例如,考虑可能包括自然断裂、人为断裂(例如水压断裂)或者自然和人为断裂的结合的页岩地层中的井。作为示例,可以针对横向延伸的储集层来钻井。在这样的示例中,在属性、压力等方面的横向变化可能存在,其中这样的变换的评估可以帮助计划、操作等以开发储集层(例如经由断裂、注入、提取等)。作为示例,装备157和/或158可以包括用于断裂、地震感测、地震数据的分析、一个或多个断裂的评估的组件、一系统、多个系统等。
如所提及的那样,系统100可以被用于执行一个或多个工作流。工作流可以是包括多个工作步骤的处理。工作步骤可以操作数据,例如以创建新的数据、更新现有数据等。作为示例,例如可以基于一个或多个算法来操作一个或多个输入并且创建一个或多个结果。作为示例,系统可以包括工作流编辑器,以用于工作流的创建、编辑、执行等。在这样的示例中,工作流编辑器可以提供对一个或多个预先定义的工作步骤、一个或多个定制工作步骤等的选择。作为示例,工作流可以是在软件中可实现的工作流,例如,其可以操作地震数据、地震属性等。作为示例,工作流可以是可在框架中实现的处理。作为示例,工作流可以包括访问诸如插件(例如外部可执行的代码等)这样的模块的一个或多个工作步骤。
图2示出地层201的示例、钻孔的示例、倾角的规范(convention)215的示例、数据获取处理220的示例以及系统250的示例。
如所示那样,地层201包括水平表面和各种子表面层。作为示例,钻孔可以是垂直的。作为另外的示例,钻孔可以是偏离的。在图2的示例中,钻孔210可以被视为垂直钻孔,例如,其中z轴向下延伸正交于地层201的水平表面。
关于倾角的规范215,平面的三个维度的朝向可以通过其倾角和走向来定义。倾角是平面在特定方向上以垂直平面测量的距水平平面(例如,影像平面)的平面的倾斜的角度。倾角可以通过量值(例如,也被称为角度或量)和方位角(例如,也被称为方向)来定义。如在图2的规范215中所示,各种角度φ指示例如距影像水平平面(例如扁平的上表面)的向下的倾斜的角度;而方位角指倾斜平面倾斜沿着的方向(例如,其可以参照度数、指南针方向等给出)。在图2的规范中示出的另外的特征是走向(strike),其是由倾斜平面和水平平面(例如,考虑平坦的上表面作为影像水平平面)相交所创建的线的朝向。
与倾角和走向有关的一些另外的术语可以应用于分析,例如取决于境况、所收集的数据的朝向等。一个术语是“真正倾角(true dip)”(例如参见图2的规范215中的DipT)。真正倾角是垂直于走向直接测量的平面的倾角(例如参见方向朝北并且被标记为“走向”和角度α90的线),并且也是倾角量值的最大可能值。另一个术语是“视在倾角(apparent dip)”(例如参见图2的规范215中的DipA)。视在倾角可以是在除了真正倾角的方向之外的任何其他方向上测量的平面的倾角(例如参见φA,作为角度α的DipA);然而,视在倾角可能等于真正倾角(例如参见φ,作为相对于走向的角度α90的DipA=DipT)。换句话说,在使用(例如,在方法、分析、算法等中)术语视在倾角的情况下,对于具体的倾斜平面,“视在倾角”的值可以等于该具体倾斜平面的真正倾角。
如图2中的规范215所示,如在恰好垂直于走向的横截面中所见到的平面的倾角是真正倾角(例如,具有φ作为相对于走向的角度α90的DipA=DipT的表面)。如所示那样,在任何其他方向上在横截面中观察到的倾角为视在倾角(例如,参见被标记有DipA的表面)。另外,如图2的规范215所示,视在倾角可以近似地为0度(例如,平行于水平表面,其中切割平面的边缘沿着走向方向行进)。
就井孔中的观察倾角而言,在垂直钻孔的井中观察到真正倾角。在任何其他方向(或偏离)上钻孔的井中,所观察到的倾角是视在倾角(例如,其有时被称为相对倾角)。为了确定在这样的钻孔中观察到的平面的真正倾角值,作为示例,可以对一个或多个视在倾角值应用向量计算(例如,基于钻孔偏离)。
如所提及的那样,在根据钻孔图像沉积学解释使用的另一术语是“相对倾角”(例如DipR)。从位于非常平静的环境中沉积的岩石中的钻孔图像测量的真正倾角的值可以从砂岩体中的倾角减去(例如使用向量减法)。在这样的示例中,所得到的倾角被称为相对倾角,并且在解释砂岩体朝向中使用。
关于分析、解释、属性等(例如参见图1的系统100的各个框)可以使用诸如规范215这样的规范。作为示例,可以部分地通过倾角来描述各种类型的特征(例如沉积性层理、断层和断裂、单斜脊、火成岩墙和岩床、变质叶理等)。
地震解释可以旨在至少部分地基于一个或多个倾角参数(例如角度或幅度、方位角等)对一个或多个地表下边界进行标识和分类。作为示例,可以至少部分地通过角度、至少部分地通过方位角等来描述各种类型的特征(例如沉积性层理、断层和断裂、单斜脊、火成炎墙和岩床、变质叶理等)。
如在图2的图表220中所示,地质体(geobody)225可以存在于地质环境中。例如,地质体225可以是盐丘。盐丘可以是由盐构成的蘑菇形状或插座形状的底辟(diapir),并且可以具有上覆盖岩(例如,或者冠岩)。盐丘可以在被埋在其他类型的沉淀物下面时由于盐的相对浮力而形成。在盐丘处或者附近,由于因在相关蒸发矿物密封(association evaporite mineral sealing)中的盐移动而形成闭圈(trap),可以发现碳氢化合物。浮力差可以使盐开始垂直地流动(例如作为盐枕),其可能造成断层。在图表220中,地质点225与多个层相遇,其每一个可以通过倾角角度φ来标识。
作为示例,可以针对区域以轨迹的形式获取地震数据。在图2的示例中,图标220示出获取装备222,其发射来自来源(例如传送器)的能量并且经由沿着主测线(inline)方向排成一列的一个或多个传感器(例如接收器)接收所反射的能量。因为该区域包括层223和地质体225,由获取装备的传送器发射的能量可以反射远离层223和地质体225。这样的反射的证据可以在所获取的轨迹中找到。关于轨迹226的一部分,所接收的能量可以通过以一采样率工作的模数转换器来离散化。例如,获取装置222可以将传感器Q所感测到的能量信号转换成以近似每4ms一个样本的速率的数字样本。给出在介质或多个介质中的声音的速度,可以将样本率转换成近似的距离。例如,声音在岩石中的速度可以是每秒大约5km的量级。因此,近似4ms的样本时间间隔将对应于大约10米的样本“深度”间隔(例如,假设从源到边界以及边界到传感器的路径长度)。作为示例,轨迹可以是大约4秒的持续时间;因此,对于大约4ms时间间隔一个样本的采样率,这样的轨迹将包括大约1000个样本,其中较晚获取的样本对应于较深的反射边界。如果将前述示例的4秒轨迹持续时间除以2(例如,以考虑反射),则对于垂直对齐的源和传感器,最深的边界深度可以被估计为大约10km(例如,假设每秒大约5km的声音速度)。
在图2的示例中,系统250包括一个或多个信息存储设备252、一个或多个计算机254、一个或多个网络260以及一个或多个模块270。关于一个或多个计算机254,每个计算机可以包括一个或多个处理器(例如处理内核)256以及用于存储指令的存储器258(例如模块),指令例如可以由一个或多个处理器中的至少一个执行。作为示例,计算机可以包括一个或多个网络接口(例如有线的或无线的)、一个或多个图形卡、显示接口(例如有线的或无线的)等。
在图2的示例中,所述一个或多个存储设备252可以存储在长度和宽度上跨越几千米并且在深度上例如大约10km的地质环境的地震数据。地震数据可以按照表面网格(例如根据主测线和十字交叉线方向定义)来获取。例如,给定大约40米乘大线40米的网格块,40km乘40km的领域可以包括大约一百万个轨迹。这样的轨迹可以被视为3D地震数据,其中时间近似于深度。作为示例,计算机可以包括网络接口,其用于经由网络访问存储在一个或多个存储设备252中的地震数据。进而,计算机可以经由指令来处理所访问的地震数据,指令可以是以一个或多个模块的形式。
作为示例,一个或多个属性模块可以提供用于处理地震数据。作为示例,属性可以包括几何属性(例如倾角角度、方位角、连续性、地震轨迹等)。这样的属性可以是结构属性库(例如图1的属性组件130)的一部分。结构属性可以帮助边缘检测、地震反射器的局部朝向和倾角、地震事件的连续性(例如,平行于所估计的层理朝向)等。作为示例,边缘可以被定义为地震数据内的水平幅度连续性中的不连续性,并且可以对应于断层、断裂等。几何属性可以是空间属性,并且依赖于多个轨迹。
如所提及的那样,作为示例,一区域的地震数据可以包括一百万个轨迹,其中每个轨迹包括一千个样本,总共十亿个样本。按照如今的标准以实时的方式处理这样的地震数据中所涉及的资源可能也是相当可观的。作为示例,倾角扫描方法可以应用于地震数据,其涉及参照离散平面(例如由离散平面界定的体积)来处理地震数据。取决于地震数据的大小,这样的方法可能涉及实时地处理的相当可观的资源。这样的方法可以查看轨迹与它们的幅度之间的局部相干性,并且因此可以以“视在倾角”的类别来分类。
作为示例,可以处理诸如表面影像(例如,卫星的、地质的、地球物理学的,等等)这样的影像。作为示例,方法可以使用拆分技术来分析影像,并且可以包括生成凸包的集合,凸包例如可以表示一个或多个结构(例如断裂的部分或断层的部分等)。作为示例,框架可以访问表面影像,并且可以访问地表下地震数据并且生成表面结构和地表下结构(其可以经由内插处理或者其他处理来结合)的三维表示(例如用于可视化)。例如,可以通过地震学并用通过卫星影像和至少部分地基于对地震数据和表面影像的分析而构造的模型来捕获潜在的结构。
作为示例,可以执行蚂蚁跟踪(ant-tracking)作为工作流的一部分。作为示例,蚂蚁跟踪可以生成蚂蚁跟踪属性、蚂蚁跟踪表面、蚂蚁跟踪体积(例如,或者立方体)等。蚂蚁跟踪可以包括使用通过模拟来涉及“蚂蚁”找到蚂蚁的巢穴和蚂蚁的食物来源之间最短路径(例如,通过使用信息素来吸引其他蚂蚁进行通信)的算法。在这样的示例中,最短路径比长的路径变成用更多的信息素来标记,使得随后的蚂蚁更有可能选取最短路径,以此类推。
在特征例如由于噪声、获取印迹等而可能是潜在(例如潜在的结构)的情况下,在蚂蚁跟踪之前执行分析可以增强跟踪潜在特征的能力,特别是在特征具有一定量的连续性(例如在表面、体积等内是连续的)的情况下。例如,由断裂处理(例如考虑水压断裂)所生成的断裂可以倾向于相对较小(例如相比于断层)并且是连续的。
作为示例,可以使用一个或多个边缘检测算法来处理数据。各种边缘检测算法可以包括确定梯度(例如数据集中的值的空间偏离)。应用边缘检测算法可以帮助揭示并强调地下环境中的结构,例如用于特定目的。例如,工作流可以旨在确定某些结构是否存在于地下环境中以及那些结构是否以相对于其他结构的某种关系统而存在。
在各种示例实施例中,可以对诸如地震数据、从地震数据推导的数据等数据应用一个或多个分析。作为示例,方法可以包括执行一个或多个分析以检测诸如例如断裂、其他潜在结构等特征。作为示例,可以在框架中将分析实现为模块、模块的集合等,例如以检测断层、断裂以及潜在的反射。作为示例,可能执行一个或多个分析以帮助检测油气探测和生产(E&P)中感兴趣的一个或多个特征。例如,来自分析的结果可以帮助良好的布局、地质建模、岩床分析、断裂区域和断裂走廊的检测以及针对非常规资源和碳酸盐领域(例如考虑页岩领域)的E&P。
作为示例,断裂走廊或纤细断层可以引起地震信号,其可能在所获取的地震数据中呈现为例如横截面中的小幅度的自不相干特征以及切片或地震表面上的轮廓(lineament)。对这样的特征的检测可以包括处理地震信号、地震数据或者两者以生成一个或多个边缘检测属性,例如,其中属性可以被视为地震数据的可测量的“特性”(例如,考虑幅度、倾角、频率、相位、极性等)。例如,属性可以是从地震信号、地震数据等推导的值或一组值,并且可以参照坐标系(例如一维的、二维的、三维的、四维的或者甚至更高维度的)来定义。作为示例,维度可以是空间维度、时间维度、频率维度等。作为示例,考虑将地震数据提供为“立方体”,其中立方体中的每个体元(体积元素)具有值。在这样的示例中,边缘检测算法可以处理立方体中的值以生成新的值,其中新的值一起被称为边缘检测属性(例如属性立方体)。
作为示例,可以处理地震立方体(例如地震体积或体积的地震数据)以生成属性立方体(例如属性体积或体积的属性值)。作为另一个示例,可以处理地震表面以生成属性表面。作为又一个示例,可以处理地震线以生成属性线。作为示例,可以处理地震点以生成属性点。
属性可以例如在一个时刻、针对多个时刻、在一时间窗口上等推导、测量等,并且例如可以在单个轨迹上、在一组轨迹上、在从地震数据解释出的表面上等测量。属性分析可以包括评估例如关于储集层的各种参数,考虑从幅度随偏移变化(AVO)分析(例如,或者幅度随角度(AVA)等)推导的碳氢化合物指标。
作为示例,模型可以是地震不连续性平面模型(例如“SDP”模型),其可以从地震数据(例如原始数据、一个或多个属性等)推导。作为示例,方法可以包括使用地震不连续性平面(SDP)来执行离散裂隙网络(DFN)转换以生成DFN模型。作为示例,可以实现一种方法,其中例如可以通过点的列表(例如,不同深度的)来表示复杂的断裂表面。
作为示例,可以从地震属性立方体或多个立方体提取一个或多个SDP。在这样的示例中,SDP可以“物理地”(例如,在空间上)跟随从地震数据推导的包括信号信息的信息,信号信息可能成Z字形,可能是波状的,等等。因此,SDP可以成Z字形,可以是波状的,等等。关于DFN模型,其可以近似自然和/或所引入的断裂的解释。例如,DFN模型可以包括符合一个或多个数字建模技术的一族凸平面多边形。DFN模型可以包括表示断裂的“补片”,其中个体断裂可以通过一个或多个凸平面多边形来表示。作为示例,方法可以包括将SDP模型(例如具有波状表面等)变换成DFN模型(例如具有凸多边形表面)。
作为示例,变换可以包括接收多面性(multi-faceted)SDP模型,其中这样的面可以朝向不同。在这样的示例中,方法可以包括逐个面地拆分SDP模型,例如至少部分地基于面朝向。作为示例,方法可以旨在生成拆分的多面的集合或多个集合,其维持断裂密度和朝向与被拆分的SDP模型的密度和朝向相对一致。这样的处理可以例如出于质量控制的目的来实现。例如,分析可以将拆分的各面与一个或多个SDP相比较(例如使用一个或多个空间标准)。在这样的示例中,可以调节一个或多个拆分参数以实现一个或多个SDP与从中拆分的各面之间的期望级别的一致性。这样的处理例如可以允许剪裁至少部分地通过从SDP模型拆分出的各面形成的DFN模型。
作为示例,DFN模型可以包括符合一个或多个连接性标准的补片。例如,考虑规定在DFN模型中的邻近补片应彼此相连(例如,其中它们在拆分之前属于同一SDP的两个邻近的部分)的连接性标准。
作为示例,可以分析裂隙网络以生成可以形成一个或多个SDP模型多个SDP。作为示例,SDP模型可以包括数千个SDP的量级(例如在裂隙网络的SDP模型中可以包括可能十万或者更多的SDP)。关于将SDP模型变换成DFN模型,作为示例,考虑包括要被拆分的形成DFN模型的补片的100,000个SDP的模型。在这样的示例中,如果个体SDP花费大约一秒来处理,则100,000个SDP可能总共花费100,000秒(例如多于大约10天)。
作为示例,方法可以包括将SDP模型变换成适合于作为DFN模型补片使用的凸多边形。这样的方法可以包括可以串行和/或并行地实现的一个或多个算法。作为示例,方法可以包括一个或多个分支,例如,以使用一种或多种技术,分析SDP或其一部分,然后决定是否处理SDP或其一部分。在这样的示例中,SDP或其一部分可以被施以一个或多个拆分处理,其可以根据计算资源的命令而不同。
作为示例,方法可以决定对SDP实现最佳似合平面处理以生成适合于用作DFN模型中的补片的凸平面多边形。作为示例,方法可以决定针对SDP或其一部分实现不同的处理,其中例如,储集层可以是断裂的、非常规的储集层(例如,使得多个多边形近似于SDP)。
作为示例,方法可以包括多边形标准,例如针对每个SDP考虑最大多边形数量,针对每个SDP考虑最大多边形密度,等等。例如,针对每个SDP的最大多边形数量可以是针对每个SDP的十个或数十个多边形的量级。作为示例,方法可以包括可以被应用于从将SDP拆分成多边形的处理得出的多边形的玫瑰图(rose diagram)分析。例如,玫瑰图分析可以示出(例如裂隙网络的)倾角和方位角分布。作为示例,可以至少部分地基于玫瑰图分析来调节一个或多个多边形标准。
图3示出方法310的示例,方法310包括:接收块314,接收定义三维地下结构的信息;拆分块318,将地下结构拆分成各部分;生成块322,生成各部分的凸包;生成块326,至少部分地基于凸包来生成离散裂隙网络(DFN);以及模拟块330,至少部分地基于包括离散裂隙网络(DFN)的模型来模拟与地下地层相关联的现象。
在图3中与各个计算机可读介质(CRM)块315、319、323、327和331相关联地示出的方法310。这样的块一般包括适合于由一个或多个处理器(或处理器内核)执行以指示计算设备或系统执行一个或多个动作的指令。虽然示出各种块,但是可以允许至少部分地执行方法300的各个动作的指令来配置单个介质。作为示例,计算机可读介质(CRM)可以是计算机可读存储介质。作为示例,计算机可读介质可以是非临时性的,并且不是载波。
作为示例,方法310或其一部分可以可选地实现为框架的插件。例如,考虑可以生成凸包(例如针对图3的生成块322)的框架的插件。在这样的示例中,框架可以包括可以至少部分地基于这样的凸包来建立至少一部分DFN的一个或多个模块。
作为示例,方法可以包括将复杂的断裂转换成离散裂隙网络(DFN),其可以是以适合于断层或断裂建模(例如非常规储集层或者其他储集层)的数据格式。
作为示例,断层和断裂模型(例如由DFN或多个DFN表示)可以帮助探测和开采储集层或多个储集层。作为示例,地震模拟框架可以被配置为输入DFN或多个DFN。
作为示例,储集层可以包括多个断裂。例如,储集层可以包括数以万计的断裂的量级。作为示列,包括拆分的方法可以被实现为将断裂表面拆分成各个部分,其中每个部分例如可以通过最佳拟合平面来近似。作为示例,拆分算法可以包括可调节的角度参数,其可以具有默认值,可以经由学习算法来确定,可以经由用户界面(例如图形用户界面等)来设置。
作为示例,拆分方法假设可以基于角度来拆分复杂表面。在这样的示例中,在所选择的方向(例如水平方向)上大于该角度的倾向改变可以要求表面的拆分。作为示例,对于相对平滑的曲面,从拆分得出的各部分可以保持彼此相连。作为示例,可以采用格雷厄姆扫描算法来确定每个部分的凸包(例如,基于拆分的表面的每个部分的最佳拟合平面)。可以提供凸包的集合作为形成DFN模型的基础。因此,作为示例,方法可以包括将复杂断裂表面或多个表面变换成可以被用于形成DFN模型的至少一部分的凸平面多边形。
图4示出方法410的示例以及信息411的近似例示,信息411定义地下结构、空间415中的点、至少部分地经由空间中的点来定义的区域419以及至少部分经由与区域419相关联的空间中的点来定义的凸包423。如所示那样,方法410包括:接收块414,接收定义三维地下结构的信息(例如考虑可以是SDP等的信息411);拆分块418,将三维地下结构的至少一部分拆分成区域(例如考虑区域419等);以及生成块422,生成区域的凸包(例如,考虑区域419的凸包423)。作为示例,凸包可以被视为一种类型的代理(proxy)(例如精细代理)。作为示例,代理可以通过其顶点来定义。例如,作为凸五边形的代理可以通过空间中的五个点来定义。因此,通过诸如点415这样的空间中的多个点定义的区域可以由减少数量的信息来表示(例如,凸多边形的空间中的顶点等)。
在图4的示例中,可以使用可以至少部分地定义一个或多个区域的分区处理来拆分所接收的信息。作为示例,可以将代理定义为表示地下结构(例如,可以是裂隙网络的一部分的断裂)的一部分的区域。这样的代理可以是诸如例如地震数据、属性数据等信息的“轻量级”表示(例如来自一个或多个数据获取技术的空间数据填充物)。作为示例,可以对地震数进行处理、解释等以形成例如空间中的点的集合。作为示例,代理可以在空间上(例如在多维空间中)通过一起定义凸包(例如凸平面多边形)的点的集合来表示。
作为示例,代理可以至少部分地基于空间中的点的集合。在这样示例中,点的集合可以投射到区域内,例如可以将点的集合投射到平面(例如平面区域等)上。作为示例,可以生成凸包的生成处理可能处理与区域(例如平面区域等)相关联的点的集合。作为示例,方法可以包括例如至少部分地基于空间中的点(例如经由对点进行分组等)来近似各个区域的自动处理信息。例如,在图4中,点415可以被用于定义一个或多个区域(例如区域419及其邻域等)。作为示例,代理可以是这样的区域,或者例如代理可能使用到一组点的最佳拟合平面来确定(例如经由最小平方或者其他类型的拟合技术)。关于如何定义代理,作为示例,方法可以包括评估对空间中的点的区域的拟合(例如参照一个或多个标准)。作为示例,评估可以考虑参照一个或多个邻近区域和/或代理(例如关于连接性等)的拟合。
作为示例,方法可以包括接收SDP以及经由区域分区处理对SDP进行分区。这样的方法还可以包括在区域上进行循环,其中例如针对个体区域,处理可以精细调谐区域,其可以被视为代理。例如,考虑包括将在区域内的点投射到平面上以及找到在平面内的凸多边形的处理。这样的处理例如可以包括帮助确保连接性的分析。例如,这样的分析可以旨在连接邻近的凸多边形(例如,以坚持DFN模型的补片标准)。
图5示出方法510的示例,方法510包括:接收块512,接收信息(例如SDP、SDP模型等);标识块516,至少部分地基于一部分信息来标识结构;以及拆分块520,拆分所标识的结构。在图5的示例中,拆分块520可以在一个或多个方向上执行拆分。例如,水平拆分块522可以经由评估相对于水平线、多个水平线、参考水平线或多个线等的结构的分析来提供拆分;而垂直拆分块524可以经由评估相对于垂直的井、垂直的结构、参考垂直线或多个线等的结构的分析来提供拆分。
在图5的示例中,方法510包括:开始块530,开始可以处理从拆分处理(例如按照拆分块520)得出的多个个体区域的循环。如所示那样,方法510包括:决定块534,决定是否已经处理了(例如从拆分结构得出的N个区域的)最后的区域。在决定块534决定已经处理了最后的区域的情况下,方法510可以继续到终止块542(例如或者其他块、工作流等)。在决定块534决定尚未处理最后的区域的情况下,方法510继续到确定块538,以便确定关于区域(例如正在被处理的当前区域)的代理,其可以包括根据可以将代理连接到邻近区域的代理的一个或多个连接性标准来确定代理。
在图5的示例中,方法510包括决定块542,以便决定代理是否实现关于作为所标识的结构(例如按照标识块516)的一部分的代理的充分拟合。例如,出于评估代理多好地拟合结构的目的,可以提供一个或多个标准。这样的一个或多个标准可以包括空间标准,可选地至少部分地基于一个或多个统计地计算出的值(例如均值、标准离差等)。作为示例,决定块542可以包括按照一个或多个连接性标准来评估拟合。例如,在先前的代理邻近于当前的代理而存在的情况下,决定块542可以确定当前的代理是否“连接”到先前的代理。这样的确定可以包括参照当前代理的边缘(例如关于相符、重叠等)来分析先前代理的边缘。
如所示那样,在决定块542决定代理拟合“很好”的情况下,方法510可继续到确定块550,以便确定当前代理上的凸包。作为示例,可以经由能够标识在其内可能存在其他点的“最外面的”点的格雷厄姆扫描技术来分析点。在这样的示例中,最外面的点可以定义凸包(例如平面的凸多边形)。如所示那样,在决定块542决定代理拟合不是“很好”的情况下,方法510可继续到确定块546,其可以确定“最佳”拟合平面(例如当前代理的最佳拟合平面)的凸包并且还可以包括填充的插入,例如用于符合一个或多个连接性标准的目的。
在图5的示例中,在已经确定出凸包或者凸包和一个或多个插入填充的情况下,方法510可以继续到添加块554,以便将至少一个凸包添加到结果(例如作为第一凸包或者作为随后的凸包,可选地具有一个或多个填充)。方法510然后可以通过选择另一个区域来处理而迭代。例如,如所示那样,可以包括选择块558以便选择要处理的下一个区域。方法510因此可以包括具有一个或多个路径的循环,例如,取决于是否满足一个或多个拟合标准。在图5的示例中,一个路径可以比另一个路径有更多的计算需求。作为示例,对一个路径或者另一个路径的代理的比率可以取决于接收块512所接收的信息、按照标识块516标识的结构、标识块516所应用的标识技术、按照拆分块520的拆分技术、一个或多个拟合标准等。作为示例,诸如多边形界限参数、多边形密度参数等参数(例如拆分处理等的)可以以可以对被取得用作代理的路径具有效果的方式来选择和/或调节。
在图5的示例中,接收块512可以包括接收SDP。作为示例,方法510可以包括输出(例如在按照终止块542终止时)凸平面多边形,其可以被用作例如DFN模型的补片。方法510可以包括逐个地自动拆分SDP集合(例如SDP模型)中的一个或多个SDP。作为示例,方法可以至少部分并行地实现。例如,在SDP模型包括多个SDP的情况下,可以并行地处理至少一部分SDP。作为示例,在“补片”被连接的情况下,可以串行地处理结构,例如以相继生成被添加到一起(可选地通过填充以帮助满足一个或多个连接性标准)的凸平面多边形的方式。
图6示出包括凸包的集合的数据630和经处理的数据650的示例。作为示例,凸包的集合可以被用于至少部分地形成DFN模型,例如其中裂隙网络的各个部分被表示为凸平面多边形。作为示例,方法可以包括接收数据以及处理数据以生成DFN(例如DFN模型)。
图7示出方法710的示例,方法710包括:提供平面(例如I、J、K坐标或指数系统的I、J平面)中的表面,其可以表示为曲线;以及基于角度标准(参见块714)来拆分曲线。拆分可以包括例如如块718所示那样地形成线段。在这样的示例中,平面中的表面的曲线可以通过分段来表示,其中每个分段包括端点。
作为示例,拆分可以在一个或多个其他平面(例如,I、J、K坐标和指数系统的)中进行。例如,拆分可以包括水平拆分和垂直拆分(例如在I、J平面中的拆分以及在K坐标方向上延伸的平面中的折分)。在拆分在多个平面中进行的情况下,作为示例,平面可以是正交的。拆分的结果可以是由点、线等定义的区域(例如在多维度空间中)。
图7还示出具有点和拆分角度(θ)的三个分段的示例性图720。如所示那样,近邻计数可以是关于可以参照诸如例如SDP的信息在空间上定义的点的参数或度量。关于拆分角度,其可以按照例如可以是对点的拟合的线来定义。作为示例,拆分角度可以相对于参考来定义。例如,给定一地点,第一条线可以定义对该地点的第一正切,而第二条线可以定义对该地点的的第二正切。在这样的示例中,第一和第二正切可以定义它们之间可以与拆分角度进行比较的角度。在这样的示例中,在角度小于拆分角度的情况下,第一和第二线可以被视为区域的公共线的一部分;而在角度满足或超出拆分角度的情况下,第一和第二线可以被视为是两个不同区域(例如两个邻近的区域)的部分。
关于对曲线的选择,因为断裂的交叉线的形状和在深度上的切片可能在结构上是相似的,所以方法可以包括选择最长的切片来表示断裂的形状(例如,如块714和718中那样的顶视图)。在这样的示例中,最长的切片的每个分段可以通过最佳拟合线分段Si来估计,其可以被写成以下形式:
Si:y=(tanαi)x+b,xi<x<xi+1,i=0,1,2...(1)
在这样的示例列中,可以给出拆分角度θ(例如由用户等),并且可以根据下面的等式将最长的切片拆分成不同的线段:
|αi-αj|>θ (2)
在这样的示例中,可以将断裂拆分成若干区域(例如从顶视图)。作为示例,可以存储用bi来标记的每个二等分平面(其可以是两个区域为一对的边界)。
作为示例,一个或多个垂直间隙存在于区域内的网孔的不同部分之间,网孔可以被进一步拆分成若干隔离的网孔。
作为示例,作为DFN变换处理的一部分,方法可以开始于多个区域中的一个区域(例如Area_0),并且逐个地处理每个区域。作为示例,对于Area_0,如果网孔可以通过单个最佳拟合平面c0来近似,则可能存在b0和c0的交叉线,其可以被称为连接线并且被标记为coni。否则,该区域内的网孔可以分别通过它们自己的最佳拟合平面来近似。作为示例,在确定区域中的网孔的代理之后,每个断裂的边界上的点可以被确定并输入给格雷厄姆扫描算法以生成相应的凸包。在这样的示例中,每个断裂网孔可以通过其自己的代理上的凸多边形pi,j来表示。
作为示例,在考虑Areai(i>0)之前,如果存在con0,则方法可以帮助确保最终的凸多边形的两个点在con0上。作为示例,如果Areai内的网孔可以通过单个最佳拟合平面ci来近似,则可以形成bi和ci的交叉线coni。例如,考虑可以在coni-1上找到两个点,并且可以在coni上找到另外的两个点。关于其他的边界点,也可以使用格雷厄姆扫描算法来找到凸包。在这样的示例中,因为在二等分平面的两侧上的网孔的代理凸多边形({pi-1,j}和{pi,j})包括同一连接线上的点,所以拆分的网孔可以通过连接线来连接(例如,图6的经处理的数据650和图9的包954和955)。如果Areai内的网孔未通过单个最佳拟合平面良好地近似,则该区域内的网孔可以分别通过它们自己的最佳拟合平面来近似。作为示例,每个区域中的网孔可以以这样的方式来变换。
作为示例,在处理区域之后,方法可以输出凸多边形集合{pi,j},其可以被用于至少部分地形成DFN模型。这样的模型可以被用于一个或多个目的。例如,模拟器可以接收DFN模型作为输入,以用于模拟DFN模型的一个或多个断裂中的流体流的目的(例如,这可以通过补片或多个补片来表示)。
作为示例,所接收的信息可以包括点。例如,SDP可以根据包括点的数据结构(例如可以被编号的点的列表等)来指定。这样的点可以使用适当的坐标系(例如笛卡尔坐标系等)在空间中定义。作为示例,方法可以包括水平地和/或垂直地扫描一个或多个点列表以检测一个或多个拆分的点。如图7所示,(例如SDP的)点列表中的点可以处于近似的深度。在扫描点列表期间,方法可以采用多个参数,例如考虑拆分角度参数和近邻计数参数(例如参见图7的图表720)。作为示例,拆分的角度可以可选地为可以经由图形用户接口、批处理文件等设置和/或选择的参数。作为示例,拆分处理可以确定在倾向的改变大于由拆分角度参数所给出的角度值时进行拆分。在这样的示例中,到拆分的任何一侧的点(例如拆分位置或拆分点)可以被分到两个不同的区域中。作为示例,近邻计数参数值可以取决于所接收的信息。例如,近邻计数可以至少部分地取决于SDP的一个或多个特性。作为示例,在包括多个SPD的SDP模型中,近邻计数可以针对SPD中的至少一个而不同。作为示例,算法可以自适应地计算近邻计数参数值,其例如可以取决于SDP大小、SDP的标准离差等之中的一个或多个。
图8示出对应于格雷厄姆扫描技术的图示810,例如给定点的集合,可以采用格雷厄姆扫描技术来确定凸包。在图8的示例中,例示被编号为从0至12的点的集合,在那些点之中,点0、1、3、10和12可以在空间上定义凸包的顶点(例如凸多边形)。在图8的示例中,凸包包括五个顶点和五个边(例如五边形),注意针对给定的点的集合,可以确定其他类型的包。
作为示例,可以使用格雷厄姆扫描技术来找到“2D”面上的有限点集的凸包。作为示例,n维度的点的集合S的凸包可能被定义为包括S的凸起集合的交集。作为示例,可以通过属性定义凸包。例如,考虑以下:凸包是凸多边形;凸包平面;以及凸包包括在代理上的投射点(例如,对应于诸如SDP数据等数据的代理)。
作为示例,可以根据可以确定凸包的算法来采用一个或多个连接性标准。例如,标准可以起作用以包括作为与当前的代理(例如当前的区域)的交叉线上的顶点的(例如代理或区域的)先前确定的凸包的点。在这样的示例中,算法可以开始以作为部分定义要确定的凸包的集合的一部分的那些点。对于第一代理或区域,这种标准不可以应用(例如,除非提供交叉线等)。作为示例,方法可以采用格雷厄姆扫描技术,其中至少部分地基于要邻近于要确定的凸包的凸包的边界来预先确定定义凸包的至少一个点。在这样的示例中,可以确保存在的凸包与“要确定的”凸包之间的至少一定量的连接性。例如,在图8的例示810中,点P3和P10是预先存在的凸多边形的顶点。
图9示出包括数据块930和经处理的数据块950的方法910。如在数据块930中所示,例如,可以使用一个或多个角度标准以各种深度(例如沿着深度坐标)来处理数据。这样的方法可以将结构拆分成分段,例如,其中每个分段包括端点。在这样的示例中,在每个深度处,可以经由拆分平面(例如,水平面)中的结构的线、曲线等来生成分段。在图9的示例中,可以形成区域(例如参见图7的块714),并且可以形成分段(例如参见图7的块718)。
如图9中的块950所示,在块930中生成的信息可以被用于按照格雷厄姆扫描或者可以找出凸包的其他类型的算法来确定极值点。如所示那样,由块930的数据所指示的表面可以被拆分成8个凸包951、952、953、954、955、956、957和958。这8个凸包一起可以被用于形成DFN模型的一部分。在这样的示例中,DFN模型可以遵循一个或多个规则,例如关于多边形。
图10示出包括具有间隙1034的多边形1031和多边形1033的多边形集合1030以及包括至少部分地位于间隙1032中的多边形1052的多边形集合1050。在这样的示例中,多边形1052可以被插入以满足关于多边形1033相对于多边形1031的一个或多个连接性标准。如参照图5的方法510的提及地,一个或多个拟合标准(例如,其可以包括一个或多个连接性标准)可能不被满足,例如,其中代理未充分地“匹配”所接收的数据集(例如SDP数据)的区域。在这样的示例中,间隙可以存在于两个邻近的多边形之间(例如,由于一定量的失配)。
作为示例,在代理未足够良好地拟合SDP的情况下,方法可以包括丢弃该代理以及例如确定最佳拟合平面。例如,在多边形集合1030中,多边形1031可以是代理,并且多边形1033可以是最佳拟合平面。在这样的示例中,由多边形1033所表示的最佳拟合平面未充分地连接到多边形1031。例如,在多边形1031与多边形1033之间未共用连接线(例如,连接性可能丧失)。如参照多边形集合1050所例示的那样,可以在多边形1031与多边形1033之间插入多边形1052,使得存在连接性。作为示例,所插入的多边形可以是具有适当尺度和朝向的三角形以提供连接性。
作为示例,标准可以是到代理的点的经归一化的平均距离,其可以被定义为参数“d”。在这样的示例列中,可以设置阈值,其中,例如,如果d大于阈值,则方法可以考虑代理拟合不充分。在这样的示例中,方法可以确定最佳拟合平面,并且使用该最佳拟合平面作为对代理的替代。因为最佳拟合平面可能不确保到相邻的多边形的连接性,所以方法可以包括插入一个或多个三角形作为填充。作为示例,阈值可以由用户例如经由图形用户接口、批处理文件来定义。作为示例,阈值可以取决于地质信息、SDP、SDP模型、DFN、DFN的一部分等之中的一个或多个特性。
图11示出图表1110、1130、1152和1154的一系列示例。作为示例,可以分析多边形,例如以评估地下区域。例如,关于倾角角度,可以分析经由诸如图4的方法410、图5的方法510等方法生成的多边形。例如,图表1110示出关于倾角角度的样本的数量(例如多边形的数量)。如所示那样,存在样本的数量一般随着倾角角度接近90度而增加的分布。
作为示例,多边形可以参照走向来分析,并且例如可以经由玫瑰图来进行绘制。在图11中,图表1130示出超过40,000个样本的走向信息,其中,如可以经由图表1110所确保的那样,样本的一大部分跨越玫瑰图图表1130、从大约90度至大约270度地。
图11的图表1152和1154分别对应于近似地南北延伸的断裂和近似地东西延伸的断裂。如所示那样,南北断裂倾向于在长度上短于东西断裂。
作为示例,方法可以包括至少部分地基于多边形数据(例如凸包数据)来执行一个或多个统计分析。作为示例,统计可以针对断裂倾角、朝向和断裂长度中的一个或多个来计算。如在图11中所示,断裂可以被分类为其中一部分具有近似地东西朝向的高角度。作为示例,方法可以至少部分地基于多边形来计算断裂的统计分布,其可以揭示诸如存在一部分交替定向的断裂这样的附加信息。例如,在来自现场的观察可以指示近似东西断裂的同时,统计分析可以揭示背斜或其他特征。例如,背斜可以通过在规模上可能小于近似东西断裂的一部分近似南北断裂来指示(例如,参见图表1152和1154)。
图12示出方法1200的示例,方法1200包括:接收块1210,接收地震和/或属性信息;以及生成块1230,生成地震不连续性平面(SDP)模型。在这样的示例中,方法1200可以可选地包括:接收块1240,接收诸如例如可以至少部分地由生成块1230用于生成SDP模型的日志信息这样的一个或多个其他类型的信息。
如图12的示例所示,方法1200包括:确定块1250,至少部分地基于SDP模型来确定凸包。例如,SDP模型可以包括空间中的点,其中这样的点可以被分析以确定凸包。例如,考虑图4的方法410,其与空间中的点的近似图形表示一起被例示。作为示例,确定块1250可以包括填充,例如用一个或多个多边形来填充凸包之间的一个或多个区域,以帮助在适当时确保邻近的凸包的情况下的连接性。
方法1200还包括:生成模块1270,至少部分地基于凸包的一部分(例如,可选地,基于填充,如果存在的话)来生成离散裂隙网络(DFN)模型。在这样的示例中,DFN模型可以被输入给模拟器,例如按照模拟块1290,以便模拟地质环境中的物理现象。作为示例,这样的模拟器可以是能够求解可以指示一个或多个流体如何在地质环境中流动的压力值、饱和度值等的储集层模拟器。
作为示例,方法可以包括确定地下结构的倾角。作为示例,地下结构可以是断层、多个断层、断裂、多个断裂或者其他结构。作为示例,结构可以是或者包括裂隙网络(例如包括交叉的多个断裂)。
作为示例,倾角可以被定义为距离水平面的平面的倾斜的幅度。作为示例,可以垂直于走向来测量真正(或最大)倾角。作为示例,可以在垂直于走向之外的方向上测量视在倾角。作为示例,倾角可以被定义为诸如沉积床、断层、断裂等平面特征与平面(例如水平面)之间的角度。
作为示例,倾角(例如倾角信息)可以便于分析地下地层信息。例如,倾角信息可以便于断层和/或断裂分析。在结构可以表示为平面表面的情况下,倾角可以基于该平面表面相对于另外的表面(例如水平表面)的角度。
作为示例,在结构包括曲线或者一个或多个弯曲部分(例如,或者其他形状等)的情况下,倾角确定可能更具有挑战性。作为示例,考虑被分类为非常规储集层的储集层。在这样的示例中,可以进行获取地震数据(例如,可以表示为地震数据集中的值的地震信号)的地震调查或多个调查。在这样的示例中,从可以表示断裂的地震数据中提取出的表面的形状可能比例如单个平面表面更复杂。
作为示例,分析可以包括在玫瑰图上绘制信息(例如,参见图11的玫瑰图表1130)。例如,玫瑰图可以例示一个或多个断裂走向朝向,并且可以可视地指示可能存在多于一个的断裂集合。作为示例,玫瑰图可以例示倾角信息(例如,考虑倾角方位角方向;注意,在走向信息和倾角信息之间可能存在各种关系)。
作为示例,方法可以包括以可以减少对可能导致低倾角估计的数据的影响的方式来确定复杂断裂的倾角。例如,断裂可能被期望具有相对高的倾角值。例如,断裂可能被期望具有不接近于0的角度(例如,其中0表示水平面)。
作为示例,方法可以包括至少部分地基于表示复杂表面(例如,或者复杂表面的一部分)的局部三角形网孔中的平均单位法线向量来确定复杂表面的倾角。例如,这样的方法可以包括复杂表面的三角剖分。作为示例,三角剖分处理可以是包括使用可能在形状不同以填充区域(例如表面)的三角形的镶嵌处理。在这样的示例中,可以计算出三角形网孔中的每个三角形的单位法线向量。作为示例,可以例如使用每个三角形的面积作为权重,对单位法线向量进行求和。方法然后可以基于平均单位法线向量来确定倾角(例如,复杂表面的倾角可以根据平均单位法线来估计)。
作为示例,镶嵌处理可以包括使用一个或多个类型的多边形。在这样的示例中,每个多边形可以包括相关联的信息,从中可以确定单位法线和面积。在这样的示例中,在可以确定平均法线(例如平均单位法线)的处理中,可以应用面积作为法线的权重。
作为示例,采用镶嵌用于表面的多边形填充以及法线和面积计算的方法对于可能在表面中存在的变化可以是健壮的。例如,这样的方法可以在一些部分可能不同于表面的适合的、直实倾角的情况下是健壮的。
图13示出方法1310的示例,方法1310包括:接收块1314,接收定义三维地下结构的信息;镶嵌块1318,将地下结构镶嵌到三角形中(例如,至少部分地基于所提供的信息);计算块1322,计算每个三角形的法线;计算块1326,计算每个三角形的面积;以及确定块1330,基于三角形的法线和面积来确定地下结构的倾角。作为示例,这样的方法可以包括至少部分地基于倾角来模拟与地下地层相关联的现象。
在图13中示出了与各种计算机可读介质(CRM)块1315、1319、1323、1327和1331相关联的方法1310。这样的块一般包括适合于由一个或多个处理器(或处理器内核)执行以指示计算设备或系统执行一个或多个动作的指令。虽然示出各种块,但是可以通过允许至少部分地执行方法1310的各个动作的指令来配置单个介质。作为示例,计算机可读介质(CRM)可以是计算机可计存储介质。
作为示例,方法1310或其一部分可以可选地实现为框架的插件。例如,考虑可以确定倾角信息(例如,按照图13的确定块1330)的框架的插件。
作为示例,诸如图13的方法1310这样的方法可以作为可以包括诸如图3的方法310这样的方法的至少一部分的工作流的一部分来执行。例如,在按照镶嵌块1316执行镶嵌时,可以经由分区处理和/或拆分处理(例如,参见图3的方法310的拆分块318)来分析所镶嵌的形状(例如多边形)。例如,可以组合多个所镶嵌的形状以定义区域,其可以被处理以定义代理。作为示例,可以分析来自诸如图3的方法310这样的方法的信息以及来自诸如图13的方法1310这样的方法的信息。例如,考虑图11的图表1110、1130、1152和1154。具体地,考虑示出样本数量与倾角角度的图表1110。作为示例,可以针对倾角角度来分析凸包,并且可以针对倾角角度来分析所镶嵌的多边形。作为示例,例如可以使用玫瑰图(例如,参见图111的玫瑰图1130),针对凸包和/或所镶嵌的多边形来绘制朝向。
作为示例列,方法可以包括执行工作流,其中,工作流包括复杂表面的三角剖分、每个三角形的单位法线的计算、每个三角形的面积的计算、以面积作为权重的单位法线的求和、计算垂直于加权的平均法线的平面的倾角以及在随后的处理(例如绘制玫瑰图、模拟现象等)中使用所计算的倾角。
作为示例,可以将复杂表面组织成若干点列表(例如,以不同深度)。在这样的示例中,方法可以包括在两个邻近的点列表中的点上进行循环以形成三角形,其中,在对复杂表面的各个邻近点列表执行这样的操作之后,可以完成对复杂表面的三角剖分。
作为示例,方法可以包括计算三角形网孔中的各个三角形的法线(例如,单位法线)。作为示例,给定三角形,其边缘中的两个可以被标记为和并且三角形的法线可以例如使用和的叉积来确定:
在这样的示例中,法线可以被归一化成单位向量。在这样的示例中,单位法线向量可以具有例如可以参照坐标系定义的方向。
作为示例,方法可以包括计算三角形网孔中的各个三角形的面积。例如,给定三角形的三个顶点,可以确定其边缘的长度,其可以被表记为a、b和c。在这样的示例中,可以使用Heron公式来计算三角形的面积S:
其中p=(a+b+c)/2 (2)
在前述公式中,p可以被称为半周长。作为示例,三角形的面积可以被计算为其内径(内切圆的半径)与其半周长的乘积。作为示例,取决于三角菜的角度,等式可以包括对长度a、b和c进行排序以及使用包括基于所排序的这些值的项的公式(例如可以添加数字稳定性的公式,其中可以存在小内角)。作为示例,可以应用(例如可选地取决于镶嵌形状等)诸如Brahmagupta公式、Bretschneider公式这样的各种公式。
作为示例,例如在求和以得到区域的单位法线(例如单位法线向量)时,方法可以包括使用三角形网孔中的三角形的各个面积作为权重。
作为示例,方法可以包括带权重的求和。例如,假设三角形网孔的三角形的面积{Si}以及对应单位法线向量在这样的示例中,可以使用下面的公式来确定最终的向量其可以表示复杂表面的朝向:
在这样的示例中,复杂表面的倾角可以使用来估计。作为示例,倾角可以用在一个或多个随后的处理中,例如玫瑰花图绘制、模拟等。
图14示出包括例如如在块1414中所示那样地提供数据的方法1410的示例。这样的数据可以包括表示表面1416的数据,其中表面可以包括例如可以根据玫瑰花图1417(例如,还参见图11的玫瑰花图图表1130)表示的朝向或多个朝向。作为示例,方法1410可以包括例如使用可能在形状上(例如参照坐标系)变化的三角形来镶嵌表面1416的至少一部分。
图15示出包括例如如在块1514中所示那样地提供数据的方法1510的示例。这样的数据可以包括顶点(例如可以定义三角形的顶点)的坐标。在这样的示例中,方法1510可以包括计算各个三角形的值。例如,如在块1518中所示,方法1510可以包括计算法线和计算面积。
图16示出包括例如如在块1614中所示那样地接收数据(例如,还参见图6的数据630)的方法1610的示例。这样的数据可以包括表示三维空间中的表面1616的数据。作为示例,表面1616可以被镶嵌以形成网孔,例如可以定义诸如三角形这样的多边形的节点或顶点的集合。在这样的示例中,如在块1618中所示,对于各个多边形,方法1610可以包括计算法线和面积以及确定由各个多边形表示的区域的总法线。例如,在多边形是三角形的情况下,从1到N的多个三角形可以被考虑,其中计算出法线和面积,其进而被用于确定由三角形表示的区域的总法线。在这样的示例中,总法线可以是例如使用各个单位法线的加权和所确定的平均单位法线,其中权重是对应于各个单位法线的面积。
图17示出被表示为表面1710、1720和1730的结构的三个示例。例如,表面1710可以是薄且长的表面,表面1720可以是弯曲的表面,而表面1730可以是成Z字形的表面。作为示例,表面1710、1720和1730可以是表示断裂的复杂表面。
作为示例,诸如图13的方法1310这样的方法可以被应用于这样的表面中的一个或多个,例如以确定一个或多个对应的倾角值。作为示例,这样的方法可以提供表示断裂、断层等的倾角或倾角值。作为示例,使用最佳拟合平面来拟合表面1710的方法可能是次优的,因为表面薄且长并且可能导致垂直平面而不是水平平面。作为示例,方法对于表面1730可能是次优的,因为成Z字形的表面1730的多数点的倾角倾向于近似45度,使得平均可能不会得到用户可能期望的高倾角值。作为示例,诸如图13的方法1310这样的方法可以应用于表面1710、1720和1730中的一个或多个,并且提供在地质学上可能更有意义的对应的倾角(例如表示引起表面的结构)。
作为示例,方法可以包括接收以模型的形式的信息(例如模型信息)。作为示例,模型可以是地震不连续性平面模型(例如“SDP”模型),其可以从地震数据(例如原始数据、一个或多个属性等)推导。作为示例,方法可以包括使用地震不连续性平面(SDP)模型用于倾角计算。作为示例,方法可以被用于3D空间中的点集的倾角计算。作为示例,可以实现例如通过点集来表示复杂断裂表面的方法。
作为示例,系统可以包括一个或多个模块,其可以被用于分析数据、控制处理、执行任务、执行工作步骤、执行工作流等。
作为示例,方法可以包括:接收定义三维地下结构的信息;将地下结构镶嵌到三角形中;计算每个三角形的法线;计算每个三角形的面积;以及基于三角形的法线和面积来确定地下结构的倾角。作为示例,这样的方法可以包括至少部分地基于倾角来模拟与地下地层相关联的现象。
作为示例,地下结构可以是或者包括断裂。作为示例,断裂可以包括按照至少一个量级的幅度(例如10或0.1)而不同于统一的纵横比。作为示例,断裂可以是或者包括弯曲的断裂。作为示例,断裂可以是或者包括成Z字形的断裂。
作为示例列,方法可以包括通过计算两个向量的叉积来计算法线。作为示例,方法可以包括通过应用Heron公式来计算面积。
作为示例,方法可以包括通过计算法线的加权和来确定倾角。在这样的示例中,方法可以包括使用面积作为权重来对对应的法线进行加权。作为示例,方法可以包括对法线进行归一化(例如,对多个法线进行归一化)。
作为示例,方法可以包括至少部分地基于所确定的倾角来绘制玫瑰花图。
作为示例,计算面积可以包括计算半周长、计算内径或者计算一个或多个其他参数。
作为示例,系统可以包括处理器、操作地耦合到处理器的存储器、被存储在存储器中并且包括计算机可执行的指令的一个或多个模块,所述指令指示系统:接收定义三维地下结构的信息;将地下结构镶嵌到三角形中;计算每个三角形的法线;计算每个三角形的面积;以及至少部分地基于三角形的法线和面积来确定地下结构的倾角。作为示例,这样的系统可以包括至少部分地基于包括倾角的模型来模拟与地下地层相关联的现象的指令。作为示例,系统可以包括通过计算法线的加权和(例如使用面积作为权重)来确定倾角的指令。
作为示例,一个或多个计算机可读存储介质可以包括处理器可执行的指令,所述指令指示计算机:接收定义三维地下结构的信息;将地下结构镶嵌到三角形中;计算每个三角形的法线;计算每个三角形的面积;以及基于三角形的法线和面积来确定地下结构的倾角。
作为示例,一个或多个计算机可读存储介质可以包括计算机可执行的指令,所述指令指示计算机:至少部分地基于包括所确定的倾角的模型来模拟与地下地层相关联的现象。
作为示例,一个或多个计算机可读存储介质可以包括计算机可执行的指令,所述指令指示计算机:计算法线的加权和(例如使用面积作为权重)。
图18示出系统1800的示例、各个模块1810的示例以及裂隙网络1880(例如作为DFN模型的一部分等)的近似例示的示例。在图18的示例中,系统1800包括操作地耦合到存储器1804的一个或多个处理器1802。作为示例,存储器1804可以存储诸如模块1810中的一个或多个这样的模块,其可以提供对地下环境中的储藏、流等进行建模。在图18的示例中,模块1810包括流体储集层模块1812、干燥储集层模块1814、存在的井的模块1822、预期的井的模块1824、自然断裂模块1842、人为断裂模块1844以及一个或多个求解器模块1860。在图18的示例列中,模块1810可以包括适合于由一个或多个处理器(例如处理器内核)执行以指示计算设备或系统执行一个或多个动作的指令。例如,系统1800可以由一个或多个模块1810的指令来指示。
作为示例,方法可以包括实现一个或多个模块1810来表示诸如裂隙网络1880这样的网络。在图18的示例中,裂隙网络1880可以包括自然断裂和人为断裂。作为示例,水压断裂的创建可以受到一个或多个自然断裂的影响。例如,水压断裂成长可以在东北至西南的方向上进行,其重新激活趋势在另外的方向或多个方向(参见,例如,箭头指示水压断裂的可能的传播方向)上的自然断裂(虚线)。
作为示例,方法可以包括:例如关于储集层或多个储集层,使用模型并且使用模型来模拟环境的行为(例如储藏、流动等)而对环境中的自然断裂进行建模。进而,可以针对预期的人为断裂来分析解。这样的分析例如可以包括定位一个或多个井,以便参照一个或多个自然断裂来创建一个或多个预期的人为断裂,以生成用以重新激活自然断裂作为流体流动的管道的网络。作为示例,这样的分析可以旨在避免一定的自然断裂并且重新激活(例如,利用)其他自然断裂。在这样的示例中,对自然断裂位置、特性等的精细化可以使用模型(例如DFN模型等)、可选地结合对一个或多个储集层进行建模的3D刚性模型来进行。
作为示例,模型可以考虑可能与断裂有关的压力或者一个或多个其他因素。作为示例,自然断裂模型可以数学地链接到3D环境的压力模型。作为示例,模型可以考虑化学过程(例如,酸化)。作为示例,自然断裂模型可以数学地链接到化学反应模型,以便对化学过程进行建模(例如参照一个或多个断裂特性)。在至少部分地基于自然断裂模型的解针对流执行历史匹配的情况下,对自然断裂模型的精细化可以用以更新与压力相关联的一个或多个参数(例如方向等)。
作为示例,系统可以包括用于处理信息的一个或多个处理器、操作地耦合到一个或多个处理器的存储器以及包括可在存储器中存储并且可由一个或多个处理器中的至少一个执行的指令的模块。这样的模块可以包括:储集层模拟,经由三维网格模型对地下三维环境中的储集层进行建模;自然断裂模块,经由模型(例如,DFN)对自然断裂进行建模;井模块,经由井模型对井进行建模;以及一个或多个求解器模块,至少部分地基于经由模型对自然断裂进行建模,来求解裂隙网络中的流体流的值。作为示例,系统可以包括:人为断裂模块,经由模型对人为断裂进行建模。作为示例,系统可以包括:求解器模块,求解包括至少一个自然断裂和至少一个人为断裂的裂隙网络中的流体流的值。
如所提及的那样,可以在对自然断裂、多个自然断裂等进行建模的模型的一个或多个部分上定义(例如,施加)边界条件。图19示出包括各种地层、井眼和自然断裂的环境1910的示例。如所指示的那样,地层包括诸如油、气和/或水这样的流体,其可以定义各种区域。关于边界条件,自然断裂可以包括自然断裂到自然断裂的边界条件、自然断裂到充满油的地层的边界条件、自然断裂到井眼的边界条件、自然断裂到充满气的地层的边界条件、自然断裂到充满水的地层的边界条件等。作为示例,自然断裂可以包括多个边界条件,例如针对井眼和充满流体的地层两者。
作为示例,地层可以被视为充满流体或者是空的(例如,“干的”),取决于流体的类型。例如,相对于要产生油的目标而言,充满气的地层可以被视为关于油是空的。如图19的示例性环境1910所指示的那样,油和水可以共存于地层内,并且策略可以被阐明为产生具有最少的水容量的油。作为示例,这样的策略可以通过使用参照环境对一个或多个自然断裂进行建模的模型而改善(hone)(例如,以避免激活可能导致油中的水容量增加的自然断裂)。
作为示例,方法可以包括:接收定义三维地下结构的信息;将地下结构(例如断裂、多个断裂、断层、多个断层等)拆分成各部分;生成所述部分的凸包;以及至少部分地基于凸包来生成离散裂隙网络。这样的方法还可以包括:至少部分地基于包括离散裂隙网络的模型,来模拟与地下地层相关联的现象。
作为示例,方法可以包括:将地下结构的至少一部分镶嵌到多边形中;计算每个多边形的法线;以及基于多边形的法线来确定一个或多个倾角值。在这样的示例中,多边形可以是三角形,并且方法还可以包括:计算每个三角形的面积。作为示例,方法可以包括:至少部分地基于所镶嵌的多边形的法线,并且例如至少部分地基于所镶嵌的多边形的面积,来确定地下结构的倾角(例如倾角值)。
作为示例,方法可以包括:通过应用格雷厄姆扫描算法来生成一个或多个凸包。作为示例,方法可以包括:通过实现角度拆分参数来进行拆分。作为示例,拆分可以包括生成被地下结构交叉的至少一个平面中的区域。作为示例,拆分可以生成线段,例如,其中每个线段包括端点。
作为示例,方法可以包括拆分曲线,其中通过结构与平面的交叉来确定曲线。作为示例,方法可以包括针对沿着坐标方向的平面的多个位置来重复拆分。
作为示例,拆分可以包括:用平面交叉地下结构以生成形成曲线的交叉点;然后将曲线拆分成各部分。在这样的示例中,拆分可以包括用线段表示每个部分,其中例如每个线段具有相关联的区域。
作为示例,方法可以包括通过提供最佳拟合平面来生成凸包。作为示例,方法可以包括输出离散裂隙网络作为至少部分地通过凸包的集合表示的模型。
作为示例,系统可以包括处理器、操作地耦合到处理器的存储器、被存储在存储器中并且包括处理器可执行的指令的一个或多个模块,所述指令指示系统:接收定义三维地下结构的信息;将地下结构拆分成各部分;生成所述部分的凸包;以及至少部分地基于凸包来生成离散裂隙网络。
作为示例,这样的系统可以包括至少部分地基于包括离散裂隙网络的模型来模拟与地下地层相关联的现象的指令。作为示例,生成至少一个部分的凸包的指令可以包括应用格雷厄姆扫描算法的指令。
作为示例,一个或多个计算机可读存储介质可以包括处理器可执行的指令,所述指令指示计算机:接收定义三维地下结构的信息;将地下结构拆分成各部分;生成所述部分的凸包;以及至少部分地基于凸包来生成离散裂隙网络。
作为示例,一个或多个计算机可读存储介质可以包括处理器可执行的指令,所述指令指示计算机:至少部分地基于包括离散裂隙网络的模型来模拟与地下地层相关联的现象。作为示例,一个或多个计算机可读存储介质可以包括处理器可执行的指令,所述指令指示计算机:例如通过应用格雷厄姆扫描算法来生成拆分的结构的至少一个部分的凸包。
图20示出计算系统2000的示例以及连网的系统2010的示例的组件。系统2000包括一个或多个处理器2002、存储器和/或存储组件2004、一个或多个输入和/或输出设备2006和总线2008。在示例性实施例中,指令可以被存储在一个或多个计算机可读介质(例如存储器和/或存储组件2004)中。这样的指令可以由一个或多个处理器(例如处理器2002)经由可以是有线的或无线的通信总线(例如总线2008)来读取。一个或多个处理器可以执行这样的指令以实现(整体地或者部分地)一个或多个属性(例如作为方法的一部分)。用户可以经由I/O设备(例如设备2006)查看来自处理的输出以及与处理进行交互。在示例性实施例中,计算机可读介质可以是存储组件,诸如物理存储器存储设备,例如芯片、封装芯片、存储卡等(例如计算机可读存储介质)。
在示例性实施例中,组件可以是分布式的,诸如在网络系统2010中那样。网络系统2010包括组件2022-1、2022-2、2022-3……2022-N。例如,组件2022-1可以包括(多个)处理器2002,而(多个)组件2022-3可以包括(多个)处理器2002可访问的存储器。另外,(多个)组件2002-2可以包括I/O设备以便显示以及可选地与方法进行交互。网络可以是或者包括因特网、内联网、蜂窝网络、卫星网络等。
作为示例,设备可以是包括用于信息的通信的一个或多个网络接口的移动设备。例如,移动设备可以包括无线网络接口(例如,可经由IEEE 802.11、ETSI GSM、卫星等操作)。作为示例,移动设备可以包括诸如主处理器、存储器、显示器、显示图形电路(例如,可选地,包括触摸和手势电路)、SIM插槽、音频/视频电路、运动处理电路(例如加速度仪、陀螺仪)、无线LAN电路、智能卡电路、传送器电路、GPS电路以及电池这样的组件。作为示例,移动设备可以被配置为手机、平板等。作为示例,方法可以使用移动设备来实现(例如整体地或者部分地)。作为示例,系统可以包括一个或多个移动设备。
作为示例,系统可以是分布式环境,例如,所谓的“云”环境,其中各种设备、组件等出于数据存储、通信、计算等目的而交互。作为示例,设备或系统可以包括一个或多个组件用于经由因特网(例如,其中通信经由一个或多个因特网协议进行)、蜂窝网络、卫星网络等之中的一个或多个进行信息通信。作为示例,方法可以在分布式环境中实现(例如,整体地或者部分地作为基于云的服务)。
作为示例,信息可以从显示器(例如考虑触摸屏)输入、向显示器输出或者这两者。作为示例,信息可以输出给投影仪、激光设备、打印机等,使得信息可以被查看。作为示例,信息可以以立体画或全息术的方式输出。关于打印机,考虑2D或3D打印机。作为示例,3D打印机可以包括可以输出以进行构成3D对象的一个或多个物质。例如,数据可以提供给3D打印机以构造地下地层的3D表示。作为示例,层可以以3D(例如水平面等)、以3D构造的地质体等来构造。作为示例,洞、断裂等可以以3D构造(例如,作为正结构、作为负结构等)。
虽然在上文仅详细描述了一些示例性实施例,但是本领域技术人员将容易意识到,在示例性实施例中可能有很多修改。因此,所有这样的修改旨在被包括在由所附的权利要求所限定的本公开的范围内。在权利要求中,装置加功能的句子旨在涵盖在本文中被描述为执行所述功能的结构,并且不仅是结构等效物,而且是等效的结构。因此,虽然钉子和螺丝钉可能不是结构等效的,其中钉子采用圆柱表面以将木制部分安装在一起,而螺丝钉采用螺旋状表面,但是在紧固木制部件的环境中,钉子和螺丝钉可以是等效的结构。申请人的明确意图是启动35U.S.C.§112第6段用于对本文中的任何权利要求的任何限制,除了在权利要求书中明确地使用措词“用于…的部件”以及相关联的功能的那些之外。
Claims (20)
1.一种方法(310),包含:
接收定义三维地下结构的信息(314);
将三维地下结构拆分成各部分(318);
生成各部分的凸包(322);以及
至少部分地基于凸包来生成离散裂隙网络(326)。
2.根据权利要求1所述的方法,还包含:至少部分地基于包含离散裂隙网络的模型来模拟与地下地层相关联的现象。
3.根据权利要求1所述的方法,其中,生成凸包包含应用格雷厄姆扫描算法。
4.根据权利要求1所述的方法,其中,拆分包含:实现角拆分参数。
5.根据权利要求1所述的方法,其中,拆分生成线段。
6.根据权利要求5所述的方法,其中,每个线段包含端点。
7.根据权利要求1所述的方法,其中,拆分包含:拆分曲线,其中曲线通过地下结构与平面的交叉来确定。
8.根据权利要求7所述的方法,还包含:沿着坐标方向重复针对平面的多个位置的拆分。
9.根据权利要求1所述的方法,其中,拆分包含:用一平面交叉三维地下结构来生成用于形成曲线的交叉点,然后将曲线拆分成各部分。
10.根据权利要求9所述的方法,其中,拆分包含:用线段表示每个部分。
11.根据权利要求10所述的方法,其中,每个线段包含相关联的区域。
12.根据权利要求1所述的方法,其中,生成凸包包含:确定最佳拟合平面。
13.根据权利要求1所述的方法,其中,输出处理输出离散裂隙网络(DFN)作为至少部分地由凸包的集合表示的模型。
14.根据权利要求1所述的方法,还包含:将至少一部分地下结构镶嵌到多边形中;计算每个多边形的法线;以及基于多边形的法线来确定一个或多个倾角值。
15.一种系统(250),包含:
处理器(256);
存储器(258),操作地耦合到处理器;
一个或多个模块(270),被存储在存储器中并且包含指示系统的处理器可执行指令,其中,所述指令包含用于以下的指令:
接收定义三维地下结构的信息(315);
将三维地下结构拆分成各部分(319);
生成各部分的凸包(323);以及
至少部分地基于凸包来生成离散裂隙网络(327)。
16.根据权利要求15所述的系统,其中,被存储在存储器中并且包含指示系统的处理器可执行指令的一个或多个模块包含至少部分地基于包含离散裂隙网络的模型来模拟与地下地层相关联的现象的指令。
17.根据权利要求15所述的系统,其中,针对至少一个部分生成凸包的指令包含应用格雷厄姆扫描算法的指令。
18.一个或多个计算机可读存储介质,包含处理器可执行指令以指示计算机:
接收定义三维地下结构的信息;
将三维地下结构拆分成各部分;
生成各部分的凸包;以及
至少部分地基于凸包来生成离散裂隙网络。
19.根据权利要求18所述的一个或多个计算机可读存储介质,包含指示计算机至少部分地基于包含离散裂隙网络的模型来模拟与地下地层相关联的现象的处理器可执行指令。
20.根据权利要求18所述的一个或多个计算机可读存储介质,包含指示计算机通过应用格雷厄姆扫描算法来生成至少一个部分的凸包的处理器可执行指令。
Applications Claiming Priority (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361908480P | 2013-11-25 | 2013-11-25 | |
US201361908469P | 2013-11-25 | 2013-11-25 | |
US61/908,469 | 2013-11-25 | ||
US61/908,480 | 2013-11-25 | ||
US14/551,667 | 2014-11-24 | ||
US14/551,667 US11125912B2 (en) | 2013-11-25 | 2014-11-24 | Geologic feature splitting |
PCT/US2014/067391 WO2015077768A1 (en) | 2013-11-25 | 2014-11-25 | Geologic feature splitting |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105934689A true CN105934689A (zh) | 2016-09-07 |
Family
ID=53180297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480074060.6A Pending CN105934689A (zh) | 2013-11-25 | 2014-11-25 | 地质特征拆分 |
Country Status (4)
Country | Link |
---|---|
US (1) | US11125912B2 (zh) |
EP (1) | EP3074957B1 (zh) |
CN (1) | CN105934689A (zh) |
WO (1) | WO2015077768A1 (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11250182B2 (en) | 2015-11-10 | 2022-02-15 | Landmark Graphics Corporation | Target object simulation using orbit propagation |
US10388065B2 (en) * | 2015-11-10 | 2019-08-20 | Landmark Graphics Corporation | Fracture network triangle mesh adjustment |
EP3475854A4 (en) * | 2016-06-22 | 2020-01-22 | Services Petroliers Schlumberger | DISPLAYS OF TANK SIMULATIONS WITH FRACTURE NETWORKS |
WO2018022030A1 (en) * | 2016-07-27 | 2018-02-01 | Schlumberger Technology Corporation | System automation tools |
CA3037276A1 (en) | 2016-10-05 | 2018-04-12 | Magic Leap, Inc. | Surface modeling systems and methods |
WO2018117890A1 (en) * | 2016-12-21 | 2018-06-28 | Schlumberger Technology Corporation | A method and a cognitive system for predicting a hydraulic fracture performance |
US10241903B1 (en) | 2017-11-15 | 2019-03-26 | Accenture Global Solutions Limited | Parallel testing and reporting system |
US10409553B2 (en) * | 2017-11-15 | 2019-09-10 | Accenture Global Solutions Limited | Optimized construction of a sample imprint for selecting a sample dataset for comparison testing |
US10514890B2 (en) | 2017-11-15 | 2019-12-24 | Accenture Global Solutions Limited | Test case and data selection using a sampling methodology |
CA3156233A1 (en) | 2019-10-28 | 2021-05-06 | Conocophillips Company | Integrated machine learning framework for optimizing unconventional resource development |
CN111611632A (zh) * | 2020-04-17 | 2020-09-01 | 长沙矿山研究院有限责任公司 | 基于光固化3d打印的离散裂隙网络岩体模型的建模方法 |
EP3933339B1 (en) * | 2020-06-30 | 2024-03-20 | Mitutoyo Corporation | Method and computer program product for filtering a measurement data set usable for specifying and/or verifying an internal feature of a workpiece |
US11434759B2 (en) | 2020-07-07 | 2022-09-06 | Saudi Arabian Oil Company | Optimization of discrete fracture network (DFN) using streamlines and machine learning |
EP4256177A1 (en) * | 2020-12-07 | 2023-10-11 | Services Pétroliers Schlumberger | Fluid production network leak detection system |
CN113933909B (zh) * | 2021-08-30 | 2024-03-19 | 核工业北京地质研究院 | 一种快速确定裂谷盆地边界的方法 |
US11921250B2 (en) | 2022-03-09 | 2024-03-05 | Saudi Arabian Oil Company | Geo-mechanical based determination of sweet spot intervals for hydraulic fracturing stimulation |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1166207A (zh) * | 1995-10-06 | 1997-11-26 | 阿莫科公司 | 地震信号处理和勘察的方法和装置 |
US20030112704A1 (en) * | 2001-12-14 | 2003-06-19 | Goff Douglas Francis | Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume |
US20070168169A1 (en) * | 2005-10-26 | 2007-07-19 | Roxar Software Solutions As | Analysis and characterization of fault networks |
CN102667866A (zh) * | 2010-02-26 | 2012-09-12 | 雪佛龙美国公司 | 关心的地质体积的地球模型内的表面平滑 |
Family Cites Families (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6014343A (en) * | 1996-10-31 | 2000-01-11 | Geoquest | Automatic non-artificially extended fault surface based horizon modeling system |
FR2757947B1 (fr) * | 1996-12-30 | 1999-01-29 | Inst Francais Du Petrole | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches |
US7777606B2 (en) | 2007-01-09 | 2010-08-17 | Westerngeco L.L.C. | Fracture cluster mapping |
US20100191516A1 (en) * | 2007-09-07 | 2010-07-29 | Benish Timothy G | Well Performance Modeling In A Collaborative Well Planning Environment |
US8103493B2 (en) * | 2007-09-29 | 2012-01-24 | Schlumberger Technology Corporation | System and method for performing oilfield operations |
GB2466412B (en) * | 2007-10-22 | 2012-03-21 | Schlumberger Holdings | Formation modeling while drilling for enhanced high angle or horizontal well placement |
CA2705340C (en) * | 2007-12-21 | 2016-09-27 | Exxonmobil Upstream Research Company | Method and apparatus for analyzing three-dimensional data |
US9228415B2 (en) * | 2008-10-06 | 2016-01-05 | Schlumberger Technology Corporation | Multidimensional data repository for modeling oilfield operations |
EP2411912A2 (en) * | 2009-03-24 | 2012-02-01 | Chevron U.S.A. Inc. | A system and method for characterizing fractures in a subsurface reservoir |
WO2011077227A2 (en) * | 2009-12-21 | 2011-06-30 | Schlumberger Technology B.V. | Identification of reservoir geometry from microseismic event clouds |
GB201017898D0 (en) * | 2010-10-22 | 2010-12-01 | Internat Res Inst Of Stavanger | Earth model |
GB2500332B (en) * | 2010-12-30 | 2018-10-24 | Schlumberger Holdings | System and method for performing downhole stimulation operations |
FR2981475B1 (fr) * | 2011-10-12 | 2013-11-01 | IFP Energies Nouvelles | Methode pour construire un maillage d'un reservoir fracture avec un nombre limite de noeuds dans le milieu matrice |
WO2013119245A1 (en) * | 2012-02-10 | 2013-08-15 | Landmark Graphics Corporation | Systems and methods for selecting facies model realizations |
US9377546B2 (en) * | 2012-05-06 | 2016-06-28 | Schlumberger Technology Corporation | Automatic extraction and characterization of fault and fracture populations |
US20140058713A1 (en) * | 2012-08-21 | 2014-02-27 | Schlumberger Technology Corporation | Seismic modeling system and method |
US9523781B2 (en) * | 2012-12-27 | 2016-12-20 | Schlumberger Technology Corporation | Normalization seismic attribute |
US20140358510A1 (en) * | 2013-05-29 | 2014-12-04 | Chevron U.S.A. Inc. | System and method for characterizing uncertainty in subterranean reservoir fracture networks |
US10309216B2 (en) * | 2014-11-04 | 2019-06-04 | Shell Oil Company | Method of upscaling a discrete fracture network model |
US10408955B2 (en) * | 2014-11-19 | 2019-09-10 | Halliburton Energy Services, Inc. | Filtering microseismic events for updating and calibrating a fracture model |
FR3046479B1 (fr) * | 2016-01-04 | 2018-07-20 | Services Petroliers Schlumberger | Extrapolation de la permeabilite effective d'un reseau de fractures discretes |
-
2014
- 2014-11-24 US US14/551,667 patent/US11125912B2/en active Active
- 2014-11-25 EP EP14863440.5A patent/EP3074957B1/en active Active
- 2014-11-25 WO PCT/US2014/067391 patent/WO2015077768A1/en active Application Filing
- 2014-11-25 CN CN201480074060.6A patent/CN105934689A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1166207A (zh) * | 1995-10-06 | 1997-11-26 | 阿莫科公司 | 地震信号处理和勘察的方法和装置 |
US20030112704A1 (en) * | 2001-12-14 | 2003-06-19 | Goff Douglas Francis | Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume |
US20070168169A1 (en) * | 2005-10-26 | 2007-07-19 | Roxar Software Solutions As | Analysis and characterization of fault networks |
CN102667866A (zh) * | 2010-02-26 | 2012-09-12 | 雪佛龙美国公司 | 关心的地质体积的地球模型内的表面平滑 |
Non-Patent Citations (5)
Title |
---|
MANDO GUIDO BLOCHER ET AL.: "Three dimensional modelling of fractured and faulted reservoirs: Framework and implementation", 《CHEMIE DER ERDE-GEOCHEMISTRY》 * |
PAUL R. LAPOINTE ET AL.: "Reservoir Compartmentalization- Fractured Reservoir Discrete Feature Network Technologies", 《COMPARTMENTALIZATION RESEARCH REPORT》 * |
徐维胜等: "基于离散裂缝网络模型的储层裂缝建模", 《大庆石油学院学报》 * |
王建华: "DFN模型裂缝建模新技术", 《断块油气田》 * |
王晖等: "渤海C油田潜山裂缝型储集层随机离散裂缝网络模型的实现与优选方法", 《岩性油气藏》 * |
Also Published As
Publication number | Publication date |
---|---|
EP3074957A1 (en) | 2016-10-05 |
EP3074957B1 (en) | 2022-02-02 |
EP3074957A4 (en) | 2016-12-14 |
US11125912B2 (en) | 2021-09-21 |
US20150149142A1 (en) | 2015-05-28 |
WO2015077768A1 (en) | 2015-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105934689A (zh) | 地质特征拆分 | |
USRE49507E1 (en) | Faulted geological structures having unconformities | |
US11209560B2 (en) | Assignment of systems tracts | |
EP3535607B1 (en) | Seismic data processing artificial intelligence | |
CA2920499C (en) | Stratigraphic function | |
US11042676B2 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
EP3298437B1 (en) | Inversion for tectonic stress | |
US20140222403A1 (en) | Geologic model via implicit function | |
US20120265510A1 (en) | Grid from depositional space | |
GB2523869A (en) | Representation of unstructured grids | |
WO2016183043A1 (en) | Geologic stratigraphy via implicit and jump functions | |
WO2014004109A1 (en) | Seismic orthogonal decomposition attribute | |
US20190243017A1 (en) | Geologic Structural Model Generation | |
US10380793B2 (en) | Geobody surface reconstruction | |
EP3785233B1 (en) | Fault radiation based grid compartmentalization | |
WO2023219936A1 (en) | Adaptive multicolor reordering simulation system | |
WO2024155670A1 (en) | Volumetric geometry quantification of near-surface carbon containers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160907 |