频率—空间域有限差分正演迭代法实现

频率域正演是频率域全波形反演的基础,此外它更适用于在处理介质参数与频率有关的问题。为了更高效地实现频率域正演,本文利用有限差分算子离散频率域—空间域波动方程,频率域正演转化为大型稀疏矩阵线性方程组的求解问题,以广义最小残差法实现频率—空间域有限差分正演迭代。与频率—空间域有限差分正演直接法相比,频率—空间域有限差分正演迭代法所需的储存空间小的特点,而两者的模拟精度和耗费时间相当。随着频率域正演的发展,数据维数越来越大,计算所需的储存空间面临极大的挑战,而迭代法能很好地解决这一问题。