油藏数值模拟原理_黑油模型培训PPT
本帖最后由 自由飞翔 于 2011-11-2 20:07 编辑黑油模型简介
黑油模型是指油质较重性质的油藏类型。
黑油模型是最完善、最成熟，也是应用最为广泛的模型。
是油藏数值模拟的基础，其它模型大都是黑油模型的扩展。
目的：软件这样一个“黑箱”是如何操作的，以便指导我们更好地去应用软件。
黑油模型简介
黑油模型的基本假设
单相渗流黑油模型基本微分方程
三维三相黑油模型基本微分方程
初边值条件
5. 网格系统
6. 黑油模型的差分方程
7. 差分方程的线性化
8. 线性代数方程组的求解
9. 黑油模型的主要数据流
1 、黑油模型基本假设
（1） 油藏中的渗流是等温渗流。
（2） 油藏中最多只有油、气、水三相，每一相均遵守达西定律。
（3）油藏烃类只含有油、气两个组分。在油藏状态下，油气两组分可能形成油气两相，油组分完全存在于油相内，气组分则可以以自由气的方式存在于气相中，也可以以溶解气的方式存在于油相中，所以地层内油相为油组分和气组分的某种组合。在常规油田中，一般不考虑油组分向气组分挥发的现象。
（4）油藏中气体的溶解和逸出是瞬间完成的，即认为油藏中油气两相瞬时达到相平衡状态。
（5）油水之间不互溶；天然气也假定不溶于水。
2.单相渗流基本微分方程
2. 单相渗流基本微分方程
渗流速度是一个空间向量，用分量表示为：
── 渗流速度矢量；──压力； ──势函数
── 绝对渗透率； ──重力加速度常数；
──粘度； ── 垂向坐标； ── 密度；
Reservoir Simulation III
2.单相渗流基本微分方程
根据物质守恒原理，在单位时间内
流入单元内的流体质量-流出单元的流体质量
=单元内流体质量的变化
取渗流场中一个微小六面体体积单元来研究
2.单相渗流基本微分方程
该六面体单元的中心点坐标为（x，y，z），其
每一侧面的质量流速均以其侧面的中心点的质量
来代替。
方向流入质量=
方向流出质量=
方向流入质量=
方向流出质量=
方向流入质量=
方向流出质量=
单元体内流体质量的变化=
2.单相渗流基本微分方程
根据物质守恒定律：
将上式两边同除以 ，
2.单相渗流基本微分方程
并令 ， ，
把源汇项（井）加入到方程中，得到
单相渗流连续性方程 ：
2.单相渗流基本微分方程
用Hamilton算子 表示，并将达西公式代入，得
到，单相渗流微分方程（流动方程）：
生产井
注入井
关 井
2.单相渗流基本微分方程
这就是地下流体所遵循的渗流规律。
总之，渗流基本微分方程包含三部分内容:
流动项、源汇项以及累积项
渗流基本微分方程由两个定律控制：
达西定律与物质守恒定律
压力场和饱和度场的分布是流体流动的结果，流
动规律由达西定律来确定，流量由物质守恒来确
定，这两个规律是数值模拟的关键所在。
3. 三维三相渗流基本微分方程
将单相渗流达西定律推广到多相渗流中，得到
三相流动达西公式：
式中，w，o，g——水、油、气诸相下标。
── 相对渗透率；
3. 三维三相渗流基本微分方程
对油、水、气三相，每一相运用物质守恒定律，可以建立以下连续性方程
3. 三维三相渗流基本微分方程
把达西公式代入连续性方程，得到
三维三相渗流基本微分方程
3. 三维三相渗流基本微分方程
令
3. 三维三相渗流基本微分方程
上述T的物理意义是反映相邻网格间的流动性，又叫传导率（transmissibility)
3. 三维三相渗流基本微分方程
3. 三维三相渗流基本微分方程
还要补充一些方程
饱和度 ：
油水毛管力：
油气毛管力：
相对渗透率和毛管力是饱和度函数：
3. 三维三相渗流基本微分方程
Stone公式：
式中，—— 束缚水饱和度
——油水两相流动时油相相对渗透率
——油气两相流动时油相相对渗透率
3. 三维三相渗流基本微分方程
流体PVT性质
水相：
气相： （Z)
油相：随压力变化，相态会发生变化。
三相状态：
两相状态：
岩石的高压物性
相渗曲线以及流体PVT性质，以表格形式给出，
称为静态数据。
3. 三维三相渗流基本微分方程
实际求解变量
初始化输入数据ljfinit.txt
3. 三维三相渗流基本微分方程
产量项（作源汇项处理）（动态数据）
3. 三维三相渗流基本微分方程
其中，
—— 在井基准面深度 上的井底压力。
—— 节点等效泄油半径。
—— 井筒半径。
△X，△Y，△Z—— 网格节点x，y和z方向步长
PI叫采油指数， 通常写成KH,可以直接输入
4 .初边值条件
4 .初边值条件
外边界
封闭边界:尖灭、断层、圈闭
流动边界：边、底水
内边界，即井的生产条件
定产边界条件（rate constraint):给定产液量、产油量、产气量、注水量、注气量
定压边界条件(pressure constraint):给定井底流压
4 .初边值条件
5.网格系统
5.网格系统
在一般情况下流动方程求不出解析解，只能用离散化方法求数值解。
离散化方法，首先要把求解区域按一定的网格系统进行剖分。网格剖分方法有很多，像直角坐标网格、径向网格、角点网格、PEBI网格等。目前，商用软件中大都采用矩形网格系统。离散化方法主要有：有限差分（Finite Difference)、有限元(Finite Element)以及流线法(Stream Line)。比较成熟的是有限差分方法。有限差分方法大多空间变量采用网格块中心（Block Center)差分, 时间变量按向前差分方法。
5.网格系统:矩形网格
5.网格系统:径向网格系统
5.网格系统:角点网格
块中心网格与角点网格的比较
5.网格系统:块中心网格的数据结构
模型被分为20×5×10个网格区域。
网格从（1，1，1）开始向下倾斜。
各网格块在X方向上是300ft，在Y方向上是1000ft。
层厚度从上至下分别是32，22，20，4，32，4，26，26，4，28ft
5.网格系统:块中心网格的数据结构
--前20个TOPS定义网格 (1, 1, 1) 到 (20, 1, 1)
6855.000 6865.000 6875.000 6885.000 6895.000
6905.000 6915.000 6925.000 6935.000 6945.000
7005.000 7015.000 7025.000 7035.000 7045.000
7055.000 7065.000 7075.000 7085.000 7095.000
--下面20个TOPS定义网格 (1, 2, 1) 到 (20, 2, 1)
6930.000 6940.000 6950.000 6960.000 6970.000
6980.000 6990.000 7000.000 7010.000 7020.000
7080.000 7090.000 7100.000 7110.000 7120.000
7130.000 7140.000 7150.000 7160.000 7170.000
--下面20个TOPS定义网格 (1, 3, 1) 到 (20, 3, 1)
7030.000 7040.000 7050.000 7060.000 7070.000
7080.000 7090.000 7100.000 7110.000 7120.000
7180.000 7190.000 7200.000 7210.000 7220.000
7230.000 7240.000 7250.000 7260.000 7270.000
--下面20个TOPS定义网格 (1, 4, 1) 到 (20, 4, 1)
7130.000 7140.000 7150.000 7160.000 7170.000
7180.000 7190.000 7200.000 7210.000 7220.000
7280.000 7290.000 7300.000 7310.000 7320.000
7330.000 7340.000 7350.000 7360.000 7370.000
--下面20个TOPS定义网格 (1, 5, 1) 到 (20, 5, 1)
7205.000 7215.000 7225.000 7235.000 7245.000
7255.000 7265.000 7275.000 7285.000 7295.000
7355.000 7365.000 7375.000 7385.000 7395.000
7405.000 7415.000 7425.000 7435.000 7445.000
/ 这完成了第一层TOPS的定义
角点网格的数据结构
--COORD
--此关键字定义坐标线
--X1Y1 Z1 X2 Y2 Z2
0.0. 6825.000 0. 0. 7023.000
300.0000 0. 6835.000 300.0000 0. 7033.000
600.00000. 6845.000 600.0000 0. 7043.000
900.00000. 6855.000 900.0000 0. 7053.000
1200.0000. 6865.000 1200.000 0. 7063.000
..........
/对于20 * 10 的网格模型, 需要21*11 条坐标线，即 231条。-每一条坐标线都用6个数来定义，所以关键字COORD包含有1386个数。
ZCORN
--这个关键字定义网格拐点的深度，数据按照X (或 R) 方向循环最快，然后是Y (或 THETA) 方向，最后是Z方向的顺序排列。
6825.000 6835.000 6835.000 6845.000 6845.000
6855.000 6855.000 6865.000 6865.000 6875.000
6875.000 6885.000 6885.000 6895.000 6895.000
6905.000 6905.000 6915.000 6915.000 6925.000
6975.000 6985.000 6985.000 6995.000 6995.000
7005.000 7005.000 7015.000 7015.000 7025.000
7025.000 7035.000 7035.000 7045.000 7045.000
7055.000 7055.000 7065.000 7065.000 7075.000
.......................
/ 对于1000个网格，需要 8000 个ZCORN值
--这些只是角点网格模型的数据体中很少的一部分
块中心网格与角点网格的比较
非结构网格技术
传统的结构网格--即规则矩形剖分的网格,具有如下不足:
存在较严重的网格取向效应
不能保证每口井都在网格中心
不能很好地描述油藏形态
粗细网格间的过渡衔接性差
研究新的网格技术势在必行.
非结构网格技术
非结构网格技术(Unstructured Grids)
常用的非结构网格
各网格中心之间连线构成三角网
PEBI网格（又叫垂直平分网格Perpendicular Bisectors):各三角形边的垂直平分线构成网格边界
CVFE网格--控制体有限元网格(Control volume finite element):各三角形边的中点与重心的连线构成网格边界
待解决的问题
非结构网格条件下代数方程组的快速求解
5.网格系统: PEBI网格
unstructured PEBI grids
Advanced Fault Modeling
非结构网格技术 Channels can be modelled
2. 非结构网格技术 (complex) faults
好东西~~~谢谢分享~下了 感谢楼主分享 这就是传说中的黑油模型啊！ 需要好好看看
页:
[1]