应用快速不动点延续算法的地震数据重建

• 作者: 付丽华 张雪敏 彭佳明
• 来源期刊: 石油地球物理勘探
• 关键词: 地震数据重建 不动点延续算法 计算复杂度 近似奇异值分解
• 年，卷（期）: 2019,（6）
• 所属期刊栏目: 处理技术
• 页码范围: 1195-1205
• 页数: 11页
• 分类号: P631
• 语种: 中文
• DOI: 10.13810/j.cnki.issn.1000-7210.2019.06.003

摘要

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