数学物理方程--有限差分法.doc

数学物理方法课程报告 题 目： 学生姓名：xxx 学 号：xxx 学 院：地球科学与技术学院 2016年 4月12日 声波有限差分法数值模拟 Xxx （地球科学与技术学院研15级 学号：xxx） 摘要：， 正文 引言 二维各向同性介质声波方程数值模拟 使用规则网格差分对二阶方程进行求解。 具体过程： 在x方向上，关于对称分布的2N个网格节点的坐标分别为， ，……，，，……，。其 中，表示节点间的最小间距；表示任意正整数。2N个网格节点所对应的函 数值已知，分别为，，……，， ……，，。利用Taylor级数展开求解 在点处的一阶导数近似值。 其中，i=1,2，…，N 将上述两式相加，省略式中的误差项，得到 (1) 将相减后得到的式子整理成矩阵形式，有 (2) 为了简化矩阵，可以记作 ， 同时，构造两个简单矩阵，辅助计算 整理的， 假设存在，使得，也可得；即为的逆，得到 。式子两边右乘向量就可得 (3) 由式（2）可得 (4) 同时，假设 （5） 将带入式（4），得 （6） 整理得 可结合式（3）和式（5），可得到矩阵计算式： （7） 当的值确定后，可根据式（7）来求解的值，从而计算出的值。 利用式（7）可以求得对称任意节点间距的一阶导数差分系数。其中，当取值为 ,则式（7）可表示为 （8） 此时，所求得的就是等节点间距的一阶导数的规则网格不同差分精度的差分系数（表1所示）。 表1 二阶导数规则网格各阶精度的差分系数 2L 2 -2.00000 1.00000 4 -2.50000 1.33333 -8.33333E-2 6 -2.72222 1.50000 -1.50000E-1 1.11111E-2 8 -2.84722 1.60000 -2.00000E-1 2.53968E-2 -1.78571E-3 10 -2.92722 1.66667 -2.38095E-1 3.96825E-2 -4.96032E-3 3.17460E-4 12 -2.98278 1.71429 -2.67857E-1 5.29101E-2 -8.92857E-3 1.03896E-3 -6.01251E-5 二维声波方程的形式可表示为： （9） 时间导数采用2阶，空间导数采用2N阶近似，即 带入式（9）中，可得到在固定网格下，差分格式为 （10） 3、模型测试： 震源选取： 正演模拟过程中采用雷克子波作为震源子波，雷克子波的表达式为 Source (it) =((1-2π fm (t-t0)2 )e-2π fm (t-t0)2 模型建立： 建立了一个两层介质模拟，其上层纵波速度为v=2000m/s，下层纵波速度为v=3000m/s。模型大小为200×200，空间采样间隔为dx=dz=10m。采用30Hz的雷克子波作为震源子波，震源位于模型(70,100)处，时间采样间隔为1ms。 结果分析： it=50 it=100 it=150 it=200 it=250 it=300 it=350 it=400 图2 不同时刻波场快照 图中可以看出，在未遇到界面前，地震波在均匀介质中的波前面一个圆。当遇到地层界面之后，在界面处发生了反射、透射和折射现象。沿测线方向的地震记录如图2所示。记录中存在两条直线状的同相轴和两条近似双曲线的同相轴。由于直达波的时距曲线是直线，因此两条直线同相轴对应直达波