摘要
有限差分法地震波场外推通常基于笛卡尔坐标系,当地下介质速度差异很大时,会导致波场的局部过采样现象,增加计算量。为此,将伪深度域的思想引入到弹性波有限差分数值模拟。首先通过坐标转换的链式法则,推导了深度域到伪深度域的坐标映射关系,获得了伪深度域下的弹性波方程(PDD-EWE),在保证精度的前提下,在伪深度域进行波场外推能减少大量内存占用;引入变长度差分算子,即在一定的误差范围内,根据速度的不同设计了不同长度的差分算子。模型测试结果表明,与常规深度域方法相比,伪深度域变长度差分算子弹性波正演模拟和成像算法在保证精度的前提下,既能减少内存占用又能提高计算效率。
- 单位