弹性模拟及其在AVA解释中的应用

为正确解释非均匀复杂构造上采集的地震资料,通常必须计及弹性效应。用弹性波波动方程模拟有前题的地质构造的地震响应是建立在观测地震资料与地层不均匀性之间的关系,和验证某项地震资料解释的有效手段。可以用好几种方法对描述弹性波传播的偏微分方程进行数值积分。伪潜法(傅里叶方法)就是最主要的数值积分技术,其最主要的优点是精度高和对向量和平行结构计算机的适应性,缺点是计算费用高。然而,对于给定的精度,伪谱法需要的网格数量要比有限差分法所需的小,这样就弥补了计算成本高的缺点。本文介绍一种在IBM 3090VF向量多处理机上实现的二维伪谱法弹性模型。这种算法已经适当的修改,从而能够充分地利用计算机结构,发挥它的最佳性能。