CN105659294A - 用于离散网络网格化的3d视距算法 - Google Patents
用于离散网络网格化的3d视距算法 Download PDFInfo
- Publication number
- CN105659294A CN105659294A CN201380077145.5A CN201380077145A CN105659294A CN 105659294 A CN105659294 A CN 105659294A CN 201380077145 A CN201380077145 A CN 201380077145A CN 105659294 A CN105659294 A CN 105659294A
- Authority
- CN
- China
- Prior art keywords
- facet
- sighting distance
- sighting
- computer
- polygon
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Remote Sensing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Nonwoven Fabrics (AREA)
- Processing Or Creating Images (AREA)
Abstract
所公开的实施方案包括用于对三维物体,例如但不限于裂缝,进行建模以进行储层模拟的方法、设备和计算机程序产品。例如,一个所公开的实施方案包括一种方法,所述方法接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;围绕相交的刻面产生闭合的3D视距表面;使用各种形状的单元来离散化所有3D视距表面;在所述闭合的3D视距表面中的每一者内产生3D单元;以及给每一单元指派适合于数值模拟程序(例如但不限于储层模拟软件)的储层特性。
Description
技术领域
本发明一般来说涉及用于产生可用于建构地下储层的模拟模型的网格的系统和方法,更明确地说,涉及被配置成用于对地质裂缝进行建模的系统和方法。
发明背景
在油气产业中,储层建模涉及建构石油储层的计算机模型以便提高储量的估计并作出关于油田开发的决策。例如,可以建立地质模型以在开采之前提供对储层的静态描述。相反地,可以建立储层模拟模型以模拟在储层的开采寿命内流体在储层内的流动。
储层模拟模型的一个挑战是储层内裂缝的建模,所述建模要求透彻地理解矩阵流特性、裂缝网络连接性和裂缝矩阵相互作用。裂缝可以被描述为地层内的张开的裂隙或空隙,并且可以天然地形成或由井筒人工地产生。裂缝的恰当建模是重要的,因为裂缝的特性如空间分布、孔隙、长度、高度、传导性和连接性显著地影响储层流体至井筒的流动。
因此,所公开的实施方案提供了用于围绕复杂且离散的裂缝产生3D混合计算网格以便进行储层模拟的系统、方法和计算机程序产品。
附图简述
下文参看附图详细地描述了本发明的说明性实施方案,所述附图以引用方式并入本文中并且其中:
图1示出了根据所公开的实施方案建模的三维裂缝的图像;
图2是示出了根据所公开的实施方案的用于对三维裂缝进行建模的方法的流程图;
图3示出了根据所公开的实施方案的具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝的实例;
图4提供了根据所公开的实施方案的用于围绕多边形刻面以指定半径产生一组3D视距的过程的图示;
图5示出了根据所公开的实施方案围绕多边形刻面产生的一组3D视距的内部视图的实例;
图5A示出了根据所公开的实施方案围绕另一多边形刻面产生的一组3D视距的内部视图的另一实例;
图6示出了根据所公开的实施方案围绕两个相邻的多边形刻面产生的一组相交的3D视距的内部视图的实例;
图6A示出了根据所公开的实施方案的图6所示的该组相交的3D视距的3D外部视图的实例;
图7示出了根据所公开的实施方案的围绕两个相邻的多边形刻面产生的一组相交的3D视距的内部视图的另一实例;
图7A示出了根据所公开的实施方案的图7所示的该组相交的3D视距的3D外部视图的实例;以及
图8是示出了用于实现所公开的实施方案的系统的一个实施方案的框图。
具体实施方式
所公开的实施方案包括用于对三维(3D)物体(例如但不限于地质裂缝)进行建模的系统和方法。参看图式中的图1至图8可以最好地理解所公开的实施方案以及其优点,相同的元件符号用于各个图中的相同且对应的零部件。在阅读了以下图式和详细描述之后,所公开的实施方案的其它特征和优点将是或将变成是本领域的普通技术人员显而易见的。希望所有此类另外的特征和优点包括在所公开的实施方案的范围内。另外,所示出的图式仅为示例性的并且不希望确定或暗示关于不同的实施方案可以在其中实现的环境、架构、设计或过程的任何限制。
图1示出了根据所公开的实施方案进行建模的三维裂缝的图像。如图像100中可见,地球地层的层包括地层内的裂缝。如上所述,这些裂缝可以被描述为地层内的张开的裂隙或空隙,并且可以天然地形成或由井筒人工地产生。理解这些裂缝的适当特性并对这些裂缝的适当特性进行建模是重要的,因为裂缝实现并影响储层流体至井筒的流动。可以使用图像测井仪来获得或产生如图像100等图像。图像测井仪使用旋转传感器来测量整个钻孔壁上的声阻抗以识别岩石裂缝的存在和方向以及理解地层的倾斜方向。
图2是示出了根据所公开的实施方案的用于对三维裂缝进行建模的方法/过程200的流程图。在所描绘的实施方案中,所述方法从接收3D域开始,所述3D域包括具有已通过多边形刻面的集合在2D流形中离散化的几何形状的一组3D裂缝表面(步骤202)。在替代实施方案中,过程200可以从将该组3D裂缝表面离散化成具有多边形刻面的集合的2D流形以产生3D域开始。
在步骤204处,所述方法围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距。在一个实施方案中,每一3D视距由平行于所述刻面的两个平坦的面、沿着每一条边的半圆柱体以及每一顶点处的球形截面组成以完全围封原始刻面,其中从每一表面至所述刻面的距离是恒定半径。在某些实施方案中,3D视距的数目和/或围绕多边形刻面的集合中的每一刻面的指定半径可以是用户可修改的参数。
在步骤206处,所述方法围绕相交的刻面产生闭合的3D视距表面。在一个实施方案中,步骤206包括计算针对每一指定半径的所有3D视距表面的交点(步骤206A)、识别完全被其它刻面的3D视距含有的3D视距表面区段(步骤206B)以及丢弃所述所含有的表面区段以产生围绕所述刻面的闭合表面(步骤206C)。
所述过程接着使用各种形状的单元(例如但不限于三角形、四边形和一般多边形)对所有3D视距表面进行离散化(步骤208)。接着用3D单元填充2D元素的所得多组流形(步骤210)。由此,所述过程可以将储层特性(例如但不限于孔隙率和渗透性)指派给所述单元中的每一者以便对储层的流体流进行建模(步骤212)。这些特性值可以由用户手动输入或者可以从测井或从含有相关地质信息的数据库中自动地提取。
最后,在一个实施方案中,所述过程将三维单元模型输入至模拟程序(例如但不限于储层模拟软件)中,用于进行数值模拟和用于评估流体流(步骤214),其后过程200结束。
图3示出了根据所公开的实施方案的具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝300的实例。在所描绘的实施方案中,该组3D裂缝300包括三个相交的裂缝302、304和306。明确地说,裂缝302被裂缝304相交并且还被裂缝306相交。裂缝302、304和306中的每一者是通过多边形刻面的集合通过2D表面离散化。虽然示出了三角形刻面,但是所公开的实施方案可以使用任何多边形形状的刻面。
图4提供了根据所公开的实施方案的用于围绕多边形刻面402以指定半径产生一组3D视距的过程的图示。在所描绘的实施方案中,产生以指定半径平行于刻面402的两个平坦的面404。接下来,在每一顶点处产生沿着每一条边的半圆柱体406和球形截面408以完全围封原始刻面402。从每一表面至刻面402的距离是恒定的半径。
接着以类似方式围绕第一3D视距以距离第一3D视距恒定的半径建构第二3D视距以完全围封第一3D视距和刻面402。在某些实施方案中,第二3D视距与第一3D视距之间的半径与第一3D视距与刻面402之间的半径相同。或者,在一些实施方案中,第二3D视距与第一3D视距之间的半径可以不同于第一3D视距与刻面402之间的半径。接着以相同的方式产生另外的3D视距以围封第二3D视距,等等。如先前所述,在某些实施方案中,3D视距的数目和/或围绕多边形刻面的集合中的每一刻面的指定半径可以是用户可修改的参数。
图5示出了根据所公开的实施方案的围绕三角形形状的多边形刻面产生的一组3D视距的内部俯视图的实例。图5A示出了根据所公开的实施方案的围绕梯形形状的多边形刻面产生的一组3D视距的内部俯视图的另一实例。如从所示实施方案中可见,多边形形状的刻面中的每一者被一组类似形状的3D视距环绕。
由于所述过程围绕多边形刻面产生所述多组3D视距,因此用于邻近/相邻多边形刻面的多组3D视距将彼此相交。例如,图6示出了根据所公开的实施方案的围绕彼此相交的两个相邻多边形刻面产生的两组3D视距的内部视图的实例。如过程200中所描述,在一个实施方案中,所公开的方法通过计算针对每一指定半径的所有3D视距表面的交点以识别被其它刻面的3D视距完全含有的3D视距表面区段并且将所含有的表面区段丢弃以产生围绕所述刻面的闭合表面来围绕这些相交刻面产生闭合3D视距表面。例如,图6A示出了图6所示的所述多组相交3D视距的所产生的3D外部视图的实例,其中相交的被完全含有的表面区段已被丢弃。
图7示出了根据所公开的实施方案的围绕彼此相交的两个相邻多边形刻面产生的两组3D视距的内部视图的另一实例。图7A示出了根据所公开的实施方案的围绕图7所示的两组相交的3D视距的刻面的闭合表面的3D外部视图的实例。
图8是示出了用于实现所公开的实施方案的特征和功能的系统800的一个实施方案的框图。除了其它组件外,系统800还包括处理器800、主存储器802、辅助存储单元804、输入/输出接口模块806和通信接口模块808。处理器800可以是能够执行用于执行所公开的实施方案的特征和功能的指令的任何类型或任何数目的单核或多核处理器。
输入/输出接口模块806使得系统800能够接收用户输入(例如,从键盘和鼠标)以及将信息输出至一个或多个装置,例如但不限于打印机、外部数据存储装置和音频扬声器。系统800可以任选地包括使得信息能够在集成或外部显示装置上显示的单独的显示模块810。例如,显示模块810可以包括用于提供与一个或多个显示装置相关联的增强的图形、触摸屏和/或多触点功能性的指令或硬件(例如,图形卡或芯片)。
主存储器802是存储目前执行的指令/数据或为了执行而预先取出的指令/数据的易失性存储器。辅助存储单元804是用于存储持久性数据的非易失性存储器。辅助存储单元804可以是或包括任何类型的数据存储部件,例如硬盘、闪存盘或存储卡。在一个实施方案中,辅助存储单元804存储用于使用户能够执行所公开的实施方案的特征和功能的计算机可执行代码/指令和其它相关数据。
例如,根据所公开的实施方案,辅助存储单元804可以永久地存储用于对三维(3D)物体(例如但不限于地质裂缝)进行建模的上述3D视距算法820的可执行代码/指令。接着在由如图8中所示的处理器800执行期间,将与3D视距算法820相关联的指令从辅助存储单元804载入至主存储器802。
通信接口模块808使系统800能够与通信网络830通信。例如,网络接口模块808可以包括用于使系统800能够通过通信网络830和/或直接地与其它装置发送和接收数据的网络接口卡和/或无线收发器。
通信网络830可以是任何类型的网络,包括以下网络中的一者或多者的组合:广域网、局域网、一个或多个专用网络、因特网、电话网(例如公共交换电话网(PSTN))、一个或多个蜂窝式网络、和无线数据网络。通信网络830可以包括多个网络节点(未描绘),例如路由器、网络接入点/网关、交换机、DNS服务器、代理服务器以及用于帮助在装置之间路由数据/通信的其它网络节点。
例如,在一个实施方案中,系统800可以与一个或多个服务器834或数据库832交互以执行本发明的特征。例如,系统800可以查询数据库832以获得地质信息,以便将储层特性指派给单元以执行模拟。系统800可以查询数据库832以获得测井信息,以便确定裂缝取向或密度以使得能够根据所公开的实施方案对裂缝进行建模。另外,在某些实施方案中,系统800可以充当用于一个或多个客户端装置的服务器系统或用于与一个或多个装置的对等通信或并行处理的对等系统。
因此,如上所述,所公开的实施方案的优点包括但不限于使得围绕复杂几何形状快速地产生具有结构化元素的未结构化网格。另外,不需要很多的专业知识就使用户能够利用所公开的实施方案来产生适合于许多数值模拟器的高质量的网格单元。例如,所公开的实施方案实现了使非专业人员将先进的数值建模技术用于复杂几何形状的工作流程,这在之前需要用户进行粗略近似和/或需要来自数值建模专家的使用前帮助。
虽然已描述了关于上述实施方案的具体细节,但是上述硬件和软件描述希望仅作为示例实施方案并且不希望限制所公开的实施方案的结构或实现方式。例如,虽然未示出系统800的许多其它内部组件,但是本领域的普通技术人员将了解此类组件以及其互连是熟知的。
另外,所公开的实施方案的某些方面(如上文所概述)可以体现为使用一个或多个处理单元/组件执行的软件。技术的程序方面可以被认为是通常呈载运在一类机器可读介质上或包括在一类机器可读介质中的可执行代码和/或相关联数据的形式的“产品”或“制品”。有形非暂时性“存储”类型介质包括用于计算机、处理器或类似者的任何或所有存储器或其它存储装置或可以在任何时候为软件编程提供存储的其相关联模块,例如各种半导体存储器、磁带驱动器、磁盘驱动器、光盘或磁盘和类似者。
另外,图式中的流程图和框图示出了根据本发明的各种实施方案的系统、方法和计算机程序产品的可能实现方式的架构、功能性和操作。还应指出,在一些替代实现方式中,框中所指出的功能可以不按图式中所指出的次序来发生。例如,接连示出的两个框实际上可以基本上同时执行,或者所述框有时可以按相反次序来执行,这取决于所涉及的功能性。还将指出,框图和/或流程图图示的每一框以及框图和/或流程图图示中的框的组合可以通过专用的基于硬件的系统来实现,所述系统执行所指定的功能或动作或者专用硬件和计算机指令的组合。
所公开的实施方案包括用于围绕复杂且离散的断裂产生3D混合计算网格以便进行储层模拟的方法、设备和计算机程序产品。例如,一个所公开的实施方案是用于对三维(3D)地质裂缝进行建模的计算机实现的方法。所述方法包括接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝的步骤。所述方法围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距。所述方法接着围绕相交的刻面产生闭合的3D视距表面。所述方法使用各种形状的单元来离散化所有3D视距表面以及在所述闭合3D视距表面中的每一者内产生3D单元。所述方法还可以给每一单元指派适合于数值模拟程序(例如但不限于储层模拟软件)的储层特性。
在一些实施方案中,围绕相交的刻面产生闭合的3D视距表面包括计算所有3D视距表面的交点、识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段以及丢弃所含有的3D视距表面区段而产生围绕相交刻面的闭合3D视距表面。各种形状的单元可以包括三角形、四边形和一般多边形。每一3D视距可以包括平行于刻面的两个平坦的面、沿着每一条边的半圆柱体和在每一顶点处的球形截面以完全围封所述刻面,并且从每一表面至所述刻面的距离可以是恒定半径。在另一实施方案中,所述计算机实现的方法还可以包括给每一单元指派适合于数值模拟程序的储层特性,并且所述数值模拟程序可以是储层模拟软件。
在又一实施方案中,非暂时性计算机可读介质包括用于对三维(“3D”)结构进行建模的计算机可执行指令。所述计算机可执行指令在被执行时使一个或多个机器执行多个操作,所述操作包括接收具有通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝。围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距。围绕相交的刻面产生闭合的3D视距表面。使用各种形状的单元对所有3D视距表面进行离散化,并且在闭合的3D视距表面中的每一者内产生3D单元。
在一些实施方案中,围绕相交的刻面产生闭合的3D视距表面可以包括计算所有3D视距表面的交点、识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段以及丢弃所含有的3D视距表面区段而产生围绕相交刻面的闭合的3D视距表面。各种形状的单元可以包括三角形、四边形和一般多边形。每一3D视距可以包括平行于刻面的两个平坦的面、沿着每一条边的半圆柱体和在每一顶点处的球形截面以完全围封所述刻面,并且从每一表面至所述刻面的距离可以是恒定半径。在另一实施方案中,所述计算机可读介质还可以包括用于给每一单元指派适合于数值模拟程序的储层特性的计算机可执行指令,并且所述数值模拟程序可以是储层模拟软件。
在又一实施方案中,一种系统包括至少一个处理器和耦接至所述至少一个处理器并存储指令的至少一个存储器。当执行所述指令时,所述至少一个处理器执行操作,所述操作包括:接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;围绕相交的刻面产生闭合的3D视距表面;使用各种形状的单元来离散化所有3D视距表面;以及在所述闭合的3D视距表面中的每一者内产生3D单元。
在另一实施方案中,围绕相交的刻面产生闭合的3D视距表面可以包括计算所有3D视距表面的交点、识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段以及丢弃所含有的3D视距表面区段而产生围绕相交刻面的闭合的3D视距表面。各种形状的单元可以包括三角形、四边形和一般多边形。每一3D视距可以包括平行于刻面的两个平坦的面、沿着每一条边的半圆柱体和在每一顶点处的球形截面以完全围封所述刻面,并且从每一表面至所述刻面的距离可以是恒定半径。在又一实施方案中,所述系统还可以包括用于给每一单元指派适合于数值模拟程序的储层特性的计算机可执行指令,并且所述数值模拟程序可以是储层模拟软件。
所公开的实施方案的一个优点是它使得能够围绕复杂几何形状快速地产生具有结构化元素的未结构化网格。
本文中使用的术语是仅用于描述特定实施方案而不希望限制本发明。如本文中使用,单数形式“一”、“一个”和“所述”希望也包括复数形式,除非上下文另外清楚指明。将进一步理解,术语“包括”(comprise/comprising)在本说明书和/或权利要求书中使用时指明了所述特征、整体、步骤、操作、元件和/或组件的存在,但不排除一个或多个其它特征、整体、步骤、操作、元件、组件和/或其群组的存在或增添。以下的权利要求书中的所有装置或步骤加功能元件的对应结构、材料、动作和等效物希望包括用于结合如明确地要求权利保护的其它要求权利保护的元件来执行功能的任何结构、材料或动作。已为了进行说明和描述而呈现本发明的描述,但本发明的描述不意欲为详尽的或将本发明限于所公开的形式。在不脱离本发明的范围和精神的情况下,许多修改和变化将是本领域的普通技术人员显而易见的。实施方案经过选择和描述以阐释本发明的原理和实际应用,并且使本领域的普通技术人员能够理解本发明以获得具有适合于所预期的特定使用的各种修改的各种实施方案。权利要求书的范围意欲广泛地涵盖所公开的实施方案和任何此类修改。
权利要求书(按照条约第19条的修改)
1.一种用于对三维(3D)地质裂缝进行建模的计算机实现的方法,所述方法包括:
接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;
围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;
围绕相交的刻面产生闭合的3D视距表面;
使用各种形状的单元来离散化所有3D视距表面;以及
在所述闭合的3D视距表面中的每一者内产生3D单元。
2.如权利要求1所述的计算机实现的方法,其中围绕相交的刻面产生闭合的3D视距表面包括:
计算所有3D视距表面的交点;
识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段;以及
丢弃所述所含有的3D视距表面区段以产生围绕相交刻面的所述闭合的3D视距表面。
3.如权利要求1所述的计算机实现的方法,其中所述各种形状的单元包括三角形、四边形和一般多边形。
4.如权利要求1所述的计算机实现的方法,其中每一3D视距由平行于所述刻面的两个平坦的面、沿着每一条边的半圆柱体以及每一顶点处的球形截面组成以完全围封所述刻面,并且其中从每一表面至所述刻面的距离是恒定半径。
5.如权利要求1所述的计算机实现的方法,还包括给每一单元指派适合于数值模拟程序的储层特性。
6.如权利要求5所述的计算机实现的方法,其中所述数值模拟程序是用于对储层地层中的流体流进行建模的储层模拟程序。
7.一种包括用于对三维(“3D”)结构进行建模的计算机可执行指令的非暂时性计算机可读介质,所述计算机可执行指令在被执行时使一个或多个机器执行操作,所述操作包括:
接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;
围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;
围绕相交的刻面产生闭合的3D视距表面;
使用各种形状的单元来离散化所有3D视距表面;以及
在所述闭合的3D视距表面中的每一者内产生3D单元。
8.如权利要求7所述的计算机可读介质,其中围绕相交的刻面产生闭合的3D视距表面包括:
计算所有3D视距表面的交点;
识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段;以及
丢弃所述所含有的3D视距表面区段而产生围绕相交刻面的所述闭合的3D视距表面。
9.如权利要求7所述的计算机可读介质,其中所述各种形状的单元包括三角形、四边形和一般多边形。
10.如权利要求7所述的计算机可读介质,其中每一3D视距由平行于所述刻面的两个平坦的面、沿着每一条边的半圆柱体以及每一顶点处的球形截面组成以完全围封所述刻面,并且其中从每一表面至所述刻面的距离是恒定半径。
11.如权利要求7所述的计算机可读介质,还包括用于给每一单元指派适合于数值模拟程序的储层特性的计算机可执行指令。
12.如权利要求11所述的计算机可读介质,其中所述数值模拟程序是用于对储层地层中的流体流进行建模的储层模拟程序。
13.一种系统,包括:
至少一个处理器;以及
至少一个存储器,所述至少一个存储器耦接至所述至少一个处理器并存储在由所述至少一个处理器执行时执行操作的指令,所述操作包括:
接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;
围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;
围绕相交的刻面产生闭合的3D视距表面;
使用各种形状的单元来离散化所有3D视距表面;以及
在所述闭合的3D视距表面中的每一者内产生3D单元。
14.如权利要求13所述的系统,其中围绕相交的刻面产生闭合的3D视距表面包括:
计算所有3D视距表面的交点;
识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段;以及
丢弃所述所含有的3D视距表面区段而产生围绕相交刻面的所述闭合的3D视距表面。
15.如权利要求13所述的系统,其中所述各种形状的单元包括三角形、四边形和一般多边形。
16.如权利要求13所述的系统,其中每一3D视距由平行于所述刻面的两个平坦的面、沿着每一条边的半圆柱体以及每一顶点处的球形截面组成以完全围封所述刻面,并且其中从每一表面至所述刻面的距离是恒定半径。
17.如权利要求13所述的系统,还包括用于给每一单元指派适合于数值模拟程序的储层特性的计算机可执行指令。
18.如权利要求17所述的系统,其中所述数值模拟程序是用于对储层地层中的流体流进行建模的储层模拟程序。
Claims (18)
1.一种用于对三维(3D)地质裂缝进行建模的计算机实现的方法,所述方法包括:
接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;
围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;
围绕相交的刻面产生闭合的3D视距表面;
使用各种形状的单元来离散化所有3D视距表面;以及
在所述闭合的3D视距表面中的每一者内产生3D单元。
2.如权利要求1所述的计算机实现的方法,其中围绕相交的刻面产生闭合的3D视距表面包括:
计算所有3D视距表面的交点;
识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段;以及
丢弃所述所含有的3D视距表面区段以产生围绕相交刻面的所述闭合的3D视距表面。
3.如权利要求1所述的计算机实现的方法,其中所述各种形状的单元包括三角形、四边形和一般多边形。
4.如权利要求1所述的计算机实现的方法,其中每一3D视距由平行于所述刻面的两个平坦的面、沿着每一条边的半圆柱体以及每一顶点处的球形截面组成以完全围封所述刻面,并且其中从每一表面至所述刻面的距离是恒定半径。
5.如权利要求1所述的计算机实现的方法,还包括给每一单元指派适合于数值模拟程序的储层特性。
6.如权利要求5所述的计算机实现的方法,其中所述数值模拟程序是储层模拟软件。
7.一种包括用于对三维(“3D”)结构进行建模的计算机可执行指令的非暂时性计算机可读介质,所述计算机可执行指令在被执行时使一个或多个机器执行操作,所述操作包括:
接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;
围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;
围绕相交的刻面产生闭合的3D视距表面;
使用各种形状的单元来离散化所有3D视距表面;以及
在所述闭合的3D视距表面中的每一者内产生3D单元。
8.如权利要求7所述的计算机可读介质,其中围绕相交的刻面产生闭合的3D视距表面包括:
计算所有3D视距表面的交点;
识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段;以及
丢弃所述所含有的3D视距表面区段而产生围绕相交刻面的所述闭合的3D视距表面。
9.如权利要求7所述的计算机可读介质,其中所述各种形状的单元包括三角形、四边形和一般多边形。
10.如权利要求7所述的计算机可读介质,其中每一3D视距由平行于所述刻面的两个平坦的面、沿着每一条边的半圆柱体以及每一顶点处的球形截面组成以完全围封所述刻面,并且其中从每一表面至所述刻面的距离是恒定半径。
11.如权利要求7所述的计算机可读介质,还包括用于给每一单元指派适合于数值模拟程序的储层特性的计算机可执行指令。
12.如权利要求11所述的计算机可读介质,其中所述数值模拟程序是储层模拟软件。
13.一种系统,包括:
至少一个处理器;以及
至少一个存储器,所述至少一个存储器耦接至所述至少一个处理器并存储在由所述至少一个处理器执行时执行操作的指令,所述操作包括:
接收具有已通过多边形刻面的集合通过2D表面离散化的几何形状的一组3D裂缝;
围绕所述多边形刻面的集合中的每一刻面以指定半径产生一组3D视距;
围绕相交的刻面产生闭合的3D视距表面;
使用各种形状的单元来离散化所有3D视距表面;以及
在所述闭合的3D视距表面中的每一者内产生3D单元。
14.如权利要求13所述的系统,其中围绕相交的刻面产生闭合的3D视距表面包括:
计算所有3D视距表面的交点;
识别被其它刻面的3D视距完全含有的所含有的3D视距表面区段;以及
丢弃所述所含有的3D视距表面区段而产生围绕相交刻面的所述闭合的3D视距表面。
15.如权利要求13所述的系统,其中所述各种形状的单元包括三角形、四边形和一般多边形。
16.如权利要求13所述的系统,其中每一3D视距由平行于所述刻面的两个平坦的面、沿着每一条边的半圆柱体以及每一顶点处的球形截面组成以完全围封所述刻面,并且其中从每一表面至所述刻面的距离是恒定半径。
17.如权利要求13所述的系统,还包括用于给每一单元指派适合于数值模拟程序的储层特性的计算机可执行指令。
18.如权利要求17所述的系统,其中所述数值模拟程序是储层模拟软件。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2013/049145 WO2015002642A1 (en) | 2013-07-02 | 2013-07-02 | 3d stadia algorithm for discrete network meshing |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105659294A true CN105659294A (zh) | 2016-06-08 |
Family
ID=52144087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380077145.5A Pending CN105659294A (zh) | 2013-07-02 | 2013-07-02 | 用于离散网络网格化的3d视距算法 |
Country Status (9)
Country | Link |
---|---|
US (1) | US9715762B2 (zh) |
CN (1) | CN105659294A (zh) |
AU (1) | AU2013393345B2 (zh) |
CA (1) | CA2913898C (zh) |
DE (1) | DE112013007213T5 (zh) |
GB (1) | GB2529958B (zh) |
MX (1) | MX2015016558A (zh) |
SG (1) | SG11201510115UA (zh) |
WO (1) | WO2015002642A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017082870A1 (en) * | 2015-11-10 | 2017-05-18 | Landmark Graphics Corporation | Fracture network triangle mesh adjustment |
US11077587B2 (en) * | 2017-10-13 | 2021-08-03 | Cornell University | Stretchable surfaces with programmable texture |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080091396A1 (en) * | 2006-10-13 | 2008-04-17 | Kennon Stephen R | Method and system for modeling and predicting hydraulic fracture performance in hydrocarbon reservoirs |
CN102612682A (zh) * | 2009-11-12 | 2012-07-25 | 埃克森美孚上游研究公司 | 用于储层建模和模拟的方法和设备 |
US20120267104A1 (en) * | 2011-04-19 | 2012-10-25 | Halliburton Energy Services, Inc. | System and Method for Improved Propped Fracture Geometry for High Permeability Reservoirs |
US20130124169A1 (en) * | 2010-06-18 | 2013-05-16 | Industry-University Cooperation Foundation Hanyang University | Method for oil prediction in fractured reservoirs and recorded medium therefor |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6549879B1 (en) | 1999-09-21 | 2003-04-15 | Mobil Oil Corporation | Determining optimal well locations from a 3D reservoir model |
US6980940B1 (en) * | 2000-02-22 | 2005-12-27 | Schlumberger Technology Corp. | Intergrated reservoir optimization |
WO2005033739A2 (en) | 2003-09-30 | 2005-04-14 | Exxonmobil Upstream Research Company Corp-Urc-Sw348 | Characterizing connectivity in reservoir models using paths of least resistance |
RU2321064C2 (ru) | 2004-06-03 | 2008-03-27 | Мурманский государственный технический университет | Способ построения обратимой трехмерной гидродинамической модели земли, калибруемой в реальном времени в процессе бурения |
EA201170472A1 (ru) * | 2008-09-19 | 2011-10-31 | Шеврон Ю.Эс.Эй. Инк. | Реализованные посредством компьютера системы и способы для использования при моделировании геомеханической системы коллектора |
RU2563862C2 (ru) | 2010-10-06 | 2015-09-20 | Лоджинд Б.В. | Атрибут напряжения в горных породах |
-
2013
- 2013-07-02 WO PCT/US2013/049145 patent/WO2015002642A1/en active Application Filing
- 2013-07-02 MX MX2015016558A patent/MX2015016558A/es unknown
- 2013-07-02 AU AU2013393345A patent/AU2013393345B2/en not_active Ceased
- 2013-07-02 CA CA2913898A patent/CA2913898C/en not_active Expired - Fee Related
- 2013-07-02 CN CN201380077145.5A patent/CN105659294A/zh active Pending
- 2013-07-02 DE DE112013007213.0T patent/DE112013007213T5/de not_active Withdrawn
- 2013-07-02 US US14/889,110 patent/US9715762B2/en not_active Expired - Fee Related
- 2013-07-02 GB GB1520860.6A patent/GB2529958B/en not_active Expired - Fee Related
- 2013-07-02 SG SG11201510115UA patent/SG11201510115UA/en unknown
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080091396A1 (en) * | 2006-10-13 | 2008-04-17 | Kennon Stephen R | Method and system for modeling and predicting hydraulic fracture performance in hydrocarbon reservoirs |
CN102612682A (zh) * | 2009-11-12 | 2012-07-25 | 埃克森美孚上游研究公司 | 用于储层建模和模拟的方法和设备 |
US20130124169A1 (en) * | 2010-06-18 | 2013-05-16 | Industry-University Cooperation Foundation Hanyang University | Method for oil prediction in fractured reservoirs and recorded medium therefor |
US20120267104A1 (en) * | 2011-04-19 | 2012-10-25 | Halliburton Energy Services, Inc. | System and Method for Improved Propped Fracture Geometry for High Permeability Reservoirs |
Also Published As
Publication number | Publication date |
---|---|
US9715762B2 (en) | 2017-07-25 |
GB2529958B (en) | 2019-11-13 |
SG11201510115UA (en) | 2016-01-28 |
AU2013393345B2 (en) | 2016-02-11 |
GB201520860D0 (en) | 2016-01-13 |
DE112013007213T5 (de) | 2016-04-28 |
US20160104320A1 (en) | 2016-04-14 |
GB2529958A (en) | 2016-03-09 |
MX2015016558A (es) | 2016-05-16 |
WO2015002642A1 (en) | 2015-01-08 |
CA2913898A1 (en) | 2015-01-08 |
AU2013393345A1 (en) | 2015-12-10 |
CA2913898C (en) | 2017-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9922142B2 (en) | Systems and methods for subsurface reservoir simulation | |
US10254440B2 (en) | Process for constructing a volume mesh for modeling geological structures | |
US8463586B2 (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids | |
CA2941020C (en) | Modeling intersecting faults and complex wellbores in reservoir simulation | |
US20140136171A1 (en) | Unstructured Grids For Modeling Reservoirs | |
US20130218539A1 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
US10529131B2 (en) | Simulating fractured reservoirs using multiple meshes | |
EP3018502A2 (en) | Modeling fluid-conducting fractures in reservoir simulation grids | |
CN103097657B (zh) | 用于模拟由储集层的网格界定的多个油气井的生产特性的机器及计算机实施方法 | |
CN104533370A (zh) | 压裂水平井油藏、裂缝、井筒全耦合模拟方法 | |
CN102612682A (zh) | 用于储层建模和模拟的方法和设备 | |
CN105934689A (zh) | 地质特征拆分 | |
EP2631685A2 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
CN102609982A (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
US11500123B2 (en) | Hybrid 3D geocellular representation of selected natural fracture network subsets | |
AU2011271175A1 (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids | |
US20160101474A1 (en) | Spiral toolpaths for high-speed machining of polygonal pockets | |
CN105659294A (zh) | 用于离散网络网格化的3d视距算法 | |
CN105378800A (zh) | 2.5d围场网格化 | |
CN105068141A (zh) | 确定断缝系统的岩溶储层分布的方法和装置 | |
Yanbing et al. | On 3D GIS spatial modeling | |
Min et al. | Selection of geologic models based on pareto-optimality using surface deformation and CO2 injection data for the in Salah gas sequestration project | |
US11802989B2 (en) | Systems and methods for generating vertical and lateral heterogeneity indices of reservoirs | |
CN105474272A (zh) | 2.75d网格划分算法 | |
Bachi et al. | An Efficient Hydraulic Fracture Geometry Calibration Workflow Using Microseismic Data |
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: 20160608 |