摘要
贾晓峰,胡天跃,王润秋.复杂介质中地震波模拟的无单元法.石油地球物理勘探,2006,41(1):43~ 48
复杂介质中地震波的传播多是通过求取单程或双程波动方程的数值解进行模拟的。波动方程
的数值解法多种多样,本文将工程和岩石力学领域中广泛应用的一种无网格算法——无单元
法应用于地震模拟。由于无单元法的节点完全独立存在,因此在计算成本上与有限元相比优
势明显;再采用滑动最小二乘的拟合方法使得无单元法具有精度高、独立变量解高次连续等
优点。本文在介绍无单元法原理的基础上,通过薄膜震动算例表明该方法与有限差分法相比
有精度优势;然后将无单元法用于叠前地震模拟,并讨论了几种吸收边界方法在地震模拟中
与无单元法的结合。模型计算的结果表明,辅以有效的吸收边界,无单元法能够获得波场完 整、走时精确的理论地震图,其精度和稳定性皆令人满意。
|