最优化广义离散牦犺犪状状狅状奇异核褶积微分算子地震波场模拟.pdf

年 月 第 卷 第 期 ２０１３ ６ ４８ ３ 　 　 　 ·地震模拟 · 文章编号： （ ） １０００７２１０２０１３０３０４１００７ 最优化广义离散犛犺犪狀狀狅狀奇异核 褶积微分算子地震波场模拟 ①② ① ①③ ①② ①② 刘少林 刘有山 汪文帅 李小凡 张美根 　 　 　 　 （ 中国科学院地质与地球物理研究所，北京 １０００２９； 中国科学院地球深部重点实验室，北京 １０００２９； ① ② ③宁夏大学数学计算机学院，宁夏银川 ７５００２１） 刘少林，刘有山，汪文帅，李小凡，张美根 最优化广义离散 奇异核褶积微分算子地震波场模拟 石油地 ． 犛犺犪狀狀狅狀 ． 球物理勘探， ， （）： ２０１３４８３ ４１０ ４１６ ～ 摘要 本文运用广义离散 奇异核褶积微分算子（ ）计算波动方程的空间导数，推导了一阶和二 Ｓｈａｎｎｏｎ ＧＤＳＣＤ 　 阶ＧＤＳＣＤ的具体形式，并提出了系数优化方法，即在频率域逼近平面波的真实导数，得到不同半径和采样下限 的最优权系数。通过计算滤波响应分析算子精度，与多种数值方法进行了对比。模型试算结果表明，本文构造 的最优化ＧＤＳＣＤ方法模拟地震波具有高精度、高效率的特点。 关键词 地震波场 广义Ｓｈａｎｎｏｎ奇异核 数值模拟 　 　 　 中图分类号： 文献标识码： Ｐ６３１ Ａ 　　 ［ ］ １０ 有较高的计算效率 ，但算子无法覆盖高波数域， ［ ， ］ １ 引言 对地震波高频部分模拟效果并不理想 １１１２ 。基于 　 快速傅里叶变换的伪谱法虽能在波数域对波数有效 地震波数值模拟是研究地震波传播规律的有效 覆盖，但伪谱法缺乏因果律，有悖于微分的局部属 ［ ］ ［ ］ 手段。常用的地震波数值模拟方法有：几何射线 １３