- 1、本文档共19页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
差分方程组的建立
主讲教师:李亚军
《油藏数值模拟基础》第三章
1
基本渗流微分方程的一般式:
非线性偏微分方程:
分析地层内流体压力、饱和度随时间的变化状况
研究流体在油层中的渗透率特征
解析求解 (Analytical Solution)
数值求解 (Numerical Solution)
回顾复习
有限差分法(FDM)
有限元法(FEM)
边界元法(BEM)
连续问题
离散化
2
第一节 基本有限差分(离散化概念 差商)
第二节 差分方程组的建立
第三节 网格排列方式及其系数矩阵
第四节 差分方程组的稳定性分析
第五节 边界条件的处理
第三章 差分方程组的建立
3
第一节 基本有限差分(离散化概念 差商)
第二节 差分方程组的建立
第三节 网格排列方式及其系数矩阵
第四节 差分方程组的稳定性分析
第五节 边界条件的处理
第三章 差分方程组的建立
4
连续问题
偏微分方程
偏导数
解析解
第一节 基本有限差分
离散问题
差分方程
差商
数值解
数学模型
数值模型
有限差分法
连续问题离散化
数值求解:将实际油藏的连续求解区域划分成有限个离散点,通过求解离散点上的压力、饱和度来近似实际油藏问题的解。
离散化:将连续的问题离散成有限个离散点并可进行数值计算的方法。
5
几何形状:边界不规则、非等厚、构造起伏等
储层参数:渗透率非均质与各向异性、孔隙度非均匀等
流体分布:流体性质差异、初始油气水分布等
开发过程:井网不规则、注采强度差异、非等温渗流等
离散化的背景
一、离散化的概念
离散化分类
空间离散
时间离散
6
1、空间离散
空间离散的影响因素:
油藏的空间形状(一维、两维、三维)
油藏的空间大小(网格划分的大小)
研究问题的需要(精度、油藏非均质性、局部加密)
计算机的速度和容量
定义:将连续的油藏空间研究区域按一定的网格系统剖分成有限个单元或网格。
一、离散化的概念
7
(1)一维问题
1、空间离散
一、离散化的概念
单管模型实验
一维油藏模型空间离散网格
一维垂向重力驱替
Dz
8
(2)二维问题
1、空间离散
一、离散化的概念
①平面问题(x-y)
等距离网格
不等距网格
求解下列问题:
(a) 计算水驱或气驱的平面采收率
(b) 计算平面上剩余油的分布情况
(c) 对比各种开发方案及调整方案
9
油、水井附近网格要密些,油藏边缘可稀些
井之间一定要间隔一个网格以上
研究开发初期的可行性方案时,网格可粗些
当历史拟合、动态预测时网格要密些
(2)二维问题
1、空间离散
一、离散化的概念
①平面问题(x-y)
10
(2)二维问题
1、空间离散
一、离散化的概念
②剖面问题(x-z)
(a) 层状油藏的垂向波及效率
(b) 剖面上的水驱效率
(c) 模拟气驱的过程、重力的影响等
11
研究对象:厚油层,平面、垂向上不均质,可模拟油田规模的水驱、气驱
一般选三维笛卡尔网格系统,在x,y,z三个方向划分网格
(3)三维问题
1、空间离散
一、离散化的概念
12
定义:在所研究的时间范围内把时间离散成一定数量的 时间段(时间步长△t),在时间步长内进行求解。
影响时间步长大小的因素:
所研究的油藏问题所经历的时间长短
开发时间长, △t可大些,反之可小些
所研究油藏问题的性质
初期可行性方案, △t可大些
历史拟合时, △t可取不等距的
解的稳定性及准确性要求
2、时间离散
一、离散化的概念
13
块中心网格
点中心网格
二、网格系统
离散化
空间离散 → 单元(网格)
时间离散 → 时间步长
笛卡尔坐标系下的记号
14
1、块中心网格系统
i
j
I=5
J=4
实质:将研究区域剖分成小块(各块的大小可以不同),以块的几何中心作为节点。
1 2 3 4 5
4
3
2
1
二、网格系统
适用于定流量边界条件
15
① 节点位置
本节点(i,j)
邻节点:左(i-1,j)
右(i+1,j)
上(i,j+1)
下(i,j-1)
节点(i,j)的边界
左 i-1/2
右 i+1/2
上 j+1/2
下 j-1/2
二、网格系统
1、块中心网格系统——网格位置表示方法
16
二、网格系统
1、块中心网格系统——网格位置表示方法
② 坐标位置
本节点(xi,yj)
邻节点:左(xi-1,yj)
右(xi+1,yj)
上(xi,yj+1)
下(xi,yj-1)
节点(xi,yj)的边界
文档评论(0)