二阶声波正演c语言程序,声波波动方程正演模拟程序总结 - 图文

此处，为了成图完整，我用的是t2，而不是 t1，也就是把雷克子波向右移动了一段距离，使主要部分都显示出来。(频率采用的是30hz)

从图中可以看出程序是正确的，符合理论上雷克子波的波形。 第二部分：主程序，编写声波正演模拟算子。

首先定义了各种变量，然后指定震源位置，选择权系数，给速度赋值，然后是差分算子的编写，这是主要部分，最后再进行时间转换，即把n-1时刻的值给n时刻，把n时刻的值给n+1时刻。 此处，我编写的是均匀介质声波方程规则网格的正演模拟程序，时间导数采用二阶中心差分、空间导数为2N阶差分精度，网格大小为 200*200，总时间为400。

第三部分：这一部分就是记录文件。

首先记录Un文件，然后记录record文件。

模型构建与试算：

1、我首先建立了一个均匀介质模型，首先利用不同时间，进行了数值模拟，得到波场快照如图所示：

100ms

200ms 300ms

此处，纵波速度为v=3000m/s。模型大小为200×200，空间采样间隔为dx=dz=10m。采用30Hz的雷克子波作为震源子波，时间采样间隔为1ms，图中可以看出，波场快照中的同相轴是圆形的，说明在均匀各向同性介质中，点源激发的波前面是一个圆，这与理论也是吻合的。并且随着时间的增大，波前面的面积逐渐增大，说明地震波从震源中心向外传播。

2、我在建立的均匀模型的基础上，改变差分算子的精度ÿ