上机习题20101014.pptVIP

* 隐式差分格式 n n＋1 出发点在（i,n ＋1)，取关于t的一阶向后差商和关于x的二阶差商。 取 i＝1……Nx 误差 实验1：三对角系数矩阵的解法 实验1：三对角系数矩阵的解法 对于一维渗流问题的隐式差分方程组的系数矩阵为三对角矩阵，追赶法（THOMAS）就是用来求解三对角矩阵方程组的一种比较简单、 应用也极为广泛的解法。 它的基本思路是将三对角矩阵A分解成两个特定形式的三对角矩阵的乘积。 实验1：三对角系数矩阵的解法 ．　．　． ．　． ． ． 求解步骤 LU分解及求Y 求P 实验1：三对角系数矩阵的解法 实际计算时，没有必要存储 ，算法如下 (1) 令 (2) 对i=2，3，4，……，n，计算： （3）令 （4）对i=n-1，n-2，……，1计算 实验1：三对角系数矩阵的解法 求解三对角方程组： 实验1：三对角系数矩阵的解法 实验2:单相流数值模拟方法 已知一维均质等厚单相渗流，已知边界条件定压、初始压力分布，求任意时刻的压力分布。渗流微分方程为： 网格系统为均匀的块中心网格，求t=任意时刻的压力分布（隐式差分格式）。 数据准备 输入p,Q(或pwf) ttmax Yes 计算△x 开始 系数矩阵和右端项 输入tmax，△t t＝t＋△t Mi， t＝0 解压力方程组 打印t,r,p 结束 一维单相流数值模拟程序框图 实验2:单相流数值模拟方法 2 4 5 6 10 6 7 8 3 4 5 4 5 6 7 10 初始压力水平 设某油藏由5×4的排列组成，边界压力为常数，x，y步长为2，求第一个时间段的压力分布： 二维问题： 0 1 2 3 4 0 1 2 3 实验3：