黏弹介质波动方程有限差分解的稳定性研究  

Study on the stability of finite difference solution of visco-elastic wave equations

作　　者：孙成禹[1] 肖云飞[1] 印兴耀[1] 彭洪超 

机构地区：[1]中国石油大学(华东)地球资源与信息学院,中国山东青岛266555 [2]中国石油集团东方地球物理公司,中国河北涿州072751

出　　处：《地震学报》2010年第2期147-156,共10页Acta Seismologica Sinica

基　　金：国家重点基础研究发展计划(2007CB209605);国家高技术研究发展计划(2007AA060502)资助

摘　　要：稳定性问题是地震波数值模拟的一个重要问题.基于地震波传播理论,从黏弹介质本构方程出发,对矩形网格下不同黏弹模型波动方程有限差分解的稳定性进行了理论分析,导出了Kelvin-Voigt黏弹模型和Maxwell黏弹模型在任意空间差分精度下稳定性条件的表达式;给出了品质因子Q≥5时的简化式,并通过数值算例验证了理论研究结论的正确性;总结了地震波速度、频率、空间网格大小、差分系数以及品质因子与稳定性条件的关系;通过误差分析给出了近似公式的使用条件.The stability problem is a very important aspect in seismic wave numerical simulation.Based on the theory of seismic waves and constitutive relations of visco-elastic models,the stability problem of finite difference scheme with rectangular grids for two visco-elastic models is analyzed.Expressions of stability condition under arbitrary spatial accuracy for Kelvin-Voigt and Maxwell models are derived.With an approximation of medium quality factor Q≥5,a simplified expression is developed and some digital models are given to show the validity of theoretical results.Influence of velocity,frequency,size of grid,differential coefficients,as well as quality factor,on the stability condition is summarized.Finally,the working condition of the simplified stability expression is given by means of its error analysis.

关 键 词：有限差分 稳定性 黏弹模型 波动方程 品质因子 

分 类 号：P315.31[天文地球—地震学]