CN105659113A - 用于模拟天然裂缝网络调节器井筒观察的3d几何形状的地质统计程序 - Google Patents
用于模拟天然裂缝网络调节器井筒观察的3d几何形状的地质统计程序 Download PDFInfo
- Publication number
- CN105659113A CN105659113A CN201380078588.6A CN201380078588A CN105659113A CN 105659113 A CN105659113 A CN 105659113A CN 201380078588 A CN201380078588 A CN 201380078588A CN 105659113 A CN105659113 A CN 105659113A
- Authority
- CN
- China
- Prior art keywords
- crack
- length
- trilateral
- target
- computer implemented
- 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 title claims abstract description 81
- 238000004088 simulation Methods 0.000 title claims abstract description 28
- 238000003384 imaging method Methods 0.000 claims description 32
- 238000009826 distribution Methods 0.000 claims description 22
- 230000004044 response Effects 0.000 claims description 15
- 238000007405 data analysis Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000002372 labelling Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 abstract description 4
- 230000001143 conditioned effect Effects 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 20
- 238000003860 storage Methods 0.000 description 13
- 238000004891 communication Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000013500 data storage Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Processing Or Creating Images (AREA)
- User Interface Of Digital Computer (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
公开的实施方案包括用于提供用于通过井筒观察调节的天然裂缝网络的3D几何形状的模拟的地质统计程序的方法、设备和计算机程序产品。
Description
发明背景
1.发明领域
本发明一般来说涉及计算机化储层建模的领域,且更具体地说,涉及被配置用来提供用于天然裂缝网络的三维(3D)几何形状的条件模拟的地质统计程序的系统和方法。
2.现有技术的论述
在石油和天然气行业中,储层建模涉及建构石油储层的计算机模型以用于改进对储量的估计和关于该领域的发展做出决定的目的。举例来说,可创建地质模型以在生产之前提供对储层的静态描述。相比之下,可创建储层模拟模型以模拟在储层生产寿命内在储层内的流体的流动。
储层建模的一项挑战是对储层内的裂缝建模,这需要对流动特性、裂缝网络连接度和裂缝-基质交互的透彻理解。裂缝的正确建模是重要的,因为裂缝的性质(诸如空间分布、孔径、长度、高度、传导率和连接度)会显著影响储层流体到井筒的流动。
附图简述
下文参看附图详细描述本发明的说明性实施方案,附图以引用的方式并入本文中且其中:
图1A至图1C是说明根据公开的实施方案的用于天然裂缝网络的3D几何形状的条件模拟的方法的流程图;
图2是说明根据公开的实施方案的用于根据数据分析和地质类比来确定模拟的裂缝参数的算法的图,数据分析和地质类比用于执行天然裂缝网络的3D几何形状的条件模拟;
图3是说明根据公开的实施方案的用于根据井筒中的成像测井来确定实际裂缝参数的算法的图,井筒中的成像测井用于执行天然裂缝网络的3D几何形状的条件模拟;以及
图4是说明用于实现公开的实施方案的系统的一个实施方案的框图。
详细描述
通过参看附图的图1至图4来最佳地理解公开的实施方案和其优点,相同数字用于各图的相同和对应部件。在解释以下附图和详细描述后,公开的实施方案的其它特征和优点对本领域技术人员将为或将变得显而易见。希望所有此类额外特征和优点包括在公开的实施方案的范围内。另外,说明的附图仅为示例性的,且无意声明或暗示关于可实现不同实施方案的环境、架构、设计或过程的任何限制。
图1A至图1C是说明根据公开的实施方案的用于提供天然裂缝网络的3D几何形状的条件模拟的方法100的流程图。方法100在步骤102处以设置根据数据分析和地质类比确定的模拟的裂缝参数开始。举例来说,图2说明根据公开的实施方案的用于根据数据分析和地质类比来确定模拟的裂缝参数的算法200,数据分析和地质类比用于执行天然裂缝网络的3D几何形状的条件模拟。如图2中所说明,在一个实施方案中,用于设置模拟的裂缝参数的算法200包括确定将被模拟的裂缝的最小尺寸(M)以及也确定裂缝集合的数目(NS)。
对于裂缝集合(i)中的每一者,算法200确定集合中的裂缝的走向长度(l)的累积频率分布(fl,i)。另外,算法200确定集合中的裂缝的走向长度与倾斜长度的比(r)的概率分布(fr,i)。算法200还确定集合中的裂缝的走向角(θ)的概率分布(fθ,i)和集合中的裂缝的倾斜角(δ)的概率分布(fδ,i)。另外,对于裂缝集合(i)中的每一者,算法200确定集合中的裂缝的裂缝聚类程度(Ci)和集合中的裂缝的裂缝表面的平滑程度(Si)。在一个实施方案中,裂缝的聚类程度被指派范围从0到1中的值,其中1指示集合内的裂缝的最高聚类程度。类似地,在一个实施方案中,裂缝表面的平滑程度可被指派范围从0到1中的值,其中1指示集合内的裂缝表面的最高平滑程度。
另外,对于每一对裂缝集合(i,j),算法200确定集合i中的裂缝对集合j中的裂缝进行截断的概率(Prob{Tij})。或者,在某些实施方案(未描绘)中,算法200可确定网格化辅助属性与局部裂缝密度之间的相关系数rs。
返回参看图1,在步骤102之后,方法100在步骤104处根据井筒中的成像测井来确定实际裂缝参数。举例来说,如图3中所说明,在一个实施方案中,方法100可执行用于确定实际裂缝参数的算法300。算法300可将从一个或多个井筒中取得的一个或多个成像测井作为输入。算法300根据成像测井确定与一个或多个井筒相关联的观察到的裂缝的数目。对于每一观察到的裂缝,算法确定裂缝与井筒交叉的位置坐标(xi,yi,zi),且还确定在裂缝与井筒的交叉点处裂缝的走向长度(θi)和倾斜长度(δi)。另外,算法300包括用于确定成像测井经以提供数据的成像间隔的数目的指令。对于每一图像间隔,算法300确定成像间隔的端点的坐标(x1,i,y1,i,z1,i)和(x2,i,y2,i,z2,i)。
一旦确定了所有参数,方法100在步骤106处确定裂缝的种子点的位置。举例来说,在一个实施方案中,方法以观察到的裂缝的位置开始,且添加额外的未观察到的裂缝的可能选择的位置以实现每一裂缝集合中的目标频率FL,i(M)。在一个实施方案中,位置的选择由参数Ci(裂缝集合内的裂缝的聚类程度)控制。在某些实施方案中,位置的选择也可由辅助数据的网格(如果提供的话)控制。
在步骤108处,方法100将目标走向长度(SL)和目标倾斜长度(DL)指派给每一种子点。在一个实施方案中,方法100通过从集合中的裂缝的走向长度(l)的累积频率分布(Fl,i)和从集合中的裂缝的走向长度与倾斜长度的比(r)的概率分布(fr,i)中抽取随机值来指派目标走向长度和目标倾斜长度。方法100在步骤110处使用观察到的位置处的已知定向,且从集合中的裂缝的走向角的概率分布(fθ)以及集合中的裂缝的倾斜角的概率分布(fδ)中随机地抽取来将初始走向(θ)和倾斜(δ)指派给每一种子点。
在步骤112处,方法100确定任何未观察到的裂缝是否与任何成像间隔不一致地交叉。如果方法100确定未观察到的裂缝与成像间隔不一致地交叉,那么方法100在步骤114处更改SL、DL、θ和δ的指派的值。在一个实施方案中,值中的每一者可用从相应概率分布中随机抽取的不同值替换。或者,在一些实施方案中,方法100可确定仅更改某些值。举例来说,在一个实施方案中,方法100可确定仅更改初始倾斜(δ)。而且,在一些实施方案中,方法100可执行算法以确定对值中的一个或多个进行调整,而不是用随机选择的值来替换一个或多个值。方法100接着重复步骤112且基于更改的值确定任何未观察到的裂缝是否与任何成像间隔不一致地交叉。
如果方法100在步骤112处确定没有未观察到的裂缝与成像间隔中的任一者不一致地交叉,那么方法100如图1B中所说明在步骤120处在具有指定定向的每一种子点处定位尺寸为M(即,将被模拟的裂缝的最小尺寸)的三角形。方法100在步骤122处将所有三角形边缘标记为开放的以用于增长。在步骤124处,方法100随机地选择开放边缘且将尺寸为M的新三角形添加至所选择的开放边缘。在一个实施方案中,方法100使用克里格(kriging)方法以使用邻近三角形作为调节数据来模拟新三角形的定向。
当执行克里格方法时,方法100检查可能发生的三种可能情形中的一者。第一种,如果在步骤128处裂缝的走向长度达到目标走向长度(SL),那么方法100在步骤136处将所有水平边缘标记为对增长闭合。第二种,如果在步骤132处倾斜长度达到目标倾斜长度(DL),那么方法100在步骤138处将所有垂直边缘标记为对增长闭合。第三种情形是在步骤134处三角形是否触碰不同的裂缝。如果此情形发生,那么在一个实施方案中,方法100使用当前裂缝集合中的裂缝对另一裂缝集合进行截断的概率(Prob{Tij})来决定是否应截断扩展的裂缝。如果方法100确定应截断扩展的裂缝,那么方法100闭合达到裂缝的所有边缘,裂缝在边缘处终止。
一旦方法100基于以上情形中的一者的发生而闭合开放边缘,那么方法100在步骤142处确定任何开放边缘是否有剩余,且如果有,那么方法100对剩余的开放边缘重复步骤124处的过程。一旦闭合了所有边缘,方法100将三角测量写入至输出文件。方法100可在步骤152处执行额外后处理以添加孔径φ和k。
在步骤154处,方法100确定是否执行另一实现。如果需要额外实现,那么方法100重复在步骤106处开始的以上过程。如果不需要额外实现,那么方法100终止。
图4是说明用于实现公开的实施方案的特征和功能的系统400的一个实施方案的框图。系统400除了其它组件之外包括处理器400、主存储器402、辅助存储单元404、输入/输出接口模块406和通信接口模块408。处理器400可以是任何类型或任何数目的单核或多核处理器,其能够执行用于执行公开的实施方案的特征和功能的指令。
输入/输出接口模块406使得系统400能够(例如,从键盘和鼠标)接收用户输入且将信息输出至一个或多个装置,诸如但不限于打印机、外部数据存储装置和音频扬声器。系统400可任选地包括单独的显示模块410,其使得信息能够显示在集成的或外部显示装置上。举例来说,显示模块410可包括用于提供与一个或多个显示装置相关联的增强的图形、触摸屏和/或多点触控功能的指令或硬件(例如,图形卡或芯片)。举例来说,在一个实施方案中,显示模块410是QuadroFX类型的图形卡,其使得能够观看和操纵三维物体。
主存储器402是易失性存储器,其存储当前执行的指令/数据或被预取用于执行的指令/数据。辅助存储单元404是用于存储永久数据的非易失性存储器。辅助存储单元404可以是或可包括任何类型的数据存储组件,诸如硬盘驱动器、闪存驱动器或存储卡。在一个实施方案中,辅助存储单元404存储用于使得用户能够执行公开的实施方案的特征和功能的计算机可执行代码/指令和其它相关数据。
举例来说,根据公开的实施方案,辅助存储单元404可永久地存储算法的可执行代码/指令420,算法用于提供如上文所描述的通过井筒观察调节的天然裂缝网络的3D几何形状的条件模拟的地质统计程序。接着在处理器400的执行期间将与算法420相关联的指令从辅助存储单元404加载至主存储器402以用于执行公开的实施方案。另外,辅助存储单元1104可存储其它可执行代码/指令和数据422,诸如但不限于用于与公开的实施方案一起使用的储层模拟应用。
通信接口模块408使得系统400能够与通信网络430进行通信。举例来说,网络接口模块408可包括网络接口卡和/或无线收发器,其用于使得系统400能够通过通信网络430和/或直接与其它装置发送和接收数据。
通信网络430可以是任何类型的网络,包括以下网络中的一个或多个的组合:广域网、局域网、一个或多个专用网络、因特网、电话网络,诸如公共交换电话网络(PSTN)、一个或多个蜂窝式网络,和无线数据网络。通信网络430可包括多个网络节点(未描绘),诸如路由器、网络接入点/网关、交换器、DNS服务器、代理服务器,和用于辅助装置之间的数据/通信的路由的其它网络节点。
举例来说,在一个实施方案中,系统400可与用于执行本发明的特征的一个或多个服务器434或数据库432交互。举例来说,系统400可根据公开的实施方案向数据库432询问钻井日志信息。在一个实施方案中,数据库432可利用可从LandmarkGraphicsCorporation购得的Open软件来有效地管理、评估和分析单个数据库中的广泛的一系列油田项目数据。另外,在某些实施方案中,系统400可充当一个或多个客户端装置的服务器系统或与一个或多个装置/计算系统(例如,群集器、网格)的对等通信或并行处理的对等系统。
尽管已描述关于以上实施方案的特定细节,但以上硬件和软件描述仅希望作为实例实施方案且无意限制公开的实施方案的结构或实现方式。举例来说,尽管未示出系统400的许多其它内部组件,但本领域技术人员将了解,此类组件和其互连是众所周知的。
另外,如上文概述的公开的实施方案的某些方面可体现在使用一个或多个处理单元/组件执行的软件中。技术的程序方面可被看作通常呈可执行代码和/或相关联的数据的形式的“产品”或“制造物品”,可执行代码和/或相关联的数据在机器可读媒体的类型上执行或体现。有形的非暂时性“存储”类型媒体(即,计算机程序产品)包括用于计算机、处理器或类似物,或其相关联的模块的存储器或其它存储装置中的任一者或全部,诸如各种半导体存储器、磁带驱动器、磁盘驱动器、光盘或磁盘和类似物,其可在任何时间提供存储以用于软件编程。
另外,附图中的流程图和框图说明根据本发明的各种实施方案的系统、方法和计算机程序产品的可能的实现方式的架构、功能性和操作。还应注意,在一些替代实现方式中,框图或所说明的伪代码中指出的功能、指令或代码可不按附图中指出的顺序发生。举例来说,实际上,连续示出的两个方框可基本上同时执行,或所述方框有时可按颠倒顺序执行,这取决于涉及的功能性。还将指出,框图和/或流程图说明中的每一方框,和框图和/或流程图说明中的方框的组合可由执行指定功能或动作的基于专用硬件的系统,或专用硬件与计算机指令的组合来实现。
因此,公开的实施方案提供用于提供用于通过井筒观察调节的天然裂缝网络的3D几何形状的条件模拟的地质统计程序的系统、计算机程序产品和方法。除了上文描述的实施方案之外,特定组合的许多实例在本公开的范围内,下文对其中一些进行详细描述。
一个实例实施方案是一种用于产生天然裂缝网络的三维(3D)几何形状的模拟的计算机实现的方法,所述方法包括:确定用于执行天然裂缝网络的3D几何形状的模拟的参数;确定裂缝的种子点的位置。方法将目标走向长度和目标倾斜长度指派给每一种子点。方法还将初始走向长度和初始倾斜长度指派给每一种子点。方法在每一种子点处定位具有将被模拟的裂缝的最小尺寸的三角形且将三角形的所有边缘标记为开放的以用于增长。当三角形具有至少一个开放边缘时,方法选择三角形的开放边缘且将具有将被模拟的裂缝的最小尺寸的新三角形添加至开放边缘,直到满足条件集合中的一个条件为止。
在一个实施方案中,条件集合包括裂缝的走向长度达到目标走向长度,裂缝的倾斜长度达到目标倾斜长度,和新三角形触碰不同的裂缝。响应于满足一个条件,方法闭合三角形的至少一个开放边缘。举例来说,在一个实施方案中,计算机实现的方法响应于确定条件集合中的被满足的一个条件是裂缝的走向长度达到目标走向长度,闭合三角形的所有水平边缘。在同一实施方案中,方法可响应于确定条件集合中的被满足的一个条件是裂缝的倾斜长度达到目标倾斜长度,闭合三角形的所有垂直边缘。
在以上实例实施方案的某些实施方案中,根据数据分析和地质类比来确定参数。在一些实施方案中,参数可包括以下中的一个或多个:将被模拟的裂缝的最小尺寸,裂缝集合的数目,每一裂缝集合中的裂缝的走向长度的累积频率分布,裂缝集合中的裂缝的走向长度与倾斜长度的比的概率分布,裂缝集合中的裂缝的走向角的概率分布和倾斜角的概率分布,裂缝集合中的裂缝的聚类程度,以及裂缝集合中的裂缝的平滑程度。在一些实施方案中,以上实例计算机实现的方法实施方案可进一步包括确定第一集合中的裂缝对第二集合中的裂缝进行截断的概率。
另外,在以上实例计算机实现的方法的一些实施方案中,根据至少一个井筒的至少一个成像测井来确定参数。另外,以上计算机实现的方法可包括根据至少一个井筒的至少一个成像测井来确定观察到的裂缝的数目。对于每一观察到的裂缝,计算机实现的方法确定观察到的裂缝与井筒交叉的位置坐标,还确定在观察到的裂缝与井筒的交叉点处观察到的裂缝的走向长度和倾斜长度。而且,在某些实施方案中,计算机实现的方法可包括以下步骤:确定成像测井经以提供数据的成像间隔的数目;以及对于每一图像间隔,确定对应于成像间隔的端点的坐标。
作为确定裂缝的种子点的位置的过程的一部分,在某些实施方案中,计算机实现的方法从观察到的裂缝的位置开始,且添加额外的未观察到的裂缝的可能选择的位置,直到对于裂缝集合达到目标频率为止。
在一些实施方案中,计算机实现的方法通过从裂缝集合中的裂缝的走向长度的累积频率分布和从裂缝集合中的裂缝的走向长度与倾斜长度的比的概率分布中分别抽取随机值来执行将目标走向长度和目标倾斜长度指派给每一种子点。
以上公开的计算机实现的方法的另一实施方案可包括:确定任何未观察到的裂缝是否与任何成像间隔不一致地交叉;以及响应于确定未观察到的裂缝与成像间隔不一致地交叉而更改目标走向长度、目标倾斜长度、初始走向长度和初始倾斜长度中的至少一者。
基于以上公开内容的第二实例实施方案是一种系统,其包括:至少一个处理器;以及至少一个存储器,其耦合至至少一个处理器且存储用于产生天然裂缝网络的三维(3D)几何形状的模拟的计算机可执行指令,计算机可执行指令包括用于执行以下操作的指令:确定用于执行天然裂缝网络的3D几何形状的模拟的参数;确定裂缝的种子点的位置;将目标走向长度和目标倾斜长度指派给每一种子点;将初始走向长度和初始倾斜长度指派给每一种子点;在每一种子点处定位具有将被模拟的裂缝的最小尺寸的三角形;将三角形的所有边缘标记为开放的以用于增长;以及当三角形具有至少一个开放边缘时,选择三角形的开放边缘,重复将具有将被模拟的裂缝的最小尺寸的新三角形添加至开放边缘的步骤,直到满足条件集合中的一个条件为止,条件集合包括裂缝的走向长度达到目标走向长度,裂缝的倾斜长度达到目标倾斜长度,和新三角形触碰不同的裂缝,且响应于满足一个条件,闭合三角形的至少一个边缘。
而且,另一实例是一种非暂时性计算机可读媒体,其包括用于产生天然裂缝网络的三维(3D)几何形状的模拟的计算机可执行指令,计算机可执行指令在被执行时使一个或多个机器执行包括以下的操作:确定用于执行天然裂缝网络的3D几何形状的模拟的参数;确定裂缝的种子点的位置;将目标走向长度和目标倾斜长度指派给每一种子点;将初始走向长度和初始倾斜长度指派给每一种子点;在每一种子点处定位具有将被模拟的裂缝的最小尺寸的三角形;将三角形的所有边缘标记为开放的以用于增长;以及当三角形具有至少一个开放边缘时,选择三角形的开放边缘,重复将具有将被模拟的裂缝的最小尺寸的新三角形添加至开放边缘的步骤,直到满足条件集合中的一个条件为止,条件集合包括裂缝的走向长度达到目标走向长度,裂缝的倾斜长度达到目标倾斜长度,和新三角形触碰不同的裂缝,且响应于满足一个条件,闭合三角形的至少一个边缘。
以上第二和第三实例实施方案可类似地在上文关于第一实例实施方案描述的各种实施方案中进行修改。然而,以上特定实例实施方案和修改无意限制权利要求书的范围。举例来说,实例实施方案可通过包括、排除或组合本公开中描述的一个或多个特征或功能来修改。
如本文中所使用,单数形式“一”和“所述”希望也包括复数形式,除非上下文另外清楚地指示。应进一步理解,在本说明书和/或权利要求书中使用时,术语“包括(comprise和/或comprising)”指定所陈述特征、整数、步骤、操作、元件和/或组件的存在,但不排除一个或多个其它特征、整数、步骤、操作、元件、组件和/或其群组的存在或添加。所附权利要求书中的所有构件或步骤加功能元件的对应结构、材料、动作和等效物希望包括用于结合具体地要求的其它要求的元件执行功能的任何结构、材料或动作。本发明的描述已被呈现用于说明和描述的目的,但无意为详尽的或以公开的形式限制本发明。许多修改和变化对于本领域技术人员来说将为显而易见的,而不脱离本发明的范围和精神。选择和描述实施方案以解释本发明和实践应用的原理,且使本领域技术人员能够理解本发明的各种实施方案具有各种修改,所述修改适合于预期的特定用途。权利要求书的范围希望广泛地涵盖公开的实施方案和任何此修改。
Claims (20)
1.一种用于产生天然裂缝网络的三维(3D)几何形状的模拟的计算机实现的方法,所述方法包括:
确定用于执行所述天然裂缝网络的所述3D几何形状的所述模拟的参数;
确定裂缝的种子点的位置;
将目标走向长度和目标倾斜长度指派给每一种子点;
将初始走向长度和初始倾斜长度指派给每一种子点;
在每一种子点处定位具有将被模拟的所述裂缝的最小尺寸的三角形;
将所述三角形的所有边缘标记为开放的以用于增长;以及
当所述三角形具有至少一个开放边缘时,选择所述三角形的开放边缘,重复将具有将被模拟的所述裂缝的所述最小尺寸的新三角形添加至所述开放边缘的步骤,直到满足条件集合中的一个条件为止,所述条件集合包括所述裂缝的走向长度达到所述目标走向长度,所述裂缝的倾斜长度达到所述目标倾斜长度,和所述新三角形触碰不同的裂缝,且响应于满足所述一个条件而闭合所述三角形的至少一个边缘。
2.如权利要求1所述的计算机实现的方法,其中根据数据分析和地质类比来确定所述参数。
3.如权利要求2所述的计算机实现的方法,其中所述参数包括将被模拟的所述裂缝的所述最小尺寸和裂缝集合的数目。
4.如权利要求3所述的计算机实现的方法,其进一步包括确定每一裂缝集合中的所述裂缝的所述走向长度的累积频率分布和所述裂缝集合中的所述裂缝的所述走向长度与所述倾斜长度的比的概率分布。
5.如权利要求4所述的计算机实现的方法,其进一步包括确定所述裂缝集合中的所述裂缝的走向角的概率分布和倾斜角的概率分布。
6.如权利要求5所述的计算机实现的方法,其进一步包括确定所述裂缝集合中的所述裂缝的聚类程度和平滑程度。
7.如权利要求6所述的计算机实现的方法,其进一步包括确定第一集合中的裂缝对第二集合中的裂缝进行截断的概率。
8.如权利要求2所述的计算机实现的方法,其中根据至少一个井筒的至少一个成像测井来进一步确定所述参数。
9.如权利要求8所述的计算机实现的方法,其进一步包括:
根据所述至少一个井筒的所述至少一个成像测井来确定观察到的裂缝的数目;以及
对于每一观察到的裂缝,确定观察到的裂缝与井筒交叉的位置坐标,确定在所述观察到的裂缝与所述井筒的交叉点处所述观察到的裂缝的所述走向长度和所述倾斜长度。
10.如权利要求9所述的计算机实现的方法,其进一步包括:
确定成像测井经以提供数据的成像间隔的数目;以及
对于每一图像间隔,确定对应于所述成像间隔的端点的坐标。
11.如权利要求1所述的计算机实现的方法,其中确定所述裂缝的种子点的所述位置包括从观察到的裂缝的位置开始以及添加额外未观察到的裂缝的可能选择的位置直到对于裂缝集合达到目标频率为止。
12.如权利要求1所述的计算机实现的方法,其中通过从裂缝集合中的所述裂缝的所述走向长度的累积频率分布和从所述裂缝集合中的所述裂缝的所述走向长度与所述倾斜长度的比的概率分布中抽取随机值来执行将所述目标走向长度和所述目标倾斜长度指派给每一种子点。
13.如权利要求1所述的计算机实现的方法,其进一步包括确定任何未观察到的裂缝是否与任何成像间隔不一致地交叉;以及响应于确定未观察到的裂缝与成像间隔不一致地交叉而更改所述目标走向长度、所述目标倾斜长度、所述初始走向长度和所述初始倾斜长度。
14.如权利要求1所述的计算机实现的方法,其中响应于满足所述一个条件而闭合所述三角形的至少一个边缘包括响应于确定所述条件集合中的被满足的所述一个条件是所述裂缝的所述走向长度达到所述目标走向长度,闭合所述三角形的所有水平边缘。
15.如权利要求1所述的计算机实现的方法,其中响应于满足所述一个条件而闭合所述三角形的至少一个边缘包括响应于确定所述条件集合中的被满足的所述一个条件是所述裂缝的所述倾斜长度达到所述目标倾斜长度,闭合所述三角形的所有垂直边缘。
16.一种系统,其包括:
至少一个处理器;以及
至少一个存储器,其耦合至所述至少一个处理器且存储用于产生天然裂缝网络的三维(3D)几何形状的模拟的计算机可执行指令,所述计算机可执行指令包括用于执行以下操作的指令:
确定用于执行所述天然裂缝网络的所述3D几何形状的所述模拟的参数;
确定裂缝的种子点的位置;
将目标走向长度和目标倾斜长度指派给每一种子点;
将初始走向长度和初始倾斜长度指派给每一种子点;
在每一种子点处定位具有将被模拟的所述裂缝的最小尺寸的三角形;
将所述三角形的所有边缘标记为开放的以用于增长;以及
当所述三角形具有至少一个开放边缘时,选择所述三角形的开放边缘,重复将具有将被模拟的所述裂缝的所述最小尺寸的新三角形添加至所述开放边缘的步骤,直到满足条件集合中的一个条件为止,所述条件集合包括所述裂缝的走向长度达到所述目标走向长度,所述裂缝的倾斜长度达到所述目标倾斜长度,和所述新三角形触碰不同的裂缝,且响应于满足所述一个条件而闭合所述三角形的至少一个边缘。
17.如权利要求16所述的系统,其中所述参数是根据至少一个井筒的至少一个成像测井确定的,且其中所述计算机可执行指令进一步包括用于执行以下操作的指令:
根据所述至少一个井筒的所述至少一个成像测井来确定观察到的裂缝的数目;以及
对于每一观察到的裂缝,确定观察到的裂缝与井筒交叉的位置坐标,确定在所述观察到的裂缝与所述井筒的交叉点处所述观察到的裂缝的所述走向长度和所述倾斜长度。
18.如权利要求17所述的系统,其中所述计算机可执行指令进一步包括用于执行以下操作的指令:
确定成像测井经以提供数据的成像间隔的数目;以及
对于每一图像间隔,确定对应于所述成像间隔的端点的坐标。
19.如权利要求16所述的系统,其中所述计算机可执行指令进一步包括用于执行以下操作的指令:确定任何未观察到的裂缝是否与任何成像间隔不一致地交叉;以及响应于确定未观察到的裂缝与成像间隔不一致地交叉而更改所述目标走向长度、所述目标倾斜长度、所述初始走向长度和所述初始倾斜长度。
20.一种非暂时性计算机可读媒体,其包括用于产生天然裂缝网络的三维(3D)几何形状的模拟的计算机可执行指令,所述计算机可执行指令在被执行时使一个或多个机器执行包括以下的操作:
确定用于执行所述天然裂缝网络的所述3D几何形状的所述模拟的参数;
确定裂缝的种子点的位置;
将目标走向长度和目标倾斜长度指派给每一种子点;
将初始走向长度和初始倾斜长度指派给每一种子点;
在每一种子点处定位具有将被模拟的所述裂缝的最小尺寸的三角形;
将所述三角形的所有边缘标记为开放的以用于增长;以及
当所述三角形具有至少一个开放边缘时,选择所述三角形的开放边缘,重复将具有将被模拟的所述裂缝的所述最小尺寸的新三角形添加至所述开放边缘的步骤,直到满足条件集合中的一个条件为止,所述条件集合包括所述裂缝的走向长度达到所述目标走向长度,所述裂缝的倾斜长度达到所述目标倾斜长度,和所述新三角形触碰不同的裂缝,且响应于满足所述一个条件而闭合所述三角形的至少一个边缘。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2013/057639 WO2015030809A1 (en) | 2013-08-30 | 2013-08-30 | A geostatistical procedure for simulation of the 3d geometry of a natural fracture network conditioner well bore observations |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105659113A true CN105659113A (zh) | 2016-06-08 |
Family
ID=52587155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380078588.6A Pending CN105659113A (zh) | 2013-08-30 | 2013-08-30 | 用于模拟天然裂缝网络调节器井筒观察的3d几何形状的地质统计程序 |
Country Status (11)
Country | Link |
---|---|
US (1) | US9645281B2 (zh) |
CN (1) | CN105659113A (zh) |
AR (1) | AR097499A1 (zh) |
AU (1) | AU2013399054B2 (zh) |
BR (1) | BR112016001586A2 (zh) |
CA (1) | CA2919871C (zh) |
DE (1) | DE112013007368T5 (zh) |
GB (1) | GB2537004B (zh) |
MX (1) | MX2016000791A (zh) |
SG (1) | SG11201600642XA (zh) |
WO (1) | WO2015030809A1 (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017082871A1 (en) * | 2015-11-10 | 2017-05-18 | Landmark Graphics Corporation | Target object simulation using orbit propagation |
WO2017082870A1 (en) | 2015-11-10 | 2017-05-18 | Landmark Graphics Corporation | Fracture network triangle mesh adjustment |
CN109648693A (zh) * | 2018-12-20 | 2019-04-19 | 昆明理工大学 | 一种基于3d打印技术的充填型裂隙网络岩石试件制备方法 |
WO2020167282A1 (en) | 2019-02-13 | 2020-08-20 | Landmark Graphics Corporation | Planning a well configuration using geomechanical parameters |
CN113343423B (zh) * | 2021-05-08 | 2022-04-29 | 武汉大学 | 基于强度空间变异性的随机裂隙网络生成方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050015231A1 (en) * | 1998-05-04 | 2005-01-20 | Schlumberger Technology Corporation | Near wellbore modeling method and apparatus |
US20090194273A1 (en) * | 2005-12-01 | 2009-08-06 | Surjaatmadja Jim B | Method and Apparatus for Orchestration of Fracture Placement From a Centralized Well Fluid Treatment Center |
US20100066560A1 (en) * | 2007-09-12 | 2010-03-18 | Hexion Specialty Chemicals, Inc. | Wellbore casing mounted device for determination of fracture geometry and method for using same |
US8437962B2 (en) * | 2009-11-25 | 2013-05-07 | Halliburton Energy Services, Inc. | Generating probabilistic information on subterranean fractures |
US20130124166A1 (en) * | 2011-11-15 | 2013-05-16 | Jack Gammill Clemens | Modeling operation of a tool in a wellbore |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2002239619A1 (en) * | 2000-12-08 | 2002-06-18 | Peter J. Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
US7188092B2 (en) * | 2002-07-12 | 2007-03-06 | Chroma Energy, Inc. | Pattern recognition template application applied to oil exploration and production |
FR2870621B1 (fr) * | 2004-05-21 | 2006-10-27 | Inst Francais Du Petrole | Methode pour generer un maillage hybride conforme en trois dimensions d'une formation heterogene traversee par une ou plusieurs discontinuites geometriques dans le but de realiser des simulations |
US7480205B2 (en) * | 2005-04-20 | 2009-01-20 | Landmark Graphics Corporation | 3D fast fault restoration |
GB2465504C (en) * | 2008-06-27 | 2019-12-25 | Rasheed Wajid | Expansion and sensing tool |
US20110029293A1 (en) * | 2009-08-03 | 2011-02-03 | Susan Petty | Method For Modeling Fracture Network, And Fracture Network Growth During Stimulation In Subsurface Formations |
US8392165B2 (en) * | 2009-11-25 | 2013-03-05 | Halliburton Energy Services, Inc. | Probabilistic earth model for subterranean fracture simulation |
CN102918423B (zh) * | 2010-05-28 | 2016-09-07 | 埃克森美孚上游研究公司 | 用于地震烃体系分析的方法 |
US9442205B2 (en) * | 2011-03-23 | 2016-09-13 | Global Ambient Seismic, Inc. | Method for assessing the effectiveness of modifying transmissive networks of natural reservoirs |
-
2013
- 2013-08-30 SG SG11201600642XA patent/SG11201600642XA/en unknown
- 2013-08-30 CN CN201380078588.6A patent/CN105659113A/zh active Pending
- 2013-08-30 DE DE112013007368.4T patent/DE112013007368T5/de not_active Withdrawn
- 2013-08-30 BR BR112016001586A patent/BR112016001586A2/pt not_active IP Right Cessation
- 2013-08-30 US US14/894,320 patent/US9645281B2/en active Active
- 2013-08-30 GB GB1601200.7A patent/GB2537004B/en active Active
- 2013-08-30 CA CA2919871A patent/CA2919871C/en active Active
- 2013-08-30 MX MX2016000791A patent/MX2016000791A/es unknown
- 2013-08-30 WO PCT/US2013/057639 patent/WO2015030809A1/en active Application Filing
- 2013-08-30 AU AU2013399054A patent/AU2013399054B2/en not_active Ceased
-
2014
- 2014-08-29 AR ARP140103243A patent/AR097499A1/es unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050015231A1 (en) * | 1998-05-04 | 2005-01-20 | Schlumberger Technology Corporation | Near wellbore modeling method and apparatus |
US20090194273A1 (en) * | 2005-12-01 | 2009-08-06 | Surjaatmadja Jim B | Method and Apparatus for Orchestration of Fracture Placement From a Centralized Well Fluid Treatment Center |
US20100066560A1 (en) * | 2007-09-12 | 2010-03-18 | Hexion Specialty Chemicals, Inc. | Wellbore casing mounted device for determination of fracture geometry and method for using same |
US8437962B2 (en) * | 2009-11-25 | 2013-05-07 | Halliburton Energy Services, Inc. | Generating probabilistic information on subterranean fractures |
US20130124166A1 (en) * | 2011-11-15 | 2013-05-16 | Jack Gammill Clemens | Modeling operation of a tool in a wellbore |
Non-Patent Citations (5)
Title |
---|
张淑娟等: "蚂蚁追踪技术在潜山油藏裂缝预测中的应用", 《断块油气田》 * |
王浩明等: "一种基于三角形区域生长的图像篡改检测方法", 《武汉大学学报(理学版)》 * |
秦巍等: "成像测井井壁图像裂缝自动识别", 《测井技术》 * |
艾廷华等: "基于三角网的"种子法"多边形生成", 《武汉大学学报信息科学版》 * |
赵军等: "基于图像分割的成像测井资料目标拾取与计算", 《地球物理学进展》 * |
Also Published As
Publication number | Publication date |
---|---|
AU2013399054B2 (en) | 2017-08-31 |
US9645281B2 (en) | 2017-05-09 |
US20160209545A1 (en) | 2016-07-21 |
CA2919871C (en) | 2023-03-21 |
WO2015030809A8 (en) | 2016-02-18 |
CA2919871A1 (en) | 2015-03-05 |
MX2016000791A (es) | 2016-08-05 |
GB201601200D0 (en) | 2016-03-09 |
GB2537004B (en) | 2020-12-02 |
GB2537004A (en) | 2016-10-05 |
AR097499A1 (es) | 2016-03-16 |
WO2015030809A1 (en) | 2015-03-05 |
SG11201600642XA (en) | 2016-02-26 |
DE112013007368T5 (de) | 2016-05-12 |
BR112016001586A2 (pt) | 2017-08-29 |
AU2013399054A1 (en) | 2016-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2008357712B2 (en) | Distribution of properties in a 3D volumetric model using a maximum continuity field | |
AU2017202784B2 (en) | Gridless simulation of a fluvio-deltaic environment | |
CN105659113A (zh) | 用于模拟天然裂缝网络调节器井筒观察的3d几何形状的地质统计程序 | |
CN105684047A (zh) | 动态更新表示一个或多个地质构造的区划 | |
JP2017174393A (ja) | マルチ解像度画像システム | |
US10260318B2 (en) | Three-dimensional interactive wellbore model simulation system | |
US9378587B2 (en) | Systems and methods for creating a three-dimensional texture atlas | |
AU2013398344B2 (en) | Local updating of 3D geocellular model | |
RU2685005C1 (ru) | Способ и компьютерная система для проектирования размещения кустовых площадок на месторождении | |
GB2519630A (en) | Interactive visualization of reservoir simulation data sets | |
US9342923B2 (en) | Systems and methods for creating a three-dimensional texture atlas | |
CN105408941A (zh) | 通过优化去断层法在逆向断层区域中的全局网格构建 | |
WO2017213631A1 (en) | Systems and methods for unfaulting point clouds | |
CN107038176B (zh) | 网络图页面渲染方法、装置及设备 | |
Khishvand et al. | Developing an algorithm for reconstruction of blocky systems in discontinuous media: two-dimensional analysis | |
Van Aalsburg et al. | Interactive editing of digital fault models | |
KR101769028B1 (ko) | 3차원 지리 공간상에서의 객체 표현 방법 | |
WO2017053080A1 (en) | Subsurface volume evaluation |
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 |
Application publication date: 20160608 |
|
WD01 | Invention patent application deemed withdrawn after publication |