CN102239507B - 用于离散裂隙建模的网格生成***和方法 - Google Patents
用于离散裂隙建模的网格生成***和方法 Download PDFInfo
- Publication number
- CN102239507B CN102239507B CN200980148461.0A CN200980148461A CN102239507B CN 102239507 B CN102239507 B CN 102239507B CN 200980148461 A CN200980148461 A CN 200980148461A CN 102239507 B CN102239507 B CN 102239507B
- Authority
- CN
- China
- Prior art keywords
- crack
- grid
- polyhedral elements
- refinement
- interface
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000005516 engineering process Methods 0.000 claims description 17
- 239000012530 fluid Substances 0.000 claims description 17
- 239000011435 rock Substances 0.000 claims description 12
- 238000002347 injection Methods 0.000 claims description 8
- 239000007924 injection Substances 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 12
- 238000010586 diagram Methods 0.000 description 12
- 229930195733 hydrocarbon Natural products 0.000 description 12
- 150000002430 hydrocarbons Chemical class 0.000 description 12
- 238000003860 storage Methods 0.000 description 11
- 238000004088 simulation Methods 0.000 description 9
- 238000013507 mapping Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 241001415846 Procellariidae Species 0.000 description 1
- 208000037656 Respiratory Sounds Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000013065 commercial product Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- HJUFTIJOISQSKQ-UHFFFAOYSA-N fenoxycarb Chemical compound C1=CC(OCCNC(=O)OCC)=CC=C1OC1=CC=CC=C1 HJUFTIJOISQSKQ-UHFFFAOYSA-N 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Acoustics & Sound (AREA)
- Geochemistry & Mineralogy (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了生成精确网格的计算机实现***和方法,该精确网格可以用于模拟像有裂缝地下储层那样的高度复杂的地下储层。提供有裂缝地下储层和被表示在有裂缝地下储层内的裂缝的表示。有裂缝地下储层的表示被分解成多面体单元。选择多面体单元之间的界面来近似被显性表示在有裂缝地下储层内的离散裂缝。所选界面用于生成可以用于构建有裂缝地下储层的模拟模型的网格。
Description
相关申请的交叉参考
本专利申请要求2008年12月3日提交、通过引用全文并入本文中的美国临时专利申请第61/119,604号的利益。
技术领域
本发明一般涉及生成可以用于构建地下储层的模拟模型(simulationmodel)的网格的***和方法,更具体地说,涉及生成可以用于构建高度复杂的地下储层(例如有裂缝地下储层)的模拟模型的网格的***和方法。
背景技术
烃生产领域旨在取得陷在地下储层中的烃。通常,这些储层由各自通过不同的沉积学和流体性质表征的岩石和流体物质的平行层组成。烃积聚在非多孔或低渗透率岩石层下面或之间,形成储层。这些烃可以通过将井钻到储层中来开采。于是,烃能够从储层流到井中并向上流到地面。烃流到井中的生产率对于石油工业来说是至关重要的,因此,人们把大量心血花在开发技术上,以便更好地预测地下储层的流体流动和地质力学特性。这些技术之一涉及本文后面将更详细讨论的储层的网格化。
高度复杂的地质地下储层(例如具有裂缝网络的储层)对储层模拟提出了独特和专业挑战。具有裂缝网络的地下储层通常具有低渗透率岩石基体(rockmatrix),使烃难以通过地层。裂缝可以被描述成地层内的裂口和空洞,可以是自然产生的或从井筒中人工生成的。因此,裂缝的存在可能在使流体流过地层到达井方面起重要作用。例如,取决于油井是否遇到大裂缝,井的烃生产率往往差别很大。有时,将诸如水、化学物品、气体或它们的组合物的流体注入储层中,以帮助烃流入生产井。在裂缝允许生产井与流体注入井之间的直接连通的情况下,注入流体可能流过裂缝,绕过注入流体本来要帮助生产的地层内的大部分烃。因此,期望表征烃储层中裂缝的范围和取向,以便适当预测通过地下地层的地质力学和流体流动特性。为了计算这些特性,必须首先应用网格化技术。
储层网格化技术可以描述成将3D储层体(3Dreservoirvolume)分解成通常是凸3D体的多个更小和更简单的3D体的过程。于是,这些技术将连续模拟域分解成随后可以用于通过离散化描述流体流动、热传递、地质力学性质或它们的组合的控制方程来构建模拟模型的离散对应物。在储层模拟领域内,取决于正在使用的离散化和模拟技术,离散体通常被称为单元(cell)、有限体(finitevolume)、控制体(controlvolume)、或有限元。
对于有裂缝地下储层,由于裂缝网络的几何复杂性和随机性质,网格化提出了独特的挑战。例如,传统储层模拟的网格化策略通常未设计用来管理像裂缝那样的大量内部几何特征。许多网格化策略只适用于通常称为“不漏水”几何结构的、连通良好但不包括裂纹或重叠的几何结构。例如,这些网格化策略通常预先计算所有几何特征的相交部分。其它网格化策略在捕捉诸如稍微穿过另一个裂缝的平面的一个裂缝、相互接近但没有相交的两个裂缝、或以小夹角相交的两个裂缝的特定几何复杂体的时候不能取得良好网格质量。一般说来,包括内部特征的网格必须平衡精确地近似(approximate)特征和保持良好质量的相反目标。如果储层网格化技术不能在保持良好网格质量的同时适应这些细节,模型的模拟就可以导致精度下降、运行时间延长、收敛问题、或兼而有之。
现有技术试图通过预处理裂缝组以消除最有问题的构造来提高网格质量。但是本领域技术人员应该理解,由于消除一个问题往往引起另一个问题,所以这不是直截了当的手段。另外,由于在计算裂缝相交部分时需要考虑浮点运算,所以这些过去的尝试的实现是不容易的。
发明内容
根据本发明的一个方面,公开了一种生成可以用于构建有裂缝地下储层的模拟模型的网格的方法。该方法包括提供有裂缝地下储层(fracturedsubsurfacereservoir)和被表示在有裂缝地下储层内的裂缝的表示。构建具有在多面体单元的边界上定义单元间界面(cell-to-cellinterface)的多个多面体单元的网格。选择与裂缝相邻的界面,以形成近似裂缝的几何形状的连续界面路径。将连续界面路径与预定裂缝性质相关联。输出具有与预定裂缝性质相关联的连续界面路径的网格,以便构建有裂缝地下储层的模拟模型。
在一个或多个实施例中,提供一个井,并选择与该井相邻的单元间界面,以形成近似该井的几何形状的界面井区域。将该界面井区域与预定井性质相关联,使得输出网格包括与预定井性质相关联的界面井区域。
在一个或多个实施例中,在选择与裂缝相邻的多个单元间界面,以形成近似裂缝的几何形状的连续界面路径之前,细化至少一个多面体单元。在一个或多个实施例中,使用自适应网格细化技术来细化至少一个多面体单元。在一个或多个实施例中,使用迫近函数细化至少一个多面体单元。在一个或多个实施例中,该迫近函数确定每个多面体单元与诸如裂缝、注入井、和生产井的与每个多面体单元相邻的特征之间的距离。在一个或多个实施例中,多面体单元是具有最长边和与最长边相对的边的四面体。至少一个多面体单元通过用在最长边的中点和与最长边相对的边之间延伸的切面(planarcut)来分割四面体来细化。
在一个或多个实施例中,在输出网格之前合并两个或更多个多面体单元。
在一个或多个实施例中,优化节点的位置。
在一个或多个实施例中,通过在输出网格上执行控制方程的离散化构建模拟模型,并且利用储层模拟器模拟该模拟模型,以便直观显示有裂缝地下储层内的流体流动。
在一个或多个实施例中,预定裂缝性质定义网格内的裂缝-裂缝、裂缝-基体、和基体-基体连接。
本发明的另一个方面包括生成可以用于构建有裂缝地下储层的模拟模型的网格的***。该***包括用户控制接口、数据库、计算机处理器、和软件程序。该用户控制接口用于向***输入例如有裂缝地下储层的表示的信息。该数据库为***存储例如通过用户控制接口输入的有裂缝地下储层的表示的信息。该处理器为***执行操作。该软件程序包括多个模块,包括分解模块和界面选择模块。该分解模块将有裂缝地下储层的表示分解成在多面体单元之间的边界上定义单元间界面的多面体单元。该界面选择模块选择与裂缝相邻的界面,以形成近似被表示在有裂缝地下储层内的裂缝的连续界面路径。将连续界面路径与预定裂缝性质相关联。输出具有与预定裂缝性质相关联的连续界面路径的网格,以便构建有裂缝地下储层的模拟模型。
在一个或多个实施例中,该软件程序包括使用迫近函数自适应地细化至少一个多面体单元的细化模块。
在一个或多个实施例中,该软件程序包括合并两个或更多个多面体单元的合并模块。
在一个或多个实施例中,该软件程序包括优化网格内节点的位置的网格调整模块。
在一个或多个实施例中,该界面选择模块选择与井相邻的界面,以形成近似井的几何形状的界面井区域。
本发明的另一个方面包括存储在处理器可读媒体中,用于生成可以用于构建有裂缝地下储层的模拟模型的网格的软件程序。该软件程序被配置成构建具有在多面体单元的边界上定义单元间界面的多个多面体单元的网格,和选择与裂缝相邻的界面,以形成近似裂缝的几何形状的连续界面路径。将连续界面路径与预定裂缝性质相关联。输出具有与预定裂缝性质相关联的连续界面路径的网格,以便构建有裂缝地下储层的模拟模型。
附图说明
图1是例示根据本发明生成可以用于构建有裂缝地下储层的模拟模型的网格的方法的步骤的流程图。
图2是例示根据本发明生成可以用于构建有裂缝地下储层的模拟模型的网格的方法的步骤的流程图。
图3A-C是例示根据本发明生成可以用于构建有裂缝地下储层的模拟模型的网格的如图1和2所示的方法的步骤的区域的示意图。
图4是例示根据本发明的四面体的一系列细化的示意图。
图5是例示根据本发明生成可以用于构建有裂缝地下储层的模拟模型的网格的***的示意图。
图6A是根据本发明的有裂缝地下储层的区域的示意图。
图6B是根据本发明的近似例示在图6A中的区域的裂缝的多面体单元的界面的示意图。
图7A和7B是例示在图6A中的区域的各种分辨率的网格的示意图。
图8A是显示在图7A中的网格的映射图。
图8B是显示在图6A中的网格的映射图。
图8C是显示在图7B中的网格的映射图。
图9A是根据本发明的多面体单元的界面近似离散裂缝的区域的示意图。
图9B是根据本发明的例示在图9A中的区域的网格的示意图。
图10是根据本发明的包括生产和注入井的有裂缝地下储层的区域的示意图。
图11是根据本发明的例示在图10中的区域的网格的示意图。
具体实施方式
全球石油资源的很大部分储藏在有裂缝地下储层中。为了最佳地提取这些烃资源,必须利用与储层模拟有关的各种技术,例如储层网格化。例如,在各种储层模拟工作流中可以利用包括本发明的一些方面的储层网格化***和方法。
图1例示了根据本发明的一些方面的概述与储层网格化有关的步骤的方法10。具体地,采用这些步骤生成可以用于构建有裂缝地下储层的模拟模型的网格。该模拟模型可以表示整个有裂缝地下储层或有裂缝地下储层的特定部分。在步骤11中提供有裂缝地下储层和被表示在有裂缝地下储层内的显性离散裂缝的表示。该表示可以表示整个有裂缝地下储层或有裂缝地下储层的特定部分。在步骤13中将有裂缝地下储层的表示分解成多面体单元,以便创建多面体单元的网格。例如,可以将有裂缝地下储层分解成均匀、结构化四面体栅格。在多面体单元之间的边界上定义单元间界面,并且在步骤15中,选择近似被表示在有裂缝地下储层内的显性离散裂缝的多面体单元之间的界面。具体地,选择与裂缝相邻的单元间界面,以形成近似裂缝的几何形状的连续界面路径。正如本文所使用的那样,相邻意味着位于裂缝附近或非常接近裂缝,包括与裂缝毗邻或相交的单元间界面,以及不与裂缝相交而作为裂缝穿过的多面体单元或邻接多面体单元的边界的那些单元间界面。因此,所生成的网格(包括裂缝近似)与预定裂缝性质相关联,使得该网格定义网格单元之间的裂缝-裂缝、裂缝-基体、和基体-基体连接。所生成的网格(包括通过多面体单元界面定义的裂缝的近似)可以用于构建有裂缝地下储层的模拟模型。
图2例示了根据本发明的一些方面的概述与储层网格化有关的步骤的方法20。具体地,图2按照本发明的一些方面,描绘了可以用于生成网格的用虚线显示的可选步骤。在步骤21中提供有裂缝地下储层和被表示在有裂缝地下储层内的显性离散裂缝的表示。该表示可以表示整个有裂缝地下储层或有裂缝地下储层的特定部分。在步骤23中将有裂缝地下储层的表示分解成定义四面体单元的均匀、结构化四面体栅格。在步骤25中,可选地细化在显性离散裂缝附近或与显性离散裂缝相邻的四面体单元。本领域技术人员应该理解,也可以细化围绕诸如注入井和生产井的有裂缝地下储层的其它内部特征或与之相邻的四面体单元,本文后面将对此作更详细描述。在步骤27中,选择近似被表示在有裂缝地下储层内的显性离散裂缝的四面体单元之间的界面,以便生成可以用于构建有裂缝地下储层的模拟模型的网格。这两种方法10,20都可以称为裂缝的小尺度处理。
可以有选择地执行许多步骤,以便提高网格的质量,使网格更适用于数值计算,优化显性离散裂缝的近似,或其组合。例如,可以在步骤29中调整网格节点的位置,以便在保持网格质量的同时改进显性离散裂缝的近似。在步骤31中,可以修改网格布局,以便在保持显性离散裂缝的近似的同时提高网格质量。在现有技术中已知的各种技术,例如执行边缘翻转,可以用在步骤31中。在步骤33中可以合并四面体单元,以便在保持网格分辨率的同时减少网格单元的数量。例如,可以将四面体单元合并以形成四棱锥、三棱柱和六面体。在步骤35中,可以以适当格式输出生成的网格,以便对控制方程进行离散化来构建模拟模型。
图3A-C例示了区域40的映射图,例示了根据本发明的一些方面生成可以用于构建有裂缝地下储层的模拟模型的网格的如图1和2所示的方法的步骤。在图3A中描绘了有裂缝地下储层的区域40,它例示了有裂缝地下储层39和显性离散裂缝41的表示。具体地,在方法10的步骤11和方法20的步骤21中规定了这种有裂缝地下储层39和显性离散裂缝41的表示。如图3B所示,在方法10的步骤13和方法20的步骤23中用多面体单元分解有裂缝地下储层的区域40。在图3B中,像在方法20的步骤25中规定的那样,在显性离散裂缝41附近细化多面体单元,以便提高显性离散裂缝41附近的分辨率。网格节点43、单元界面45、和多面体单元47各自在图3B中都用标号指示。网格节点43是两个或更多个多面体单元47共享的顶点,单元界面45是两个或更多个多面体单元47共享的单元边缘或表面。在方法10的步骤15和方法20的步骤27中,选择近似被表示在有裂缝地下储层内的显性离散裂缝41的多面体单元47之间的界面,以便生成可以用于构建有裂缝地下储层的模拟模型的网格。图3C示出了近似被表示在有裂缝地下储层内的显性离散裂缝41的多面体单元界面49。
在方法10的步骤11和方法20的步骤21中提供的表示定义包括有裂缝地下储层39和被表示在有裂缝地下储层内的显性离散裂缝41的表示的区域。例如,这些步骤通常包括代表岩石基体的区域,以及描述裂缝网络的裂缝表示。可以将通常随机生成的裂缝表示输入地质模型中,以构建有裂缝地下储层的更真实表征。地质模型通常包括具有诸如渗透率和孔隙率分布的岩石性质,以及诸如初始流体饱和度分布的流体性质的储层的结构和地层框架。这些性质或参数可以从包括地震图像、岩芯、生产记录、井下测量值、钻探信息、和露头的各种来源中获得。
存在许多可用于构建地质模型或表示的商用产品,例如,总部在荷兰阿姆斯特丹的ParadigmGeotechnologyBV发布的EarthDecisionSuite(GOCADTM支持)、和来自总部在得克萨斯休斯敦的SchlumbergerLimited的PetrelTM。还存在许多用于构建裂缝表示的商用产品,例如,由总部在佐治亚州亚特兰大的GolderAssociatesInc.发布的FracManTM。本领域技术人员应该理解,可以作出相同地质情况的许多实现,导致各自通过准随机变化成形的相似但不相同的储层模型。储层模型的实现之间的这些差异是固有的,因为只能从地下储层中提取一定数量的确定信息,并且通常依赖于与所获取地质数据结合地应用概率方法来获取真实储层模型。
一旦提供了有裂缝地下储层39和显性离散裂缝41的表示,则在方法10的步骤13和方法20的步骤23中将区域分解成多面体单元47,例如均匀、结构化四面体栅格。连续储层模型的分解生成定义随后可以离散化的模型内的特定控制体的网眼或网格。有裂缝地下储层的区域可以以各种单元尺寸和形状分解,并且通常以能充分捕捉构建模拟模型所需的精细尺度细节的分辨率网格化。本领域技术人员应该理解,分解步骤以能够在随后步骤中容易地细化或合并多面体单元的方式执行。另外,所得到的网格被实现为使得它可以通过当前在现有技术中实施的方法(例如有限体离散化)来离散化。另外,取决于要研究的区域,生成的网格可以覆盖有裂缝地下储层的特定部分或整个有裂缝地下储层。
在方法10的步骤15和方法20的步骤27中,选择近似被表示在有裂缝地下储层内的显性离散裂缝的多面体单元之间的界面,以便生成可以用于构建有裂缝地下储层的模拟模型的网格。该界面通常是平面的或***面的;但是,在一些情况下,该界面可以是非平面。为了选择最佳地近似裂缝的多面体单元的界面,使用目标函数。例如,目标函数可以确定最接近裂缝的多面体单元,以便可以通过那些多面体单元之间的界面近似裂缝几何形状。本领域技术人员应该理解,可以将不同准则用作目标函数的输入。例如,该输入可以只允许来自多面体单元的一个边缘被用来定义裂缝的几何形状,只允许被选择来定义裂缝的几何形状的多面体单元的两个相邻边缘之间的关联角的预定范围,根据与每个裂缝的端点最接近的节点之间的加权最短路径选择多面体单元的边缘,或它们的组合。而且,接近度可以通过裂缝与界面之间的欧几里得距离来确定,或者,也可以通过非米制或非欧几里得测量值形式来确定。在一些情况下,可以预处理通过裂缝的所有水平片(slice),以便找出包括最佳地近似裂缝几何形状的多面体单元之间的界面的示范性水平片。可以使示范性水平片中多面体单元之间的界面垂直延伸,以将裂缝定义为水平层内的平坦长方形结构。另外,有裂缝地下储层的区域可以包括诸如注入井或生产井的一个或多个井,并且可以类似地通过多面体单元的界面来近似这些井。
一旦为有裂缝地下储层生成网格,就可以将它离散化,并且可以使用储层模拟器进行储层内的流体流动的模拟。存在许多用于进行储层模拟的商用产品,例如,Chevron的专有CHEARSTM模拟软件包,或Schlumberger的ECLIPSETM储层模拟器。另外,Chevron和Schlumberger共同拥有的INTERSECTTM是可以用于模拟储层内的流体流动的另一种专有储层模拟软件包。
为有裂缝地下储层生成网格最好基本上是自动的。另外,网格分辨率可以是空间自适应的。例如,可以在裂缝和井附近使用细化单元,而在其它地方使用更大的多面体单元以便提高效率。如果单元是四面体状的,使得它们具有最长边和与最长边相对的边,那么可以通过用在最长边的中点和与最长边相对的边之间延伸的切面分割四面体单元来进行细化。
图4描绘了四面体的一系列细化,如图所示,四面体可以描述成定义每个立方体六个四面体单元的叠加均匀网格的一个单元。如图所示,通过用在最长边59的中点57和与最长边59相对的边61之间延伸的切面分割四面体51,可以将四面体51细化成四面体53和四面体55。类似地,通过使用这种“最长边”细化技术可以将四面体53进一步细化成四面体63和四面体65。四面体63可以被进一步细化成四面体67和四面体69。
本领域技术人员应该理解,可以一直细化到裂缝或井附近的多面体单元被细化成预定数量。而且,在构建有裂缝地下储层的模拟模型之前,可以合并两个或更多个多面体单元形成更大的多面体单元,以减少网格单元的数量。在一些实施例中,可以使用单个输入参数调整网格生成***和方法,以便在网格单元质量与裂缝表示的近似之间取得适当折衷。按照本发明一些方面的网格生成***和方法也不需要预处理裂缝组或显性计算裂缝相交部分。
图5例示了根据本发明的一些方面为有裂缝地下储层生成网格的***100。***100包括用户接口110,以便操作人员可以主动输入信息和观察***100的运行。用户接口110可以是诸如键盘、鼠标、触摸屏显示器、或包括个人数字助理(PDA)在内的手持图形用户接口(GUI)的使得人们能够与***100交互的任何手段。通过用户接口110进入***100的输入可以存储在数据库120中。另外,***100生成的信息也可以存储在数据库120中。例如,数据库120可以存储有裂缝地下储层的表示121,以及例如网格123的***生成信息。
***100包括执行多种操作的网格化软件130。正如本文更详细讨论的那样,网格化软件130包括分解模块131、界面选择模块133、细化模块135、合并模块137、和网格调整模块139。本领域技术人员应该理解,***100的运行不需要所有这些模块。例如,在一些情况下,网格化软件130可以只利用分解模块131和界面选择模块133来生成网格。处理器140为***100解释执行操作的指令,以及响应预定条件生成运行网格化软件130的自动指令。来自用户接口110和网格化软件130两者的指令被处理器140处理,供***100运行用。
在某些实施例中,***100可以包括向操作人员或向其它***(未示出)提供信息的报告单元150。例如,报告单元150可以是打印机、显示屏、或数据存储设备。但是,应该明白,***100无需包括报告单元150,可替代地,用户接口110可以用于向操作人员报告***100的信息。例如,可以使用监视器或像包括个人数字助理(PDA)在内的手持图形用户接口(GUI)那样的用户接口设备向用户直观显示输出。
诸如用户接口110、数据库120、网格化软件130、处理器140和报告单元150的***100的任何部件之间的通信最好在通信网络160上传送。通信网络160可以是实现信息传送的任何手段。这种通信网络160的例子目前包括计算机内的交换器、个人区域网(PAN)、局域网(LAN)、广域网(WAN)、和全球范围网络(GAN),但不局限于这些。通信网络160也可以包括诸如光纤或无线射频的用于连接网络中的各个设备的任何硬件技术。
在***100运行时,操作人员通过用户接口110输入数据(例如描述有裂缝地下储层的有裂缝地下储层的表示121),以便存储在数据库120中。然后,启动网格化软件130为模拟模型生成网格123。网格化软件130使用分解模块131,利用多面体单元分解有裂缝地下储层的表示121。网格化软件130使用界面选择模块133选择近似被表示在有裂缝地下储层内的显性离散裂缝的多面体单元之间的界面,以便生成可以用于构建有裂缝地下储层的模拟模型的网格。
细化模块135、合并模块137、和网格调整模块139可以用于在网格单元质量与裂缝表示的近似之间取得更适当折衷。细化模块135能够细化裂缝和井附近的多面体单元,使得界面选择模块133可以更好地用多面体单元的界面来近似裂缝和井的几何形状。合并模块137可以组合两个或更多个多面体单元,以便通过减少模拟模型所需的数值计算量来提高网格的效率。网格调整模块139能够在平衡网格的质量和效率的同时优化裂缝或井的近似。例如,网格调整模块139可以扫描网格内的节点,并且移动它们以便更好地近似裂缝和井的几何形状。描述裂缝或井的节点的移动可以响应惩罚函数,以保证移动节点不会很大地负面影响网格的质量。例如,如果网格调整模块139确定应该移动节点来改进裂缝的近似,则惩罚函数能够根据对网格质量的负面影响约束节点的移动。也可以移动其它附近节点,以便提高新调整的多面体单元周围的网格的质量,由此使形变扩大到更大区域上。网格调整模块139也可以修改网格布局,以便在保持显性离散裂缝的近似的同时提高网格质量。例如,网格调整模块139可以进行边缘翻转,以便变更最差关联多面体单元。
例子
图6A是例示描述具有两个水平层(未示出)和裂缝表示201的有裂缝地下储层的区域200的例子。区域200包括40个用2D多边形表示的离散裂缝203。在本例中,裂缝表示201内的裂缝被表示成在水平层内延伸的平坦长方形。穿过区域200的两个水平层的裂缝从区域200的底部延伸到区域200的顶部。裂缝205代表由顶部水平层界定的裂缝,裂缝207代表由底部水平层界定的裂缝,裂缝209穿过两个水平层,以便从区域200的底部延伸到区域200的顶部。区域200是可以在方法10的步骤11和方法20的步骤21中提供的区域的例子。
图6B示出了按照本发明的一些方面,被选择来近似区域200的裂缝表示201的多面体单元211的界面。具体地,与在方法10的步骤13和方法20的步骤23中进行的操作类似,利用四面体单元分解区域200。像在方法10的步骤15和方法20的步骤27中进行的那样,通过四面体单元之间的界面来近似离散裂缝。如方法20的步骤29和31所述,对网格进行网格调整,以便提高对裂缝近似的质量。四面体单元的界面或边缘描绘了可以用于构建有裂缝地下储层的模拟模型的网格的一部分。
本领域技术人员应该理解,可以生成各种分辨率的网格,以便在单元数量、网格单元质量与裂缝表示的近似之间取得适当折衷。例如,图7示出了使用各种分辨率为区域200构建网格的比较。具体地,图7A是使用约54,000个四面体单元生成的,而图7B是使用约209,000个四面体单元生成的。如方法20的步骤29和31所述,对显示在图7A和7B中的网格进行网格调整,以便提高对裂缝近似的质量。
图8通过描绘例示被选择来近似裂缝的四面体单元的网格的映射图比较了以图7A和7B的分辨率生成的网格。具体地,图8A描绘了使用约54,000个四面体单元生成网格的显示在图7A中的网格的映射图,而图8C描绘了使用约209,000个四面体单元生成网格的显示在图7B中的网格的映射图。作为强调与如何在不同分辨率上通过每个网格中的四面体单元来近似裂缝有关的差异的参考,在图8B中的映射图中示出了描绘在图6A中的裂缝表征的表示。区域213,215和217强调了图8C的较精细分辨率如何能够更好地捕捉特定的几何复杂体,诸如相互接近但没有相交的两个裂缝(213)、以小夹角相交的两个裂缝(215)、或几乎穿过另一个裂缝的平面的一个裂缝(217)。
图9描绘了根据本发明的一些方面的有裂缝地下储层的大模拟模型,其中生成网格。具体地,该模拟模型包括在物理尺度上约为17mi×1.4mi×1.1mi并包括1021个离散裂缝的区域。模拟模型的内部体积被分解成四面体单元,最精细的四面体单元具有60米分辨率。四面体单元是使用自适应网格细化技术细化的。具体地,四面体单元是通过用在其最长边的中点和与最长边相对的边之间延伸的切面来分割单元而细化的。通过四面体单元近似裂缝,以便通过四面体单元的三角形界面表示它们。具体地,通过显示在图9A中的70,000个三角形界面或边缘近似裂缝。单元的界面描绘了有裂缝地下储层的网格的一部分。在图9A中还示出了网格外部的一部分。在图9B中示出了网格的其余外部部分。整个网格被分解成约一百万个四面体单元。
图10描绘了按照本发明的一些方面的显示在图9中的大模拟模型的另一种表示,其中生成另一种网格。模拟模型300包括用标号310表示的1230个离散裂缝、和47个当前或未来计划井(8个流体注入井320和39个生产井330)。
参照图11,采用显示在图1和2中的网格化技术来创建具有与离散裂缝310和井320,330相符的自适应空间分辨率的高质量网格。为了在保持网格分辨率的同时减少网格单元的数量,如方法20的步骤33所述,合并多面体单元。网格340将模型300分解成197,802个节点和344,805个精细尺度多面体单元350,以便捕捉裂缝网络和井周围的细节。具体地,305,751个基体和39,054个裂缝控制体定义了可以用于在储层模拟器中模拟模型的893,987个流动连接。虽然图9-11执行公开在图1和2中的网格化技术来构建大模拟模型,但本领域技术人员应该理解,这些技术可以应用于包括这些较大地下储层的一个分区的各种尺度的地下储层。
虽然在前面的说明中,结合本发明的一些优选实施例对本发明作了描述,并且为了例示的目的给出了许多细节,但对于本领域技术人员来说,显而易见,可以容易地对本发明加以更改,并且可以不偏离本发明的基本原理地对本文所述的一些其它细节作较大改变。
Claims (20)
1.一种生成用于构建有裂缝地下储层的模拟模型的网格的计算机实现的方法,该方法包括:
提供有裂缝地下储层、被表示在该有裂缝地下储层内的裂缝网络、被表示在该有裂缝地下储层内的岩石基体、以及每个裂缝的几何形状的表示;
构建包括代表岩石基体的多个多面体单元的均匀、结构化的网格,所述多面体单元定义多面体单元之间的边界上的单元间界面;
细化所述均匀、结构化的网格的表示岩石基体的至少一个多面体单元,其中所细化的多面体单元在裂缝网络附近,并且其中细化发生在构建之后;
对于每个裂缝来选择与裂缝相邻的在表示岩石基体的多面体单元之间的边界上的多个单元间界面,以形成近似裂缝的几何形状的连续界面路径,其中选择发生在细化之后;
将所形成的连续界面路径中的每一个与预定裂缝性质相关联,其中关联发生在选择之后;以及
输出包括与预定裂缝性质相关联的每个连续界面路径的网格,以便用于构建有裂缝地下储层的模拟模型。
2.如权利要求1所述的方法,进一步包括:
提供一个井和该井的几何形状;
选择与该井相邻的多个单元间界面,以形成近似该井的几何形状的界面井区域;以及
将该界面井区域与预定井性质相关联,使得输出的网格包括与预定井性质相关联的界面井区域。
3.如权利要求1所述的方法,其中:
所述多面体单元是四面体,每个四面体具有最长边和与最长边相对的边;以及
所述细化包括用在最长边的中点和与最长边相对的边之间延伸的切面来分割四面体。
4.如权利要求1所述的方法,其中:
所述细化包括使用自适应网格细化技术。
5.如权利要求1所述的方法,其中:
细化使得迫近函数确定每个多面体单元与从由裂缝、注入井和生产井构成的组中选择的和每个多面体单元相邻的特征之一之间的距离。
6.如权利要求1所述的方法,进一步包括:
输出网格之前合并两个或更多个多面体单元。
7.如权利要求1所述的方法,进一步包括:
识别节点的位置,所述节点通过由两个或更多个多面体单元共享的顶点来表征;以及
输出网格之前优化节点的位置。
8.如权利要求1所述的方法,其中:
所述均匀、结构化的网格的所述多面体单元是四面体。
9.如权利要求1所述的方法,其中:
通过输出的网格上执行控制方程的离散化来构建模拟模型;以及
利用储层模拟器模拟所述模拟模型,以便直观显示有裂缝地下储层内的流体流动。
10.如权利要求1所述的方法,其中:
输出的网格定义网格内的裂缝-裂缝、裂缝-基体、和基体-基体连接。
11.一种生成用于构建有裂缝地下储层的模拟模型的网格的***,该***包括:
用户接口,用于提供有裂缝地下储层、被表示在有裂缝地下储层内的裂缝网络、被表示在该有裂缝地下储层内的岩石基体、以及每个裂缝的几何形状的表示;
分解模块,用于构建包括多个代表岩石基体的多面体单元的均匀、结构化的网格,所述多面体单元定义多面体单元之间的边界上的单元间界面;
细化模块,用于细化所述均匀、结构化的网格的表示岩石基体的至少一个多面体单元,其中所细化的多面体单元在裂缝网络附近,并且其中细化发生在构建之后;
界面选择模块,用于对于每个裂缝来选择与裂缝相邻的在表示岩石基体的多面体单元之间的边界上的多个单元间界面,以形成近似裂缝的几何形状的连续界面路径,其中选择发生在细化之后;
关联模块,用于将所形成的连续界面路径中的每一个与预定裂缝性质相关联,其中关联发生在选择之后;以及
输出模块,用于输出包括与预定裂缝性质相关联的每个连续界面路径的网格,以便构建有裂缝地下储层的模拟模型。
12.如权利要求11所述的***,进一步包括:
用于提供一个井和该井的几何形状的模块;
用于选择与该井相邻的多个单元间界面,以形成近似该井的几何形状的界面井区域的模块;以及
用于将该界面井区域与预定井性质相关联,使得由输出模块输出的网格包括与预定井性质相关联的界面井区域的模块。
13.如权利要求11所述的***,其中:
所述多面体单元是四面体,每个四面体具有最长边和与最长边相对的边;以及
由细化模块进行的所述细化包括用在最长边的中点和与最长边相对的边之间延伸的切面来分割四面体。
14.如权利要求11所述的***,其中:
由细化模块进行的所述细化包括使用自适应网格细化技术。
15.如权利要求11所述的***,其中:
在由细化模块进行的细化使得迫近函数确定每个多面体单元与从由裂缝、注入井和生产井构成的组中选择的和每个多面体单元相邻的特征之一之间的距离。
16.如权利要求11所述的***,进一步包括:
用于在由输出模块输出网格之前合并两个或更多个多面体单元的模块。
17.如权利要求11所述的***,进一步包括:
用于识别节点的位置的模块,所述节点通过由两个或更多个多面体单元共享的顶点来表征;以及
用于在输出模块输出网格之前优化节点的位置的模块。
18.如权利要求11所述的***,其中:
所述均匀、结构化的网格的所述多面体单元是四面体。
19.如权利要求11所述的***,其中:
通过在由输出模块输出的网格上执行控制方程的离散化来构建模拟模型;以及
利用储层模拟器模拟所述模拟模型,以便直观显示有裂缝地下储层内的流体流动。
20.如权利要求11所述的***,其中:
由输出模块输出的网格定义网格内的裂缝-裂缝、裂缝-基体、和基体-基体连接。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11960408P | 2008-12-03 | 2008-12-03 | |
US61/119,604 | 2008-12-03 | ||
PCT/US2009/066632 WO2010065769A2 (en) | 2008-12-03 | 2009-12-03 | System and method of grid generation for discrete fracture modeling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102239507A CN102239507A (zh) | 2011-11-09 |
CN102239507B true CN102239507B (zh) | 2016-02-24 |
Family
ID=42223611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980148461.0A Expired - Fee Related CN102239507B (zh) | 2008-12-03 | 2009-12-03 | 用于离散裂隙建模的网格生成***和方法 |
Country Status (9)
Country | Link |
---|---|
US (1) | US9026416B2 (zh) |
EP (1) | EP2370957B1 (zh) |
CN (1) | CN102239507B (zh) |
AU (1) | AU2009322303A1 (zh) |
BR (1) | BRPI0922214A2 (zh) |
CA (1) | CA2745316A1 (zh) |
EA (1) | EA201170749A1 (zh) |
MX (1) | MX2011005663A (zh) |
WO (1) | WO2010065769A2 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526104A (zh) * | 2017-08-28 | 2017-12-29 | 电子科技大学 | 基于多机多核的裂缝介质地震波场数值模拟方法 |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2477901A (en) * | 2008-12-03 | 2011-08-17 | Chevron Usa Inc | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs |
US8965745B2 (en) | 2011-04-14 | 2015-02-24 | Schlumberger Technology Corporation | Grid from depositional space |
EP2712440A4 (en) * | 2011-05-17 | 2016-05-25 | Exxonmobil Upstream Res Co | METHOD FOR DIVIDING PARALLEL-BASED SIMULATIONS OF VESSELS IN THE PRESENCE OF WELLS |
US20130151215A1 (en) * | 2011-12-12 | 2013-06-13 | Schlumberger Technology Corporation | Relaxed constraint delaunay method for discretizing fractured media |
EP2761529B1 (en) * | 2011-12-16 | 2023-01-18 | Landmark Graphics Corporation | System and method for flexible and efficient simulation of varying fracture density in a reservoir simulator |
US20140136171A1 (en) * | 2012-11-13 | 2014-05-15 | Chevron U.S.A. Inc. | Unstructured Grids For Modeling Reservoirs |
GB2512372B (en) * | 2013-03-28 | 2020-07-29 | Total Sa | Method of modelling a subsurface volume |
CN103279991A (zh) * | 2013-04-16 | 2013-09-04 | 西南石油大学 | 一种利用离散裂缝端点变形数模提高油藏开发效果的方法 |
CN103399970B (zh) * | 2013-04-16 | 2016-02-03 | 西南石油大学 | 用离散裂缝线处理进行数模测算油藏流动状况的方法 |
MX2016000151A (es) * | 2013-07-02 | 2016-03-01 | Landmark Graphics Corp | Mallado de estadios de 2 , 5d. |
CN103400413B (zh) * | 2013-07-19 | 2016-02-10 | 河海大学 | 组独有边界单元面的无厚度3节点衬砌单元智能构建方法 |
MX2016000415A (es) * | 2013-08-16 | 2016-08-03 | Landmark Graphics Corp | Creacion de una grilla de espesor para determinar estimaciones de reserva de un yacimiento. |
US10208577B2 (en) | 2013-10-09 | 2019-02-19 | Chevron U.S.A. Inc. | Method for efficient dynamic gridding |
CN104656128B (zh) * | 2013-11-19 | 2018-11-30 | 中国石油天然气集团公司 | 一种应用于石油勘探中地震仪器排列的监控*** |
WO2015102516A1 (en) * | 2013-12-31 | 2015-07-09 | Schlumberger Canada Limited | System and methodology for evaluating fracture networks |
US10808501B2 (en) * | 2014-03-17 | 2020-10-20 | Saudi Arabian Oil Company | Modeling intersecting faults and complex wellbores in reservoir simulation |
GB2545608B (en) * | 2014-11-12 | 2020-06-17 | Halliburton Energy Services Inc | Reservoir mesh creation using extended anisotropic, geometry-adaptive refinement of polyhedra |
US10344204B2 (en) | 2015-04-09 | 2019-07-09 | Diversion Technologies, LLC | Gas diverter for well and reservoir stimulation |
US10012064B2 (en) | 2015-04-09 | 2018-07-03 | Highlands Natural Resources, Plc | Gas diverter for well and reservoir stimulation |
CN104899358B (zh) * | 2015-05-13 | 2016-01-06 | 山东科技大学 | 奥灰岩溶裂隙水网络横向分布的预测方法 |
US10846445B2 (en) * | 2015-06-18 | 2020-11-24 | M-I L.L.C. | Method of waste distribution and formation pressure profile determination |
CN106326517B (zh) * | 2015-07-03 | 2019-07-30 | 中国石油化工股份有限公司 | 层状裂缝-基质杂交网格建模方法和装置 |
CA2994093A1 (en) * | 2015-09-08 | 2017-03-16 | Halliburton Energy Services, Inc. | Domain-adaptive hydraulic fracture simulators and methods |
AU2015409650A1 (en) * | 2015-09-24 | 2018-03-08 | Halliburton Energy Services Inc. | Simulating fractured reservoirs using multiple meshes |
EP3374596B1 (en) * | 2015-11-10 | 2023-06-07 | Landmark Graphics Corporation | Fracture network triangle mesh adjustment |
GB2562916A (en) * | 2016-03-01 | 2018-11-28 | Halliburton Energy Services Inc | Method for flexible structured gridding using nested locally refined grids |
US10982520B2 (en) | 2016-04-27 | 2021-04-20 | Highland Natural Resources, PLC | Gas diverter for well and reservoir stimulation |
EP3246831A1 (en) * | 2016-05-20 | 2017-11-22 | Dassault Systemes Simulia Corp. | Scalable finite element simulation of additive manufacturing |
WO2018022114A1 (en) * | 2016-07-29 | 2018-02-01 | Halliburton Energy Services, Inc. | Time-dependent spatial distribution of multiple proppant types or sizes in a fracture network |
WO2018022115A1 (en) * | 2016-07-29 | 2018-02-01 | Halliburton Energy Services, Inc. | Time-dependent spatial distribution of proppant effects in a discrete fracture network |
CN106407503B (zh) * | 2016-08-22 | 2019-10-29 | 南京特雷西能源科技有限公司 | 储层裂缝预测方法及装置 |
US11530600B2 (en) * | 2017-05-03 | 2022-12-20 | Schlumberger Technology Corporation | Fractured reservoir simulation |
CN109507723A (zh) * | 2017-09-15 | 2019-03-22 | 中国石油化工股份有限公司 | 微地震压裂裂缝模型的压裂体积计算方法及*** |
US11343775B1 (en) * | 2018-11-29 | 2022-05-24 | United Services Automobile Association (Usaa) | Method and system for managing power consumption of a mobile device in an emergency |
US10983233B2 (en) | 2019-03-12 | 2021-04-20 | Saudi Arabian Oil Company | Method for dynamic calibration and simultaneous closed-loop inversion of simulation models of fractured reservoirs |
CN110021220B (zh) * | 2019-04-28 | 2024-03-26 | 中国华能集团清洁能源技术研究院有限公司 | 一种积木式地热尾水回灌分析***及使用方法 |
CN110259443A (zh) * | 2019-07-12 | 2019-09-20 | 西南石油大学 | 一种基于3dec离散元的煤层井壁稳定性预测方法 |
US20210182460A1 (en) * | 2019-12-11 | 2021-06-17 | Exxonmobil Upstream Research Company | Semi-Elimination Methodology for Simulating High Flow Features in a Reservoir |
US11320555B2 (en) * | 2020-06-08 | 2022-05-03 | Sim Tech Llc | Systems and methods for calibration of indeterministic subsurface discrete fracture network models |
CN114001712B (zh) * | 2021-12-24 | 2022-03-22 | 四川公路工程咨询监理有限公司 | 危岩体变形监测装置及预警*** |
Family Cites Families (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5740342A (en) | 1995-04-05 | 1998-04-14 | Western Atlas International, Inc. | Method for generating a three-dimensional, locally-unstructured hybrid grid for sloping faults |
US6106561A (en) | 1997-06-23 | 2000-08-22 | Schlumberger Technology Corporation | Simulation gridding method and apparatus including a structured areal gridder adapted for use by a reservoir simulator |
US6266062B1 (en) | 1997-10-08 | 2001-07-24 | Maria-Cecilia Rivara | Longest-edge refinement and derefinement system and method for automatic mesh generation |
GB2352036B (en) | 1998-05-04 | 2002-11-27 | Schlumberger Evaluation & Prod | Near wellbore modelling method and apparatus |
FR2787219B1 (fr) | 1998-12-11 | 2001-01-12 | Inst Francais Du Petrole | Methode pour modeliser les flux de fluides dans un milieu poreux multi-couches fissure et les interactions correlatives dans un puits de production |
US6826520B1 (en) | 1999-06-24 | 2004-11-30 | Exxonmobil Upstream Research Company | Method of upscaling permeability for unstructured grids |
FR2801710B1 (fr) * | 1999-11-29 | 2002-05-03 | Inst Francais Du Petrole | Methode pour generer un maillage hybride permettant de modeliser une formation heterogene traversee par un ou plusieurs puits |
FR2802324B1 (fr) | 1999-12-10 | 2004-07-23 | Inst Francais Du Petrole | Methode pour generer un maillage sur une formation heterogene traversee par une ou plusieurs discontinuites geometriques dans le but de realiser des simulations |
US6490526B2 (en) | 2000-03-20 | 2002-12-03 | Exxonmobil Upstream Research Company | Method for characterization of multi-scale geometric attributes |
FR2809494B1 (fr) | 2000-05-26 | 2002-07-12 | Inst Francais Du Petrole | Methode pour modeliser des ecoulements dans un milieu fracture traverse par de grandes fractures |
AU2001270265A1 (en) | 2000-06-29 | 2002-01-14 | Object Reservoir, Inc. | Method and system for coordinate transformation to model radial flow near a singularity |
US7369973B2 (en) | 2000-06-29 | 2008-05-06 | Object Reservoir, Inc. | Method and system for representing reservoir systems |
GB2387000B (en) * | 2002-03-20 | 2005-06-01 | Inst Francais Du Petrole | Method for modelling fluid flows in a multilayer porous medium crossed by an unevenly distributed fracture network |
JP2004309998A (ja) * | 2003-02-18 | 2004-11-04 | Nec Corp | 確率分布推定装置および異常行動検出装置,ならびにその確率分布推定方法および異常行動検出方法 |
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 |
FR2891383B1 (fr) * | 2005-09-26 | 2008-07-11 | Inst Francais Du Petrole | Methode pour simuler des ecoulements de fluides au sein d'un milieu discretise par un maillage hybride |
US7925482B2 (en) | 2006-10-13 | 2011-04-12 | Object Reservoir, Inc. | Method and system for modeling and predicting hydraulic fracture performance in hydrocarbon reservoirs |
US7565278B2 (en) * | 2006-12-04 | 2009-07-21 | Chevron U.S.A. Inc. | Method, system and apparatus for simulating fluid flow in a fractured reservoir utilizing a combination of discrete fracture networks and homogenization of small fractures |
US20090150124A1 (en) | 2007-12-07 | 2009-06-11 | Schlumberger Technology Corporation | Model based workflow for interpreting deep-reading electromagnetic data |
CA2720638C (en) * | 2008-06-16 | 2016-01-26 | Schlumberger Canada Limited | Streamline flow simulation of a model that provides a representation of fracture corridors |
GB2477901A (en) | 2008-12-03 | 2011-08-17 | Chevron Usa Inc | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs |
US8274859B2 (en) | 2010-02-22 | 2012-09-25 | Landmark Graphics Corporation | Systems and methods for modeling 3D geological structures |
CN102750725B (zh) | 2012-06-05 | 2016-04-06 | 南京安讯网络服务有限公司 | 基于图形处理器加速的地形渲染方法 |
-
2009
- 2009-12-03 BR BRPI0922214A patent/BRPI0922214A2/pt not_active IP Right Cessation
- 2009-12-03 CN CN200980148461.0A patent/CN102239507B/zh not_active Expired - Fee Related
- 2009-12-03 CA CA2745316A patent/CA2745316A1/en not_active Abandoned
- 2009-12-03 US US12/630,734 patent/US9026416B2/en active Active
- 2009-12-03 EA EA201170749A patent/EA201170749A1/ru unknown
- 2009-12-03 AU AU2009322303A patent/AU2009322303A1/en not_active Abandoned
- 2009-12-03 WO PCT/US2009/066632 patent/WO2010065769A2/en active Application Filing
- 2009-12-03 EP EP09831131.9A patent/EP2370957B1/en active Active
- 2009-12-03 MX MX2011005663A patent/MX2011005663A/es active IP Right Grant
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526104A (zh) * | 2017-08-28 | 2017-12-29 | 电子科技大学 | 基于多机多核的裂缝介质地震波场数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
US9026416B2 (en) | 2015-05-05 |
US20100138202A1 (en) | 2010-06-03 |
EA201170749A1 (ru) | 2011-12-30 |
WO2010065769A2 (en) | 2010-06-10 |
EP2370957A2 (en) | 2011-10-05 |
WO2010065769A3 (en) | 2010-09-10 |
EP2370957A4 (en) | 2017-05-31 |
MX2011005663A (es) | 2011-06-16 |
BRPI0922214A2 (pt) | 2015-12-29 |
EP2370957B1 (en) | 2018-11-14 |
AU2009322303A1 (en) | 2011-06-23 |
CA2745316A1 (en) | 2010-06-10 |
CN102239507A (zh) | 2011-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102239507B (zh) | 用于离散裂隙建模的网格生成***和方法 | |
US9068448B2 (en) | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs | |
EP3497302B1 (en) | Machine learning training set generation | |
CN102362262B (zh) | 用于表征地下储油层中的裂面的***和方法 | |
US8818780B2 (en) | Forming a model of a subsurface region | |
US8140310B2 (en) | Reservoir fracture simulation | |
US8931580B2 (en) | Method for using dynamic target region for well path/drill center optimization | |
US20130218539A1 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
NO346116B1 (en) | Multiphase flow in a wellbore and associated hydraulic fracturing | |
Kalantari-Dahaghi et al. | Numerical simulation and multiple realizations for sensitivity study of shale gas reservoir | |
EP2631685A2 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
EP3488073A1 (en) | Modeling of oil and gas fields for appraisal and early development | |
US20230111179A1 (en) | Predicting oil and gas reservoir production | |
CA3012429C (en) | Hybrid 3d geocellular representation of selected natural fracture network subsets | |
EP3912060A1 (en) | Reservoir performance system | |
US20240037819A1 (en) | Predictive geological drawing system and method | |
CN112292714B (zh) | 基于断层辐射的网格隔分 | |
Sabatino et al. | Structural uncertainty integration within reservoir risk analysis and history matching | |
US20230152485A1 (en) | Method for modeling the damage zone of faults in fractured reservoirs | |
Baig et al. | Digitally Augmented Subsurface History Match (DASH-Sim): A New Frontier in Reservoir Simulation History Matching |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160224 Termination date: 20161203 |
|
CF01 | Termination of patent right due to non-payment of annual fee |