偏微分方程的离散化方法课件.ppt

偏微分方程的离散化方法 第一页，编辑于星期四：十三点 三十二分。 一、离散化的概念 油藏是非均质的，岩石和流体性质伴随时间常常是发生变化的，建立的偏微分方程一般是非线性的，求解偏微分方程的解析解比较困难，常用数值求解。 目前工程上应用的离散化方法有：有限差分法、有限元法、边界元法、变分法等。 离散化的核心是把整体分成若干单元来处理，而每个小单元的形状是规则的，并可以认为是均质的，从而把形状不规则的非均质的问题转化为形状规则的均质的问题——非线性问题线性化。 计算过程中可以控制精度。要求的精度越高，则需要划分的单元就越多，计算工作量相应就越大，反之，单元划分得少些，计算工作量就小，但精度变差些。 微分方程离散化，主要在空间和时间两方面被离散化 第二页，编辑于星期四：十三点 三十二分。 （1）离散空间：把所研究的空间划分成某种类型的网格，大的空间转化为若干小单元组成，网格之间动态连接，通常采用矩形网格（正方体）。 （2）离散时间：把研究的时间域分成若干小的时间段，在每个时间段内，对问题求解，时间段之间有机连接。步长大小取决于所要解决的实际问题。 第三页，编辑于星期四：十三点 三十二分。 离散空间 第四页，编辑于星期四：十三点 三十二分。 t P 离散时间 第五页，编辑于星期四：十三点 三十二分。 1、网格系统 它有x，y两个自变量，在平面上用平行线分割成许多网格，如考虑时间，则。编号：x→i，y→j，t→n。为步长（对三维z→k）。 节点：网格的交点叫网格节点。取一些与边界s接近的网格节点，把他们连成折线Sh，Sh所围成的区域记为Dh，Dh内的节点为内部节点、边界上的节点为边界节点。 2、 等距网格就是指建立差分网格时，所采用的步长都是相等的，反之称为不等距网格。 第六页，编辑于星期四：十三点 三十二分。 3、网格类型 常规网格系统： （1）块中心网格：用网格小块的几何中心来表示小块的坐标 （2）点中心网格：用节点的坐标来表示小块的坐标 块中心网格和点中心网格的离散点数不同，但最终形成一样的差分方程，只有在处理边界条件时各有方便之处，块中心网格比较容易处理定流量边界，点中心网格比较容易处理定压边界。 非常规网格系统： （1）局部网格加密 （2）混合网格 （3）多边形网格 第七页，编辑于星期四：十三点 三十二分。 无效网格 有效网格 点中心网格 块中心网格 x y 第八页，编辑于星期四：十三点 三十二分。 y z x 局部网格加密 第九页，编辑于星期四：十三点 三十二分。 第十页，编辑于星期四：十三点 三十二分。 模拟区网格图（井位、边界、断层） 第十一页，编辑于星期四：十三点 三十二分。 五点法注水开发5年后XW3层含水饱和度分布图 第十二页，编辑于星期四：十三点 三十二分。 五点法注水开发20年后XW3层含水饱和度分布图 第十三页，编辑于星期四：十三点 三十二分。 r z 混合网格 第十四页，编辑于星期四：十三点 三十二分。 二、有限差分法----导数的差商逼近 P x 第十五页，编辑于星期四：十三点 三十二分。 第十六页，编辑于星期四：十三点 三十二分。 第十七页，编辑于星期四：十三点 三十二分。 第十八页，编辑于星期四：十三点 三十二分。 Δx1 Δx2 Δx 第十九页，编辑于星期四：十三点 三十二分。 三、有限差分方程的建立 第二十页，编辑于星期四：十三点 三十二分。 第二十一页，编辑于星期四：十三点 三十二分。