在非结构化栅格上对地下过程进行建模的制作方法

专利名称：在非结构化栅格上对地下过程进行建模的制作方法
技术领域：
该申请一般地涉及计算机建模，具体地说，涉及用非结构化栅格对地下过程进行建模。
背景技术：
在地质勘探中，期望获得关于地球表面下存在的各种岩层和结构的信息。这类信 息可以包括测定的地质层、密度、孔隙度、组成成分等。接着，这个信息用来对地下盆地进行 建模，即使用获得的数据预测碳氢化合物储藏的位置并且帮助碳氢化合物的提取。非结构化栅格对于对复杂的地质结构(例如地下盆地)中的物理过程进行建模具 有许多有吸引力的特性。这种栅格也可以用于其他工业中，例如航空工业或汽车工业。所 关心的盆地或域可以被建模或表示为堆叠在一起的不同厚度的层的集合。地质层可沿着垂 直表面或倾斜表面断裂并且退化，生成所谓的尖灭(Pinch-out)。尖灭被定义为具有接近零 厚度的地质层的部分。栅格应该考虑这种复杂性以产生地质层的良好模型。非结构化栅格 提供比结构化栅格更好的模型。非结构化栅格可以包括由它们的顶点定义的多面体元素或 单元的集合并且具有完全任意的拓扑图。例如，栅格的顶点可以属于许多个单元并且每一 单元可以具有任何数目的边或面。可以通过对流_扩散类型的数学方程来描述许多物理地下过程。这类过程的示例 可以是多孔介质中的流体流、温度分布和/或压力分布。用于石油勘探的重要过程是温度 分布或热建模。热建模包含从地壳下的岩浆并且穿过沉积层和源岩运动的热量。源岩是包 含在石油和其他碳氢化合物材料的地层中的岩石。石油和/或其他碳氢化合物材料从源岩 中被排出并且迁移到其他地方。碳氢化合物的质量由源岩和它们周围区域上承受的温度和 压力条件确定。质量也受源岩和其当前位置之间的迁移路径的温度和压力条件影响。因此， 盆地在其整个历史中的压力和温度条件是重要的。为了更精确地对过程进行建模，重要的是不仅对主变量(例如压力或温度)进行 建模，而且对在任何给定的表面之上的主变量的通量或能量、流体的流速等进行建模。有多 种已知的方法用于对这些过程进行建模，例如有限差分、有限体积或有限元方法。在这些方 法中，在考虑物理过程的地方，域被栅格覆盖。接着，通过以下方法来在栅格上逼近域在栅 格单元的指定位置处引入称为自由度的一组未知量并且为每一位置导出代数方程，所述方 程将在那个位置的自由度与其他自由度连接起来。对于上面所提到的不同方法，导出这类 方程的方式和自由度的位置是不同的，但是所有这些方法具有共同的特征，即，它们仅包含 主变量，例如温度或压力。
为了计算通量，感兴趣的研究人员将首先使用上面所描述方法中的一个来计算所 期望的主变量，接着使用数值微分来计算主变量的通量。在规则栅格(例如矩形或平行六 面体栅格)上精确的所有现有数值微分方法在非结构化栅格上是不精确的并且在计算上 非常耗时，尤其是如果考虑物理过程的域是高度不均勻的或异质的。此外，使用有限差分方 法求解对流-扩散问题的方法需要笛卡尔(Cartesian)栅格，因此不可应用于不得不使用 非结构化栅格的许多地下应用。能够对复杂几何形状进行建模的有限元方法不具有局部守 恒属性并且不能够被应用在许多地下过程中。相反地，有限体积方法是局部守恒的并且能 够被应用在局部正交的非结构化栅格的子集上。然而，当非结构化栅格不具有局部正交属 性时，有限体积方法提供不精确的解。因此，从上面所提到的所有三种方法中，没有一个可 应用于描述以非结构化栅格建模的盆地中的地下对流-扩散过程。有另一数学方法来同时逼近主未知量和它们的通量，叫做混合有限元方法，其在 以下文献中描述F. Brezzi 和 M. Fortin 的“Mixed and hybrid finite elementmethods”， Springer Verlag，柏林，1991年。这种方法被证明在存在不均勻介质的情况下是局部质量 守恒的、精确的，并且提供对主未知量和通量两者的精确逼近。直到最近，混合有限元方法 才能够被直接应用到由非结构化多面体栅格覆盖的域，所述非结构化多面体栅格对于地下 应用是非常常见的。在以下文献中提出了针对在任意多面体栅格上的扩散类型方程的混合 有限元方法的新变型Yu. Kuznetsov 禾口 S. R印in 的"New mixed finite element method on polygonal andpolyhedral meshes", Russian Journal of Numerical Analysis and MathematicalModeling，第 18 卷，261-278 页，2003 年。

发明内容
本发明涉及对在经过地质时期形成的盆地中的能量传递和/或压力分布提供精 确模型的系统和方法。在任何给定的时间，盆地被表示为堆叠在一起的不同厚度的层的集 合。在盆地中的某些位置，层的厚度退化到零，形成尖灭。本发明的实施例使用棱柱体网 格和混合有限元分析来对盆地中的各种过程进行建模，这些过程包括能量传递(例如热能 量)和压力。因此，本发明的实施例求解主未知量(例如温度或压力)和从未知量(例如温 度通量或压力通量)两者。下述方面中的一个或更多个可以用来提供经过地质时期形成的 盆地中的能量传递和/或压力分布(例如物理过程)的精确模型。该模型可以用来解释现 代储油层，并且继而依赖于该模型，从而基于模型的模拟结果来控制碳氢化合物开采活动。 基于从由下述方面中的一个或多个产生的模拟盆地模型中解释的结果，可以控制碳氢化合 物的开采，例如，可以控制来自地面设施的开采速率，可以有策略地打井，和/或一般性地 表征储油层。在一个一般方面，一种用于在计算机上对物理区域进行建模的方法，其中所述物 理区域包括多个地层，所述方法包括接收定义所述物理区域的至少一个物理特性的数据； 在所述物理区域的模型的平面上提供三角形网格，其中所述网格包括多个单元；以非均勻 方式对所述三角形网格进行粗化，使得较不满意的单元较大；以及将粗化的三角形单元在 与所述物理区域中的所述平面正交的方向上进行投影以形成棱柱体栅格，其中根据所述地 层，粗化的三角形网格的每一单元被分割成子单元。这个方面的实现可以包括下述特征中的一个或更多个。例如，粗化可以包括将所
6述数据投影在平面上以及使用所述数据来确定哪些单元是较不满意的。所述模型可以包 括对所述物理区域中的物理特征建模的建模特征，以及其中使用可以包括将优先值分配给 每一单元，其中基于每一单元是否接近建模特征和该建模特征的类型来确定该值。提供三 角形网格可以包括在所述平面上提供矩形网格；以及沿着至少一个对角线分割所述矩形网 格的每一单元。粗化可以包括通过消除两个相邻三角形的公共边来合并所述两个相邻三角 形。所述棱柱体栅格可以包括多个棱柱体单元、多个棱锥体单元和多个四面体单元。所述 方法可以用来对所述物理区域中的物理过程的至少一个通量进行建模，所述方法还包括为 每一子单元中的通量分配多个自由度；将混合有限元分析应用到每一子单元以产生矩阵； 以及求解该矩阵以确定所述区域中的物理过程的通量。分配可以包括对于每一单元，为所述物理过程分配一个自由度；以及为所述单元 的每一面分配另一自由度。应用可以包括使用恒定划分法来形成有限元空间。所述物理过 程可以是对流-扩散过程。所述物理过程可以是温度和压力中的一个并且所述物理区域是 地下地质盆地。所述物理过程可以包含碳氢化合物材料的形成。所述物理过程可以包含碳 氢化合物材料的运动。可以从来自传感器的信息导出所述数据，所述传感器测量所述物理 区域的至少一个物理特性。在另一一般方面，一种用于在计算机上对物理过程和该物理过程的通量进行建模 的方法，其包括形成对物理区域建模的非结构化棱柱体栅格，其中所述物理过程在所述物 理区域内发生并且所述棱柱体栅格包括多个单元；对于每一单元，为所述物理过程和所述 通量分配多个自由度；将混合有限元分析应用到每一单元以产生矩阵；以及求解该矩阵以 确定所述区域中的所述物理过程和所述通量。这个方面的实现可以包括下述特征中的一个或更多个。例如，形成可包括在所述 物理区域的模型的平面上提供三角形网格，其中所述网格包括多个单元；以非均勻方式对 所述三角形网格进行粗化，使得较不满意的单元较大；以及将粗化的三角形单元在与所述 物理区域中的所述平面正交的方向上进行投影以形成所述棱柱体栅格。所述棱柱体栅格可 以包括多个棱柱体单元、多个棱锥体单元以及多个四面体单元。分配可以包括对于每一单 元，为所述物理过程分配一个自由度；以及对于每一单元，为所述单元的每一面分配另一自 由度。应用可以包括使用恒定划分法来形成有限元空间。所确定的物理过程和通量可以用 来影响所述物理区域中的改变。所述物理过程可以是温度和压力中的一个并且所述物理区 域是地下地质盆地。在另一一般方面，一种具有计算机可读介质的计算机程序产品，所述计算机可读 介质具有在其上记载的计算机程序逻辑，所述计算机程序逻辑用于对物理区域中的物理过 程和该物理过程的通量进行建模，所述计算机程序产品包括用于形成对所述物理区域进行 建模的非结构化棱柱体栅格的代码；用于将混合有限元分析应用到所述棱柱体栅格以产生 矩阵的代码；以及用于求解该矩阵从而确定所述区域中的物理过程和通量的代码。这个方面的实现可以包括下述特征中的一个或更多个。例如，用于形成的代码可 以包括这样的代码用于在所述物理区域的模型的平面上提供三角形网格的代码，其中所 述网格包括多个单元；以非均勻方式对所述三角形网格进行粗化，使得较不满意的单元是 较大的；以及将粗化的三角形网格在与所述物理区域中的所述平面正交的方向上进行投影 以形成所述棱柱体栅格。所述棱柱体栅格可以包括多个单元，以及所述用于应用的代码可以包括对于每一单元，为所述物理过程分配一个自由度；对于每一单元，为所述单元的每一 面分配另一自由度；以及使用恒定划分法来形成有限元空间。所述用于求解的代码可以包 括使用预调的共轭梯度分析来求解所述矩阵。本发明的实施例通过将某些或大多数地质和几何特征，例如尖灭边界投影到水平 平面中来工作。注意，投影可以是非正交的或倾斜的。接着，本发明的实施例在那个平面 上生成分解所有期望的特征的非结构化栅格。注意到，栅格可以由多边形、四边形、三角形 或其组合组成。接着，本发明的实施例将获得的栅格投影回到所有层的所有边界表面上，从 而构造棱柱体栅格。所述棱柱体栅格可以包括多个单元，这些单元可以是棱柱体、四面体形 状、棱锥体或其组合。注意，非结构化棱柱体栅格逼近所有层的边界表面。接着本发明的实施例可以通过以下步骤来工作对于主未知量，在单元中心处为 每一单元关联一个自由度，并且对于通量的法向分量，在面中心处为所述单元的每一面关 联一个自由度。接着本发明的实施例使用混合有限元方法，例如Yu. Kuznetsov和S. Repin 的方法来使问题离散化。空间离散化产生稀疏矩阵方程。接着本发明的实施例可以求解该 矩阵方程，以得到主未知量和在所述单元的各面处的通量的法向分量两者。因此，本发明的 实施例在没有极大地扩展需要求解的未知量的数目的情况下，提供更精确的建模。本发明的实施例可以提供覆盖所述物理区域的水平平面的三角形网格来形成所 述棱柱体栅格。接着可以以非均勻方式对所述网格进行粗化，使得较不满意的单元是较大 的，而使满意的单元保留更细的形式。接着将粗化的网格在所述物理区域中的垂直方向上 进行投影以形成所述棱柱体网格。为了可以更好地理解后面的本发明的详细描述，前面已经相当广泛地概述了本发 明的特征和技术优势。在下文中将描述本发明的另外特征和优势，其构成本发明的权利要 求的主题。本领域技术人员应意识到，公开的概念和具体的实施例可以容易地用作修改或 设计其他结构的基础，其他结构用于实现本发明的相同目的。本领域技术人员也应该认识 到，这类等效的构造不偏离随附的权利要求中陈述的本发明的精神和范围。当结合附图考 虑时，从下述描述中将更好地理解被认为是本发明的特性的新颖特征(关于其组织架构和 操作方法两者)以及其他目标和优势。然而，应清楚地明白，提供每一附图是仅为了解释和 描述的目的，无意作为本发明的限制的定义。


为了更彻底地理解本发明，现在结合附图参考下述描述，在附图中图1描绘根据本发明的实施例的被分割成层的域。图2A和2B描绘根据本发明的实施例的域和用矩形网格覆盖的域。图3描绘根据本发明的实施例的非均勻粗化的三角形栅格。图4描绘根据本发明的实施例形成的不同类型的单元。图5描绘根据本发明的实施例形成的3D棱柱体栅格。图6描绘本发明的实施例使用的四面体单元。图7描绘本发明的实施例使用的棱锥体单元。图8A-8D根据本发明的实施例描绘本发明的实施例使用的、被分成三个四面体的 棱柱体单元。
图9描绘根据本发明的实施例的分割相邻单元的独立面；图10描绘根据本发明的实施例的形成棱柱体栅格的方法；图11描绘根据本发明的实施例的求解矩阵的方法；以及图12描绘适合使用本发明的计算机系统的框图。
具体实施例方式本发明的实施例对于对地下油田进行建模是有用的。在此所描述的实施例的示例 可以参考这类油田。然而，实施例可以用来对包含其他材料和/或过程的的其他域进行建 模。例如，可以包含其他碳氢化合物材料，例如煤炭。本发明的实施例对于采矿或掘坑道是 有用的。本发明的实施例可以用于其他域类型，例如大气，并且对于对天气、温度和/或污 染进行建模是有用的。另一域可以是海洋，并且实施例可以用来测量声音、温度、含盐度和 /或污染度。可以使用本发明的实施例来对任何类型的分层域进行建模。可以使用本发明 的实施例来对通过对流-扩散过程而运动的任何类型的材料进行建模。可以使用本发明的 实施例来对域或材料中存在的任何类型的通量进行建模。如前面所述的，可以将本发明的实施例应用到任何对流_扩散过程。下面是3D对 流-扩散类型方程的示例-V.(/▽/ ) + ￡ ·； = /在 Ω 中(1.1)这里ρ是未知函数(称为压力)，K = K(X)是扩散张量，c是非负函数，f是源函 数，以及Ωci 3是有界计算域。假设κ是均勻正定矩阵，并且域Ω的边界an被分割成两个 不重叠的集合rD和rN。用边界条件补充方程(1. 1)ρ = gD 在 Γ D 上(/:Vp).n + iT.p = gA^rN±(1.2)这里η是。的外向单位法向向量，σ是非负函数，并且gD和&是给定函数。假 设方程(1.1)-(1.2)具有唯一解。可以由等价一阶系统代替偏微分方程(1. 1)-(1. 2)U+ ArVp = O在 Ω 中▽·!! + ￡·./ = /在 Ω 中(1.3)ρ = gD 在 Γ D 上-U · η+ σ · ρ = gN 在 Γ N 上(1.4)方程(1.3)-(1.4)是方程(1. 1)-(1.2)的混合公式表示。注意，以此方式，可以同 时逼近主未知量P及其通量U。如前面所述的，本发明的实施例可以在不同的域工作。因此，设G是R2中具有规则 成形的边界5Ω的域，例如分段光滑并且段之间的角度大于0。设计算域Ω被定义为如下Ω = {(X，y, ζ) e R3 (χ, y) E G, Zmin (x，y) ^ ζ ^ Zmax (x，y)}这里Zmin(χ, y)和 Zmax (χ, y)是光滑表面。设m是正整数并且ζ = Zi (X，y)，i = 0，. . .，m，是定义在G上的单值连续函数，使 得
Z0 (x, y) ^ Zmin(χ, y) 在G中Zi^1 (x, y) ^ Zi (χ, y)在G中=Zm (X，y) = Zmax (X，y)在G 中这些函数定义地质层之间的界面。换句话说，计算域Ω可以被分成m个子域(带 或层)，其被定义如下对于所有i = l，...，m，Qi = {(χ, y, ζ) e Ω (χ, y) E G, Zi^1 (χ，y) ^ ζ ^ Zi (χ, y)}图1描绘了将计算域Ω 100分割成多个子域或层101-107的示例。注意，图1以 分解图描绘了不同的层，然而对于计算，不需要将层隔开。还注意到，假设子域Qi满足锥 条件，这里子域的边界不具有奇点(零角度等)，并且此外，所有集合Pi = {(x,^) e G Zi^x, y) = Zt(x,y)}由有限数目的多边形组成。用表示对应的集合PiW边界。图1描绘了盆地的不同地层。可以通过各种技术例如地层分析和/或地震反演， 使用传感器来测量盆地的各种特性来确定用来形成图1的不同层的数据。可以以下述方式将集合Pi的边界投影到平坦平面。对于来自閉的任何给定的点 (x，y，z)，投影的点具有坐标(x，y，0)。所有这类点组成像图2A那样的闭合线的集合，其用 来生成平面三角测量，在图2B和图3中示出了这种示例。微分方程(1.3)_(1.4)的变分混合公式可以被写成如下找到《ε^^Ω)， P e L2(Q)和 λ e L2(ΓN)，使得对于所有νe ^ftv(Q)，q e L2(Q)和 μ e L2(Tn) 这里Hdiv(Q) = {ν: ν e [I2(Ω)]3,V· ν e L2(Q), J|vn|2^< }。
an注意，λ是压力函数ρ = ρ(χ)在ΓΝ上的约束。在这个公式中，所有边界条件是 自然边界条件。在σ =0的情况下，可以以不同的形式将变分混合公式写为如下找到 Iie^iiv(Q),U ·η = -&(在ΓΝ上)和ρ G L2 (Ω)，使得对于所有ν e点咖…)，ν ·η = 0 (在 1\上)和 q e L2(Q) 在下述分析中，考虑方程(1.5)，然而也可以将结论应用到方程(1.6)而不失一般性。
10
本发明的实施例使用棱柱体栅格，其在对对流-扩散地下过程进行建模方面提供 许多有吸引力的特性。在许多情况下，域可以被表示为堆叠在一起的不同厚度的层的集合。 棱柱体栅格满意地表示了堆叠层中的分片地质结构和非结构化几何特征。通过地质统计信 息的后处理来提供2D几何数据。通常，数据是与2D细矩形栅格的节点或单元关联的材料 属性。注意，可能有数百万个节点。节点中存在材料数据意味着其是计算节点，而不存在材 料数据则意味着其是局外节点(node-outlier)。计算节点的集合定义了计算域。地质层之间的界面是地质统计数据，并且可能互相相交，导致拓扑上不正确的情 形。底部界面是最低的地质层的底部边界，并且也由地质统计数据表示。地质层可以沿着 垂直表面断裂并且退化。尖灭被定义为厚度模不大于用户定义的阈值δ >0的地质层的 部分。断层折线被定义为底部地质界面和断层的相交处。为了简化算法的描述，在本发明的实施例中要使用的栅格应该满足某些非常自然 的要求。目标棱柱体栅格应该是ID栅格和2D三角测量的逻辑积，所述ID栅格和2D三角 测量的节点形成不具有局外节点的地质统计矩形栅格的子集，其中节点密度等于所关心的 用户指定区域中原始栅格中的节点密度。除此之外，在用户定义的断层和井以及自动检测 到的尖灭附近必须精炼(改进)三角测量。此外，2D栅格中三角形的数目不应该大于用户 定义的数目。关于棱柱体栅格，棱柱体的侧面不得不逼近地质界面并且形成形状规则的2D 三角测量。下面所描述的过程是可以用来构造棱柱体栅格的过程的示例。注意，可以使用其 他过程。此外，那种类型的棱柱体栅格不限制混合有限元方法。首先是为了将尖灭投影在 表示井的底部地质界面、断层折线和点上，对所述2D规则三角测量的产生进行改进。接下 来，将2D三角测量投影在由函数Zi (X，y)，i = 1，. . .，m所定义的表面上，以形成得到的3D 棱柱体栅格。在下面的段落中进一步描述这个过程。为了形成2D规则三角形栅格，示例性过程可以开始于矩形栅格。给定节点在χ-和 y-方向上的节点坐标，产生覆盖域G的矩形一致网格Gh，其由具有带材料数据的四个节点 中的至少一个的单元组成。例如，图2A描绘域200，图2B描绘覆盖该域200的矩形栅格 201。不失一般性，假设网格Gh由正方形组成，即在χ-和y-方向上的网格尺寸相等，hx = hy = h0接下来，根据本发明的实施例，每一矩形单元被其对角线分割成两个三角形。下面 描述可以用来形成三角形的一个过程，注意，可以使用其他过程。可以根据下述规则来进行 两个可能的对角线之间的选择。设每一矩形单元被分配一整数，该整数等于其节点的最小 χ-和y_索引的和。对于具有偶数的单元，分割对角线具有带最小χ-和y_索引的单元的节 点。对于具有奇数的单元，选择另一对角线。对于具有给定的材料数据分布的给定的节点 集合，上述过程唯一地指定三角测量。改变对角线的方向减少了栅格朝向的问题。作为形 成三角形的替换过程，可以通过使用两个对角线来将每一矩形单元分割成四个三角形。将产生的三角测量投影在底部地质界面上，如公式(1. 2)那一段中所描述的。注 意，在域G中可以定义感兴趣的区域ω”这里栅格的修改是不必要的或不期望的。设Pi11表示Gh的矩形元素的最大子集，其属于Pi。如果没有属于Pi的Gh的元素， 但是有属于Pi的Gh的顶点，则这个顶点被认为属于ΡΛ那么集合被定义为
这被称为“尖灭”投影。根据该定义，“尖灭”投影是属于Gh的边和顶点的子集。接下来，根据本发明的实施例，将不同的优先级分配给三角形。将优先级或整数标 记分配给细栅格的每一三角形。优先级的值控制粗化过程。在开始时，将零优先级分配到所有三角形。对于其闭合线与断层相交的三角形，它 们的优先级被改变为1。为了找到与断层相交的三角形，可以使用下述方法。首先，从断层 三角测量中提取三角形，所述断层三角测量与最底部地质界面相交。其次，检查每一提取的 三角形是否与细栅格三角形相交。对于其闭合线与尖灭相交的三角形，它们的优先级被改 变为2。这些三角形被定义为违反下述条件的三角形在所有三角形节点中，地质层的厚度 或大于δ或小于-δ。对于其闭合线包含井点的三角形，它们的优先级被改变为3。对于 属于用户定义的感兴趣的区域Qi的三角形，它们的优先级被改变为4。对于满足若干条件 的三角形，它们的优先级被改变为最大优先级值。在已经分配优先级之后，可以对三角形栅格进行非均勻粗化。栅格的细部分可以 具有大量的同样小的三角形。这些区域是更期望的，因为它们包含更多的信息，包括感兴趣 的地质特征，例如井、断层、尖灭和/或用户指示为期望的特征。栅格的粗化部分与栅格的 较细部分不是同样满意的。栅格可以具有粗化范围，这里最粗化的区域指示具有较不满意 或不满意质量的部分，并且不具有粗化的区域指示最期望的区域。在最粗化的区域和不粗 化的区域之间的粗化区域指示具有某些期望方面的区域。粗化是一系列合并三角形的过程。例如，可以通过消除它们的公共边来将两个三 角形合成一个。这个过程包括两个阶段。首先，标记特定三角形用于进行粗化。其次，对它 们进行粗化。应该注意到，栅格一致性可能引起对未标记三角形的粗化。每一粗三角形继 承了两个合并的三角形的最大优先级。除了优先级外，每一三角形还被分配表示级别的另 一整数。初始细栅格的任何三角形具有级别1。可以将粗化应用到相同级别j的一对三角 形，并且得到具有级别j+ι的更粗的三角形。下面是根据本发明的实施例的粗化过程的一个示例。注意，可以使用其他过程。粗化过程可以被描述为下述循环1)设置 k=l。2)由具有零优先级的三角形形成集合M，零优先级三角形的粗化不会引起对具有 非零优先级的其他三角形进行粗化。3)如果M为空，则到6。4)对来自M的三角形进行粗化。5)到 2。6)如果新栅格中三角形的数目不大于用户定义的阈值Nt_，则停止。7)由非零优先级不大于k的三角形形成集合M，其粗化不会引起对优先级大于k 的其他三角形进行粗化。8)对来自M的三角形进行粗化。9)如果k≤3，则设置k = k+Ι，否则设置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，，并且用ει^Λ = 1，2，3表示它们与表面ζ = 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。换句话说，棱锥体是一个边消失的棱柱体，而四面体是两个边消失的棱柱体。对《中的所有三角形执行上述操作将提供域Ω的分区Qh。在特定的单元中，通 过分段三角形表面ζ = ΖΛ) (x,y)来逼近每一表面ζ = Zi (X，y)，所述分段三角形表面由棱 柱体的顶部(底部)三角形以及棱锥体和四面体的特定面组成。图5描绘了 3D棱柱体栅 格500的示例。在完成3D棱柱体栅格之后，可以对栅格进行混合有限元分析。在前面的部分中 指出，栅格括元素{Ek}，其或是垂直棱柱体，或是棱锥体，或是四面体。为了针对方程 (1.5)用公式表示混合有限元(MFE)方法，应该定义空间".(Ω)、L2(Q) *L2(rN)的有限 元子空间。有限元空间、Ci2(Q)包括函数ph，其在每一栅格单元尽上是恒定的。有限 元空间CjL2OTw)包括函数Ah，其在Qh中的栅格单元Ek与边界部分Γ,的每一交界处上 是恒定的。这些交界处可以是四边形或是三角形。混合有限元方法中的一个问题是设计空间^iv(Q)的有限元子空间Vh。为了计算 效率，应该仅考虑那些有限元向量函数，它们在相邻单元1)之间的界面rkl 上以及在单元Ek与边界Γ ,的交界处上具有恒定的法向分量。空间&的有限元子 空间的维数等于不同界面{rkl}和^1fov丨的总数目。可以基于在以下文献中描述的“恒定划 分(div-constant) ”法来构造这个有限元空间Yu. Kuznetsov 和 S. Repin 的 “New mixed finite element method on polygonal andpolyhedral meshes", Russian Journal of Numerical Analysis and MathematicalModeling，H 18 ·，261—278 M，2003 $。对于四面体单元T，有限元空间Vh I τ与经典的最低阶Ravi art-Thomas有限元空间 RT0(T)—致(参见F. Brezzi 禾ΡΜ· Fortin 的"Mixed and hybrid finite elementmethods，，， Springer Verlag，柏林，1991年)。有限元向量值函数wh e RTtl (T)具有四个自由度(DOF)， 即 这里ΦΑχ)是与四面体T的面Y i (i = 1，2，3，4)关联的基向量函数。用γ j表示的四面体T的面与顶点Aj相对，即，Y工是面A2A3A4，γ 2是面A3A4A1，γ 3 是面A4A1A2，Y4是面A1A2Ay设η」是在面Y」上的外向单位向量，并且h」是从顶点Aj到面 Y j上的垂直线长度，j = 1，2，3，4。在图6中示出了这种四面体单元600。在四面体T上最低阶Raviart-Thomas元素的空间可以被定义为RT0(T) = 8ρΒπ{Φ1, Φ2, Φ3，ΦJ这里基向量函数Cti满足条件论并且δ。是kronecker符号，如果i = j，则\」等于1，否则其等于0。简单明了的计算显示了基函数可以被明确地定义为 这里χ⑴是顶点Ai的坐标，i = 1,2,3,4.对于棱锥体单元P，当单元P e Qh是四边形棱锥体时，使用“恒定划分”法来构造 有限元空间Vh ι P。用γ」，j = 1，2，3，4，5 来表示棱锥体 P 的面，SP，Y1 是面 A2A3A5，Y2 是面 A1A3A4, Y3 是面A1A2A4A5，γ 4是面A1A2A3，以及Y 5是面A3A4A5。在图7中示出了这种棱锥体单元700。为了描述“恒定划分”法，将棱锥体P分成两个四面体T1和1~2。可以以两种不同的 方式来完成，并且两种方式都提供下面所描述的可工作的算法。因此，不失一般性，假设棱 锥体P被划分成两个四面体T1 = A1A2A3A4和T2 = A2A3A4A50设Iii表示到面、,的单位外向 法向向量，i = 1，2，3，4，5，并且Ii6表示到四面体T1和T2之间的界面γ 6 = A2A3A4的、从T1 指向T2的单位法向向量。在每一四面体上，构造向量函数的经典的最低阶Raviart-Thomas空间，确定基函 数MwA1 (k = 1,2)的集合，并且P中的向量场Uh被定义为如下 这个表达式显示出向量函数Uh在四面体T1和T2中的每一个上是线性的，其属于 空间Hdiv(P)，并且满足所要求的条件 为了确定在界面Y6上通量的法向分量的未知值U6，下述条件应该是可操作的▽·、=COnst (常数)在 P 中(2.2)从Uh的定义(2. 1)中获得TjPT2中Uh的散度的表达式。将散度算子应用于(2. 1) 的两边，得到 以及

的值
由于V·
'A|r2 = .Φ + 5V^,(2) +^-^+ ,V·^
(2.4)
T- = V-UaI7-，因此得到
ra J ▼ -Λ ηivt^'J
二 t/sV ·) + "3V·软、+ ,V·^ -U2V ■ φ^ -^3VUiV ·φ^(2
6 V.^+V·^v .)
可以以不同的方式得到在假设(2.2)下U6的值。第一，应用Stokes公式
udx =
■nds
PSP
由下式确定▽· 的值
▽ — "！ h 1 + ^2 Kl + 3 1^1+ 4 hi + 5 VsI
AW
这里I YiI是对应的面Yi的面积，|p|是棱锥体的体积。第二，从(2.4)确定U6
—U5V.戎2) + M3V · φψ 土 U1V. ff - V · "6=
(2.6) 在这里，重新构造在棱锥体P上的基函数{<}二 ,其满足条件於I,/ ’ i，
(i) 6 ‘
1j2J3J4J5o由于在棱锥体P的面上的基函数的法向分量的值是已知的，因此其法向分量Ul i = 1，2，3，4，5在内部面Y6上的值可以从公式(2.5)和(2.6)中得到。因此，由以下给出
基函数的明确的表达式
在乃中
Φ2
φ^-Ui6iY^ 在T2 中 Φ^+Η(62ψ 在乃中
= <
在T2中
Φ
对于具有三角形底的棱柱体单元，设Π e Qho使用“恒定划分”法来构造有限元 空间VhIp用Yj.来表示棱柱体Π的面，即，Y1是面A2A3A5A6, Y 2是面A1A3A4A6，Y3是面 A1A2A4A5，γ 4是底面A1A2A3，以及γ 5是顶面A4A5A60在图8Α中示出了这种棱柱体单元800。为了根据实施例应用“恒定划分”法，将棱柱体Π分成三个四面体1\、T2和Τ3。在 图8B-8D中示出了分割的四面体801、802、803。不失一般性，假设棱柱体Π被分成四面体 T1 = A1A2A3A6^T2 = A1A4A5A6 和 T3 = A1A2A5A6q 在图 8B-8D 中示出了分割的四面体 801、802、 803。设ni，i = 1，2，3，4，5表示到面Y i的外向单位法向向量，n6表示到四面体T1和T3 之间的界面Y 6 = A1A2A6的、从T1指向T3的单位法向向量，并且最后，n7表示到四面体T2和 T3之间的界面Y7 = A1A5A6的、从T2指向T3的单位法向向量。在每一四面体上构造经典的 最低阶Raviart-Thomas有限元空间，确定基函数，k = 1，2，3的集合，并且Π中的向
量场％被定义为如下 根据这个表达式，向量函数Uh在四面体1\、T2和T3中的每一个上是线性的，属于 空间Hdiv ( Π )，并且满足所要求的条件 对于分别在附属界面^和Y7上未知法向向量116和117的自然选择来自于以下条 件在Π上Uh的散度是恒定的，即uA =COnst在 π 上(2.8)以类似于针对棱锥体单元所讨论的方式得到未知值U6和1!7。首先，使用Stokes公 式来确定▽ ^在！!上的值 因此 这里I Y i I是对应的面γ i的面积并且I Π I是棱柱体的体积。其次，U6和U7的 值计算如下 在这里，重新构造在棱柱体π上的基函数M.t ,其满足下述条件由于在棱柱体Π的面上的基函数的法向分量的值是已知的，因此可以通过公式 (2. 10)和(2. 11)来分别确定法向分量u6(i)和u7(i)，i = 1，2，3，4，5在内部面^6和”上
的值。因此，由下面给出基函数的明确的表达式
强调“恒定划分”法的下述区别是合情合理的。将栅格单元分割成四面体不依赖于其相邻单元的分割，例如，两个相邻单元的公共的四边形面可以以不同的方式从不同的 边进行划分。例如，图9描绘了两个相邻单元901、902，其具有以不同的方式分割的邻接的 面903、904。与四面体栅格相比，这个特征简化了棱柱体栅格的产生和其上的边界值离散化 问题。对于一般的棱柱体单元，上面所描述的针对具有三角形底的棱柱体单元的方法可 以被推广到覆盖一般的棱柱体Π。每一棱柱体可以被独立地分成四面体并且使用在前面的 部分中描述的“恒定划分”法来构造有限元空间VhI π。上面的段落描述了对于通量的法向分量，单元的每一个面在面中心处具有一个自 由度。对于每一主未知量，每一单元在单元中心也被分配关联的一个自由度。例如，(多个) 主未知量可以包括周围的温度和/或流体压力。下述段落描述了将混合有限元(MFE)分析的杂化变型用在3D棱柱体栅格的单元 上，其中引入了与棱柱体的面关联的额外的自由度。这些自由度允许将一个单元的通量与 另一单元的通量隔开。这些自由度在数学文献中被称为拉格朗日(LaGrange)乘子。由于 原始的变量(例如一个棱柱体单元的单元中心的温度(或压力)和在边界上的通量)脱离 任何其他单元的温度和通量，引入额外的自由度允许简化数值问题的结构。因此，可以消除 那些未知量，这允许减少未知量的数目。根据上面提供的用于自由度的定义，可以引入MFE方法如下找到Uh e Vh,Ph e Lh 以及弋e Ah，使得对于所有ν e Vh，q G Lh以及// e Aa 其中鞍点矩阵
可以显示出，系(3. 2)具有唯一解。 很好地开发出了用于具有对称鞍点矩阵的代数系统的迭代方法。然而，起因于多 面体栅格上MFE离散化的、用于鞍点矩阵的高效预调技术仍然是关注点。对于高效预调，对 称正定矩阵是更好的对象。可以通过使用混合有限元问题(3. 1)的杂化来将系(3. 2)变换 为具有对称正定矩阵的等价系。在接下来段落中，这个方法被描述为求解问题(3.1)的优
有限元问题(3. 1)得到线性代数方程系
USD
AP=7JSn选方式。对于混合有限元分析的杂化，设Ek是Ω h中的栅格单元，并且Vh(k)和Lh(k)分别是有 限元空间Vh和Lh在Ek上的约束。此外，生成新的有限元空间Ah，其是在栅格单元之间的 界面rkl上以及在栅格单元与ΓΝ的边界部分交界处上定义的函数Xh= Ah(x)的空间。 在每一界面上，函数Xh e Ah等于常量。为了引入混合杂化有限元(MHFE)问题，不得不定义两个额外的有限元空间以及 许多双线性形式和线性泛函。两个新的有限元空间是 这里η是Qh中单元Ek的数目。注意到空间Vh(k)中的任何一个的维数最多是5， 并且每一 Lh(k)的维数等于1。对于元素U，ν e Vh，p，q G Lh以及λ，μ e Ah，引入下述双线性形式 这里 ， 曰丫广’么和&是口！^^/^的值，毛，^是入，μ e Ah的值，nkl是到 1\1的、从&指向&&<1)的单位法向向量，并且 此外，定义下述线性泛函 这里、GνΑα) >qk e《)和 μ e Λ h。借助上述定义，有限元问题(3.1)的等价混合杂化公式变为如下找到！！^^义，
Ph e 4禾口入h e A h，使得对于所有的ν e\hrqe 4禾口 μ e八h a(uh, v)+b (ν, ph)+c (ν, Xh) =Id (ν)b(uh，q) -d(ph, q) = If (q) (3. 4)
19 是具有对称正定子矩阵Mk，k = 1，. . .，η的块对角矩阵，D是对角正定或半正定矩 阵，以及Σ是对角半正定矩阵。通过(3.4)中的线性泛函来定义右手边子向量;Ji^ J、~iN 的分量。矩阵A具有非常有用的表达式W =，这里
是用于单元Ei的局部鞍点矩阵以及Ni是对应的组合矩阵(assemble matrix) 0 相应地，系(3.5)的右手边可以被写为如下 重要的是，发现可以通过将局部混合有限元离散化应用于下述问题来获得矩阵Ai 以及子向量i0厂和L，,·: Μ|. = O 在广上，如果朋,ηΓο=0这里r^k= l，...Si是单元Ei的面。(2. 6)中的矩阵Mi, B” C” Di和乙i的重要属性可以被展现在内部单元Ei上，即 dEt HdQfl = 0。设单元Ei具有Si个面，那么e Ww'是对称正定矩阵，则 这里VEi是单元Ei的体积并且
对于内部单元矩阵Σ i = 0。
设i； = (1 1…l)r e Rs'。于是4 = -CiSi。这个属性适用于来自(3· 6)的任何Ai。相关地，注意主变量二和;^，i = 1，...，η仅在单个单元内被连接。所以，可以容 易排除这些未知量 以及 由于矩阵M、B和D的结构，矩阵BH+D是对角矩阵，因此它是可逆的。使用关系 (3. 7)和(3. 8)来将系(3. 5)变换为系Sl= ξ(3.9)这里 以及 矩阵S被称为“浓缩矩阵(condensed matrix)”。这个矩阵是对称的和正定的， 除了在Neumann边界条件的情况下，S是半正定的，但是其具有简单的核恒定(kernel constant)向量。这个矩阵在本质上是全局的，其将所有的节点或单元连接在一起。可以同 时求解大的线性方程系。可以将任何迭代方法应用到求解具有该矩阵的线性方程系(3.9)。一个方法是预 处理的共轭梯度方法(PCG)，然而可以使用其他的方法。注意，在半正定性的情况下，应该 在与核正交的子空间中执行PCG。在求解系(3.9)之后，可以分别使用方程(3.8)和(3.7) 来逐个元素地局部恢复主未知量 和口。矩阵S也可以被表示为
，这里 以及次是对应的组合矩阵。(3.9)的右手边具有类似的表达式。注意，本发明的实施例可以以单个主未知量操作，例如温度或压力，及其相关的通 量。其他实施例可以以多于一个的主未知量操作。根据本发明的各种实施例，上面概述的各种过程和方法可以被组合在一个或更多 个不同的方法中，用在一个或更多个不同的系中，用在一个或更多个不同的计算机程序产 品中。例如，一个示例性方法1000可以用于形成棱柱体栅格，如图10中所示。使用正交 投影将感兴趣的地质特征和几何特征(例如尖灭边界、断层线或井的位置)投影到水平面 中(1001)。产生细的符合矩形的网格，该网格覆盖与在其上提供材料数据的细栅格相同大 小的投影域的所有特征(1002)。矩形栅格被分割成三角形(1003)。在栅格上的各种线和 点可以表示例如断层线和井位置。以非均勻方式对三角形进行粗化(1004)。期望保持某些
21地质或几何特征附近的细三角测量，但是远离这些特征具有较粗分辨率将允许更容易的分 析。这种栅格仅由三角形组成。将粗化的栅格垂直投影在所有层的所有边界表面上以形成 棱柱体栅格(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，其可以是raOM、EI3ROM或EEraOM。正如本领域中熟知的，RAM 1203 和ROM 1204保存用户和系统的数据和程序。总线1202也被耦合到输入/输出(1/0)控制卡1205、通信适配器卡1211、用户接 口卡1208和显示卡1209。1/0适配器卡1205将存储设备1206 (例如硬盘驱动器、⑶驱动 器、软盘驱动器、磁带驱动器中的一个或多个)连接到计算机系统。1/0适配器1205可以被 连接到打印设备，其允许系统打印信息(文件、照片、文档等)的纸质副本。注意，打印设备 可以是打印机(例如喷墨打印机、激光打印机等)、传真机或复印机。通信卡1211适于将计 算机系统1200耦合到网络1212，其可以是电话网络、局域网(LAN)和/或广域网(WAN)、以 太网和/或因特网中的一个或多个。用户接口卡1208将用户输入设备(例如键盘1213和 指针设备1207)耦合到计算机系统1200。用户接口卡1208也可以将声音经由扬声器输出 到用户。由CPU 1201驱动显示卡1209以控制在显示设备1210上的显示。虽然已经详细地描述了本发明及其优势，但是应明白，在不偏离由随附的权利要 求所限定的本发明的精神和范围的情况下，在此可以进行各种改变、替代和变换。此外，无 意将本申请的范围限制在说明书中所描述的过程、机器、制造、物质的组成成分、装置、方法 和步骤的特定实施例。正如本领域任何一个技术人员从本发明的公开中容易理解的，根据 本发明可以使用当前现有的或后来开发的、执行与在此所描述的对应的实施例实质上相同的功能或实现实质上相同的结果的过程、机器、制造、物质的组成成分、装置、方法或步骤。 因此，随附的权利要求意在将这种过程、机器、制造、物质的组成成分、装置、方法或步骤包 括在它们的范围内。
权利要求
一种用于在计算机上对物理区域进行建模的方法，其中所述物理区域包括多个地层，所述方法包括接收定义所述物理区域的至少一个物理特性的数据；在所述物理区域的模型的平面上提供三角形网格，其中所述网格包括多个单元；以非均匀方式对所述三角形网格进行粗化，使得较不满意的单元较大；以及将所粗化的三角形网格在与所述物理区域中的所述平面正交的方向上进行投影以形成棱柱体栅格，其中根据所述地层，所述粗化的三角形网格的每一单元被分割成子单元。
2.根据权利要求1所述的方法，其中粗化包括 将所述数据投影在平面上；以及使用所述数据来确定哪些单元是较不满意的。
3.根据权利要求2所述的方法，其中所述模型包括对所述物理区域中的物理特征进行 建模的建模特征，以及其中使用包括将优先级值分配到每一单元，其中基于每一单元是否接近建模特征和所述建模特征的 类型来确定所述值。
4.根据权利要求1所述的方法，其中提供三角形网格包括 在所述平面上提供矩形网格；以及沿着至少一个对角线分割所述矩形网格的每一单元。
5.根据权利要求1所述的方法，其中粗化包括通过消除两个相邻三角形的公共的边来合并所述两个相邻三角形。
6.根据权利要求1所述的方法，其中所述棱柱体栅格包括多个棱柱体单元、多个棱锥 体单元和多个四面体单元。
7.根据权利要求1所述的方法，其中所述方法用来对所述物理区域中的物理过程的至 少一个通量进行建模，所述方法还包括为每一子单元中的所述通量分配多个自由度；将混合有限元分析应用到所述子单元中的每一个以产生矩阵；以及求解所述矩阵以确定所述区域中所述物理过程的通量。
8.根据权利要求7所述的方法，其中分配包括，对于每一单元 为所述物理过程分配一个自由度；以及为所述单元的每一面分配另一自由度。
9.根据权利要求7所述的方法，其中应用包括 使用恒定划分法来形成所述有限元空间。
10.根据权利要求7所述的方法，其中所述物理过程是对流_扩散过程。
11.根据权利要求7所述的方法，其中所述物理过程是温度和压力中的一个并且所述 物理区域是地下的地质盆地。
12.根据权利要求7所述的方法，其中所述物理过程包含碳氢化合物材料的形成。
13.根据权利要求7所述的方法，其中所述物理过程包含碳氢化合物材料的运动。
14.根据权利要求1所述的方法，还包括从来自传感器的信息导出所述数据，所述传感器测量所述物理区域的所述至少一个物 理特性。
15.一种用于在计算机上对物理过程和该物理过程的通量进行建模的方法，所述方法 包括形成对物理区域进行建模的非结构化棱柱体栅格，其中所述物理过程在所述物理区域 内发生，并且所述棱柱体栅格包括多个单元；对于每一单元，为所述物理过程和所述通量分配多个自由度； 将混合有限元分析应用到所述单元中的每一个以产生矩阵；以及 求解所述矩阵以确定所述区域中的所述物理过程和所述通量。
16.根据权利要求15所述的方法，其中形成包括在所述物理区域的模型的平面上提供三角形网格，其中所述网格包括多个单元； 以非均勻方式对所述三角形网格进行粗化，使得较不满意的单元较大；以及 将所粗化的三角形网格在与所述物理区域中的所述平面正交的方向上进行投影，以形 成所述棱柱体栅格。
17.根据权利要求15所述的方法，其中所述棱柱体栅格包括多个棱柱体单元、多个棱 锥体单元和多个四面体单元。
18.根据权利要求15所述的方法，其中分配包括对于每一单元，为所述物理过程分配一个自由度；以及 对于每一单元，为所述单元的每一面分配另一自由度。
19.根据权利要求15所述的方法，其中应用包括 使用恒定划分法来形成所述有限元空间。
20.根据权利要求15所述的方法，还包括使用所确定的物理过程和通量来影响所述物理区域中的改变。
21.根据权利要求15所述的方法，其中所述物理过程是温度和压力中的一个，并且所 述物理区域是地下的地质盆地。
22.—种具有计算机可读介质的计算机程序产品，所述计算机可读介质具有记载在其 上的计算机程序逻辑，所述计算机程序逻辑用于对物理区域中的物理过程和该物理过程的 通量进行建模，所述计算机程序产品包括用于形成对所述物理区域进行建模的非结构化棱柱体栅格的代码； 用于将混合有限元分析应用到所述棱柱体栅格以产生矩阵的代码；以及 用于求解所述矩阵从而确定所述区域中的所述物理过程和所述通量的代码。
23.根据权利要求22所述的计算机程序产品，其中所述用于形成的代码包括用于在所述物理区域的模型的平面上提供三角形网格的代码，其中所述网格包括多个 单元；以非均勻方式对所述三角形网格进行粗化，使得较不满意的单元较大；以及 将所粗化的三角形网格在与所述物理区域中的所述平面正交的方向上进行投影，以形 成所述棱柱体栅格。
24.根据权利要求22所述的计算机程序产品，其中所述棱柱体栅格包括多个单元，并 且所述用于应用的代码包括对于每一单元，为所述物理过程分配一个自由度； 对于每一单元，为所述单元的每一面分配另一自由度；以及使用恒定划分法来形成所述有限元空间。
25.根据权利要求22所述的计算机程序产品，其中所述用于求解的代码包括 使用预处理的共轭梯度分析来求解所述矩阵。
全文摘要
本发明的实施例包含形成棱柱体栅格以及使用该棱柱体栅格和混合有限元分析来求解对流-扩散问题。可以通过在模型的平面上提供三角形网格来形成所述棱柱体栅格。接着对所述网格进行粗化以使较不满意的网格变大。接着将所述粗化的网格进行投影以形成所述棱柱体栅格。接着所述栅格的每一单元被分配多个自由度。所述栅格的混合有限元分析产生矩阵，接着求解该矩阵以得到所述对流-扩散问题的解。
文档编号G01V1/00GK101903803SQ200880121142
公开日2010年12月1日 申请日期2008年10月20日 优先权日2007年12月14日
发明者S·马里亚索威 申请人:埃克森美孚上游研究公司