用于弹性波方程数值模拟的有限差分系数确定方法.pdfVIP

第 卷 第 期 地 球 物 理 学 报 ， ５８ ７ Ｖｏｌ．５８ Ｎｏ．７ 　 　 　 　 　 年 月 ， ２０１５ ７ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｊｕｌ．２０１５ 　 　 　 辛维，闫子超，梁文全等 用于弹性波方程数值模拟的有限差分系数确定方法 地球物理学报， （）： ， ： ．２０１５． ． ５８７ ２４８６２４９５ｄｏｉ ／ １０．６０３８ ｊｇ ， ， ， ＸｉｎＷ ＹａｎＺＣＬｉａｎ Ｗ Ｑ ｅｔａｌ．２０１５．Ｍｅｔｈｏｄｓｔｏｄｅｔｅｒｍｉｎｅｔｈｅｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｃｏｅｆｆｉｃｉｅｎｔｓｆｏｒｅｌａｓｔｉｃｗａｖｅｅｕａｔｉｏｎ ｇ ｑ （ ）， （）： ， ： ／ ｍｏｄｅｌｉｎ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ５８７ ２４８６２４９５ｄｏｉ１０．６０３８ ｇ 狆狔 ｊｇ 用于弹性波方程数值模拟的有限差分系数确定方法 ， ， １２ １２ １ １ １ 辛维 ，闫子超 ，梁文全 ，陈雨红 ，杨长春 １中国科学院地质与地球物理所 中国科学院油气资源研究重点实验室，北京　１０００２９ ２中国科学院大学，北京 １０００４９ 　 摘要 有限差分方法是波场数值模拟的一个重要方法，交错网格差分格式比规则网格差分格式稳定性更好，但方 　 法本身都存在因网格化而形成的数值频散效应，这会降低波场模拟的精度与分辨率 为了缓解有限差分算子的数 ． 值频散效应，精确求解空间偏导数，本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限 差分系数；同时应用最大最小准则作为模拟退火（ ）优化算法求解差分系数的数值频散误差判定标准来求解有限 ＳＡ 差分系数 通过上述两种方法，分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析， ． 并与传统泰勒展开算法、最小二乘算法进行比较，验证了线性化方法和模拟退火方法都能有效压制数值频散，并比 较了各个算法的特点． 关键词 弹性波方程；数值模拟；线性化方法；模拟退火算法 　 ： ／ 中图分类号 收稿日期 ， 收修定稿 ｄｏｉ１０．６０３８ Ｐ６３１ ２０１４０５１４２０１５０４２３ ｊｇ 　 犕犲狋犺狅犱