%0 Journal Article
%T LU DECOMPOSITION WITH SPECTRAL FACTORIZATION IN SEISMIC IMAGING<br>用于波场成像的谱法LU分解
%A LUO MINGQIU
%A <br>罗明秋
%J 地球物理学报
%D 2003
%I 
%X The inversion of Helmholtz matrix, a finite-difference representation of Helmholtz operator, plays a critical role in the numerical methods of seismic modeling and 3D depth migration. The Helmholtz matrix has a Toeplitz structure with a helical boundary condition. It is easy to be inversed by an efficient LU decomposition method, a spectral factorization method (SF method). In this paper we discuss the efficiency of SF method in LU decomposition, the errors of SF method in decomposing the matrix H of several different velocity models, errors distribution and their effects on seismic imaging. The result in this paper indicates that for constant velocity model all the arrays of matrix H have same nonzero values, and the errors brought by SF method will not affect the computation of wave fields if an absorb boundary condition is used; for non-constant velocity model, each array in matrix H has different nonzero values, and the error increases with variability degree of the velocity model, which are fatal to wave fields calculation. Therefore, once a SF method is used to calculate the wave field in seismic modeling and imaging, the errors given by this method should be considered to ensure a correct wave fields, which will be discussed in another paper.
%K Seismic imaging
%K LU decomposition
%K Spectral factorization<br>波场成像
%K Lu分解
%K 谱因式分解
%K 模拟
%U http://www.alljournals.cn/get_abstract_url.aspx?pcid=E62459D214FD64A3C8082E4ED1ABABED5711027BBBDDD35B&cid=1E44AE713D8A6DE0&jid=14DC41C59CBF6770055A7D610D53AE46&aid=98524FC14367D098&yid=D43C4A19B2EE3C0A&vid=D997634CFE9B6321&iid=38B194292C032A66&sid=65C08888CCE4801E&eid=CB3428B1EFB1C133&journal_id=0001-5733&journal_name=地球物理学报&referenced_num=5&reference_num=19