基于加权平均导数的频率-空间域正演模拟及GPU实现

摘要

传统基于旋转坐标系的频率-空间域正演模拟方法仅适用于方形网格,而实际生产中矩形网格广泛存在,本文提出一种适用性广的正演差分算子,不仅适用于方形网格而且适用于矩形网格.通过综合运用平均导数法、加速项加权平均、模拟退火法压制频散和减少单个波长所需网格点数,从而提高算法精度和减少计算量.在该方法的基础上采用不完全LU分解作为求解Helmholtz方程的预条件,并利用图形处理器加速计算速度,很大程度上提高了频率域正演的效率.