《油藏数值模拟》差分方程.pdf

第三章 有限差分与差分方程 主讲人:薛永超 中国石油大学（北京） 中国石油大学 北京 油藏数值模拟研究中心 第1节 基本概念 有限差分方法是目前的油藏数值模拟研究中应用 最广泛的方法。用有限差分方法求解偏微分方程组， 将连 首先就是将偏微分方程组所表示的连续问题进行离散 续的量 化处理，即以差商来代替偏导数，然后将偏微分方程 用不连 续的点 组化为差分方程组，在求解区域内的有限个点上形成 来表示。 相应的代数方程组 （线性方程组）。 中国石油大学 北京 油藏数值模拟研究中心 第1节 基本概念 用离散化方法求解数学模型的主要步骤： 1、将渗流区域剖分成若干单元，然后将单元按一定顺序排列； 2、用网格上节点的值代替函数值； 3、在网格化基础上，从微分方程出发，建立起相邻网格节点之间的关系 式，并对其进行线性化； 4、把每个网格节点上所建立的方程合在一起，再利用定解条件，使之成为 存在唯一解的方程组； 5、解这一方程组，得到各网格节点的未知值。 中国石油大学 北京 油藏数值模拟研究中心 第2节 有限差分法 一、网格系统 求解油藏数值模拟问题时，应用有限差分方法把油藏处理为离散体 积单元（整个油气藏则是这些离散体积单元的合成 ），然后，在一定条 件下，计算每个体积单元在每个离散时间段（多个时间段 ）的变化。在 这里，通常把油藏空间体积单元叫做网格块，时间间隔叫做时间步。 网格块视为各面 都是渗透的，其内部 如经过充分搅拌溶混 的储罐。 1 2 3 4 5 6 中国石油大学 北京 油藏数值模拟研究中心 第2节 有限差分法 一、网格系统 既然数值离散要求将微分意义上连续的物理关系近似地表示成有限 个相互联系的、具有一定体积和时间单位的单元体（或节点）间的物理 关系，如何构建这些单元体就成了用数值解法来求解基本渗流方程的首 要任务。构建单元体的过程就是所谓的网格剖分或空间离散。通过网格 剖分，将当前研究的整个物理空间分为若干个小的单元体。 中国石油大学 北京 油藏数值模拟研究中心 第2节 有限差分法 一、网格系统 1、全局正交网格(Globally Orthogonal Grid) 全局正交网格单元的外边界通常和坐标轴平行，而在外边界处，为 了顺应边界的复杂形状，即可以对边界形状进行微弱的扭曲，也可以用 台阶来近似。 一维径向 网格 一维网格 中国石油大学 北京 油藏数值模拟研究中心 第2节 有限差分法 一、网格系统 1、全局正交网格(Globally Orthogonal Grid) 二维剖面网格