作者:彭佳明; 付丽华; 张雪敏地震数据重建不动点延续算法计算复杂度近似奇异值分解
摘要:地形条件或采集成本等因素往往导致现场采集到的地震数据呈现不完整分布,从而影响后续地震数据的分析与处理,因此对原始地震数据做高精度重建显得尤为必要。不动点延续算法是一种基于核范数最小化的重建方法,但该算法需进行奇异值分解(SVD,其计算复杂度为O[mnmin(m,n)],m、n为矩阵的维度),且当矩阵维度较高时运算耗时较长;传统方法是直接利用PROPACK加速包,将计算复杂度降低为O(rmn)(r为矩阵的秩),但此加速方法依然耗时较长。为此,提出一种快速不动点延续算法,通过利用块克雷洛夫迭代近似奇异值分解算法和子空间复用技术,将SVD的计算复杂度降低为O[mcmin(m,c)](c?min(m,n),c∈R~+)为复杂度常数。仿真地震数据和实际地震数据重建结果表明,在确保一定信噪比的情况下,文中提出的快速不动点延续算法的计算效率显著高于传统加速型不动点延续算法。
注:因版权方要求,不能公开全文,如需全文,请咨询杂志社