储能模量与频率关系_伪谱法频散关系与稳定性准则

参考文献^[1][2]

---

声学波动方程解法

当介质密度和地震波速度均为空间变量时，声学波动方程写作：

其中

代表压力，

代表密度，

代表波速，

代表震源项，等于体力的散度除以密度。在方程（1）中以及之后，变量上一点表示时间一阶微分，两点表示时间二阶微分。

通过对时间和空间进行离散化，用傅里叶法求解（1），得到方程（1）的一个离散近似：

其中

和

分别代表在时间

时、

处的压力和震源的值。

代表方程（1）等式左端的数值近似。

方程（2）包含一个显式二阶时间微分格式，类似于在有限差分和有限元代码中经常使用的格式。傅里叶法的出发点是计算空间导数的方式。

方程（2）中的项

是在两个不同的方向计算的，一个是对于包含在

方向求导的项，另一个是对于包含在

方向求导的项。在对

求导时，在平行于

方向的每一个网格点（称为

线）上计算

。计算沿每条

线的

，首先对

进行空间快速傅里叶变换，随后，结果与空间波数向量