摘要
讨论了频率-空间域优化有限差分算子地震波场模拟方法。推导了频率-空间域双程波高阶有限差分算子,计算了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,
出处
《石油物探》
EI
CSCD
北大核心
2009年第1期20-26,共7页
Geophysical Prospecting For Petroleum
基金
国家高技术研究发展计划(863)项目(2006AA092323)
国家重点基础研究发展计划(973)项目(2006CB202402)共同资助
关键词
频率-空间域
有限差分法
标量波方程
最佳匹配层吸收边界条件
地震波正演模拟
frequency-space domain
finite difference method
scalar wave equation
perfectly matched layer absorbing boundary conditions
seismic forward modeling