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

辛维 闫子超 梁文全 陈雨红 杨长春 地球物理学报 2015年第07期

摘要:有限差分方法是波场数值模拟的一个重要方法,交错网格差分格式比规则网格差分格式稳定性更好,但方法本身都存在因网格化而形成的数值频散效应,这会降低波场模拟的精度与分辨率.为了缓解有限差分算子的数值频散效应,精确求解空间偏导数,本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限差分系数;同时应用最大最小准则作为模拟退火(SA)优化算法求解差分系数的数值频散误差判定标准来求解有限差分系数.通过上述两种方法,分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析,并与传统泰勒展开算法、最小二乘算法进行比较,验证了线性化方法和模拟退火方法都能有效压制数值频散,并比较了各个算法的特点.

关键词:弹性波方程数值模拟线性化方法模拟退火算法

单位:中国科学院地质与地球物理所中国科学院油气资源研究重点实验室 北京100029 中国科学院大学 北京100049

注:因版权方要求,不能公开全文,如需全文,请咨询杂志社

地球物理学报

北大期刊、CSCD期刊、统计源期刊

期刊详情