CN101903803B - 在非结构化栅格上对地下过程进行建模 - Google Patents
在非结构化栅格上对地下过程进行建模 Download PDFInfo
- Publication number
- CN101903803B CN101903803B CN2008801211426A CN200880121142A CN101903803B CN 101903803 B CN101903803 B CN 101903803B CN 2008801211426 A CN2008801211426 A CN 2008801211426A CN 200880121142 A CN200880121142 A CN 200880121142A CN 101903803 B CN101903803 B CN 101903803B
- Authority
- CN
- China
- Prior art keywords
- unit
- grid
- basin
- physical process
- subsurface geology
- 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 description 151
- 230000008569 process Effects 0.000 title claims description 72
- 239000011159 matrix material Substances 0.000 claims abstract description 49
- 238000009792 diffusion process Methods 0.000 claims abstract description 13
- 230000004907 flux Effects 0.000 claims description 36
- 239000000463 material Substances 0.000 claims description 19
- 239000004215 Carbon black (E152) Substances 0.000 claims description 14
- 229930195733 hydrocarbon Natural products 0.000 claims description 14
- 150000002430 hydrocarbons Chemical class 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000013459 approach Methods 0.000 claims description 8
- 238000005520 cutting process Methods 0.000 claims description 6
- 238000000638 solvent extraction Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 10
- 230000006870 function Effects 0.000 description 31
- 230000014509 gene expression Effects 0.000 description 13
- 238000013517 stratification Methods 0.000 description 11
- 238000003325 tomography Methods 0.000 description 10
- 239000011435 rock Substances 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000000470 constituent Substances 0.000 description 4
- 238000009472 formulation Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000009396 hybridization Methods 0.000 description 2
- 238000007639 printing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 102100035971 Molybdopterin molybdenumtransferase Human genes 0.000 description 1
- 101710119577 Molybdopterin molybdenumtransferase Proteins 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000005641 tunneling Effects 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
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- 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)
- Image Generation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明的实施例包含形成棱柱体栅格以及使用该棱柱体栅格和混合有限元分析来求解对流-扩散问题。可以通过在模型的平面上提供三角形网格来形成所述棱柱体栅格。接着对所述网格进行粗化以使较不满意的网格变大。接着将所述粗化的网格进行投影以形成所述棱柱体栅格。接着所述栅格的每一单元被分配多个自由度。所述栅格的混合有限元分析产生矩阵,接着求解该矩阵以得到所述对流-扩散问题的解。
Description
相关申请的交叉引用
该申请要求2007年12月14日提交的题为“MODELING-SUBSURFACEPROCESSES ON UNSTRUCTURED GRID”的美国临时专利申请61/007,761的利益,通过引用将其整体合并到本文。
技术领域
该申请一般地涉及计算机建模,具体地说,涉及用非结构化栅格对地下过程进行建模。
背景技术
在地质勘探中,期望获得关于地球表面下存在的各种岩层和结构的信息。这类信息可以包括测定的地质层、密度、孔隙度、组成成分等。接着,这个信息用来对地下盆地进行建模,即使用获得的数据预测碳氢化合物储藏的位置并且帮助碳氢化合物的提取。
非结构化栅格对于对复杂的地质结构(例如地下盆地)中的物理过程进行建模具有许多有吸引力的特性。这种栅格也可以用于其他工业中,例如航空工业或汽车工业。所关心的盆地或域可以被建模或表示为堆叠在一起的不同厚度的层的集合。地质层可沿着垂直表面或倾斜表面断裂并且退化,生成所谓的尖灭(pinch-out)。尖灭被定义为具有接近零厚度的地质层的部分。栅格应该考虑这种复杂性以产生地质层的良好模型。非结构化栅格提供比结构化栅格更好的模型。非结构化栅格可以包括由它们的顶点定义的多面体元素或单元的集合并且具有完全任意的拓扑图。例如,栅格的顶点可以属于许多个单元并且每一单元可以具有任何数目的边或面。
可以通过对流-扩散类型的数学方程来描述许多物理地下过程。这类过程的示例可以是多孔介质中的流体流、温度分布和/或压力分布。用于石油勘探的重要过程是温度分布或热建模。热建模包含从地壳下的岩浆并且穿过沉积层和源岩运动的热量。源岩是包含在石油和其他碳氢化合物材料的地层中的岩石。石油和/或其他碳氢化合物材料从源岩中被排出并且迁移到其他地方。碳氢化合物的质量由源岩和它们周围区域上承受的温度和压力条件确定。质量也受源岩和其当前位置之间的迁移路径的温度和压力条件影响。因此,盆地在其整个历史中的压力和温度条件是重要的。
为了更精确地对过程进行建模,重要的是不仅对主变量(例如压力或温度)进行建模,而且对在任何给定的表面之上的主变量的通量或能量、流体的流速等进行建模。有多种已知的方法用于对这些过程进行建模,例如有限差分、有限体积或有限元方法。在这些方法中,在考虑物理过程的地方,域被栅格覆盖。接着,通过以下方法来在栅格上逼近域:在栅格单元的指定位置处引入称为自由度的一组未知量并且为每一位置导出代数方程,所述方程将在那个位置的自由度与其他自由度连接起来。对于上面所提到的不同方法,导出这类方程的方式和自由度的位置是不同的,但是所有这些方法具有共同的特征,即,它们仅包含主变量,例如温度或压力。
为了计算通量,感兴趣的研究人员将首先使用上面所描述方法中的一个来计算所期望的主变量,接着使用数值微分来计算主变量的通量。在规则栅格(例如矩形或平行六面体栅格)上精确的所有现有数值微分方法在非结构化栅格上是不精确的并且在计算上非常耗时,尤其是如果考虑物理过程的域是高度不均匀的或异质的。此外,使用有限差分方法求解对流-扩散问题的方法需要笛卡尔(Cartesian)栅格,因此不可应用于不得不使用非结构化栅格的许多地下应用。能够对复杂几何形状进行建模的有限元方法不具有局部守恒属性并且不能够被应用在许多地下过程中。相反地,有限体积方法是局部守恒的并且能够被应用在局部正交的非结构化栅格的子集上。然而,当非结构化栅格不具有局部正交属性时,有限体积方法提供不精确的解。因此,从上面所提到的所有三种方法中,没有一个可应用于描述以非结构化栅格建模的盆地中的地下对流-扩散过程。
有另一数学方法来同时逼近主未知量和它们的通量,叫做混合有限元方法,其在以下文献中描述:F.Brezzi和M.Fortin的“Mixed and hybrid finite elementmethods”,Springer Verlag,柏林,1991年。这种方法被证明在存在不均匀介质的情况下是局部质量守恒的、精确的,并且提供对主未知量和通量两者的精确逼近。直到最近,混合有限元方法才能够被直接应用到由非结构化多面体栅格覆盖的域,所述非结构化多面体栅格对于地下应用是非常常见的。在以下文献中提出了针对在任意多面体栅格上的扩散类型方程的混合有限元方法的新变型:Yu.Kuznetsov和S.Repin的“New mixed finite element method on polygonal andpolyhedral meshes”,Russian Journal of Numerical Analysis and MathematicalModeling,第18卷,261-278页,2003年。
发明内容
本发明涉及对在经过地质时期形成的盆地中的能量传递和/或压力分布提供精确模型的系统和方法。在任何给定的时间,盆地被表示为堆叠在一起的不同厚度的层的集合。在盆地中的某些位置,层的厚度退化到零,形成尖灭。本发明的实施例使用棱柱体网格和混合有限元分析来对盆地中的各种过程进行建模,这些过程包括能量传递(例如热能量)和压力。因此,本发明的实施例求解主未知量(例如温度或压力)和从未知量(例如温度通量或压力通量)两者。下述方面中的一个或更多个可以用来提供经过地质时期形成的盆地中的能量传递和/或压力分布(例如物理过程)的精确模型。该模型可以用来解释现代储油层,并且继而依赖于该模型,从而基于模型的模拟结果来控制碳氢化合物开采活动。基于从由下述方面中的一个或多个产生的模拟盆地模型中解释的结果,可以控制碳氢化合物的开采,例如,可以控制来自地面设施的开采速率,可以有策略地打井,和/或一般性地表征储油层。
在一个一般方面,一种用于在计算机上对物理区域进行建模的方法,其中所述物理区域包括多个地层,所述方法包括:接收定义所述物理区域的至少一个物理特性的数据;在所述物理区域的模型的平面上提供三角形网格,其中所述网格包括多个单元;以非均匀方式对所述三角形网格进行粗化,使得较不满意的单元较大;以及将粗化的三角形单元在与所述物理区域中的所述平面正交的方向上进行投影以形成棱柱体栅格,其中根据所述地层,粗化的三角形网格的每一单元被分割成子单元。
这个方面的实现可以包括下述特征中的一个或更多个。例如,粗化可以包括将所述数据投影在平面上以及使用所述数据来确定哪些单元是较不满意的。所述模型可以包括对所述物理区域中的物理特征建模的建模特征,以及其中使用可以包括将优先值分配给每一单元,其中基于每一单元是否接近建模特征和该建模特征的类型来确定该值。提供三角形网格可以包括在所述平面上提供矩形网格;以及沿着至少一个对角线分割所述矩形网格的每一单元。粗化可以包括通过消除两个相邻三角形的公共边来合并所述两个相邻三角形。所述棱柱体栅格可以包括多个棱柱体单元、多个棱锥体单元和多个四面体单元。所述方法可以用来对所述物理区域中的物理过程的至少一个通量进行建模,所述方法还包括为每一子单元中的通量分配多个自由度;将混合有限元分析应用到每一子单元以产生矩阵;以及求解该矩阵以确定所述区域中的物理过程的通量。
分配可以包括对于每一单元,为所述物理过程分配一个自由度;以及为所述单元的每一面分配另一自由度。应用可以包括使用恒定划分法来形成有限元空间。所述物理过程可以是对流-扩散过程。所述物理过程可以是温度和压力中的一个并且所述物理区域是地下地质盆地。所述物理过程可以包含碳氢化合物材料的形成。所述物理过程可以包含碳氢化合物材料的运动。可以从来自传感器的信息导出所述数据,所述传感器测量所述物理区域的至少一个物理特性。
在另一一般方面,一种用于在计算机上对物理过程和该物理过程的通量进行建模的方法,其包括形成对物理区域建模的非结构化棱柱体栅格,其中所述物理过程在所述物理区域内发生并且所述棱柱体栅格包括多个单元;对于每一单元,为所述物理过程和所述通量分配多个自由度;将混合有限元分析应用到每一单元以产生矩阵;以及求解该矩阵以确定所述区域中的所述物理过程和所述通量。
这个方面的实现可以包括下述特征中的一个或更多个。例如,形成可包括在所述物理区域的模型的平面上提供三角形网格,其中所述网格包括多个单元;以非均匀方式对所述三角形网格进行粗化,使得较不满意的单元较大;以及将粗化的三角形单元在与所述物理区域中的所述平面正交的方向上进行投影以形成所述棱柱体栅格。所述棱柱体栅格可以包括多个棱柱体单元、多个棱锥体单元以及多个四面体单元。分配可以包括对于每一单元,为所述物理过程分配一个自由度;以及对于每一单元,为所述单元的每一面分配另一自由度。应用可以包括使用恒定划分法来形成有限元空间。所确定的物理过程和通量可以用来影响所述物理区域中的改变。所述物理过程可以是温度和压力中的一个并且所述物理区域是地下地质盆地。
在另一一般方面,一种具有计算机可读介质的计算机程序产品,所述计算机可读介质具有在其上记载的计算机程序逻辑,所述计算机程序逻辑用于对物理区域中的物理过程和该物理过程的通量进行建模,所述计算机程序产品包括用于形成对所述物理区域进行建模的非结构化棱柱体栅格的代码;用于将混合有限元分析应用到所述棱柱体栅格以产生矩阵的代码;以及用于求解该矩阵从而确定所述区域中的物理过程和通量的代码。
这个方面的实现可以包括下述特征中的一个或更多个。例如,用于形成的代码可以包括这样的代码:用于在所述物理区域的模型的平面上提供三角形网格的代码,其中所述网格包括多个单元;以非均匀方式对所述三角形网格进行粗化,使得较不满意的单元是较大的;以及将粗化的三角形网格在与所述物理区域中的所述平面正交的方向上进行投影以形成所述棱柱体栅格。所述棱柱体栅格可以包括多个单元,以及所述用于应用的代码可以包括对于每一单元,为所述物理过程分配一个自由度;对于每一单元,为所述单元的每一面分配另一自由度;以及使用恒定划分法来形成有限元空间。所述用于求解的代码可以包括使用预调的共轭梯度分析来求解所述矩阵。
本发明的实施例通过将某些或大多数地质和几何特征,例如尖灭边界投影到水平平面中来工作。注意,投影可以是非正交的或倾斜的。接着,本发明的实施例在那个平面上生成分解所有期望的特征的非结构化栅格。注意到,栅格可以由多边形、四边形、三角形或其组合组成。接着,本发明的实施例将获得的栅格投影回到所有层的所有边界表面上,从而构造棱柱体栅格。所述棱柱体栅格可以包括多个单元,这些单元可以是棱柱体、四面体形状、棱锥体或其组合。注意,非结构化棱柱体栅格逼近所有层的边界表面。
接着本发明的实施例可以通过以下步骤来工作:对于主未知量,在单元中心处为每一单元关联一个自由度,并且对于通量的法向分量,在面中心处为所述单元的每一面关联一个自由度。接着本发明的实施例使用混合有限元方法,例如Yu.Kuznetsov和S.Repin的方法来使问题离散化。空间离散化产生稀疏矩阵方程。接着本发明的实施例可以求解该矩阵方程,以得到主未知量和在所述单元的各面处的通量的法向分量两者。因此,本发明的实施例在没有极大地扩展需要求解的未知量的数目的情况下,提供更精确的建模。
本发明的实施例可以提供覆盖所述物理区域的水平平面的三角形网格来形成所述棱柱体栅格。接着可以以非均匀方式对所述网格进行粗化,使得较不满意的单元是较大的,而使满意的单元保留更细的形式。接着将粗化的网格在所述物理区域中的垂直方向上进行投影以形成所述棱柱体网格。
为了可以更好地理解后面的本发明的详细描述,前面已经相当广泛地概述了本发明的特征和技术优势。在下文中将描述本发明的另外特征和优势,其构成本发明的权利要求的主题。本领域技术人员应意识到,公开的概念和具体的实施例可以容易地用作修改或设计其他结构的基础,其他结构用于实现本发明的相同目的。本领域技术人员也应该认识到,这类等效的构造不偏离随附的权利要求中陈述的本发明的精神和范围。当结合附图考虑时,从下述描述中将更好地理解被认为是本发明的特性的新颖特征(关于其组织架构和操作方法两者)以及其他目标和优势。然而,应清楚地明白,提供每一附图是仅为了解释和描述的目的,无意作为本发明的限制的定义。
附图说明
为了更彻底地理解本发明,现在结合附图参考下述描述,在附图中:
图1描绘根据本发明的实施例的被分割成层的域。
图2A和2B描绘根据本发明的实施例的域和用矩形网格覆盖的域。
图3描绘根据本发明的实施例的非均匀粗化的三角形栅格。
图4描绘根据本发明的实施例形成的不同类型的单元。
图5描绘根据本发明的实施例形成的3D棱柱体栅格。
图6描绘本发明的实施例使用的四面体单元。
图7描绘本发明的实施例使用的棱锥体单元。
图8A-8D根据本发明的实施例描绘本发明的实施例使用的、被分成三个四面体的棱柱体单元。
图9描绘根据本发明的实施例的分割相邻单元的独立面;
图10描绘根据本发明的实施例的形成棱柱体栅格的方法;
图11描绘根据本发明的实施例的求解矩阵的方法;以及
图12描绘适合使用本发明的计算机系统的框图。
具体实施方式
本发明的实施例对于对地下油田进行建模是有用的。在此所描述的实施例的示例可以参考这类油田。然而,实施例可以用来对包含其他材料和/或过程的的其他域进行建模。例如,可以包含其他碳氢化合物材料,例如煤炭。本发明的实施例对于采矿或掘坑道是有用的。本发明的实施例可以用于其他域类型,例如大气,并且对于对天气、温度和/或污染进行建模是有用的。另一域可以是海洋,并且实施例可以用来测量声音、温度、含盐度和/或污染度。可以使用本发明的实施例来对任何类型的分层域进行建模。可以使用本发明的实施例来对通过对流-扩散过程而运动的任何类型的材料进行建模。可以使用本发明的实施例来对域或材料中存在的任何类型的通量进行建模。
如前面所述的,可以将本发明的实施例应用到任何对流-扩散过程。下面是3D对流-扩散类型方程的示例:
用边界条件补充方程(1.1)
p=gD 在ΓD上
这里n是ΓN的外向单位法向向量,σ是非负函数,并且gD和gN是给定函数。假设方程(1.1)-(1.2)具有唯一解。
可以由等价一阶系统代替偏微分方程(1.1)-(1.2):
p=gD 在ΓD上
-u·n+σ·p=gN 在ΓN上 (1.4)
方程(1.3)-(1.4)是方程(1.1)-(1.2)的混合公式表示。注意,以此方式,可以同时逼近主未知量p及其通量u。
Ω={(x,y,z)∈R3:(x,y)∈G,Zmin(x,y)≤z≤Zmax(x,y)}
这里Zmin(x,y)和Zmax(x,y)是光滑表面。
设m是正整数并且z=Zi(x,y),i=0,...,m,是定义在G上的单值连续函数,使得
Zi-1(x,y)≤Zi(x,y) 在中
这些函数定义地质层之间的界面。换句话说,计算域Ω可以被分成m个子域(带或层),其被定义如下:对于所有i=1,...,m,
Ωi={(x,y,z)∈Ω:(x,y)∈G,Zi-1(x,y)≤z≤Zi(x,y)}
图1描绘了将计算域Ω100分割成多个子域或层101-107的示例。注意,图1以分解图描绘了不同的层,然而对于计算,不需要将层隔开。还注意到,假设子域Ωi满足锥条件,这里子域的边界不具有奇点(零角度等),并且此外,所有集合
图1描绘了盆地的不同地层。可以通过各种技术例如地层分析和/或地震反演,使用传感器来测量盆地的各种特性来确定用来形成图1的不同层的数据。
可以以下述方式将集合Pi的边界投影到平坦平面。对于来自的任何给定的点(x,y,z),投影的点具有坐标(x,y,0)。所有这类点组成像图2A那样的闭合线的集合,其用来生成平面三角测量,在图2B和图3中示出了这种示例。
这里
注意,λ是压力函数p=p(x)在ΓN上的约束。在这个公式中,所有边界条件是自然边界条件。
在下述分析中,考虑方程(1.5),然而也可以将结论应用到方程(1.6)而不失一般性。
本发明的实施例使用棱柱体栅格,其在对对流-扩散地下过程进行建模方面提供许多有吸引力的特性。在许多情况下,域可以被表示为堆叠在一起的不同厚度的层的集合。棱柱体栅格满意地表示了堆叠层中的分片地质结构和非结构化几何特征。通过地质统计信息的后处理来提供2D几何数据。通常,数据是与2D细矩形栅格的节点或单元关联的材料属性。注意,可能有数百万个节点。节点中存在材料数据意味着其是计算节点,而不存在材料数据则意味着其是局外节点(node-outlier)。计算节点的集合定义了计算域。
地质层之间的界面是地质统计数据,并且可能互相相交,导致拓扑上不正确的情形。底部界面是最低的地质层的底部边界,并且也由地质统计数据表示。地质层可以沿着垂直表面断裂并且退化。尖灭被定义为厚度模不大于用户定义的阈值δ≥0的地质层的部分。断层折线被定义为底部地质界面和断层的相交处。
为了简化算法的描述,在本发明的实施例中要使用的栅格应该满足某些非常自然的要求。目标棱柱体栅格应该是1D栅格和2D三角测量的逻辑积,所述1D栅格和2D三角测量的节点形成不具有局外节点的地质统计矩形栅格的子集,其中节点密度等于所关心的用户指定区域中原始栅格中的节点密度。除此之外,在用户定义的断层和井以及自动检测到的尖灭附近必须精炼(改进)三角测量。此外,2D栅格中三角形的数目不应该大于用户定义的数目。关于棱柱体栅格,棱柱体的侧面不得不逼近地质界面并且形成形状规则的2D三角测量。
下面所描述的过程是可以用来构造棱柱体栅格的过程的示例。注意,可以使用其他过程。此外,那种类型的棱柱体栅格不限制混合有限元方法。首先是为了将尖灭投影在表示井的底部地质界面、断层折线和点上,对所述2D规则三角测量的产生进行改进。接下来,将2D三角测量投影在由函数Zi(x,y),i=1,...,m所定义的表面上,以形成得到的3D棱柱体栅格。在下面的段落中进一步描述这个过程。
为了形成2D规则三角形栅格,示例性过程可以开始于矩形栅格。给定节点在x-和y-方向上的节点坐标,产生覆盖域G的矩形一致网格Gh,其由具有带材料数据的四个节点中的至少一个的单元组成。例如,图2A描绘域200,图2B描绘覆盖该域200的矩形栅格201。不失一般性,假设网格Gh由正方形组成,即在x-和y-方向上的网格尺寸相等,hx=hy=h。
接下来,根据本发明的实施例,每一矩形单元被其对角线分割成两个三角形。下面描述可以用来形成三角形的一个过程,注意,可以使用其他过程。可以根据下述规则来进行两个可能的对角线之间的选择。设每一矩形单元被分配一整数,该整数等于其节点的最小x-和y-索引的和。对于具有偶数的单元,分割对角线具有带最小x-和y-索引的单元的节点。对于具有奇数的单元,选择另一对角线。对于具有给定的材料数据分布的给定的节点集合,上述过程唯一地指定三角测量。改变对角线的方向减少了栅格朝向的问题。作为形成三角形的替换过程,可以通过使用两个对角线来将每一矩形单元分割成四个三角形。
将产生的三角测量投影在底部地质界面上,如公式(1.2)那一段中所描述的。注意,在域G中可以定义感兴趣的区域ωi,这里栅格的修改是不必要的或不期望的。
设Pi h表示Gh的矩形元素的最大子集,其属于Pi。如果没有属于Pi的Gh的元素,但是有属于Pi的Gh的顶点,则这个顶点被认为属于Pi h。那么集合被定义为
这被称为“尖灭”投影。根据该定义,“尖灭”投影是属于Gh的边和顶点的子集。
接下来,根据本发明的实施例,将不同的优先级分配给三角形。将优先级或整数标记分配给细栅格的每一三角形。优先级的值控制粗化过程。
在开始时,将零优先级分配到所有三角形。对于其闭合线与断层相交的三角形,它们的优先级被改变为1。为了找到与断层相交的三角形,可以使用下述方法。首先,从断层三角测量中提取三角形,所述断层三角测量与最底部地质界面相交。其次,检查每一提取的三角形是否与细栅格三角形相交。对于其闭合线与尖灭相交的三角形,它们的优先级被改变为2。这些三角形被定义为违反下述条件的三角形:在所有三角形节点中,地质层的厚度或大于δ或小于-δ。对于其闭合线包含井点的三角形,它们的优先级被改变为3。对于属于用户定义的感兴趣的区域ωi的三角形,它们的优先级被改变为4。对于满足若干条件的三角形,它们的优先级被改变为最大优先级值。
在已经分配优先级之后,可以对三角形栅格进行非均匀粗化。栅格的细部分可以具有大量的同样小的三角形。这些区域是更期望的,因为它们包含更多的信息,包括感兴趣的地质特征,例如井、断层、尖灭和/或用户指示为期望的特征。栅格的粗化部分与栅格的较细部分不是同样满意的。栅格可以具有粗化范围,这里最粗化的区域指示具有较不满意或不满意质量的部分,并且不具有粗化的区域指示最期望的区域。在最粗化的区域和不粗化的区域之间的粗化区域指示具有某些期望方面的区域。
粗化是一系列合并三角形的过程。例如,可以通过消除它们的公共边来将两个三角形合成一个。这个过程包括两个阶段。首先,标记特定三角形用于进行粗化。其次,对它们进行粗化。应该注意到,栅格一致性可能引起对未标记三角形的粗化。每一粗三角形继承了两个合并的三角形的最大优先级。除了优先级外,每一三角形还被分配表示级别的另一整数。初始细栅格的任何三角形具有级别1。可以将粗化应用到相同级别j的一对三角形,并且得到具有级别j+1的更粗的三角形。
下面是根据本发明的实施例的粗化过程的一个示例。注意,可以使用其他过程。
粗化过程可以被描述为下述循环:
1)设置k=1。
2)由具有零优先级的三角形形成集合M,零优先级三角形的粗化不会引起对具有非零优先级的其他三角形进行粗化。
3)如果M为空,则到6。
4)对来自M的三角形进行粗化。
5)到2。
6)如果新栅格中三角形的数目不大于用户定义的阈值Ntusr,则停止。
7)由非零优先级不大于k的三角形形成集合M,其粗化不会引起对优先级大于k的其他三角形进行粗化。
8)对来自M的三角形进行粗化。
9)如果k≤3,则设置k=k+1,否则设置k=1。
10)到1。
输出三角形栅格具有与投影的输入矩形网格的特定节点重合的节点和感兴趣区域中的细三角形,以及朝着井点、断层折线和尖灭进行改进的三角形。图3描绘了上述栅格粗化过程的示例。得到的栅格300描绘了非均匀区域,其具有最粗三角形301、粗化三角形302和细三角形303。注意,细三角形可能具有某些粗化或没有粗化。注意,在图3中,通过使用上面所描述的两对角线方法来形成三角形。
在粗化之后,可以形成3D棱柱体栅格。设eh是G的三角测量中的三角形,并且a(k)=(ax (k),ay (k)),k=1,2,3是eh的顶点。考虑在R3中的三条垂直线(x,y)=a(k),k=1,2,3,,并且用a(k,i),k=1,2,3表示它们与表面z=Zi(x,y),i=0,...,m的交界处。于是,具有位于相邻表面上的顶点的任何多面体或是垂直棱柱体(所有6个顶点是不同的)、或是棱锥体(两个顶点重合,例如对应的顶点a(k),k=1,2,3,属于)、或四面体(两对顶点重合,即集合的两个顶点a(k),k=1,2,3,属于Ph)。图4描绘了具有三种类型的棱柱体的3D棱柱形栅格的部分,所述三种类型的棱柱体即垂直棱柱体401、棱锥体402和四面体403。换句话说,棱锥体是一个边消失的棱柱体,而四面体是两个边消失的棱柱体。
对中的所有三角形执行上述操作将提供域Ω的分区Ωh。在特定的单元中,通过分段三角形表面z=Zi (h)(x,y)来逼近每一表面z=Zi(x,y),所述分段三角形表面由棱柱体的顶部(底部)三角形以及棱锥体和四面体的特定面组成。图5描绘了3D棱柱体栅格500的示例。
在完成3D棱柱体栅格之后,可以对栅格进行混合有限元分析。在前面的部分中指出,栅格Ωh包括元素{Ek},其或是垂直棱柱体,或是棱锥体,或是四面体。为了针对方程(1.5)用公式表示混合有限元(MFE)方法,应该定义空间L2(Ω)和L2(ΓN)的有限元子空间。
混合有限元方法中的一个问题是设计空间的有限元子空间Vh。为了计算效率,应该仅考虑那些有限元向量函数,它们在相邻单元Ek和El(k>l)之间的界面Γkl上以及在单元Ek与边界ΓN的交界处上具有恒定的法向分量。空间的有限元子空间的维数等于不同界面{Γkl}和的总数目。可以基于在以下文献中描述的“恒定划分(div-constant)”法来构造这个有限元空间:Yu.Kuznetsov和S.Repin的“New mixed finite element method on polygonal andpolyhedral meshes”,Russian Journal of Numerical Analysis and MathematicalModeling,第18卷,261-278页,2003年。
对于四面体单元T,有限元空间Vh|T与经典的最低阶Raviart-Thomas有限元空间RT0(T)一致(参见F.Brezzi和M.Fortin的“Mixed and hybrid finite elementmethods”,Springer Verlag,柏林,1991年)。有限元向量值函数wh∈RT0(T)具有四个自由度(DOF),即
这里φi(x)是与四面体T的面γi(i=1,2,3,4)关联的基向量函数。
用γj表示的四面体T的面与顶点Aj相对,即,γ1是面A2A3A4,γ2是面A3A4A1,γ3是面A4A1A2,γ4是面A1A2A3。设nj是在面γj上的外向单位向量,并且hj是从顶点Aj到面γj上的垂直线长度,j=1,2,3,4。在图6中示出了这种四面体单元600。
在四面体T上最低阶Raviart-Thomas元素的空间可以被定义为
RT0(T)=span{φ1,φ2,φ3,φ4}
简单明了的计算显示了基函数可以被明确地定义为
这里x(i)是顶点Ai的坐标,i=1,2,3,4。
对于棱锥体单元P,当单元P∈Ωh是四边形棱锥体时,使用“恒定划分”法来构造有限元空间Vh|P。
用γj,j=1,2,3,4,5来表示棱锥体P的面,即,γ1是面A2A3A5,γ2是面A1A3A4,γ3是面A1A2A4A5,γ4是面A1A2A3,以及γ5是面A3A4A5。在图7中示出了这种棱锥体单元700。
为了描述“恒定划分”法,将棱锥体P分成两个四面体T1和T2。可以以两种不同的方式来完成,并且两种方式都提供下面所描述的可工作的算法。因此,不失一般性,假设棱锥体P被划分成两个四面体T1=A1A2A3A4和T2=A2A3A4A5。设ni表示到面γi的单位外向法向向量,i=1,2,3,4,5,并且n6表示到四面体T1和T2之间的界面γ6=A2A3A4的、从T1指向T2的单位法向向量。
这个表达式显示出向量函数uh在四面体T1和T2中的每一个上是线性的,其属于空间Hdiv(P),并且满足所要求的条件
为了确定在界面γ6上通量的法向分量的未知值u6,下述条件应该是可操作的
从uh的定义(2.1)中获得T1和T2中uh的散度的表达式。将散度算子应用于(2.1)的两边,得到
以及
由于因此得到
可以以不同的方式得到在假设(2.2)下u6的值。第一,应用Stokes公式
这里|γi|是对应的面γi的面积,|P|是棱锥体的体积。第二,从(2.4)确定u6的值:
由于在棱锥体P的面上的基函数的法向分量的值是已知的,因此其法向分量u6 (i),i=1,2,3,4,5在内部面γ6上的值可以从公式(2.5)和(2.6)中得到。因此,由以下给出基函数的明确的表达式:
对于具有三角形底的棱柱体单元,设∏∈Ωh。使用“恒定划分”法来构造有限元空间Vh|∏。
用γj来表示棱柱体∏的面,即,γ1是面A2A3A5A6,γ2是面A1A3A4A6,γ3是面A1A2A4A5,γ4是底面A1A2A3,以及γ5是顶面A4A5A6。在图8A中示出了这种棱柱体单元800。
为了根据实施例应用“恒定划分”法,将棱柱体∏分成三个四面体T1、T2和T3。在图8B-8D中示出了分割的四面体801、802、803。不失一般性,假设棱柱体∏被分成四面体T1=A1A2A3A6、T2=A1A4A5A6和T3=A1A2A5A6。在图8B-8D中示出了分割的四面体801、802、803。
设ni,i=1,2,3,4,5表示到面γi的外向单位法向向量,n6表示到四面体T1和T3之间的界面γ6=A1A2A6的、从T1指向T3的单位法向向量,并且最后,n7表示到四面体T2和T3之间的界面γ7=A1A5A6的、从T2指向T3的单位法向向量。在每一四面体上构造经典的最低阶Raviart-Thomas有限元空间,确定基函数k=1,2,3的集合,并且∏中的向量场uh被定义为如下:
根据这个表达式,向量函数uh在四面体T1、T2和T3中的每一个上是线性的,属于空间Hdiv(∏),并且满足所要求的条件
对于分别在附属界面γ6和γ7上未知法向向量u6和u7的自然选择来自于以下条件:在∏上uh的散度是恒定的,即
因此
这里|γi|是对应的面γi的面积并且|∏|是棱柱体的体积。其次,u6和u7的值计算如下:
以及
由于在棱柱体∏的面上的基函数的法向分量的值是已知的,因此可以通过公式(2.10)和(2.11)来分别确定法向分量u6 (i)和u7 (i),i=1,2,3,4,5在内部面γ6和γ7上的值。因此,由下面给出基函数的明确的表达式:
强调“恒定划分”法的下述区别是合情合理的。将栅格单元分割成四面体不依赖于其相邻单元的分割,例如,两个相邻单元的公共的四边形面可以以不同的方式从不同的边进行划分。例如,图9描绘了两个相邻单元901、902,其具有以不同的方式分割的邻接的面903、904。与四面体栅格相比,这个特征简化了棱柱体栅格的产生和其上的边界值离散化问题。
对于一般的棱柱体单元,上面所描述的针对具有三角形底的棱柱体单元的方法可以被推广到覆盖一般的棱柱体∏。每一棱柱体可以被独立地分成四面体并且使用在前面的部分中描述的“恒定划分”法来构造有限元空间Vh|∏。
上面的段落描述了对于通量的法向分量,单元的每一个面在面中心处具有一个自由度。对于每一主未知量,每一单元在单元中心也被分配关联的一个自由度。例如,(多个)主未知量可以包括周围的温度和/或流体压力。
下述段落描述了将混合有限元(MFE)分析的杂化变型用在3D棱柱体栅格的单元上,其中引入了与棱柱体的面关联的额外的自由度。这些自由度允许将一个单元的通量与另一单元的通量隔开。这些自由度在数学文献中被称为拉格朗日(LaGrange)乘子。由于原始的变量(例如一个棱柱体单元的单元中心的温度(或压力)和在边界上的通量)脱离任何其他单元的温度和通量,引入额外的自由度允许简化数值问题的结构。因此,可以消除那些未知量,这允许减少未知量的数目。
有限元问题(3.1)得到线性代数方程系
其中鞍点矩阵
这里M=MT是正定矩阵,并且D=DT和∑=∑T或是正定矩阵或是半正定矩阵。可以显示出,系(3.2)具有唯一解。
很好地开发出了用于具有对称鞍点矩阵的代数系统的迭代方法。然而,起因于多面体栅格上MFE离散化的、用于鞍点矩阵的高效预调技术仍然是关注点。对于高效预调,对称正定矩阵是更好的对象。可以通过使用混合有限元问题(3.1)的杂化来将系(3.2)变换为具有对称正定矩阵的等价系。在接下来段落中,这个方法被描述为求解问题(3.1)的优选方式。
对于混合有限元分析的杂化,设Ek是Ωh中的栅格单元,并且Vh (k)和Lh (k)分别是有限元空间Vh和Lh在Ek上的约束。此外,生成新的有限元空间Λh,其是在栅格单元之间的界面Γkl上以及在栅格单元与ΓN的边界部分交界处上定义的函数λh=λh(x)的空间。在每一界面上,函数λh∈Λh等于常量。
为了引入混合杂化有限元(MHFE)问题,不得不定义两个额外的有限元空间以及许多双线性形式和线性泛函。两个新的有限元空间是
这里n是Ωh中单元Ek的数目。注意到空间Vh (k)中的任何一个的维数最多是5,并且每一Lh (k)的维数等于1。
对于元素u,v∈Vh,p,q∈Lh以及λ,μ∈Λh,引入下述双线性形式:
此外,定义下述线性泛函
a(uh,v)+b(v,ph)+c(v,λh)=lD(v)
b(uh,q)-d(ph,q)=lf(q) (3.4)
c(uh,μ) -σ(λh,μ)=lN(μ)
有限元问题(3.4)得到线性代数方程系
其中鞍点矩阵
是用于单元Ei的局部鞍点矩阵以及Ni是对应的组合矩阵(assemble matrix)。
相应地,系(3.5)的右手边可以被写为如下:
这里Γi k,k=1,...si是单元Ei的面。
这里VEi是单元Ei的体积并且对于内部单元矩阵∑i=0。
以及
由于矩阵M、B和D的结构,矩阵BTM-1B+D是对角矩阵,因此它是可逆的。使用关系(3.7)和(3.8)来将系(3.5)变换为系:
这里
S=CTM-1C-CTM-1CB(BTM-1B+D)-1BTM-1C+∑
以及
矩阵S被称为“浓缩矩阵(condensed matrix)”。这个矩阵是对称的和正定的,除了在Neumann边界条件的情况下,S是半正定的,但是其具有简单的核恒定(kernel constant)向量。这个矩阵在本质上是全局的,其将所有的节点或单元连接在一起。可以同时求解大的线性方程系。
可以将任何迭代方法应用到求解具有该矩阵的线性方程系(3.9)。一个方法是预处理的共轭梯度方法(PCG),然而可以使用其他的方法。注意,在半正定性的情况下,应该在与核正交的子空间中执行PCG。在求解系(3.9)之后,可以分别使用方程(3.8)和(3.7)来逐个元素地局部恢复主未知量和
注意,本发明的实施例可以以单个主未知量操作,例如温度或压力,及其相关的通量。其他实施例可以以多于一个的主未知量操作。
根据本发明的各种实施例,上面概述的各种过程和方法可以被组合在一个或更多个不同的方法中,用在一个或更多个不同的系中,用在一个或更多个不同的计算机程序产品中。
例如,一个示例性方法1000可以用于形成棱柱体栅格,如图10中所示。使用正交投影将感兴趣的地质特征和几何特征(例如尖灭边界、断层线或井的位置)投影到水平面中(1001)。产生细的符合矩形的网格,该网格覆盖与在其上提供材料数据的细栅格相同大小的投影域的所有特征(1002)。矩形栅格被分割成三角形(1003)。在栅格上的各种线和点可以表示例如断层线和井位置。以非均匀方式对三角形进行粗化(1004)。期望保持某些地质或几何特征附近的细三角测量,但是远离这些特征具有较粗分辨率将允许更容易的分析。这种栅格仅由三角形组成。将粗化的栅格垂直投影在所有层的所有边界表面上以形成棱柱体栅格(1005)。这种栅格将包含可以是三角棱柱体、四面体或棱锥体的单元。以这种方式建立的非结构化棱柱体栅格逼近所有层的边界表面。注意,在对流-扩散地下问题中,输入数据与数百万个节点关联。
另一方法1100可以用于求解针对地质盆地的对流-扩散问题,如图11中所示。形成对盆地建模的栅格(1101)。注意,可以通过图10中所示的方法1000来形成栅格,或可以使用另一方法。该方法对于主未知量在单元中心处为每一栅格单元关联一个自由度,以及对于通量的法向分量在面中心处为单元的每一面关联一个自由度(1102)。使用混合有限元法来分析具有关联的自由度的栅格(1103)。这个分析产生稀疏矩阵方程。接着该方法可以求解矩阵方程,得到(多个)主未知量和(多个)未知量的通量在单元的面处的法向分量两者(1104)。
注意到,可以以硬件、软件和/或固件和/或其任何组合来实施在此描述的功能中的任何一个。当以软件实施时,本发明的元素本质上是执行必要任务的代码段。程序或代码段可以被存储在计算机可读介质中或通过计算机数据信号发送。“计算机可读介质”可以包括能够存储或传递信息的任何介质。计算机可读介质的示例包括电子电路、半导体存储器器件、ROM、闪存、可擦除ROM(EROM)、软盘、压缩盘CD-ROM、光盘、硬盘、光纤介质等。计算机数据信号可以包括能够在传输介质(例如电子网络信道、光纤、空气、电磁、RF链接等)上传播的任何信号。可以经由计算机网络(例如因特网、内部网等)下载代码段。
图12示出了适于使用本发明的计算机系统1200。中央处理器(CPU)1201被耦合到系统总线1202。CPU 1201可以是任何通用CPU,例如英特尔奔腾(IntelPentium)处理器。然而,本发明不被CPU 1201的架构所限制,只要CPU 1201支持在此描述的发明操作。总线1202被耦合到随机存取存储器(RAM)1203,其可以是SRAM、DRAM或SDRAM。ROM 1204也被耦合到总线1202,其可以是PROM、EPROM或EEPROM。正如本领域中熟知的,RAM 1203和ROM 1204保存用户和系统的数据和程序。
总线1202也被耦合到输入/输出(I/O)控制卡1205、通信适配器卡1211、用户接口卡1208和显示卡1209。I/O适配器卡1205将存储设备1206(例如硬盘驱动器、CD驱动器、软盘驱动器、磁带驱动器中的一个或多个)连接到计算机系统。I/O适配器1205可以被连接到打印设备,其允许系统打印信息(文件、照片、文档等)的纸质副本。注意,打印设备可以是打印机(例如喷墨打印机、激光打印机等)、传真机或复印机。通信卡1211适于将计算机系统1200耦合到网络1212,其可以是电话网络、局域网(LAN)和/或广域网(WAN)、以太网和/或因特网中的一个或多个。用户接口卡1208将用户输入设备(例如键盘1213和指针设备1207)耦合到计算机系统1200。用户接口卡1208也可以将声音经由扬声器输出到用户。由CPU 1201驱动显示卡1209以控制在显示设备1210上的显示。
虽然已经详细地描述了本发明及其优势,但是应明白,在不偏离由随附的权利要求所限定的本发明的精神和范围的情况下,在此可以进行各种改变、替代和变换。此外,无意将本申请的范围限制在说明书中所描述的过程、机器、制造、物质的组成成分、装置、方法和步骤的特定实施例。正如本领域任何一个技术人员从本发明的公开中容易理解的,根据本发明可以使用当前现有的或后来开发的、执行与在此所描述的对应的实施例实质上相同的功能或实现实质上相同的结果的过程、机器、制造、物质的组成成分、装置、方法或步骤。因此,随附的权利要求意在将这种过程、机器、制造、物质的组成成分、装置、方法或步骤包括在它们的范围内。
Claims (17)
1.一种用于对地下地质盆地进行建模的方法,所述方法被用于对所述地下地质盆地中物理过程的至少一个通量建模,以便可控制所述地下地质盆地中的碳氢化合物开采,其中所述地下地质盆地包括多个地层,所述方法包括:
接收定义所述地下地质盆地的至少一个物理特性的数据;
在所述地下地质盆地的模型的平面上提供三角形网格,其中所述网格包括多个单元;
以非均匀方式对所述三角形网格进行粗化,使得较不满意的单元比较满意的单元大;以及
将所粗化的三角形网格在与所述地下地质盆地中的所述平面正交的方向上进行投影以形成棱柱体栅格,
将所述粗化的三角形网格的每一单元分割成子单元从而逼近所述地层;
为每一子单元中的所述通量分配多个自由度,其中分配包括,对于每一子单元:
为所述物理过程分配一个自由度;以及
为子单元的每一面分配另一自由度;
将混合有限元分析应用到所述子单元中的每一个以产生矩阵,其中应用包括使用恒定划分法来形成有限元空间;以及
求解所述矩阵以确定所述地下地质盆地中所述物理过程的通量。
2.根据权利要求1所述的方法,其中粗化包括:
将所述数据投影在平面上;以及
使用所述数据来确定哪些单元是较不满意的。
3.根据权利要求2所述的方法,其中所述模型包括对所述地下地质盆地中的物理特征进行建模的建模特征,以及其中使用包括:
将优先级值分配到每一单元,其中基于每一单元是否接近建模特征和所述建模特征的类型来确定所述优先级值。
4.根据权利要求1所述的方法,其中提供三角形网格包括:
在所述平面上提供矩形网格;以及
沿着至少一个对角线分割所述矩形网格的每一单元。
5.根据权利要求1所述的方法,其中粗化包括:
通过消除两个相邻三角形的公共的边来合并所述两个相邻三角形。
6.根据权利要求1所述的方法,其中所述棱柱体栅格包括多个棱柱体单元、多个棱锥体单元和不同于所述多个棱柱体单元和所述多个棱锥体单元的多个四面体单元。
7.根据权利要求1所述的方法,其中所述物理过程是对流-扩散过程。
8.根据权利要求1所述的方法,其中所述物理过程是温度和压力中的一个。
9.根据权利要求1所述的方法,其中所述物理过程包含碳氢化合物材料的形成。
10.根据权利要求1所述的方法,其中所述物理过程包含碳氢化合物材料的运动。
11.根据权利要求1所述的方法,还包括:
从来自传感器的信息导出所述数据,所述传感器测量所述地下地质盆地的所述至少一个物理特性。
12.一种用于对在地下地质盆地中的物理过程和该物理过程的通量进行建模以控制来自所述地下地质盆地的碳氢化合物开采的方法,所述方法包括:
形成对地下地质盆地进行建模的非结构化棱柱体栅格,其中所述物理过程在所述地下地质盆地内进行,并且所述棱柱体栅格包括多个单元;
对于每一单元,为所述物理过程和所述通量分配多个自由度;
将混合有限元分析应用到所述单元中的每一个以产生矩阵;以及
求解所述矩阵以确定所述地下地质盆地中的所述物理过程和所述通量;
其中形成包括:
在所述地下地质盆地的模型的平面上提供三角形网格,其中所述网格包括多个单元;
以非均匀方式对所述三角形网格进行粗化,使得较不满意的单元较大;以及
将所粗化的三角形网格在与所述地下地质盆地中的所述平面正交的方向上进行投影,以形成所述棱柱体栅格。
13.根据权利要求12所述的方法,其中所述棱柱体栅格包括多个棱柱体单元、多个棱锥体单元和不同于所述多个棱柱体单元和所述多个棱锥体单元的多个四面体单元。
14.根据权利要求12所述的方法,其中分配包括:
对于每一单元,为所述物理过程分配一个自由度;以及
对于每一单元,为所述单元的每一面分配另一自由度。
15.根据权利要求12所述的方法,其中应用包括:
使用恒定划分法来形成有限元空间。
16.根据权利要求12所述的方法,还包括:
使用所确定的物理过程和通量来影响所述地下地质盆地中的改变。
17.根据权利要求12所述的方法,其中所述物理过程是温度和压力中的一个。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US776107P | 2007-12-14 | 2007-12-14 | |
US61/007,761 | 2007-12-14 | ||
PCT/US2008/080515 WO2009079088A1 (en) | 2007-12-14 | 2008-10-20 | Modeling subsurface processes on unstructured grid |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101903803A CN101903803A (zh) | 2010-12-01 |
CN101903803B true CN101903803B (zh) | 2013-05-08 |
Family
ID=40795839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008801211426A Expired - Fee Related CN101903803B (zh) | 2007-12-14 | 2008-10-20 | 在非结构化栅格上对地下过程进行建模 |
Country Status (6)
Country | Link |
---|---|
US (1) | US8396699B2 (zh) |
EP (1) | EP2223158A4 (zh) |
CN (1) | CN101903803B (zh) |
BR (1) | BRPI0820830B1 (zh) |
CA (1) | CA2703253C (zh) |
WO (1) | WO2009079088A1 (zh) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2930350B1 (fr) * | 2008-04-17 | 2011-07-15 | Inst Francais Du Petrole | Procede pour rechercher des hydrocarbures dans un bassin geologiquement complexe,au moyen d'une modelisation de bassin |
AU2009341852B2 (en) | 2009-03-11 | 2015-08-20 | Exxonmobil Upstream Research Company | Adjoint-based conditioning of process-based geologic models |
CA2754690C (en) * | 2009-03-11 | 2017-01-24 | Exxonmobil Upstream Research Company | Gradient-based workflows for conditioning of process-based geologic models |
FR2948215B1 (fr) * | 2009-07-16 | 2011-06-24 | Inst Francais Du Petrole | Methode pour generer un maillage hexa-dominant d'un milieu souterrain faille |
US10061060B2 (en) | 2009-11-12 | 2018-08-28 | Exxonmobil Upstream Research Company | Method and apparatus for generating a three-dimensional simulation grid for a reservoir model |
CA2783977C (en) | 2010-02-12 | 2018-11-20 | Exxonmobil Upstream Research Company | Method and system for partitioning parallel simulation models |
US9754056B2 (en) | 2010-06-29 | 2017-09-05 | Exxonmobil Upstream Research Company | Method and system for parallel simulation models |
US9123161B2 (en) | 2010-08-04 | 2015-09-01 | Exxonmobil Upstream Research Company | System and method for summarizing data on an unstructured grid |
US8594985B2 (en) * | 2011-02-08 | 2013-11-26 | International Business Machines Corporation | Techniques for determining physical zones of influence |
US10175386B2 (en) | 2011-02-09 | 2019-01-08 | Saudi Arabian Oil Company | Sequential fully implicit well model with tridiagonal matrix structure for reservoir simulation |
US10113400B2 (en) | 2011-02-09 | 2018-10-30 | Saudi Arabian Oil Company | Sequential fully implicit well model with tridiagonal matrix structure for reservoir simulation |
US9164191B2 (en) | 2011-02-09 | 2015-10-20 | Saudi Arabian Oil Company | Sequential fully implicit well model for reservoir simulation |
MX2013009702A (es) * | 2011-02-22 | 2014-07-09 | Landmark Graphics Corp | Generacion de datos para modelado geomecanico. |
US8928661B2 (en) | 2011-02-23 | 2015-01-06 | Adobe Systems Incorporated | Representing a field over a triangular mesh |
CN102194252A (zh) * | 2011-05-17 | 2011-09-21 | 北京航空航天大学 | 一种基于地质层面结构的三角形格架网格生成方法 |
CN102436514B (zh) * | 2011-08-25 | 2013-05-08 | 上海现代建筑设计(集团)有限公司 | 一种流固耦合网格更新的方法 |
US10241218B2 (en) * | 2012-05-30 | 2019-03-26 | Pgs Geophysical As | Methods and systems for computing notional source signatures from near-field measurements and modeled notional signatures |
CN102722622A (zh) * | 2012-06-08 | 2012-10-10 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机机舱建模方法 |
CN102819865B (zh) * | 2012-08-09 | 2015-07-01 | 成都理工大学 | 一种大地电磁三维地质结构模型的建模方法 |
AU2013350307A1 (en) * | 2012-11-20 | 2015-07-09 | Stochastic Simulation Limited | Method and system for characterising subsurface reservoirs |
US20140257701A1 (en) * | 2013-03-11 | 2014-09-11 | Schlumberger Technology Corporation | Output compression for geological strata physical property modeling |
US10311180B2 (en) * | 2014-07-15 | 2019-06-04 | Dassault Systemes Simulia Corp. | System and method of recovering Lagrange multipliers in modal dynamic analysis |
GB2533847B (en) * | 2014-11-06 | 2017-04-05 | Logined Bv | Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory |
US10061044B2 (en) * | 2015-02-10 | 2018-08-28 | Transcend Engineering and Technology, LLC | Systems, methods, and software for detecting the presence of subterranean tunnels and tunneling activity |
CA2972149A1 (en) * | 2015-02-25 | 2016-09-01 | Exxonmobil Upstream Research Company | Method for parameterizing a 3d domain with discontinuities |
CN105205302A (zh) * | 2015-04-08 | 2015-12-30 | 辽宁达能电气股份有限公司 | 基于光纤测温主机的电缆动态流量计算方法 |
US10061878B2 (en) | 2015-12-22 | 2018-08-28 | Dassault Systemes Simulia Corp. | Effectively solving structural dynamics problems with modal damping in physical coordinates |
US10113421B2 (en) | 2016-05-09 | 2018-10-30 | Schlumberger Technology Corporation | Three-dimensional fracture abundance evaluation of subsurface formations |
US10853533B2 (en) | 2016-05-09 | 2020-12-01 | Schlumberger Technology Corporation | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof |
US10650107B2 (en) | 2016-05-09 | 2020-05-12 | Schlumberger Technology Corporation | Three-dimensional subsurface formation evaluation using projection-based area operations |
US11609354B2 (en) | 2016-06-02 | 2023-03-21 | Shell Usa, Inc. | Method of processing a geospatial dataset |
CA3043231C (en) | 2016-12-23 | 2022-06-14 | Exxonmobil Upstream Research Company | Method and system for stable and efficient reservoir simulation using stability proxies |
US10913901B2 (en) | 2017-09-12 | 2021-02-09 | Saudi Arabian Oil Company | Integrated process for mesophase pitch and petrochemical production |
CN108459989B (zh) * | 2018-03-19 | 2021-05-07 | 中国气象科学研究院 | 一种非结构网格气象数值模式计算系统 |
US20200050722A1 (en) * | 2018-08-10 | 2020-02-13 | OnScale, Inc. | Hybrid meshing method for finite element analysis |
CN109345018A (zh) * | 2018-10-09 | 2019-02-15 | 浙江海洋大学 | 一种海洋污染物的扩散预测方法 |
CN109376433B (zh) * | 2018-10-26 | 2020-06-09 | 北京市水文总站 | 基于土壤非饱和水和地下水耦合的区域水流运动模拟方法 |
US11713675B2 (en) | 2020-05-01 | 2023-08-01 | Landmark Graphics Corporation | Determining exploration potential ranking for petroleum plays |
CN115906559B (zh) * | 2022-10-31 | 2023-08-15 | 重庆大学 | 一种基于混合网格的大地电磁自适应有限元正演方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1358289A (zh) * | 1999-06-24 | 2002-07-10 | 埃克森美孚上游研究公司 | 按比例增大非结构化栅格(grid)的渗透率的方法 |
US7222126B2 (en) * | 2002-07-30 | 2007-05-22 | Abel Wolman | Geometrization for pattern recognition, data analysis, data merging, and multiple criteria decision making |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT1294492B1 (it) | 1997-09-16 | 1999-04-12 | Enel Spa | Metodo per la rappresentazione tramite griglia di calcolo semistrutturata di fenomeni fisici estendentisi in un dominio spaziale |
US7043413B2 (en) | 2000-06-29 | 2006-05-09 | Object Reservoir, Inc. | Method for modeling an arbitrary well path in a hydrocarbon reservoir using adaptive meshing |
US7369973B2 (en) | 2000-06-29 | 2008-05-06 | Object Reservoir, Inc. | Method and system for representing reservoir systems |
GB0021284D0 (en) | 2000-08-30 | 2000-10-18 | Schlumberger Evaluation & Prod | Compositional simulation using a new streamline method |
US7099805B2 (en) | 2001-08-30 | 2006-08-29 | International Business Machines Corporation | Tetrahedralization of non-conformal three-dimensional mixed element meshes |
WO2003060754A1 (en) | 2001-12-31 | 2003-07-24 | The Board Of Regents Of The University And Community College System, On Behalf Of The University Of Nevada, Reno | Multiphase physical transport modeling method and modeling system |
CN100408401C (zh) * | 2002-06-26 | 2008-08-06 | 杰里·M·罗恩 | 大量运输系统 |
JP3988925B2 (ja) | 2002-07-04 | 2007-10-10 | 学校法人慶應義塾 | 混合格子型解適合格子法を用いた数値解析装置 |
GB2396448B (en) | 2002-12-21 | 2005-03-02 | Schlumberger Holdings | System and method for representing and processing and modeling subterranean surfaces |
US7496488B2 (en) | 2003-03-06 | 2009-02-24 | Schlumberger Technology Company | Multi-scale finite-volume method for use in subsurface flow simulation |
US8148475B2 (en) | 2003-06-30 | 2012-04-03 | Lubrizol Advanced Materials, Inc. | Melt spun polyether TPU fibers having mixed polyols and process |
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 |
-
2008
- 2008-10-20 WO PCT/US2008/080515 patent/WO2009079088A1/en active Application Filing
- 2008-10-20 EP EP08861582.8A patent/EP2223158A4/en not_active Withdrawn
- 2008-10-20 BR BRPI0820830A patent/BRPI0820830B1/pt not_active IP Right Cessation
- 2008-10-20 US US12/738,950 patent/US8396699B2/en active Active
- 2008-10-20 CA CA2703253A patent/CA2703253C/en not_active Expired - Fee Related
- 2008-10-20 CN CN2008801211426A patent/CN101903803B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1358289A (zh) * | 1999-06-24 | 2002-07-10 | 埃克森美孚上游研究公司 | 按比例增大非结构化栅格(grid)的渗透率的方法 |
US7222126B2 (en) * | 2002-07-30 | 2007-05-22 | Abel Wolman | Geometrization for pattern recognition, data analysis, data merging, and multiple criteria decision making |
Also Published As
Publication number | Publication date |
---|---|
BRPI0820830A2 (pt) | 2015-06-16 |
US8396699B2 (en) | 2013-03-12 |
BRPI0820830B1 (pt) | 2019-08-13 |
WO2009079088A1 (en) | 2009-06-25 |
US20100211370A1 (en) | 2010-08-19 |
EP2223158A1 (en) | 2010-09-01 |
CA2703253C (en) | 2015-06-02 |
EP2223158A4 (en) | 2017-12-27 |
CA2703253A1 (en) | 2009-06-25 |
CN101903803A (zh) | 2010-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101903803B (zh) | 在非结构化栅格上对地下过程进行建模 | |
US11352855B2 (en) | Modeling reservoir flow and transport processes in reservoir simulation | |
US11066907B2 (en) | Sequential fully implicit well model with tridiagonal matrix structure for reservoir simulation | |
US10769326B2 (en) | Parallel solution for fully-coupled fully-implicit wellbore modeling in reservoir simulation | |
US9069102B2 (en) | Scalable simulation of multiphase flow in a fractured subterranean reservoir with multiple interacting continua by matrix solution | |
EP2599028B1 (en) | Method and system for parallel multilevel simulation | |
US9922142B2 (en) | Systems and methods for subsurface reservoir simulation | |
Mustapha et al. | An efficient method for discretizing 3D fractured media for subsurface flow and transport simulations | |
Garcia et al. | Automatic grid generation for modeling reservoir heterogeneities | |
CN102667804A (zh) | 使用均匀化混合有限元建模地质性质的方法和系统 | |
US20160161635A1 (en) | A simulation-to-seismic workflow construed from core based rock typing and enhanced by rock replacement modeling | |
US10175386B2 (en) | Sequential fully implicit well model with tridiagonal matrix structure for reservoir simulation | |
Mustapha et al. | Discretizing two‐dimensional complex fractured fields for incompressible two‐phase flow | |
Abbas et al. | Clustering analysis and flow zone indicator for electrofacies characterization in the upper shale member in Luhais Oil Field, Southern Iraq | |
US11073001B2 (en) | Sequential fully implicit horizontal well model with tridiagonal matrix structure for reservoir simulation | |
CN109072688B (zh) | 用于储层模拟的具有三对角线矩阵结构的连续的全隐式井模型 | |
US11300706B2 (en) | Designing a geological simulation grid | |
Foroud et al. | An efficient optimization process for hydrocarbon production in presence of geological uncertainty using a clustering method: A case study on Brugge field | |
Clemo | Dual permeability modeling of fractured media | |
Zhang et al. | Model errors inherent in conditioning a stochastic channel to pressure data | |
Martin et al. | Reservoir modeling of low-sinuosity channel sands: a network approach | |
US20220259959A1 (en) | Deforming Well Trajectories | |
Mohamed et al. | Complex unstructured-grid groundwater modeling using centroidal Voronoi tessellation refinement and curve fitting | |
Balaha | Modelling of multiphase fluid flow in heterogeneous reservoirs | |
Alpak | A mimetic finite-volume-discretization operator for reservoir simulation |
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: 20130508 Termination date: 20201020 |
|
CF01 | Termination of patent right due to non-payment of annual fee |