油藏数值模拟原理 .pptVIP

;黑油模型简介;黑油模型简介;1 、黑油模型基本假设;2.单相渗流基本微分方程;渗流速度是一个空间向量，用分量表示为： ── 渗流速度矢量； ──压力； ──势函数 ── 绝对渗透率； ──重力加速度常数； ──粘度； ── 垂向坐标； ── 密度； ;Reservoir Simulation III;2.单相渗流基本微分方程;9;2.单相渗流基本微分方程;根据物质守恒定律： 将上式两边同除以 ， ;并令 ， ， 把源汇项（井）加入到方程中，得到 单相渗流连续性方程 ： ;???Hamilton算子 表示，并将达西公式代入，得 到，单相渗流微分方程（流动方程）： 生产井 注入井 关 井 ;这就是地下流体所遵循的渗流规律。 总之，渗流基本微分方程包含三部分内容: 流动项、源汇项以及累积项 渗流基本微分方程由两个定律控制： 达西定律与物质守恒定律 压力场和饱和度场的分布是流体流动的结果，流 动规律由达西定律来确定，流量由物质守恒来确 定，这两个规律是数值模拟的关键所在。;将单相渗流达西定律推广到多相渗流中，得到 三相流动达西公式： 式中，w，o，g——水、油、气诸相下标。 ── 相对渗透率； ; 对油、水、气三相，每一相运用物质守恒定律，可以建立以下连续性方程 ; 把达西公式代入连续性方程，得到 三维三相渗流基本微分方程 ;3. 三维三相渗流基本微分方程;上述T的物理意义是反映相邻网格间的流动性，又叫传导率（transmissibility) ;3. 三维三相渗流基本微分方程;3. 三维三相渗流基本微分方程;Stone公式： 式中， —— 束缚水饱和度 ——油水两相流动时油相相对渗透率 ——油气两相流动时油相相对渗透率 ;3. 三维三相渗流基本微分方程;3. 三维三相渗流基本微分方程;产量项（作源汇项处理）（动态数据） ;3. 三维三相渗流基本微分方程;4 .初边值条件;4 .初边值条件;4 .初边值条件; 5. 网格系统;在一般情况下流动方程求不出解析解，只能用离散化方法求数值解。 离散化方法，首先要把求解区域按一定的网格系统进行剖分。网格剖分方法有很多，像直角坐标网格、径向网格、角点网格、PEBI网格等。目前，商用软件中大都采用矩形网格系统。离散化方法主要有：有限差分（Finite Difference)、有限元(Finite Element)以及流线法(Stream Line)。比较成熟的是有限差分方??。有限差分方法大多空间变量采用网格块中心（Block Center)差分, 时间变量按向前差分方法。;5.网格系统:矩形网格;33;5.网格系统:径向网格系统;5.网格系统:角点网格;块中心网格与角点网格的比较;5.网格系统:块中心网格的数据结构;5.网格系统:块中心网格的数据结构;角点网格的数据结构;块中心网格与角点网格的比较;;非结构网格技术;非结构网格技术;非结构网格技术(Unstructured Grids);5.网格系统: PEBI网格;unstructured PEBI grids ;New Faulted Grid;非结构网格技术 Channels can be modelled ;2. 非结构网格技术(complex) faults;1. 概述 局部网格加密(Local Grid Refinement，简称LGR)技术主要分： 静态局部网格加密(Static Local Grid Refinement)局部加密区域在整个计算过程中固定不变的，常用于断层、尖灭、裂缝以及井附近的加密。该技术已经成熟，在多家软件公司的软件产品中都有这一功能。 动态局部网格加密(Dynamic Local Grid Refinement)局部加密区域在整个计算过程中随流体的运动位置而变化，主要用于追踪驱替前缘。该方法在研究之中，最早是Heinemann将动态局部网格用于水驱油藏数值模拟 。;Productive Area 1;局部网格加密技术;随驱替前缘移动的动态局部网格加密系统;1.3 动态局部网格加密技术(注化学剂时的So）;Single Well Block;开窗技术(Windowing Technique) ;Window Grid;开窗技术;开窗技术; 差分方程 ; ; 该差分方程组无论是左端还是右端都是非线性的，是一组非线性代数方程组。非线性最强的参数是相对渗透率和毛管力。这两个参数在时间方面