大创答辩-地震波场高精度有限差分算子题库.pptVIP

波动方程详尽描述了均匀各向同性完全弹性介质中各质点在不同时刻的位移位（或者相对体变）情况和纵波在该介质中的传播规律。 公式推导 二维均匀各向同性介质的纵波波动方程： 为了利用波动方程进行数值模拟，必须将公式进行离散化，利用差分代替微分，从而方便编制程序。下面介绍如何利用差分算子来近似计算函数的微分的。 令 两式相加整理得到二阶差分算子的公式： 公式推导 二阶差分算子的推导： * 高阶差分算子的推导： 公式推导 （1） （2） （3） （4） * 其矩阵形式为： 根据待定系数法： 系数=1， 的系数=0，则可得如下方程组： 或： 将求得系数c1和c2带回（5）式，即可得到的差分格式： 对1、2式3、4式分别相加并且整理得： 公式推导 公式推导 利用相同的方法，我们逐渐推出了六阶、八阶、十阶、十二阶的差分算子公式，然后带入波动方程，得出递推公式： 公式推导 虽然，通过使用高阶的差分算子，数值频散现象已经明显改善，但是，仍然存在比较明显的数值频散，为了进一步减少数值频散，我们又对算子进行了一定的优化： 传统的二阶空间分化函数f基本上是一个在x=0被截断的泰勒级数，如公式 其中偶数N按顺序排列， 是沿x均匀的时间间隔， 是一个二项式系数公式定义的常系数 如公式 我们可以通过使用优化公式的常数系数来减小数值频散，多数以前的工作希望设计一个适当的窗 得到优化差分算子，这可概况为以下形式，如公式： 是我们已使用的优化的最终形式，我们下一步是确定bn的优化方案在本文中，我们通过检查波数域和分析波数中优化有限差分算子的绝对误差来直接优化系数。我们应用模拟退火算法优化下面的目标函数。 公式推导 事实上，我们可以直接搜索优化的系数的最终形式。也就是说，窗 与公式（3）的其他部分相结合。如公式 程序编写 for(it=0;itNt;it++) { 递推公式 加载震源 数据替换 震源子波生成 记录单炮数据 生成波场快照 } 记录文件 结果分析 利用不同阶数的高阶有限差分算法进行了数值模拟，得到相应的波场快照。图中可以看出，在阶数较低时，出现很多同相轴，说明数值频散现象严重；随着算子阶数的增加，数值频散现象逐渐减弱。对于高阶差分算子来说，算子阶数越高，压制数值频散效果越好，精度越高。 二阶 四阶 十阶 十二阶 结果分析 我们求出bn的得数后，替换程序中原有的Cn，得到优化之后的程序，同 时，得到相应的波场快照，并进行对比： 优化前 优化后 优化前 优化后 结果分析 图中可以看出，对于相同阶数的算子，优化有限差分模拟的波场快照对比高阶差分算子的结果，优化算子压制数值频散效果更好，精度更高，说明优化算子优于高阶差分算子。我们取出两种算子模拟时，得到模型中间测线的波形曲线图如图所示，通过两个图对比可以更加清楚地看到波尾附近的数值频散明显减小。 结果分析 我们建立了一个均匀介质模拟，其纵波速度为v=3000m/s。模型大小为200×200，空间采样间隔为dx=dz=10m。采用30Hz的雷克子波作为震源子波，震源位于模型中心处，时间采样间隔为1ms。利用优化有限差分算子进行数值模拟，波场快照如图3-6所示。图中可以看出，波场快照中的同相轴是圆形的，说明在均匀各向同性介质中，点源激发的波前面是一个圆，这与理论也是吻合的。并且随着时间的增大，波前面的面积逐渐增大，说明地震波从震源中心向外传播。 均匀介质模型 结果分析 (a)50ms （b）100ms (c)150ms (d)200ms (e)250ms (f)300ms 结果分析 结果分析 层状介质模型 　　我们建立了一个两层介质模拟，其上层纵波速度为v=2000m/s，