变网格有限差分弹性波方程数值模拟方法.pdfVIP

2007 12 42 6 !! 变网格有限差分弹性波方程数值模拟方法 * ∀ # ∀ ∀ ∃ 朱生旺 曲寿利 魏修成 刘春园 ( ∀ ; # ( ) ; ∃ ( ) ) , , , . . , 2007, 42( 6) : 634~ 639 研究复杂介质中地震波传播规律及地震响应特征, 要在细小网格剖分下进行弹性波方程数值模拟计 算, 而小网格下的数值模拟计算将带来巨大的计算量问题, 采用变网格计算是减少计算量的有效途径本文给 出一种变网格差分计算的实现方法, 在局部复杂介质区域采用细网格计算, 其余区域采用粗网格计算, 在两种 网格的过渡区通过改变差分算子和波场插值实现波传播的过渡衔接理论分析和数值模拟结果表明, 变网格 时在粗细网格的过渡区不会对地震波传播模拟带来影响, 从而达到了既减少计算量又保证计算精度的目的 弹性波方程 数值模拟 变网格差分 频散 ( ) , 1 , , , , , , , , Jas [ 3] , tram C , [2] , , , , , , , , , [ 1, 2] , , , , n , n2 , , 2 1m % 1m , , , , * 31 , 100083 2007 1 16 ( 973) ( 2006CB202402) 42 6 : 635 , f ( x 0 + q1 x / 2) - f ( x 0 - q 1 x / 2) 1 f ( x 0 + q2 x / 2) - f ( x 0 - q2 x / 2) [ 1, 2] = x ( 2) ,