二维弹性及粘弹性TTI介质中地震波场数值模拟：四种不同网格高阶有限差分算法研究.pdfVIP

第 28卷 第4期 地 球 物 理 学 进 展 Vo1．28，No．4 2013年 8月(页码 ：1817—1827) PROGRESS IN GEOPHYSICS Aug．，2013 孙耀充，张延腾，白超英．二维弹性及粘弹性TTI介质中地震波场数值模拟：四种不 同网格高阶有限差分算法研究．地球物理 学进展，2013，28(4)：1817—1827，doi：10．6O38／pg20l3O423． SUN Yao-chong，ZHANG Yan-teng，BAIChao-ying． Seismicwavefield simulation in 2D elasticandviscoelasticmedia： comparisonbetweenfourdifferentkindsoffinite-differencegrids．ProgressinGeophys．(inChinese)，2013，28(4)：1817—1827， doi：10．6O38／pg2O130423． 二维弹性及粘弹性 TTI介质 中地震波场数值模拟 ： 四种不 同网格高阶有限差分算法研究 孙耀充 。，张延腾 ， 白超英 。 (1．长安大学地质工程与测绘学院地球物理系，西安 710054； 2．中国科学技术大学地球和空间科学学院，合肥 230026； 3．长安大学计算地球物理研究所，西安 710054) 摘 要 本文利用交错网格、辅助网格、旋转交错网格、同位网格有限差分方法分别模拟 了二维弹性 1vrI介质和二 维黏弹性TTI介质中的地震波传播．在稳定性条件 内，选用不同的网格间距及时间间隔，通过波场快照、合成理论 地震图较为系统分析对比了这四种不同网格有限差分数值模拟在计算精度、CPU 时间、相移、频散、以及保幅方面 的优缺点．数值模拟结果表明：1)这四种不同网格有限差分算法都是很好的波场数值模拟算法；2)就CPU计算 时间而言，旋转交错网格有限差分算法的计算效率最高；3)从计算精度来看，同位 网格有限差分的计算精度最 高；4)从振幅保护方面来看，四种网格的保护振幅的能力相当；5)相移方面，当网格间距增大时，交错网格和旋 转交错网格有可能出现相移现象；6)频散方面，同位网格的频散现象不明显． 关键词 有限差分，交错网格，辅助网格，旋转交错网格，同位网格，TTI介质，黏弹性TTI介质 doi：10．6O38／pg2O13O423 中图分类号 P315 文献标识码 A Seismicwavefieldsimulation jn2D elasticandviscoelasticmedia： comparison betweenfourdifferentkinds0ffinite-diffeFencegrids SUN Yao-chong ， ZHANG Yan-teng ， BAIChao—ying，。 (1．CollegeofGeologyEngineeringandGeomatic，Chang’anUniversity，Xian710054，China； 2．SchoolofEarthandSpaceScience，Uni．ofScienceandTechnologyofChina，Hefei，230026，China 3．InstituteofCoraputingGeophysics，Chang’anUniversity，xi’an710054，China) Abstract Inthispaper，weusethestaggeredgrid，auxiliarygrid，rotated staggeredgrid and non-staggered grid finitedifferencemethodstOsimulatethewavefieldin2D elasticandviscoelasticTTImedia。respect