CN112292714A - 基于断层辐射的网格隔分 - Google Patents
基于断层辐射的网格隔分 Download PDFInfo
- Publication number
- CN112292714A CN112292714A CN201980040259.XA CN201980040259A CN112292714A CN 112292714 A CN112292714 A CN 112292714A CN 201980040259 A CN201980040259 A CN 201980040259A CN 112292714 A CN112292714 A CN 112292714A
- Authority
- CN
- China
- Prior art keywords
- fault
- radiation
- cells
- model
- cell
- 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.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 112
- 238000000638 solvent extraction Methods 0.000 title description 7
- 238000000034 method Methods 0.000 claims abstract description 233
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 13
- 239000003086 colorant Substances 0.000 claims description 9
- 230000008569 process Effects 0.000 description 56
- 239000010410 layer Substances 0.000 description 39
- 238000004088 simulation Methods 0.000 description 38
- 238000013459 approach Methods 0.000 description 25
- 239000012530 fluid Substances 0.000 description 22
- 238000005286 illumination Methods 0.000 description 17
- 239000000243 solution Substances 0.000 description 16
- 238000004422 calculation algorithm Methods 0.000 description 15
- 238000004519 manufacturing process Methods 0.000 description 15
- 238000001514 detection method Methods 0.000 description 14
- 238000012545 processing Methods 0.000 description 13
- 239000000463 material Substances 0.000 description 11
- 238000004891 communication Methods 0.000 description 10
- 230000008021 deposition Effects 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- 230000015572 biosynthetic process Effects 0.000 description 9
- 230000003993 interaction Effects 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 8
- 238000005192 partition Methods 0.000 description 8
- 230000009471 action Effects 0.000 description 7
- 229930195733 hydrocarbon Natural products 0.000 description 7
- 150000002430 hydrocarbons Chemical class 0.000 description 7
- 238000007726 management method Methods 0.000 description 7
- 238000003325 tomography Methods 0.000 description 7
- 238000013508 migration Methods 0.000 description 6
- 230000005012 migration Effects 0.000 description 6
- 239000011435 rock Substances 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000002347 injection Methods 0.000 description 5
- 239000007924 injection Substances 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 5
- 238000013507 mapping Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 230000001186 cumulative effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000002452 interceptive effect Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 3
- 230000007717 exclusion Effects 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 238000013517 stratification Methods 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 238000010796 Steam-assisted gravity drainage Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000012792 core layer Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- -1 for example Chemical class 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000003094 perturbing effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011867 re-evaluation Methods 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000013515 script Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/12—Geometric CAD characterised by design entry means specially adapted for CAD, e.g. graphical user interfaces [GUI] specially adapted for CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
- G06T2219/20—Indexing scheme for editing of 3D models
- G06T2219/2012—Colour editing, changing, or manipulating; Use of colour codes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- Geophysics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Geometry (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- Remote Sensing (AREA)
- General Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Mathematical Optimization (AREA)
- Evolutionary Computation (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Human Computer Interaction (AREA)
Abstract
一种方法可以包括访问地下区域的模型,其中该模型包括与地下区域中的断层相对应的断层。经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元;以及经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;并且至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该侧面性属性值表示单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。
Description
相关申请的交叉引用
本申请要求于2018年4月22日提交的序列号为62/661,070的美国临时申请的优先权和权益,其通过引用合并于此。
背景技术
与沉积盆地相关的现象可以使用作为机器的模拟器的各种方程来建模,该模拟器可以是分布式机器。为了应用数字技术,可以使用包含节点、单元等的网格来离散这些方程式。其中盆地包含各种类型的特征(例如,地层,断层等)、节点、单元等,其中网格可以表示或分配给这些特征。反过来,离散方程可以更好地表示盆地及其特征。
发明内容
一种方法可以包括访问地下区域的模型,其中该模型包括对应于该地下区域中的断层的断层;经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元;以及经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;并且至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该侧面性属性值表示单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。系统可以包括处理器;可操作地耦合到处理器的存储器;存储在存储器中的处理器可执行指令,以指示系统:访问地下区域的模型,其中该模型包括与该地下区域中的断层相对应的断层;以及经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元,并且经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;并至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该侧面性属性值表示单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。还公开了各种其他装置、系统、方法等。
提供该概述是为了介绍一些构思,这些构思将在下面的详细描述中进一步描述。该概述不旨在标识所要求保护的主题的关键或必要特征,也不旨在用于帮助限制所要求保护的主题的范围。
附图说明
通过参考以下结合附图的描述,可以更容易地理解所描述的实施方式的特征和优点。
图1示出了示例系统,其包括用于对地质环境进行建模的各种组件以及与地质环境相关联的各种设备;
图2示出了沉积盆地的示例、方法的示例、地层的示例、井孔的示例、井孔工具的示例、惯例的示例以及系统的示例;
图3示出了技术的示例、阶梯状网格的示例以及断层的阶梯状表示的示例;
图4示出了地质模型的示例和模型的一部分的示例;
图5示出了关于地面的网格模型的示例和投影技术的示例;
图6示出了地下环境的表示的示例;
图7示出了方法的示例;
图8示出了呈现给显示器的图形用户界面的示例;
图9A、9B和9C示出了图形用户界面的示例;
图10示出了方法的示例;
图11示出了方法的示例;
图12示出了方法的示例;
图13示出了地质环境的一部分的示例;和
图14示出了系统和联网系统的示例组件。
具体实施方式
该描述不应被理解为限制性的,而仅仅是出于描述实施方式的一般原理的目的而做出的。所描述的实施方式的范围应参考所给出的权利要求书确定。
与沉积盆地有关的现象(例如,地下区域、是否低于地面、水面等)可以使用各种方程式(例如,应力、流体流动、相等)来建模。作为示例,盆地的数值模型可发现用于理解与自然资源的勘探和开采有关的各种过程(估算储量、钻井、预测产量等)。
对于数值技术的应用,可以使用包括节点、单元等的网格来离散方程式。例如,诸如有限差分方法之类的数值技术可以包括相对于空间坐标离散化温度的一维差分热方程,以近似温度导数(例如,一阶,二阶等)。在感兴趣时间的情况下,还可以提供温度相对于时间的导数。关于空间坐标,数值技术可以依赖于包括各种节点的空间网格,其中在求解热方程时(例如,受边界条件、生成项等的影响)将为每个节点提供温度。这样的示例可以应用于空间中的多个维度(例如,其中离散化应用于多个维度)。因此,网格可以将感兴趣的体积(VOI)离散化为可以分配属性或与属性(例如,孔隙度、岩石类型等)关联或关联的基本元素(例如,单元或网格块),这可能与模拟物理过程(例如,流体流动、储层压实等)密切相关。
作为数值技术的另一示例,考虑其中空间可以由一维或多维“元素”表示的有限元方法。对于一个空间维度,一个元素可以由沿着一个空间坐标定位的两个节点表示。对于多个空间维度,一个元素可以包括任意数量的节点。此外,一些方程式可以由节点总数表示,而其他方程式可以由少于节点总数表示(例如,考虑Navier-Stokes方程式的示例,其中少于节点总数代表压力)。有限元方法可包括提供可定义三角形元素(例如3D中的四面体,多维空间中的高阶单纯形等)或四边形元素(例如3D中的六面体或金字塔等)或多边形元素(例如3D棱镜等)的节点。如由网格的相应节点所定义的这样的元素可以被称为网格单元。
数值技术的另一个示例是有限体积法。对于有限体积法,可以在网格的离散位置(例如,网格的其周围包含“有限体积”的节点)计算模型方程变量的值。有限体积方法可以应用发散定理来评估每个有限体积的表面处的通量,使得进入给定有限体积的通量等于离开一个或多个相邻有限体积的通量(例如,遵守守恒定律)。对于有限体积方法,网格的节点可以定义网格单元。
作为示例,有限体积流模拟器可以使用网格来模拟现象,其中由网格定义的网格单元可以包括可通过三个索引寻址的6个面(例如,长方体)(例如,使得可以将网格视为“结构化”网格)以及网格的几何形状遵守一个或多个条件(例如,单元未越过地质断层且单元未越过地质层位)。作为示例,为了满足地质断层条件,网格可以在一个或多个地质断层上偏移。在断层网络的拓扑复杂(例如,许多X和Y形交叉点)的域中构造这种网格可能是不平凡的,并且需要随着断层网络复杂度的增加而非线性扩展的资源。
作为示例,沉积盆地建模的方法可以包括由节点、柱和单元组成的柱网格。例如,在三维中,八个节点可以定义一个单元,该单元可以被称为网格单元(例如,柱状网格单元)。在柱状网格模型中,可以使用索引i、j和k(例如,索引坐标系或空间,其可以表示为I、J和K或其他字符、符号等。)在索引域中对网格单元进行索引。例如,可以在其最浅的左下角索引一个立方网格单元(即由八个角节点定义),网格单元的数量可以是模型的i、j和k维度的乘积。在这样的示例中,每个网格单元可以由其八个节点定义,其可以根据高度和罗盘方向(例如,底部西南(basesouthwest),顶部西南(topsouthwest),底部西北,顶部西北等)进行标记。柱状网格可以建模,例如,断层(例如,切割柱状网格的表面),层位(例如,“k”索引)、区域(例如,两个层位之间的体积)、线段(例如,由断层平面界定的网格单元的连续块)等,并可以用于指定属性(例如,土属性)。
虽然相对于柱状网格描述了索引坐标系,但是可以将索引坐标系与其他类型的栅格结合使用。例如,可以使用诸如i、j和k之类的索引来表征可以定义单元的网格,以表示三个空间维度。这样的索引可能能够例如以所谓的结构化方式(例如,与非结构化方式相反)表示网格。作为示例,结构化网格可以促进各种类型的操作,例如与矩阵有关的那些操作,例如,其中最近的相邻可以在矩阵内形成簇或带。依次地,可以使用带状求解器或其他合适的技术来处理矩阵。例如,对于非结构化网格的求解器,它可能依赖于连接性信息的输入,该连接性信息指定了网格节点如何与各个单元相关。在这样的示例中,矩阵可能不容易适应带状或其他矩阵处理技术,从而又会增加计算资源需求、计算时间等。
作为示例,包括自然的(i、j、k或I、J、K)索引系统的结构化网格可以改善存储,例如,有助于其中单元索引和连接性可能不会存储在存储器中并且可以从存储器中的记录/条目的顺序推导出的识别拓扑相邻。在这样的示例中,存储结构化网格可以比例如存储相似大小的非结构化网格使用更少的存储器。此外,作为示例,对于大型方程组的构建(例如,独立于其分辨率),这可能涉及在给定网格单元的拓扑相邻上重复迭代,当与使用非结构化网格相比时,这种方法例如可以快大约一个数量级。作为示例,可以生成结构化网格的方法可以提供与一个或多个框架的兼容性(例如,是否是当前的、遗留的等)。
如所提及的,在沉积盆地(例如,地下区域)包括各种类型的特征(例如,地层、断层等)的情况下,网格的节点、单元等可以表示或分配给,此类特征。反过来,离散方程可以更好地表示沉积盆地及其特征。例如,与非结构化网格相比,可以表示沉积盆地及其特征的结构化网格可以允许更多的模拟运行、更多的模型复杂性、更少的计算资源需求、更少的计算时间等。
作为示例,网格可以符合结构特征,例如,Y-断层、X-断层、低角度不整合面、盐体、侵入体等(例如,地质上的断面),以更全面地捕获地质模型的复杂性。作为示例,网格可以可选地符合地层学(例如,除了一个或多个地质断面之外)。关于地质断面,这些可能包括模型断面,例如一个或多个模型边界。作为示例,网格可以填充有例如通过地统计方法生成的属性字段。
作为示例,可以通过地震学来识别断面,其中例如存在地下边界或界面,在该地下边界或界面处物理量(例如地震波的传播速度)突然改变。例如,在地壳和地幔之间的莫霍维奇断面处,P波的速度从大约6.5km/s增加到大约8km/s。
地震学可以包括地震层析成像,其可以是层析扫描的形式。例如,可以根据用于地下区域的勘测几何形状(例如,获取几何形状)来利用源和接收器(例如,经由陆地和/或海洋设备)来获取其中源发出地震能并且接收器接收至少一部分发出的地震能量的数据。在这样的示例中,利用勘测几何的先验知识,至少部分地基于接收器数据的接收,可以由计算装置或计算系统来实现一种或多种重构技术以空间地表征地下区域,例如通过输出包含取决于地下区域的材料属性、结构等的值的模型。作为示例,模型可以是作为像素呈现到显示器的图像,或者例如,模型可以是具有空间信息的呈现到显示器的值,或者例如,模型可以是地下区域的至少一部分的基于网格(例如,基于单元)的模型,其中网格(例如,或单元)在某个时间点或者以其他方式在结构上(例如,层状蛋糕等)在空间上代表地下区域。这样的计算装置或计算系统可以是作为地震层析成像扫描仪的分布式“机器”的一部分。地震层析成像扫描仪接收地震数据并至少部分基于地震数据生成模型,以提供对被扫描对象(即地球)的理解。作为示例,地震能量可以通过诸如将流体泵送到地球中(例如,经由井)之类的间接过程来产生,使得地球的破裂发生在这种破裂涉及地震排放物(例如,微地震排放物)的地方。地震层析成像扫描可以包括接收经由地球的水力压裂产生的微地震排放物,例如,处理这种排放物(例如地震数据),以输出有关一个或多个裂缝的信息(例如一个或多个裂缝的空间位置)。
地震解释是一种涉及检查地震数据(例如,关于位置、时间或深度)以识别一种或多种类型的地下结构(例如,层位、断层、地质体等)的过程。解释过程可以考虑垂直地震剖面、直线和交叉线方向、称为水平时间切片的水平地震剖面等。地震数据可以可选地与其他数据(例如测井数据)一起解释。
作为示例,可以使用地震模拟软件例如地震模拟软件框架(Schlumberger Limited,得克萨斯州休斯敦)进行地震解释,该软件包括各种特征,例如以执行属性分析(例如关于3D地震立方体、2D地震线等),以分析其他数据,建立模型等。尽管提到了地震模拟软件框架,但是也可以采用其他类型的软件、框架等。
图1示出了系统100的示例,该系统100包括各种管理组件110,以管理地质环境150(例如,包括沉积盆地、储层151、一个或多个断层153-1、一个或多个地质体153-2等的环境)的各个方面。例如,管理组件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可以包括模拟器的一个或多个特征,例如,储层模拟器(Schlumberger Limited,休斯敦德克萨斯州)、储层模拟器(Schlumberger Limited,休斯敦德克萨斯州)等。作为示例,模拟组件、模拟器等可以包括用于实现一种或多种无网格技术(例如,求解一个或多个方程式等)的特征。作为示例,可以相对于一种或多种增强的采收技术(例如,考虑诸如SAGD等的热过程)来模拟一个或多个储层。
在示例实施例中,管理组件110可以包括诸如地震模拟软件框架(Schlumberger Limited,得克萨斯州休斯敦)之类的框架的特征。框架提供的组件可优化勘探和开发操作。框架包括地震模拟软件组件,这些组件可以输出信息,例如通过提高资产团队的生产率来提高储层性能。通过使用这样的框架,各种专业人员(例如,地球物理学家、地质学家和储层工程师)可以开发协作工作流程并集成操作以简化流程。这样的框架可以被认为是应用程序,并且可以被认为是数据驱动的应用程序(例如,出于建模、模拟等目的而输入数据的情况)。
在示例实施例中,管理组件110的各个方面可以包括根据框架环境的规范进行操作的附加组件(add-on)或插件。例如,诸如框架环境(Schlumberger Limited,得克萨斯州休斯敦)之类的框架环境允许将附加组件(或插件)集成到框架工作流程中。框架环境利用工具(Microsoft Corporation,华盛顿州雷德蒙德)。在示例实施例中,各种组件可以被实现为符合框架环境的规范并且根据框架环境的规范(例如,根据应用编程接口(API)规范等)运行的附加组件(或插件)。
图1还示出了框架170的示例,该框架170包括模型模拟层180以及框架服务层190、框架核心层195和模块层175。框架170可以包括框架,其中,模型模拟层180是托管框架应用程序的以模型为中心的软件包。在示例实施例中,软件可以被认为是数据驱动的应用程序。软件可以包括用于模型构建和可视化的框架。
作为示例,框架可以包括用于实现一种或多种网格生成技术的特征。例如,框架可以包括用于从地震数据的解释接收信息的输入组件、至少部分基于地震数据、测井数据、图像数据等的一个或多个属性。这样的框架可以包括格网生成组件,其处理输入信息(可选地结合其他信息)以生成格网。
在图1的示例中,模型模拟层180可以提供域对象182,充当数据源184,提供呈现186并提供各种用户界面188。呈现186可以提供图形环境,其中应用程序可以显示它们的数据而用户界面188可以为应用程序用户界面组件提供共同的外观和感觉。
作为示例,域对象182可以包括实体对象、属性对象以及可选地其他对象。实体对象可用于在几何上表示井、表面、体、储层等,而属性对象可用于提供属性值以及数据版本和显示参数。例如,实体对象可以表示井,其中属性对象提供日志信息以及版本信息和显示信息(例如,以显示井作为模型的一部分)。
在图1的示例中,数据可以存储在一个或多个数据源(或数据存储,通常为物理数据存储装置)中,该数据源可以位于相同或不同的物理站点,并且可以经由一个或多个网络访问。模型模拟层180可以被配置为对项目进行建模。这样,可以存储特定项目,其中存储的项目信息可以包括输入、模型、结果和案例。因此,在建模会话完成时,用户可以存储项目。在稍后的时间,可以使用模型模拟层180来访问和恢复该项目,该模型模拟层可以重新创建相关域对象的实例。
在图1的示例中,地质环境150可以包括包括储层151和一个或多个其他特征(例如断层153-1、地质体153-2等)的层(例如,分层)。作为示例,地质环境150可以配备有各种传感器、检测器、致动器等中的任何一个。例如,设备152可以包括用于接收和发送关于一个或多个网络155的信息的通信电路。这样的信息可以包括与井下设备154相关联的信息,该设备可以是获取信息、辅助资源回收等的设备。其他设备156可以位于远离井场的位置,并且包括感测、检测、发射或其他电路。这样的设备可以包括用于存储和传递数据、指令等的存储和通信电路。作为示例,可以提供一个或多个卫星以用于通信、数据采集等的目的。例如,图1示出了与可以被配置用于通信的网络155通信的卫星,请注意,该卫星可以附加地或替代地包括用于成像的电路(例如,空间、光谱、时间、辐射度等)。
图1还示出了地质环境150,其任选地包括与井相关的设备157和158,该井包括可以与一个或多个裂缝159相交的基本上水平的部分。例如,考虑页岩地层中的井可以包括天然裂缝、人造裂缝(例如水力裂缝)或天然裂缝和人造裂缝的组合。作为示例,可以在侧向扩展的储层上钻井。在这样的示例中,可能存在性质、应力等的横向变化,其中对这种变化的评估可以帮助计划、操作等来开发横向扩展的储层(例如,通过压裂、注入、开采等)。作为示例,设备157和/或158可以包括用于压裂、地震感测、地震数据分析、一个或多个裂缝的评估等的组件、一个或多个系统等。
如所提及的,系统100可以用于执行一个或多个工作流程。工作流程可以是包括多个工作步骤的过程。工作步骤可以对数据进行操作,例如,创建新数据、更新现有数据等。作为示例,例如,可以基于一个或多个算法对一个或多个输入进行操作并创建一个或多个结果。作为示例,系统可以包括用于创建、编辑、执行等工作流程的工作流程编辑器。在这样的示例中,工作流程编辑器可以提供一个或多个预定义工作步骤、一个或多个定制工作步骤等的选择。作为示例,工作流程可以是可在软件中实现的工作流程,例如,其可以对地震数据、地震属性进行操作。作为示例,工作流程可以是可在框架中实现的过程。作为示例,工作流程可以包括访问一组指令(例如,插件(例如,外部可执行代码等))的一个或多个工作步骤。
图2示出了沉积盆地210(例如,地质环境)的示例、用于模型构建(例如,用于模拟器等)的方法220的示例、地层230的示例、地层中井孔235的示例、惯例240的示例和系统250的示例。
作为示例,储层模拟、石油系统建模等可以用于表征各种类型的地下环境,包括诸如图1所示的环境。
在图2中,作为地质环境的沉积盆地210包括层位、断层、一个或多个地质体以及在一定地质时期内形成的相(facies)。这些特征例如相对于笛卡尔坐标系(例如,x,y和z)或其他坐标系(例如,圆柱、球形等)在空间中以二维或三维分布。如所示的,模型建立方法220包括数据采集块224和模型几何块228。一些数据可能涉及建立初始模型,并且此后,可以可选地响应于模型输出、时间变化、物理现象、附加数据等来更新模型。作为示例,用于建模的数据可以包括以下一项或多项:深度或厚度图以及断层几何形状和地震、遥感、电磁、重力、露头(outcrop)和测井数据的时间。此外,数据可以包括由假定跟随地质事件(“iso”时间)的相变化(例如,由于地震不整合)产生的深度和厚度图,并且数据可以包括横向相变化(例如,由于沉积特征的横向变化)。
为了进行地质过程的建模,可以提供数据,例如,诸如地球化学数据(例如,温度、干酪根类型、有机物富集度等)、定时数据(例如,来自古生物学、辐射测年,磁逆转、岩石和流体性质等)和边界条件数据(例如,热流历史、地表温度、古水深度等)的数据。
在盆地和石油系统建模中,可以例如通过使用一种或多种数值技术求解偏微分方程(PDE)来对诸如沉积物中的温度、压力和孔隙度分布之类的数量进行建模。建模还可以相对于时间对几何形状建模,例如,以考虑到源自地质事件的变化(例如,材料的沉积、材料的侵蚀、材料的移动等)。
诸如框架(Schlumberger Limited,得克萨斯州休斯敦)之类的建模框架包括用于输入各种类型的信息(例如地震、井、地质等)的特征,以模拟沉积盆地的演化。框架经由输入各种数据(例如地震数据、井数据和其他地质数据)来为石油系统建模,例如,以模拟沉积盆地的演化。框架可以预测是否以及如何向储层中注满了烃,例如包括地下或地面条件下烃产生的来源和时间、运移路线、数量、孔隙压力和烃类型。结合诸如框架之类的框架,可以构建工作流程以提供盆地到潜在规模的勘探解决方案。框架之间的数据交换可以促进模型的构建,数据分析(例如,使用框架功能分析的框架数据)以及工作流程的耦合。
如图2所示,地层230包括水平表面和各种地下层。作为示例,井孔可以是竖直的。作为另一个示例,井孔可以被偏离。在图2的示例中,井孔235可以被认为是竖直的井孔,例如,其中z轴法向于地层230的水平表面向下延伸。作为示例,可以将工具237定位在井孔中,例如,以获取信息。如所提及的,井孔工具可以被配置为获取电井孔图像。作为示例,全孔地层微观成像仪(FMI)工具(Schlumberger Limited,得克萨斯州休斯敦)可以获取井孔图像数据。这种工具的数据采集顺序可以包括在关闭获取垫的情况下将工具运行到井孔中,将该垫打开并压在井孔壁上,将电流输送到定义井孔的材料中,同时将工具平移到井孔中,以及远程感应电流,该电流因与材料的相互作用而改变。
作为示例,井孔可以是竖直的、偏离的和/或水平的。作为示例,可以定位工具以获取井孔的水平部分中的信息。对此类信息的分析可能会发现孔洞、溶蚀(dissolution)平面(例如,沿着层理平面的溶蚀)、与应力相关的特征、倾角事件等。作为示例,工具可以获取有助于描述裂缝性储层特征的信息,可选地在那里裂缝可以是天然的和/或人造的(例如水力压裂)。这些信息可能有助于完井、增产处理等。作为示例,可以使用诸如框架(Schlumberger Limited,德克萨斯州休斯敦)之类的框架来分析由工具获取的信息。
关于倾角的惯例240,如图所示,平面的三维取向可以由其倾角和走向(strike)来定义。倾角是在特定方向的竖直面上测量的平面相对于水平面(如假想平面)的倾斜角。倾角可以由幅度(例如,也称为角度或量)和方位角(例如,也称为方向)来定义。如图2的惯例240所示的,各种角度φ表示向下倾斜的角度,例如,从假想的水平面(例如,平坦的上表面);然而,倾角是指倾斜平面向其倾斜的方向(例如,其可以相对于角度、罗盘方向等给出)。)。图2中地惯例所示的另一个特征是走向,它是由倾斜平面和水平面(例如,将平坦的上表面视为假想的水平面)的交叉点所创建的直线的取向。
例如,根据情况、所收集数据的取向等,一些与倾角和走向有关的附加术语可能适用于分析。一个术语是“真实倾角”(例如,参见图2的惯例240中的DipT)。真实倾角是指直接垂直于走向的平面的倾角(参见,例如,朝北并标有“走向”和角度α90的线)以及倾角幅度的最大可能值。另一个术语是“表观倾角”(例如,参见图2的惯例240中的DipA)。表观倾角可以是平面上在真正倾角方向以外的任何方向上测量的平面的倾角(参见,例如,φA作为角度α的DipA);然而,表观倾角可能等于真实倾角(参见,例如,φ作为关于走向的角度α90的DipA=DipT)。换句话说,在使用术语表观倾角的情况下(例如,在方法、分析、算法等中),对于特定的倾角平面,用于“表观倾角”的值可能等同于该特定倾角平面的真实倾角。
如图2的惯例240所示的,在垂直于走向的横截面中看到的平面的倾角为真实倾角(参见,例如,φ作为关于走向的角度α90的DipA=DipT的表面)。如图所示,在横截面中沿任何其他方向观察到的倾角都是表观倾角(例如,参见标记为DipA的表面)。此外,如图2的惯例240所示,表观倾角可以是大约0度(例如,平行于切割平面的边缘沿走向方向延伸的水平表面)。
就观察井眼中的倾角而言,在竖直钻井的井中观察到真实倾角。在以任何其他取向(或偏差)钻井的井中,观察到的倾角是表观倾角(例如,其被某些人称为相对倾角)。为了确定在这样的井孔中观察到的平面的真实倾角值,作为示例,可以将矢量计算(例如,基于井孔偏差)应用于一个或多个表观倾角值。
如所提及的,在井孔图像的沉积学解释中发现使用的另一个术语是“相对倾角”(例如,DipR)。可以从砂体中的倾角减去(例如,使用矢量减法)从沉积在非常平静的环境中的岩石中的井孔图像测得的真实倾角的值。在这样的示例中,所得的倾角被称为相对倾角,并且可以发现用于解释砂体取向。
可以针对分析、解释、属性等来使用诸如惯例240之类的惯例(例如,参见图1的系统100的各个框)。作为示例,可以部分地通过倾角来描述各种类型的特征(例如,沉积层理、断层和裂缝、单斜嵴、火成岩堤和岩床、变质岩成层等)。作为示例,随着层接近地质体,倾角可能会在空间上发生变化。例如,考虑可能由于各种力(例如浮力等)而上升的盐体。在这样的示例中,当盐体向上移动时,倾角可能会趋向于上升。
地震解释可以旨在至少部分地基于一个或多个倾角参数(例如,角度或幅度、方位角等)来识别和/或分类一个或多个地下边界。作为示例,可以至少部分地通过角度、至少部分地通过方位角等描述各种类型的特征(例如,沉积层理、断层和裂缝、单斜嵴、火成岩堤和岩床、变质岩成层等)。
作为示例,可以提供用于石油驱出和迁移的方程式,其可以例如相对于一段时间建模和模拟。石油从原料的迁移(例如一次迁移或驱出)可能包括使用饱和度模型,其中迁移饱和度值控制驱出。关于石油(例如,石油或天然气)的二次迁移的确定可以包括利用流体的流体动力势并考虑促进流体流动的驱动力。这样的力可以包括浮力梯度、孔隙压力梯度和毛细管压力梯度。
如图2所示,系统250包括一个或多个信息存储装置252、一个或多个计算机254、一个或多个网络260以及一组或多组指令270。至于一个或多个计算机254,每一台计算机可以包括一个或多个处理器(例如,处理核)256和存储器258,用于存储例如可由一个或多个处理器中的至少一个执行的指令。作为示例,计算机可以包括一个或多个网络接口(例如,有线或无线)、一个或多个图形卡、显示接口(例如,有线或无线)等。作为示例,诸如表面图像(例如,卫星、地质、地球物理等)的图像可以被存储,处理,通信等。作为示例,数据可以包括SAR数据、GPS数据等,并且可以被存储在例如一个或多个存储装置252中。
作为示例,一个或多个指令集270可包括可由一个或多个处理器执行以指示系统250执行各种动作的指令(例如,存储在存储器中)。作为示例,系统250可以被配置为使得一组或多组指令270提供用于建立图1的框架170或其一部分。作为示例,可以使用一组或多组指令来执行一个或多个方法、技术等,例如,这些指令可以是图2的一组或多组指令270中的一个或多个。
图3示出了用于处理断层的技术的示例。特别地,技术310和/或技术330可以用于创建旨在符合地质层位的断层结构化网格。
关于技术310,其基于表面网格的挤压。例如,建立与断层网络的断层相符并被其切断的拓扑区域网格312,并提供一些坐标线314。给定网格312和坐标线314,可通过挤压过程创建3D网格316,该过程可包括沿线性或弯曲坐标线将基础网格重复多次并使其适应地层位几何形状,并以这种方式计算轨迹使得它就不会越过断层网络。
技术310可能引起一些问题,例如,可能无法设计从网格的底部到顶部的坐标线而不会与任何断层交叉,并且可能无法有效地最小化网格单元的失真(例如,将失真定义为到纯正交几何形状的距离的度量)。
关于技术330,其包括首先创建与层位表面一致的3D(例如,未断层)网格,然后,例如,栅格化网格内的断层。在这样的示例中,栅格化动作可以对应于找到与断层表面最佳对应的一组单元小平面(facet),并沿着这些小平面“拆缝”(例如,取消缝制)网格。这样的栅格化过程涉及计算与要栅格化的对象相交的网格的拓扑对偶的一组边缘,在图3的示例中,该边缘是断层,如3D网格332的侧视图和放大的简化侧视图334中所示的。
技术330可能引起一些问题。例如,代表断层的网格小面可能无法很好地近似初始断层表面的几何形状。并且,这样的栅格化操作可能涉及创建断层表面的阶梯状(或之字形)表示(参见例如334)。阶梯状几何形状的后果可能与岩石物理性质密切相关-从而可能影响流动模拟-这可能与断层表面(例如,与断层渗透率相关的那些)以及井相交于断层的情况有关,因为与断层相交点的位置以某种程度的不准确性表示。
应用于网格的栅格化可能会使网格几何形状受到一个或多个约束。例如,考虑一种希望在断层的两侧保持匹配的阶梯形(例如,避免出现间隙、重叠等),而另一种希望保持遵循地质表面的网格分层(例如,根据地震数据或井数据)。由于断层的位移可能会在断层的两侧之间产生不对应于网格竖直方向上的整数个单元(例如,k索引、时间或深度)的偏移,因此扭曲可能会在断层的附近发生(例如,通过竖直拉伸、挤压或合并网格单元,人为移动层位或两者兼有)。
图4示出了断层网格405的示例,该断层网格405包括代表地质环境中的断层406的阶梯状部分。图4还示出了柱网格模型410的一部分的示例,其包括节点408和连接409,其中沿着I、J、K索引参考系统的K索引的连接可以对应于柱。如图所示,柱状网格模型410包括断层406的阶梯状表示;然而,在实际的地质环境中,断层不太可能以与柱网格相对应的方式在空间上被阶梯化,特别是考虑到柱网格的分辨率。作为示例,柱状网格模型410可以提供层位的表示,其中一个或多个层位可以部分地在K索引的方向上移动,使得在断层406的阶梯状表示的一侧上,层位具有第一K索引,并且在断层406的阶梯式表示的另一侧,同一层位具有第二不同的K索引。换句话说,断层可以移动层位,使得层位在断层上是不连续的。
图5将图4的柱状网格模型410示为包括断层表面510,当与断层406的阶梯状表示相比时,其可能是物理上更准确的表示(例如,更现实)。换句话说,断层406的阶梯状表示可以是在由正交坐标(例如,I、J和K索引方向)限定的离散空间中的断层表面510的近似。作为示例,断层可以是平面的,并且可以与索引I、J和K中的一个或多个偏移的角度倾斜。这样,断层可能不会再次定位或与一个或多个坐标系方向对准,这可能是通过成阶梯状近似断层的原因。
可以采用成阶梯状保持网格特性。例如,出于计算的目的,例如为了帮助迭代求解技术、数组处理等中的一个或多个数值稳定性,可能需要保持网格单元的形状。作为示例,可以采用成阶梯状以方便使用索引,例如图5中所示的I、J和K索引,这些索引可以延续到数值求解器中的数组中,该求解器可解决相对于柱网格模型在空间上定义的方程组。但是,由于断层的性质在空间上可能与其在柱网格模型中的阶梯状表示在空间上不同,因此求解断层的阶梯状表示时的方程组可以提供解决方案(例如,关于压力、流量等),其不能充分匹配自然现象,尤其是在断层附近(例如,其他断面)。
图6示出了地质环境的各种结构模型表示。这样的模型至少部分地基于使用一种或多种类型的传感器(例如,地震层析成像传感器等)获取的数据。人眼看不见各种地下特征,因此利用仪器来探索、测量、感知等地下物质(例如岩石、流体等)。
三维立体胶片(也称为三维胶片、3D胶片或S3D胶片)是一种电影,其可增强深度感知的错觉,从而增加了第三维度。制作3D胶片的常用方法是从立体摄影获得的。在这种方法中,使用常规电影摄影机系统记录从两个角度看到的图像(例如,或计算机生成的图像在后期制作中生成两个角度),并且使用特殊的投影硬件和/或眼镜将每个图像的可见度限制在观看者的左眼或右眼。
代替常规的电影摄影机系统生成光图像数据(例如,像素、颜色模型数据等),对地下区域的勘探利用了不同类型的数据,如可以通过诸如可以定位在井孔中的井下工具和/或可以发射和/或获取地震能量(例如反射地震学等)的地面工具(例如,陆上、海洋等)来获取的数据。在3D电影中以及在对地球的探索中,最终会生成一模型,以允许增强的可视化,或者就地球而言,允许一种或多种其他类型的处理方法,这些方法可以帮助您告知进行决策、现场操作(例如钻井、压裂)等)等。
在图6中,各种模型提供了根据结构模型创建沉积网格(depogrid)(例如,一种模型类型)的描述。将地质空间中的结构模型(a)转换(例如,一对一和可逆映射)到沉积空间(deposapce;(b)),在该沉积空间中合适的层位是水平面。沉积空间由正交坐标系(u、v、w)表示,竖直坐标w定义了地质时间参数。沉积空间中的结构模型可以在u和v方向上均匀切片,每个区域在w方向(c)上被指定数量的平面划分。断层可以精确地表示为沉积空间中该网格的切割单元。沉积空间中应用于网格(c)的逆映射可恢复地质空间(d)中的沉积网格。剪切的单元独立转换到地质空间,以准确表示原始断层框架模型。例如,将沉积网格(d)与结构模型(a)进行比较。
在图6的示例中,三个断层标记为610-1、610-2和610-3。如在结构模型(a)中所见的,断层610-1和断层610-2之间的层被移位,使得断层610-2的相对侧上的至少一些层在它们的上部和下部层位不对齐。在沉积空间(b)中,显示了对齐,并且在网格的切片空间(c)中也显示了对齐;而在沉积网格(d)中,位移表示为结构模型(a)中的位移。如上所述,断层可以精确地表示为沉积空间中网格的切割单元(例如,沉积网格)。在这样的示例中,断层可以是分隔物,例如,沉积网格包括到断层的一侧的单元,以及到断层的另一相对侧的其他单元。在这样的示例中,可以为每个单元分配代表地下区域的属性的属性。这样的属性可以是或包括材料属性,其中可以通过岩性对材料进行分类(例如,材料的矿物含量、晶粒尺寸、质地、颜色等的宏观性质)。
作为示例,一种方法可以包括为具有断层的体单元网格生成隔区。隔分可以是储层建模工作流程、工程工作流程或其他工作流程的输入。隔分可以包括在储层的不同部分中分配属性,如通过网格建模。例如,考虑一种方法,该方法包括定义地质边界上的油气、油水接触和流体属性,以估算储量和/或预测产量。在这样的示例中,模拟器可以提供在一个或多个条件下的储层建模,该条件可以包括与一个或多个现场操作相关联的一个或多个条件。例如,将注入视为现场操作(例如,通过一个或多个井的流体注入),将生产视为现场操作(例如,通过一个或多个井的流体生产),将压裂视为现场操作等。在各种现场操作中,例如,响应于流体注入和生产中的一种或多种,油气和/或油水之间的边界可以在储层内空间地移动。在一个水平上,从储层中生产流体可以被视为如何最佳地生产流体(例如,从储层中去除流体)的问题。这样的问题非常复杂,并且如果执行不当,可能会混淆进一步的生产(例如,导致对生产计划、技术等进行重新评估)。
作为示例,隔分可以包括指定在感兴趣的断层的相对侧上的两个单元被分配不同的隔区号(例如,同时将隔区的数目保持在最小可能的最小值)。这种隔分可能是一个复杂的技术问题,尤其是可以预期的是,某些断层部分会“超出”其实际物理范围。
在勘探和生产框架内的地下建模中,可以以适合于地质建模、模拟、现场操作的计划、现场操作的执行、现场操作的控制等的不同比例生成结构化网格。
结构化网格的生成和使用在建模和模拟端到端工作流程的不同阶段提供了一些特殊的好处,其中系统可以利用单元的逻辑(I、J、K)邻域来推断建模网格属性以及沿层和贯穿层的单元方向时的关系。后者的好处可以在储层模拟中使用,以简化单元间传输率的计算;然而,这种简化引入了误差,因为网格单元本质上变得不再像长方体了(例如,当内角偏离90度时)。
关于图4和图5提到的一种产生结构化网格的方法是“柱状网格化”,其中该断层可以是光滑的表面,并且网格在三个网格方向上具有指定数量的单元。这种网格化方法要求简化复杂地质环境中的断层(例如,断层竖直会聚;“Y”或“λ”断层),并且可能导致大量的单元畸变和大范围的网格单元大小,后者对于网格属性建模是有问题的。
阶梯网格是一种结构化网格,其倾向于非常适合于储层模拟应用。网格单元的边和K坐标线是竖直的,并且网格层的几何形状由这些“柱”下方的角点标识。虽然阶梯状网格允许对更复杂的地质环境进行建模,但断层在横向和竖直方向上都是阶梯状穿透该网格的,这会降低地质储层建模的准确性。处理这种精度的一些方法可能会在断层位置定义裂隙单元(在地质空间中),目的是准确表示断层表面和断层-层位几何形状。
缺乏整齐的坐标和结构化的号惯例的网格可以被称为非结构化网格。建模过程的输入数据允许在地质空间中创建体积模型,然后再进行地质力学转换为沉积空间,在该沉积空间中,合适的层位是水平面,并且断层偏移已被移除。网格化可以在沉积空间中发生,其中可以使用沉积空间中的切割单元来精确表示断层。将逆变换应用于沉积空间中的网格会导致地质空间中的网格(沉积网格),其在遵守地层学的同时可以更准确地表示原始体积模型。在地质空间中,这种非结构化的沉积网格的沉积特性的精确模型取决于对沉积方向的了解,假定网格属性(特别是方向属性、例如渗透性)起源于这些沉积方向。以下工作流程涉及使用沉积映射来分析沉积网格并更准确地建模和模拟方向属性:
●结构模型的生成、沉积空间转换,在沉积空间中的网格以及在地质空间中的沉积网格;
●利用沉积-地质空间映射来估算沉积网格单元的局部轴;
●分析地质空间中每个网格单元的局部轴取向和尺寸,以了解由输入数据中的错误引起的潜在局部网格错误;和
●网格轴的正交化和沉积网格中方向属性的表示,供储层模拟器使用。
如上所述,在各种框架中对地下的建模和模拟倾向于采用结构化网格。这样的网格包括端到端工作流程的不同阶段的近似值。
如所提及的,结构化网格的替代可以是非结构化网格。作为示例,请考虑非结构化的沉积网格,以及对地质输入进行精确建模以及更准确地模拟此类沉积网格的好处,例如,通过避免或减少使用结构化网格时进行的简化量。
沉积网格在概念上生成在图6中示出为使用初始结构模型、从该地质模型到沉积空间的转换、根据(正交)u、v和w方向的沉积空间中该模型的划分以及最后转换回地质空间以定义断层的沉积网格。到沉积空间的转换可以定义结构模型(地质空间)中的点与沉积空间中相应模型中的点之间的一对一映射。沉积网格断层和层位可以精确地匹配原始结构模型的断层和层位,而地震层位之间的附加网格层则符合沉积映射的地层性质。在地质空间中观察时,w坐标线与沉积空间中强制的局部地震层位保持正交关系。
关于一对一映射,可以是,在地质空间中精确位于结构模型的断层表面上的点将对应于沉积空间中的两个点(或在断层-断层相交处更多),例如该点将分别视为位于该断层的两侧的每个上。因此,断层表面上的每个点都可以由几个并置的点表示,这些点具有唯一的沉积空间位置。
图7显示了非结构化沉积网格中两个相邻单元的图形710和730的示例,其可能会在多个公共单元面相遇。在断层和断面处(例如,以及对于不在断层或断面处的两个单元),可以保证面的这种匹配。
作为示例,一种方法可以包括对作为地下区域的模型的网格的属性进行隔分,其至少部分地基于地震数据,例如,由包括地震能量发射器(例如,震源)和地震能量接收器(例如,传感器)的获取设备的分布式组件获取的地震数据。这样的方法可以包括操作模拟器,该模拟器至少部分地基于隔分,生成可以在地下区域中发生的一种或多种类型的物理现象的结果。作为示例,模拟器可以是计算系统的一部分,其包括接口,该接口可以向一个或多个现场设备发出信号,以指示现场设备在现场执行动作,其可实现如模拟器指示的结果或实现以其他方式至少部分基于由模拟器指示的物理现象的结果。
作为示例,一种方法可以包括将单元的网格中的各个单元分配给两个类别之一,例如,在单元的网格(例如,网格)中表示的感兴趣的断层的上盘(hangingwall)侧或下盘(footwall)侧(例如,“正”侧或“负”侧)。可以将这种分配作为侧面特性(例如,对应于断层的侧面的特性)存储到存储装置。这样的方法可以包括在遵守给定的一个或多个断层截断规则的同时将侧属性组合成单个隔区属性。
作为示例,一种方法可以包括实现确定单元所属的断层的一侧的技术。作为示例,一种方法可以利用辐射源方法,其中可以在网格内人工施加辐射以照亮网格的单元。在这样的技术中,可以为单个单元生成照明值。在这样的示例中,照明值可以根据单元的断层的侧面而不同。例如,考虑照明在断层的一侧为红色,在断层的另一侧为蓝色;请注意,可以使用各种其他表示方式来区分一侧的照明与另一侧的照明(例如,通过不同的颜色或另一种比例尺)。作为示例,可以利用强度值。例如,考虑范围为0至255(例如,具有特定位深度等)的红色强度,并考虑范围0至255(例如,具有特定位深度等)的蓝色强度。在这样的示例中,可以将值呈现给显示器。
图8示出了网格的一部分的图形用户界面(GUI)800的示例,该网格包括用空心圆划定的断层,其中辐射技术用红色照亮一侧,而用蓝色照亮相对的另一侧。在这样的示例中,可以将断层视为一个或多个辐射源(例如,辐射不同类型的辐射的辐射器,第一种类型辐射到断层的一侧,第二种辐射到断层的另一相对侧)。作为示例,断层可以被认为是利用这种技术的辐射类型的结构(例如,天线、透镜、主体等)。在这样的示例中,断层的表面形状(例如,多维)可以在特定一个或多个方向上引导辐射,从而照亮单元(例如,可以为单个单元产生照明值)。作为示例,根据这样的技术,由断层发射的辐射可以随着距断层的距离而衰减。例如,可以利用基于距离的衰减函数,该函数可以在现实世界的物理中可选地具有类似物(例如,关于由于一种或多种影响(例如,体积增大、面积增大、吸收、衰减等)引起的来自源的能量的衰减)。
如图8的GUI 800中所示的,一种方法可以将光源分配给每个断层面,并使用“正”或“负”颜色的光照射其各自的体积(例如“半球”或其他体积空间)。作为示例,可以计算每个单元中心的累积照明量以确定单元的域颜色。例如,考虑图8的GUI 800的下部空心圆,在这种空心圆上,来自两侧的照明可能会落在(或落入)一个或多个单元上。在这样的单元中,可以存在两个值,红色值和蓝色值,其中可以基于将红色值的大小与蓝色值的大小进行比较来确定是否要分配特定的单元到断层的一侧或要分配给断层的另一侧来进行决策。例如,可以使用这种方法来确定一单元是否属于断层的正侧或负侧,从而为该单元创建一个“侧属性”。
这样的辐射/照明方法的简单实施可能执行得很差(例如,对于中型模型而言,大约为几个小时)。作为示例,为了提高执行辐射/照明方法的计算装置或计算系统的性能,可以实现基于细节水平的技术。在这样的示例中,可以实现运行时间的速度增加(例如,对于中型模型降低到几秒钟)。在这样的示例中,提高了计算装置或计算系统的速度,其可以是用于分析地震数据以生成地球的改进模型的总体装置或系统的一部分(例如,如通过地震层析成像扫描所成像的)。作为示例,可以将改进的地球模型作为改进的地球图像呈现给显示器(例如,通过单元、像素、线、节点、层等)。
作为示例,一种方法可以包括为辐射/照明方法实施近断层调整。这样的调整可以帮助提高关于单元分配的确定性。例如,为了解决由光近似方法引起的近断层不确定性,一种方法可以包括经由近断层溢流填充进行调整。在这样的示例中,单元距离内的单元(例如,一个单元远离、两个单元远离等)可以被“溢流(flood)”以填充可以是高值的值。这样的方法可以用一种辐射(例如红色)将单元溢流填充到断层的一侧,并用另一种辐射(例如蓝色)将单元溢流填充到断层的另一侧,其中此类单元满足一个或多个标准(例如,在断层等两个单元内)。
作为示例,一种方法可以包括生成侧属性,然后将侧属性(例如,每个断层一个)组合成单个隔区属性。这种方法可以考虑断层的一个或多个截断,其可以是例如用户定义的截断(例如,通过交互式GUI等设置)。作为示例,截断可以限制使用断层将空间划分为单独的隔区的区域。
作为示例,可以为每个单元分配与模型中感兴趣的断层一样多的侧属性。作为示例,可以通过将侧属性串为位(例如,真/假)并将结果读取为整数来产生隔区号。但是,这种方法可能会生成过多的隔区,使得它们需要根据一个或多个截断规则进行合并。作为示例,一种方法可以实现多维溢流填充以至少部分地执行侧属性分配。
作为示例,由于不一定容易获得断层截断,因此一种方法可以包括分析网格中的断层的配置以产生建议的截断列表,例如,用户可以将其用作隔区合并过程的输入。作为示例,考虑一个过程,其中收集直接附接到断层的单元。从该收集中,过程可以将断层-断层相互作用识别为各种不同类型的截断中的一种或多种(例如,正、负、交叉或中继)。然后,这样的过程可以利用侧属性来找到断层扩展与实际断层相互作用的情况。在这样的示例中,当应用一组排除规则时,该过程可以产生良好的起点(例如,对于用户,计算装置、计算系统等)。
作为示例,一种方法可以包括通过辐射/照明方法划分空间,作为现实世界边界的子集的扩展(例如,如地震和/或其他数据所示)。
作为示例,一种方法可以在地下区域的各种模型上可靠地实现隔分,其中这种模型可以至少部分地基于诸如地震数据的数据。作为示例,一种方法可以提供用户控制,例如,允许通过交互式GUI进行输入,该交互式GUI可以通过选择感兴趣的断层并定义一个或多个截断规则来控制计算机算法。如所提及的,一种方法可以减少计算需求,使得可以实时地(例如,从输入到基于其的结果呈现的秒的量级)处理经由GUI的交互。作为示例,一种方法可以提供足够低的运行时间,以允许进行交互式工作流程,其中用户可以迭代地修改一个或多个输入,重新运行计算机算法并重复执行直到获得期望的结果。
作为示例,一种方法可以包括利用断层辐射来确定单元在每个断层的哪一侧,并且可以包括利用侧面性来产生隔区特性,使得位于断层不同侧的两个单元(例如,在该区域具有影响力)分配了不同的隔区号。
关于侧属性的创建,可以利用断层辐射。当呈现到显示器上时,这种方法可以产生视觉上可评估的结果。例如,考虑图8的GUI 800,其中红色和蓝色区分断层的一部分的侧面,如由空心圆所划定的(例如,截断指示符等,其可以通过由GUI促进的输入来交互地进行调整,例如例如通过触摸、鼠标、手写笔等)。这样的方法可以促进关于一个或多个单元及其侧面性的分配和/或进行决策。例如,可能存在这样的区域,其中决定一单元属于断层的一侧还是另一侧的决策充其量是不明确的。如图8的GUI 800中的方法减轻了这种不明确性,因此通过计算技术改善了这种决策作出。
作为示例,一种方法可以包括以大约180度角(例如,在从大约150度到大约210度的范围内)将断层的表面近似为一系列定向聚光灯(每个断层面一个)和余弦衰减(例如,或其他衰减,其是用户可以选择的)。作为示例,可以确定光的强度与距离的平方下降(例如,对于距离d,辐射强度Ir,减小为Ir/(d2))。
如所解释的,每个断层可以具有两侧。作为示例,聚光灯可以分配给“正色”的正侧发射,而聚光灯可以分配给“负色”的负侧发射。可以针对每个单元(例如,对于每个单元中心或其他几何/空间度量)计算这些聚光灯的累积效果。作为示例,可以利用主要的“颜色”来确定单元落在哪一侧(例如,单元位于具有第1侧和第2侧的断层的第1侧或第2侧)。
作为示例,可以以串行方式和/或并行方式来实现一种方法。关于并行方式,复杂度可以是O(m*n),其中m是断层数,n是单元数,这可能导致较长的运行时间。如所提到的,可以实现一种减少运行时间的技术(例如,减少计算资源需求等)。作为示例,可以实现断层辐射加速方法。为了确保方法具有高性能(并因此提高了其实用性),可以使用详细程度方法来加速断层辐射处理。
再次参考图8的GUI 800,在视觉上,在网格的大多数区域中,单元的侧面性是相对明确的,并且可以容易地通过粗略近似来确定。例如,考虑一种近似方法,该方法是:将网格边界框拆分为24x 24x 12分辨率的常规扇区;确定每个聚光灯属于哪个扇区;计算每个扇区中的单个累积聚光灯作为其中每个聚光灯的总和,例如,其中累积光源发出的光量与单独光源组合的光量大致相等(在这样的示例中,其中根据一项或多项标准未能充分满足条件,则可以进行一项或多项调整);使用粗糙光源计算每个扇区中的粗糙照明;并找到要进行高分辨率计算的扇区。例如,例如,这样的扇区可以是相邻扇区具有与其自身不同的主要“侧颜色”的扇区。这样的方法可以在整个期望的网格体积中准备粗略的照明近似值,并确定要进行细粒度计算以进行其他决策的区域(例如,关于期望网格体积内至少一部分单元的侧面性)。
作为示例,一种方法可以包括处理最初落入已经被标记用于高分辨率计算的扇区的各个单元。这种方法可以包括以下过程:计算该扇区和相邻扇区中的光源对单元中心的照明;并使用基于扇区的粗略光源近似来增加来自其他扇区的照明。在这样的示例中,可以将落入未被选出用于高分辨率计算的扇区中的单元分配粗略的扇区照明。这种方法可以将运行时间从几小时减少到几秒钟。这样的方法可以减少运行时间,而在单元侧面性的分配上没有实质性(例如,视觉上可见的等)差异。
如所提及的,一种方法可以采用近断层溢流填充。作为示例,可以将断层面视为区域光源。但是,这种方法相对于网格在计算上可能并不实用。因此,作为示例,可以利用作为点光源的近似值。这样的近似可以以接近断层的方式使用,因为离断层较远,通常可能不再需要感兴趣该特定问题。具体地,在断层的相对较近附近中,可能会出现可能的错误,例如,单元中心落入附近点源(例如,特别有用的点源)的阴影中,从而有可能被错误分配,因为来自断层另一侧的光源比合适的光源对那个单元的贡献更大。
作为示例,可以使用近断层溢流填充,其中可以将与断层直接相邻的单元定义为断层距离等于零(例如1或其他适当值)的“种子”。从这样的种子,例如,可以通过直接将它们分配为属于“正”或“负”侧(例如,使侧颜色的强度最大)来对这些种子的直接相邻实施溢流填充。作为示例,在相邻单元(例如,种子)具有零值的情况下,可以以2的溢流填充深度(关于相邻种子)来实现溢流填充方法。
作为示例,可能存在具有与断层的两侧相互作用的面的单元。在这样的示例中,可以将“双重分配”的单元放在一边,而不通过溢流填充过程对其进行处理。
关于合并到隔区中的侧属性,一种方法可以提供将针对一个或多个断层的断层侧分配给各个单元的侧属性。作为示例,考虑一种利用单元绑定位数组的方法,其中每个位代表一断层侧。在这样的示例中,如果该位为“真”,则表示该单元位于相应断层的“正”侧;相反,如果该位为“假”,则表示该单元位于相应断路层的“负”侧。
这样的位数组可以被解释为产生一个隔区号的整数(例如,以二进制格式)。根据需要,可以对这些值重新号以产生顺序号的过程。
关于断层截断规则,用户可能期望某些断层具有有限的影响,并且存在断层扩展不会影响隔分的区域。这样的期望可以表示为一个或多个截断规则。截断规则的示例是:“断层-A在断层B的上盘壁(正侧)被截断”。这意味着用户期望断层-A对断层-B的下盘(负侧)没有影响。
图9A示出了GUI 910的示例,该GUI 910示出了断层-A(绿色)被断层-B(黄色)截断的情况。图9B和9C分别示出了两种情况下的GUI 920和930的示例。在GUI 920中,该情况涉及未合并的隔区,其中未应用截断规则,并且断层-A跨断层-B延伸;然而,在GUI 930中,该情况演示了预期的输出,其中应用了(一个或多个)截断规则,并且断层A对断层B的另一侧没有影响。根据截断规则,图9A和9B的两个情况对于隔区可以称为未合并(GUI 920)和合并(GUI 930)的。
图9A、9B和9C的示例以相对简单的结构(例如,断层几何形状)说明了这种情况;请注意,可能会出现更复杂的情况,并包括更多的断层,其中截断定义可能会更复杂。作为示例,可以将截断规则生成器包括为处理网格的单元以确定侧面性的框架的一部分。这样的生成器可以是框架的GUI的一部分,在此用户可以交互以选择和/或定义将应用于包含断层的地下区域模型的一个或多个部分的截断规则。作为示例,模型可以伴随有预定的截断或截断规则,可以选择性地采用该截断或截断规则。作为示例,一种方法可以包括截断检测,使得截断被自动检测和/或生成。
作为示例,可以实施隔区的合并计算算法以执行合并,其可以是多维溢流填充的一种,利用截断规则来确定要向哪些隔区分配特定的隔区号。这种方法可以包括一次扰动一位,并根据截断规则确定相应的断层在该区域是否有影响。
作为示例,用于合并的计算算法可以操作以:将现有的侧属性组合(例如,上述位数组)插入“输入集”中;从“输入集”中获取第一“位数组”,并将其插入“种子集”中;增加“隔区号”;将来自“种子集”中的元素处理为“种子”;用当前的“种子”分配当前的“隔区号”(例如,使用“位数组”到隔区号图);检查是否存在单个截断规则,该截断规则会“遮盖”相应的断层(在被扰动位的位置);如果是这样,则将“相邻”插入“种子集”;检查相邻后,从“种子集”和从“输入集”中删除当前的“种子”;如果“种子集”不为空,则算法继续将来自“种子集”中的元素处理为“种子”;如果“种子集”为空,则算法将继续获取一位数组(请参见上文);以及如果“输入集”为空,则算法可以终止。这种方法可以生成从多个现有“位数组”中的每一个到顺序的隔区号的映射。在这样的示例中,一种方法可以包括基于相应的侧属性,将隔区号与模型的至少一部分中的多个单元中的每个单元相关联。
关于截断检测,断层截断定义可以是方法的输入,使得它可以通知隔区合并计算算法在何处断层要具有影响和在何处它要被忽略。作为示例,截断配置集可能会基于针对用户的工作流程和对模型地质的更深入分析而有所不同。因此,用户可能拥有最终批准,并且可以进行判断以影响最终的隔分。如所提及的,一个或多个GUI可以被呈现到允许这种用户交互的显示器,其可以以这样的方式来执行,其中计算装置或系统响应于近实时地输入以执行一个或多个计算,并且将结果呈现给显示器(例如,参见图8、图9A、9B和9C的GUI)。
各种截断可以是明确的,并且可以以编程方式(例如,自动地)被检测以例如为用户提供良好的起点。作为示例,可以实现一种用于从网格本身产生初始截断集的方法(例如,可选地,没有附加信息)。
作为示例,考虑一种方法,该方法包括检测直接的断层-断层相互作用并将这种检测到的相互作用分类。这样的方法可以包括检测在何处断层扩展与另一断层相交。检测和分类的结果可用于剔除一些有关断层扩展的不必要检测。在这样的示例中,输出可以是用户可以根据需要(例如,经由GUI等)查看并可选地调整的两组截断。
关于面分配,过程可以包括迭代断层并收集与每个断层侧相邻的单元。在这样的示例中,可以排除发现与断层的两侧相邻的单元,并且不参与进一步的处理。作为示例,单元集可用于检测断层-断层相互作用。例如,一种计算算法可以针对其他断层测试每个断层,并且可以计数一个断层的单元集中有多少个单元可以在另一个断层的单元集中找到。这样的计数可以用于确定截断的类型。
作为示例,断层可以具有两个单元集,“正集”和“负集”。在这样的示例中,来自一个断层的正集可以与另一个断层的“正集”和“负集”相交。可以为第一断层的“负集”实现这种方法。这种方法的应用可以例如产生四个整数,指示相交单元的数量。
关于可能的配置,请考虑以下因素:正-正和正-负交叉点为非零->“正截断”;负-正和负-负的交叉点为非零->“负截断”;交叉点计数非零->“交叉断层”;一个交叉点计数为零->“中继断层”;和其他配置->“无截断”。
作为示例,为了减少计算需求,一种方法可以采用哈希集,该哈希集可以加快相交计数。这样的方法允许保持算法复杂度为0(n*m),其中n是与所考虑的断层相邻的单元的总数,而m是断层的数量。注意,通过隔区合并计算算法可以随后使用“正截断”和“负截断”。作为示例,可以忽略中继和交叉截断。
关于断层扩展截断检测,考虑一种方法,该方法包括检测断层扩展是否与断层相交。这种方法可以包括将在“面分配过程”中产生的信息与例如“断层侧特性”产生过程相结合。这样的方法可以再次针对其他感兴趣的断层测试断层。
作为示例,一种方法可以包括,对于所考虑的每个断层-断层迭代,针对与另一断层相邻的每个单元查找一个断层的“侧属性”值,并计算每侧上的单元数量。在这样的示例中,如果它们都不都是零,则检测到截断。换句话说,如果在另一个断层的两侧都找到与一个断层相邻的单元,则认为一个断层被另一个断层截断了。注意,断层扩展检测过程可能无法区分截断发生在哪一侧。在这种情况下,可以实施截断侧检测处理。
关于截断侧自动检测过程,隔区的合并过程可以利用具有限定的侧面性(例如,截断的断层在截断器断层的哪一侧结束)的截断。在各种情况下,对于用户而言,可能难以了解截断的侧面性。作为示例,例如,只要在将截断传递到隔区合并过程之前可能无法使用截断的侧面性,就可以实现自动检测该侧面性的过程。
关于侧面性检测,其可以包括将截断器的“侧属性”与截断的断层的“相邻单元”组合。这样的过程可以包括基于可以找到大多数相邻单元的位置来分配侧面性(例如,如果被截断的断层的大部分“相邻单元”发现位于截断器的“正侧”,则将截断标记为“正截断”,反之亦然)。
关于截断排除,断层扩展截断检测过程可能趋于过于包容,检测并添加了一些可能被认为多余的截断。作为示例,一种方法可以包括一过程,该过程可以通过一个或多个规则来帮助解决这种过度包含的问题,例如:如果断层扩展过程再次检测到了截断(例如,之前由面分配过程发现的“截断的截断器”对),则不会再次将其添加到结果集中。这样的方法可以认为面分配检测到截断类型,而断层扩展过程可能没有。诸如“中继”和“交叉”之类的截断类型可能会被隔区合并过程忽略。作为示例,经由断层扩展过程再次添加截断可能会迫使它们稍后被自动检测为“正”或“负”截断,例如最终被合并过程利用。
作为示例,在另一个断层的相对侧上发现两个断层的地方,可以忽略两个断层之间的一个或多个截断。这种方法可以帮助减轻多余的截断,这可能会使检测到的截断更加难以阅读。
图10示出了方法1000的示例,该方法1000包括用于访问地下区域的模型的访问块1012,其中该模型包括与地下区域中的断层相对应的断层;辐射块1014,用于经由断层的第一侧用第一类型辐射从其向外辐射该模型的单元,并经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;分配块1016,用于至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,其指示将所述单元的一部分中的每个单元分配给所述断层的第一侧或分配给所述断层的第二侧;以及模拟块1018,用于至少部分地基于至少一部分单元和为其分配的一个或多个属性,经由模拟器来模拟地下区域的物理现象。这样的方法可以改善模拟器的性能和结果。例如,在以更准确的方式将属性分配给包含断层的模型中的单元的地方,可以更准确地和/或以产生解的方式来模拟物理现象,由于收敛速度等原因,这可能会减少时间量。作为示例,可以通过图10的方法1000更快地准备模型以输入到模拟器,例如,该方法的某些部分(例如,块1014和1016)可以在工作流程期间以交互方式操作使得用户可以更准确、快速地准备模型及其属性,以输入到模拟器。
作为示例,可以通过将断层的一侧视为表面来执行向外辐射,其中该表面可以是辐射源,使得辐射从该表面发出。作为示例,辐射可以在一定角度范围内发出。作为示例,可以沿法向于代表断层的一侧的表面的方向发出辐射。作为示例,其中,断层局部包括多个法线,每个法线可以是辐射的方向。作为示例,辐射可以是静态的,因为它是存在于断层一侧的辐射模式,其形式可以是空间函数、覆盖区等,从而可以确定扇区和/或模型中的单个单元的值。作为示例,可以采用一种或多种辐射技术从断层的一侧向外辐射辐射,以使单元受到辐射,例如以确定侧面性、隔区等,这些特性可用于模拟器可以利用的属性的属性分配的目的。
如所提到的,一种方法可以包括为具有断层的体单元网格生成隔区。隔分可以是储层建模工作流程和/或工程工作流程的输入,例如,以按单元建模在储层的不同部分分配不同的属性。当定义油气和/或油水接触和/或流体性质时,可以使用这种方法,其可以跨越地质边界以用于估计储量、预测产量、对现场操作的响应等。
如所提及的,隔分的特征可以是在感兴趣的断层的相对侧上的两个单元被分配了不同的隔区号(例如,同时可选地将隔区的总数保持为最小)。可以以一种可以考虑一个或多个情况的方式执行这种隔分,在这种情况下,可以预期断层可能会“扩展”到其实际物理范围之外。作为示例,可以实现与截断相关联的一个或多个过程,其中,例如,一断层被一个或多个其他断层截断;注意这种截断可以分开断层(例如分成两个或更多部分)。如所提及的,一个或多个截断过程可以是交互的,因为可以将GUI呈现到显示器,以允许用户输入来评估,制作,调整等一个或多个截断。作为示例,GUI可以允许用户选择感兴趣的断层并评估,制作,调整等断层的一个或多个截断。
作为示例,一种方法可以包括断层辐射,用于断层辐射的分层方法(例如,以粗糙和精细层加速)和/或出于侧面性确定(例如,侧面性属性分配等)的目的而进行近断层溢流(例如,近断层溢流填充)。作为示例,一种方法可以包括合并侧面性属性信息以形成隔区(例如,隔区属性),这样的过程可以包括利用一个或多个断层截断规则和例如计算的隔区合并算法。作为示例,一种方法可以包括截断检测,该方法可以包括面分配过程、断层扩展截断检测过程、截断侧自动检测过程和/或一个或多个截断排除过程。作为示例,工作流程可以包括如上所述或本文其他方式的一种或多种方法和/或一种或多种过程。
在图10中示出了与各种计算机可读介质(CRM)块1013、1015、1017和1019相关联的方法1000。这样的块通常包括适合于由一个或多个处理器(或内核)执行以指令计算装置或系统执行一个或多个动作的指令。尽管示出了各种块,但是单个介质可以配置有指令以至少部分地允许方法700的各种动作的执行。作为示例,一个或多个块1013、1015、1017和1019可以是指令,例如图2的系统250的一组或多组指令270。
图11示出了方法1100的示例,该方法1100包括地震解释块1110,该地震解释块1110用于至少访问由感测在地下环境(例如,地下区域)中传播的地震能量的一个或多个传感器获取的地震数据,用于至少部分地基于数据的至少一部分来执行地震解释并生成信息(例如,结构的位置等)。方法1100还包括用于生成表示地下环境中的物理结构的一个或多个结构(例如,模型内的结构等)的结构框架块1120。例如,解释可以识别出层位,并且结构模型可以将该层位包括为基于数据的结构,如果地震数据发生变化(例如,迹线的变化等),则该层位实际上可能存在于地下环境中。方法1100包括处理块1130,处理块1130可以实施一种或多种技术,其可以增强如何处理属性,其中这样的属性可以包括一个或多个可以根据属性块1140分配给网格的一部分的属性。例如,将隔分考虑为一种基于诸如方法1000的部分的方法(例如,块1012、1014和1016)向多个隔区分配隔区标识符的技术。这样的方法可以利用一种或多种类型的辐射作为技术的一部分,该技术有助于识别哪些单元要分配特定属性,该属性可以是模拟器利用的属性(例如,参见块1018),以生成关于一种或多种物理现象的模拟结果。
如方法1100中所示的,过程1130可以包括多个层1132和1134,其可以包括粗糙层和精细层,其中在粗糙层1132中利用了扇区;然而,在精细层1134中,单个单元例如针对不通过测试的特定扇区中的单元逐个单元地被利用。例如,决策块1133可以提供是否决定将一个扇区传递给精细层1134。如上所述,在一个扇区的一个或多个特性与相邻扇区不同的情况下,可以根据该决策块1133来选择该扇区,用于通过层1134在逐个单元的基础上进行处理。如所示的,可以将层1132和1134的结果传递到该过程,该过程可以利用该结果进行属性分配,例如,经由隔分(例如,顺序号的隔区,其对应于可能被一个或多个断层物理隔离的区域)。
如所提及的,一种方法可以包括结构模型的生成、沉积空间变换、沉积空间中的网格以及在地质空间中的沉积网格。作为示例,可以将隔区定义为适合于输入到模拟器的沉积网格的隔区属性。
作为示例,一种流体流动模拟,其可以通过一种或多种数值方法(例如,有限元、有限差分等)进行,其中一个或多个表示物理现象的方程式(例如,Navier-Stokes、Darcy等)可以通过适当的“连接”(例如,单元间边界条件等)进行指定,以使流体流动模拟的结果比没有这样一种或多种特性表示更为可行和/或更准确,这可能说明了诸如断层之类的结构特征。隔分方法可以使模拟器对于正在模拟的给定储层更稳定,和/或使模拟器的结果更准确和/或使模拟器以更少的时间和/或更少的计算资源(例如,更少的处理器能力、更少的性能等)运行。作为示例,取决于方程的数量(例如,可能与网格大小等相关联),模拟运行可能会很漫长。作为示例,模拟运行可能需要数小时或例如数天。可以理解,在缺乏稳定性的情况下,占用数小时执行模拟运行的计算机资源可能不会产生结果或不会产生有意义的结果。
较低的运行周转时间可以增加分配给预算时间段的储层研究的收益。必然地,花费在重复运行中以对抗模型的不稳定性或时间步长可能会适得其反。各种因素都会影响运行时间,例如,运行时间等于(CPU时间/步长)×(时间步长数)乘积。对于隐式公式,第一个因子趋于较大,而第二个因子则趋于较小,而对于IMPES公式,则相反。IMPES是一种条件稳定的公式,涉及时间步长Δt<Δt*,以减少振荡和误差增长的风险,其中Δt*可定义为最大稳定时间步长。对于IMPES,条件稳定性源自对块间流项中非压力变量的明确处理。
模拟器可以利用预处理器,例如,考虑嵌套分解(NF)和不完全LU分解[ILU(n)]。术语“LU分解”是指将矩阵A分解为下三角矩阵L和上三角矩阵U的乘积,这在计算上可能是昂贵的,并且涉及高斯消除。术语“ILU(n)”表示不完整的LU分解,其中允许有限的填充,并且n是填充的顺序。
NF可以在特定条件下适当地执行,例如当与网格中的特定方向相关联的传递率在整个网格中均匀地主导其他方向上的传递率时。ILU(n)或红黑ILU(n)往往比NF对块的排序和主要传递率方向的空间变化不那么敏感。
传递率可以是地层相对于流体的电导率的量度。传递率可以由特定的模拟器以特定的方式定义。作为示例,可以相对于粘度定义传递率(例如,针对粘度进行调整等)。在本文中提到了关于传递率的一些示例。作为示例,当执行模拟运行时,模拟器可利用传递率,例如关于可操作地耦合到储层的储层和/或设备(例如,导管等)中的流体运动(例如,经由流体连通等)。
INTERSECT模拟器使用可以实现预处理的计算求解器,其中这种预处理可以涉及将方程组代数分解为子组,该子组可以基于它们的特定特性来处理以促进求解。可以通过迭代技术对所得储层方程进行数值求解,直到整个方程组达到收敛为止,这可以说明一个或多个井、一个或多个地面设施等。
INTERSECT模拟器框架可以提供关于现场操作的现场任务,该现场操作可以包括关于地面设施的操作。关于工作流程的示例,请考虑一个模拟,该模拟说明了100多口生产井,并且在非结构化网格中包含数百万个单元,其可以在结构上建模10,000以上的裂缝。作为示例,运行执行期间的时间可以包括评估残差和组装雅可比矩阵、线性求解器、非线性更新和非线性收敛测试等任务。在这样的示例中,时间可能主要花费在线性求解器上,然后评估残差和雅可比矩阵的组装。
图12示出了方法1210的示例,该方法包括用于计算孔体积、传递率、深度和NNC的计算块1220、用于初始化和计算初始饱和度、就地压力和流体的初始化和计算块1240,以及定义和时间进程块1260,用于定义一个或多个井和地表设施并例如通过各个单元的物料平衡(例如,以一个或多个井作为各个汇和/或源)在时间上前进。在这样的示例中,块1220可以包括用于生成、调整等一个或多个网格的一个或多个指令集,其可以表示一个或多个断层、断层的一部分等。
关于初始化和计算块1240,对于初始时间(例如,t0),可以设置地质环境的网格模型内的饱和度分布和地质环境的网格模型内的压力分布以表示例如平衡相对于重力的状态(例如,静态或“无流动”状态)。作为示例,为了近似平衡状态,可以执行计算。作为示例,这样的计算可以通过一组或多组指令来执行。例如,可以实现地震模拟框架、储层模拟器、一组专门的指令等中的一个或多个以执行一个或多个计算,这些计算可以旨在近似或促进近似平衡状态。作为示例,储层模拟器可以包括一组指令,用于使用数据进行初始化以计算毛细管和流体梯度,以及由此代表地质环境的网格模型的各个单元中的流体饱和密度。
作为示例,储层模拟器可以及时推进。作为示例,可以实现可以为单个时间增量(例如,时间点)生成解决方案的数值求解器。作为示例,求解器可以实现隐式解决方案和/或显式解决方案,注意隐式解决方案可以比显式方案允许更大的时间增量。可以在“时间表”中列出所需解决方案的时间。例如,时间表可以包括较早时间段的较小时间增量,然后是较大的时间增量。
求解器可以实施一种或多种技术来帮助确保稳定性、收敛性、准确性等。例如,当及时推进求解时,求解器可以实现时间的子增量,但是,增加增量数量会增加计算时间。作为示例,例如,基于一个或多个先前增量的信息,可以使用可调增量大小。
作为示例,数值求解器可以实现有限差分方法、有限元方法、有限体积方法等中的一种或多种。作为示例,储层模拟器可以实现空间近似的中心差异和时间上的前向差异。作为示例,代表网格单元和相关方程的矩阵可能是稀疏的、对角带状的和封闭的,以及包括非对角线的条目。作为示例,可以实现INTERSECT模拟器。
作为示例,求解器可以实现隐式压力、显式饱和(IMPES)方案。这样的方案可以被认为是显式和隐式技术的中间形式。在IMPES方案中,饱和度显式更新,而压力隐式求解。
关于质量守恒,可以指定网格单元模型的各个单元中的值(例如,对于水、天然气和石油而言)以求和为一定值,该值可以被认为是质量守恒的控制标准。由于黑油方程趋向于是非线性的,因此可以实施牛顿-拉夫森(Newton-Raphson)类型的技术,其包括确定导数、迭代等。例如,可以通过根据牛顿-拉夫森方案进行迭代来找到解决方案,其中这样的迭代可以被称为非线性迭代、牛顿迭代或外部迭代。在满足一个或多个错误标准的情况下,求解过程已收敛,并且找到了收敛的解决方案。因此,在牛顿迭代中,通过执行多个线性迭代来解决线性问题,这可以称为内部迭代。
作为示例,一种解决方案可以由以下伪算法表示:
作为示例,求解器可以执行多个内部迭代(例如,线性的)和多个外部迭代(例如,非线性的)。作为示例,内部迭代的数量可以是在外部迭代内的大约10到大约20的数量级,而对于单个时间增量,外部迭代的数量可以是大约十或更少。
如所提及的,可以至少部分地基于模拟结果来修改(例如,调整等)网格,该模拟结果可以可选地包括诸如收敛行为的结果。例如,在可能会改善收敛的情况下,可以对网格进行一个或多个调整,其中这样的一个或多个调整可以允许收敛、加快收敛等。
图13示出了可以观察到各种结构特征的地质环境的地质1310的示例。一个或多个这样的特征可以存在于地下环境(例如,地下区域)中,该地下环境可以通过一个或多个网格(例如,一组单元等之一)来表示。
作为示例,一种方法可以包括访问地下区域的模型,其中该模型包括与地下区域中的断层相对应的断层;经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元;以及经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;并且至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该值表示单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。
作为示例,一种方法可以包括将模型分成多个扇区,其中每个扇区包括该模型的单元的一部分,并且可以向该单元的该部分分配侧面性属性值。例如,可以为整个扇区分配共同的侧面性属性。在这样的示例中,可以为一个或多个扇区分配公共的隔区属性。作为示例,一种方法可以包括将多个隔区属性分配给多个扇区,其中至少一些隔区包括多个扇区。作为示例,隔区可以包括经由粗糙层被分配了侧面性属性的单元,并且可以包括经由精细层(例如,参见图11的层1132和1134)被分配了侧面性的单元。
作为示例,一种方法可以为对应数量的单元提供多个分配的侧面性属性值,并分配至少一个隔区属性值。在这样的示例中,至少一个隔区属性值中的每一个可以是或包括标识符,该标识符是模型内对应的空间定义的隔区的一组顺序标识符的成员。
作为示例,一种方法可以包括通过另一断层来检测断层的至少一个截断,其中,分配至少一个隔区属性值至少部分地基于该检测。例如,可以至少部分地基于检测到的另一层断层对断层的截断,为断层的一侧分配至少两个隔区属性值(例如,其中该断层具有延伸到截断的断层的一侧的部分,以及其中该断层在将其截断的断层的另一侧的一部分)。作为示例,截断可以是断层的划分,以便对断层进行分段。
作为示例,一种方法可以包括识别由具有主要侧面性属性值的一部分单元构成的第一扇区,其中主要侧面性属性值不同于由不同部分的单元组成的第二扇区的主要侧面性属性值,其中第一扇区和第二扇区中的至少一些单元是相邻单元。在这样的示例中,一种方法可以包括在第一扇区中的每个单元上执行辐射计算,以向每个单元分配侧面性属性值(例如,考虑精细层等)。
作为示例,一种方法可以包括执行分层分配,其中一层基于单元扇区来分配单元侧面性属性值,而另一层基于单个单元来分配单元侧面性属性值(例如,参见图11的层1132和1134)。
作为示例,第一类型的辐射和第二类型的辐射可以相对于距断层的距离衰减。例如,考虑衰减或其他函数的能量类型。
作为示例,一种方法可以包括向显示器呈现第一颜色的第一类型的辐射和第二颜色的第二类型的辐射。作为示例,模型可以包括多个断层,并且一种方法可以包括对于每个单独的断层,通过单独的断层的第一侧,用相应类型辐射从其向外辐射该模型的单元,并且经由单独的断层的相对的第二侧,从那里向外辐射具有相应的不同类型辐射的模型的单元。在这样的示例中,辐射的类型可以被指定为不同的颜色。这样的示例可以包括将从多个断层的它们对应的侧面辐射的不同颜色呈现给显示器。作为示例,一种方法可以包括生成至少一个包括不同类型的辐射的单元。
作为示例,模型可以包括三维单元。作为示例,扇区可以包括三维扇区(参见例如图11)。作为示例,一隔区可以是三维隔区。
作为示例,系统可以包括处理器;可操作地耦合到处理器的存储器;存储在存储器中的处理器可执行指令,以指示系统:访问地下区域的模型,其中该模型包括与该地下区域中的断层相对应的断层;以及经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元,并且经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;并至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该值表示单元中的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。在这样的示例中,可以包括处理器可执行指令以指示系统将图形用户界面呈现给显示器,其中图形用户界面包括与不同类型的辐射相对应的不同颜色。
作为示例,一个或多个计算机可读存储介质可以包括用于指示计算机的计算机可执行指令,该指令包括以下指令:访问地下区域的模型,其中该模型包括与地下区域中的断层相对应的断层;经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元,并且经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;并至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该值表示该单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。
作为示例,一个或多个示例方法可以包括各种计算机可读介质(CRM)块或与各种计算机可读介质(CRM)块相关联。这样的块通常包括适合于由一个或多个处理器(或处理器核)执行的指令,以指示计算装置或系统执行一个或多个动作。虽然示出了各种块,但是单个介质可以配置有指令以至少部分地允许一种或多种方法的各种动作的执行。作为示例,计算机可读介质(CRM)可以是非暂时性的并且不是载波的计算机可读存储介质。
作为示例,工作流程可以与各种计算机可读介质(CRM)块相关联。这样的块通常包括适合于由一个或多个处理器(或核)执行的指令,以指示计算装置或系统执行一个或多个动作。作为示例,单个介质可以配置有指令以至少部分地允许工作流程的各种动作的执行。作为示例,计算机可读介质(CRM)可以是计算机可读存储介质。作为示例,块可以被提供为一组或多组指令,例如,诸如图2的系统250的一组或多组指令270等。
图14示出了计算系统1400的示例和联网系统1410的示例的组件。系统1400包括一个或多个处理器1402、存储器和/或存储组件1404、一个或多个输入和/或输出装置1406和总线1408。在示例实施例中,指令可以存储在一个或多个计算机可读介质(例如,存储器/存储组件1404)中。这样的指令可由一个或多个处理器(例如,一个或多个处理器1402)经由通信总线(例如,总线1408)读取,该通信总线可以是有线的或无线的。一个或多个处理器可以执行这样的指令以(全部或部分地)实现一个或多个属性(例如,作为方法的一部分)。用户可以查看来自过程的输出并经由I/O装置(例如,装置1406)与该过程交互。在示例实施例中,计算机可读介质可以是诸如物理存储器、存储装置之类的存储组件,例如,芯片、封装件上的芯片、存储卡等(例如,计算机可读存储介质)。
在示例实施例中,诸如在网络系统1410中的组件可以被分布。网络系统1410包括组件1422-1、1422-2、1422-3,...1422-N。例如,组件1422-1可以包括处理器1402,而组件1422-3可以包括处理器1402可访问的存储器。此外,组件1402-2可以包括I/O装置,用于显示和与方法交互。网络1420可以是或包括因特网、内联网、蜂窝网络、卫星网络等。
作为示例,装置可以是包括用于信息通信的一个或多个网络接口的移动装置。例如,移动装置可以包括无线网络接口(例如,可以经由IEEE 802.11、ETSI GSM、卫星等进行操作)。作为示例,移动装置可以包括诸如主处理器、存储器、显示器、显示图形电路(例如,可选地包括触摸和手势电路)、SIM插槽、音频/视频电路、运动处理电路(例如,加速度计、陀螺仪)、无线局域网电路、智能卡电路、发射器电路、GPS电路和电池之类的组件。作为示例,移动装置可以被配置为手机、平板电脑等。作为示例,可以使用移动装置来实现(例如,全部或部分地)的方法。作为示例,系统可以包括一个或多个移动装置。
作为示例,系统可以是分布式环境,例如所谓的“云”环境,其中各种装置、组件等为了数据存储、通信,计算等目的而交互。装置或系统可以包括一个或多个组件,用于经由一个或多个因特网(例如,当通信经由一个或多个因特网协议发生时)、蜂窝网络、卫星网络等进行信息通信。作为示例,可以在分布式环境中(例如,全部或部分地作为基于云的服务)实现该方法。
作为示例,信息可以从显示器(例如,考虑触摸屏)输入,输出到显示器或两者。作为示例,信息可以被输出到投影仪、激光装置、打印机等,使得可以观看该信息。作为示例,信息可以立体地或全息地输出。对于打印机,请考虑使用2D或3D打印机。作为示例,3D打印机可以包含一种或多种可以输出以构造3D对象的物质。例如,可以将数据提供给3D打印机以构造地下地层的3D表示。作为示例,可以以3D(例如,层位等)构造层,以3D构造地质体等。作为示例,可以以3D(例如,作为正结构、作为负结构)等构造孔、裂缝等。
尽管以上仅详细描述了几个示例实施例,但是本领域技术人员将容易理解,在示例实施例中许多修改是可能的。因此,所有这些修改旨在被包括在如所附权利要求所限定的本公开的范围内。在权利要求中,装置加功能的条款旨在覆盖这里描述的执行所叙述的功能的结构,不仅覆盖结构上的等同物,而且还覆盖等同的结构。因此,尽管钉子和螺钉可能不是结构上的等同物,因为钉子采用圆柱形表面将木质零件固定在一起,而钉子采用螺旋表面,但是在紧固木制零件的环境中,钉子和螺钉可能是等效结构。申请人的明确意图是不援引35U.S.C.§112,第6段,对此处的任何权利要求进行任何限制,但权利要求中明确使用“用于……的装置”和相关功能的限制除外。
Claims (20)
1.一种方法,包括:
访问地下区域的模型,其中该模型包括与地下区域中的断层相对应的断层;
经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元;以及经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;和
至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该侧面性属性值指示该单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。
2.根据权利要求1所述的方法,包括将所述模型分成扇区,其中,每个扇区包括所述单元的一部分,并且其中分配侧面性属性值的单元的一部分是所述扇区之一。
3.根据权利要求1所述的方法,包括对于对应数量的单元的一定数量的已分配的侧面性属性值,分配至少一个隔区属性值。
4.根据权利要求3所述的方法,其中,至少一个隔区属性值中的每一个包括标识符,该标识符是用于模型内对应的空间定义的隔区的一组顺序标识符的成员。
5.根据权利要求3所述的方法,包括通过另一断层来检测所述断层的至少一个截断,其中,分配至少一个隔区属性值至少部分地基于该检测。
6.根据权利要求5所述的方法,其中,至少部分地基于检测到的另一断层对所述断层的截断,为所述断层的一侧分配至少两个隔区属性值。
7.根据权利要求1所述的方法,包括识别由包括主要侧面性属性值的单元的一部分组成的第一扇区,其中,所述主要侧面性属性值不同于由所述单元的不同部分组成的第二扇区的主要侧面性属性值,其中,第一扇区和第二扇区中的至少一些单元是相邻单元。
8.根据权利要求7所述的方法,包括对第一扇区中的每个单元执行辐射计算,以将侧面性属性值分配给每个单元。
9.根据权利要求1所述的方法,包括执行分层的分配,其中,一个层基于单元扇区来分配单元侧面性属性值,以及其中另一层基于单独单元来分配单元侧面性属性值。
10.根据权利要求1所述的方法,其中,第一类型辐射和第二类型辐射关于距所述断层的距离衰减。
11.根据权利要求1所述的方法,包括将第一颜色的第一类型辐射和第二颜色的第二类型辐射呈现给显示器。
12.根据权利要求1所述的方法,其中,模型包括多个断层,并且对于每个单独的断层,包括经由单独的断层的第一侧,用相应类型辐射从其向外辐射该模型的单元,以及经由单独的断层的相对的第二侧,用相应的不同类型辐射从其向外辐射该模型的单元。
13.根据权利要求12所述的方法,其中,将辐射的类型指定为不同的颜色。
14.根据权利要求13所述的方法,包括将从多个断层的它们对应的侧面辐射的不同颜色呈现给显示器。
15.根据权利要求12所述的方法,其中,至少一个单元包括不同类型的辐射。
16.根据权利要求1所述的方法,其中,所述单元包括三维单元。
17.根据权利要求2所述的方法,其中,所述扇区包括三维扇区。
18.一种系统,包括:
处理器;
存储器,可操作地耦合到处理器;和
处理器可执行指令,存储在存储器中,以指示该系统:
访问地下区域的模型,其中该模型包括与地下区域中的断层相对应的断层;
经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元,并且经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;和
至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该侧面性属性值指示该单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。
19.根据权利要求18所述的系统,包括处理器可执行指令,以指示所述系统将图形用户界面呈现给显示器,其中,所述图形用户界面包括与不同类型的辐射相对应的不同颜色。
20.一种或多种计算机可读存储介质,包括用于指示计算机的计算机可执行指令,所述指令包括进行以下处理的指令:
访问地下区域的模型,其中该模型包括与地下区域中的断层相对应的断层;
经由断层的第一侧,用第一类型辐射从其向外辐射该模型的单元,并且经由断层的相对的第二侧,用第二类型辐射从其向外辐射该模型的单元;和
至少部分地基于第一类型辐射的向外辐射和第二类型辐射的向外辐射,为模型的单元的一部分分配关于断层的侧面性属性值,该侧面性属性值指示该单元的一部分中的每个单元被分配给断层的第一侧或分配给断层的第二侧。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201862661070P | 2018-04-22 | 2018-04-22 | |
US62/661,070 | 2018-04-22 | ||
PCT/US2019/028028 WO2019209609A1 (en) | 2018-04-22 | 2019-04-18 | Fault radiation based grid compartmentalization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112292714A true CN112292714A (zh) | 2021-01-29 |
CN112292714B CN112292714B (zh) | 2024-04-16 |
Family
ID=68295097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201980040259.XA Active CN112292714B (zh) | 2018-04-22 | 2019-04-18 | 基于断层辐射的网格隔分 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11947071B2 (zh) |
EP (1) | EP3785233B1 (zh) |
CN (1) | CN112292714B (zh) |
CA (1) | CA3097743A1 (zh) |
WO (1) | WO2019209609A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2620556A (en) * | 2022-07-05 | 2024-01-17 | Equinor Energy As | Internal imaging and model construction |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2631685A2 (en) * | 2012-02-22 | 2013-08-28 | Service Pétroliers Schlumberger | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces |
CN105659292A (zh) * | 2013-11-26 | 2016-06-08 | 界标制图有限公司 | 复杂断层网络拓扑的全局网格建立无断层序列 |
US9728003B1 (en) * | 2009-10-23 | 2017-08-08 | Paradigm Sciences Ltd. | Systems and methods for coordinated editing of seismic data in dual model |
WO2018069744A1 (en) * | 2016-10-14 | 2018-04-19 | Schlumberger Technology Corporation | Geologic structural model generation |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2603877B1 (en) | 2010-08-09 | 2020-07-01 | Landmark Graphics Corporation | Systems and methods for creating a surface in a faulted space |
US9366776B2 (en) | 2010-11-30 | 2016-06-14 | Schlumberger Technology Corporation | Integrated formation modeling systems and methods |
US10534871B2 (en) * | 2011-03-09 | 2020-01-14 | Schlumberger Technology Corporation | Method and systems for reservoir modeling, evaluation and simulation |
EP3052756B1 (en) | 2013-10-04 | 2021-05-26 | Services Pétroliers Schlumberger | Inversion-based workflow for consistent interpretation of nuclear density images in horizontal wells |
WO2018063195A1 (en) | 2016-09-28 | 2018-04-05 | Halliburton Energy Services, Inc. | Electromagnetic reservoir monitoring systems and methods including earth |
-
2019
- 2019-04-18 EP EP19792664.5A patent/EP3785233B1/en active Active
- 2019-04-18 WO PCT/US2019/028028 patent/WO2019209609A1/en unknown
- 2019-04-18 CN CN201980040259.XA patent/CN112292714B/zh active Active
- 2019-04-18 CA CA3097743A patent/CA3097743A1/en active Pending
- 2019-04-18 US US15/733,775 patent/US11947071B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9728003B1 (en) * | 2009-10-23 | 2017-08-08 | Paradigm Sciences Ltd. | Systems and methods for coordinated editing of seismic data in dual model |
US10013800B1 (en) * | 2009-10-23 | 2018-07-03 | Emerson Paradigm Holding Llc | Systems and methods for coordinated editing of seismic data in dual model |
EP2631685A2 (en) * | 2012-02-22 | 2013-08-28 | Service Pétroliers Schlumberger | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces |
CN105659292A (zh) * | 2013-11-26 | 2016-06-08 | 界标制图有限公司 | 复杂断层网络拓扑的全局网格建立无断层序列 |
WO2018069744A1 (en) * | 2016-10-14 | 2018-04-19 | Schlumberger Technology Corporation | Geologic structural model generation |
Non-Patent Citations (1)
Title |
---|
DONGFANG QU.ET AL: "A method for generating volumetric fault zone grids for pillar gridded reservoir models", COMPUTERS & GEOSCIENCES, pages 28 - 37 * |
Also Published As
Publication number | Publication date |
---|---|
EP3785233A4 (en) | 2022-01-26 |
EP3785233A1 (en) | 2021-03-03 |
CN112292714B (zh) | 2024-04-16 |
US11947071B2 (en) | 2024-04-02 |
EP3785233B1 (en) | 2023-12-20 |
WO2019209609A1 (en) | 2019-10-31 |
CA3097743A1 (en) | 2019-10-31 |
US20210247543A1 (en) | 2021-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
USRE49507E1 (en) | Faulted geological structures having unconformities | |
US20150113379A1 (en) | Representation of unstructured grids | |
US11175434B2 (en) | Geologic stratigraphy via implicit and jump functions | |
US10571601B2 (en) | Geologic fault representation for flow simulation | |
US20130218539A1 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
US20150066460A1 (en) | Stratigraphic function | |
US11042676B2 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
US11249208B2 (en) | Geologic structural model generation | |
EP2631685A2 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
US20220163692A1 (en) | Modeling and simulating faults in subterranean formations | |
US20160202389A1 (en) | H-matrix preconditioner | |
NO20201045A1 (en) | Fluid simulator property representation | |
CN112292714B (zh) | 基于断层辐射的网格隔分 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |