适用于三维隐式叠前深度偏移中矩阵求逆的混合算法.pdf

See discussions, stats, and author profiles for this publication at: /publication/264606770 A Hybrid Matrix Inversion Method for 3-D Implicit Prestack Depth Migration Article · September 200 DOI: 10.1002/cjg2.417 CITATIONS READS 50 4 authors, including: Hong Liu Chinese Academy of Sciences 110 PUBLICATIONS 511 CITATIONS SEE PROFILE Some of the authors of this publication are also working on these related projects: seismic View project All content following this page was uploaded by Hong Liu on 19 August 2014. The user has requested enhancement of the downloaded file. 第 46 卷 第 5 期 地 球 物 理 学 报 Vol. 46 , No. 5 　　　　　　　　　　 　　　　　　　　　　 　2003 年 9 月 CHINESE JOURNAL OF GEOPHYSICS 　 Sep. , 20 0 3 适用于三维隐式叠前深度偏移中矩阵求逆的 混合算法 1 ,2 1 1 1 罗明秋 　朱国同 　刘　洪 　李幼铭 1 中国科学院地质与地球物理研究所, 北京　100029 2 University of California at Santa Cruz , USA 摘　要 　在深度偏移方法中 ,把二维隐式方法推广到三维 ,就会面对一个分块对角矩阵求逆问题. 通常 ,这种矩 阵的求逆将耗费大量计算时间 ,严重制约了三维隐式方法偏移在实际资料处理中的广泛应用. 在螺旋边界条件下 , 该矩阵 H 具有 Toeplitz 结构的正定厄密矩阵 ,其快速求逆可由谱法 LU 分解或直解法快速实现. 本文结合谱法 LU 分解和直接解法方法的优点 ,提出了一种混合算法. 文中采用谱分解方法建立起矩阵列元素的谱分解表 ,并采用直 解法的递推公式 ,可以快速给出矩阵的分解. 通过与谱法分解和直解法在分解精度和分解速度两方面的比较表明 , 本文方法与谱法相比 ,在非均匀介质中亥姆霍兹算子矩阵分解时的精度提高 10 倍 ;在计算速度方面 ,混合方法比