标量地震波频率-空间域有限差分法数值模拟
|
摘要
讨论了频率-空间域优化有限差分算子地震波场模拟方法。推导了频率-空间域双程波高阶有限差分算子,计算了9点有限差分系数;提出了一种节省内存的存储策略,对稀疏矩阵和稀疏矩阵LU分解后的矩阵采用了节省内存的存储策略,使内存空间的需求量大幅度降低,可以利用本方法进行具有生产规模的地震波正演模拟;推导了依赖于频率的最佳匹配层吸收边界条件,获得了很好的吸收效果。洼陷模型和Marmousi模型数值模拟试验表明,方法可行且有效,在节约计算机资源的同时,计算精度也能得到保证。
The seismic wave field simulation method for optimizing finite difference operator in frequency-space domain is discussed.The high-order finite difference operator of two-way wave in frequency- space domain is deduced and the nine-point finite-difference coeffi- cient is calculated.Strategy for saving the inner-memory is pro- posed.The strategy greatly decreases the requirement of the inner memory for the huge sparse matrixes before and after LU decompo- sition.Seismic modeling with production scale can be performed based on this work.In addition,the perfectly matched layer ab- sorbing boundary condition depending on frequency is deduced and achieves good absorbing effect.Numerical tests on a depression model and Marmousi model demonstrate that the method is feasible and effective.The computation accuracy can guaranteed with less resource cost,
The seismic wave field simulation method for optimizing finite difference operator in frequency-space domain is discussed.The high-order finite difference operator of two-way wave in frequency- space domain is deduced and the nine-point finite-difference coeffi- cient is calculated.Strategy for saving the inner-memory is pro- posed.The strategy greatly decreases the requirement of the inner memory for the huge sparse matrixes before and after LU decompo- sition.Seismic modeling with production scale can be performed based on this work.In addition,the perfectly matched layer ab- sorbing boundary condition depending on frequency is deduced and achieves good absorbing effect.Numerical tests on a depression model and Marmousi model demonstrate that the method is feasible and effective.The computation accuracy can guaranteed with less resource cost,
引文
1 Deng F,McMechan G A.True-amplitude prestack depth migration[J].Geophysics,2007,72(3): S155～S166
2 Deng F,McMechan G A.Viscoelastic true-amplitude prestack reverse-time depth migration[J].Geophysics, 2008,73(4):S143～S155
3 Zhang Y,Xu S,Bleistein N.True amplitude angle do- main common image gathers from one-way wave equa tion migrations[J].Geophysics,2007,72(1):S49～S58
4 Mulder W A,Plessix R E.A comparison between one- way and two-way wave-equation migration[J].Geo- physics,2004,69(6):1491～1504
5董良国，马在田，曹景忠．一阶弹性波方程交错网格高阶差分解法稳定性研究[J]．地球物理学报，2000，43(6)：411～419
6吴国忱，王华忠．波场模拟中的数值频散分析与校正策略[J]．地球物理学进展，2005，20(1)：58～65
7 Stekl I,Pratt R G.Accurate viscoelastic modeling by frequency-domain finite differences using rotated oper- ators[J].Geophysics,1998,63(5):1779～1794
8 Churl-Hyun J,Changsoo S,Jung H S.An optimal 9-point, finite-difference,frequency-space,2-D scalar wave ex- trapolator[J].Geophysics,1996,61(2):529～537
2 Deng F,McMechan G A.Viscoelastic true-amplitude prestack reverse-time depth migration[J].Geophysics, 2008,73(4):S143～S155
3 Zhang Y,Xu S,Bleistein N.True amplitude angle do- main common image gathers from one-way wave equa tion migrations[J].Geophysics,2007,72(1):S49～S58
4 Mulder W A,Plessix R E.A comparison between one- way and two-way wave-equation migration[J].Geo- physics,2004,69(6):1491～1504
5董良国，马在田，曹景忠．一阶弹性波方程交错网格高阶差分解法稳定性研究[J]．地球物理学报，2000，43(6)：411～419
6吴国忱，王华忠．波场模拟中的数值频散分析与校正策略[J]．地球物理学进展，2005，20(1)：58～65
7 Stekl I,Pratt R G.Accurate viscoelastic modeling by frequency-domain finite differences using rotated oper- ators[J].Geophysics,1998,63(5):1779～1794
8 Churl-Hyun J,Changsoo S,Jung H S.An optimal 9-point, finite-difference,frequency-space,2-D scalar wave ex- trapolator[J].Geophysics,1996,61(2):529～537
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心