前进方向

FWT2D:一个用于宽口径地震资料频域全波形层析成像的大规模并行程序第1部分算法。这是本系列文章的第一篇,该系列文章描述了一个大规模并行代码,用于复杂结构成像的宽孔径地震数据的二维频域全波形反演。全波形反演方法,即基于全波方程分辨率的定量地震成像方法,计算量大。因此,设计利用并行计算设施的高效算法对于评估这些方法应用于典型案例研究和进一步改进至关重要。全波形建模需要求解大型稀疏线性方程组,该方程组由大规模并行直接解算器MUMPS执行,以实现高效的多炮模拟。利用BLAS3库,提高了多点解算阶段(前/后置换)的效率。反问题依赖于用梯度法实现的经典局部优化方法。直接求解器根据LU因子分布驱动的区域分解,返回分布在处理器上的多炮波场解。利用波场解的区域分解并行计算目标函数的梯度和对角Hessian,后者提供了一个合适的梯度尺度。该算法可以测试多尺度频率反演的不同策略,从连续的单频反演到同时的多频反演。这些不同的反转策略将在下面的论文中加以说明。并行效率和代码的可伸缩性也将被量化。