基于广义正交多项式褶积微分算子的地震波场数值模拟方法

被引:14
作者
程冰洁
李小凡
龙桂华
机构
[1] 中国科学院地质与地球物理研究所
基金
国家自然科学基金重点项目; 国家杰出青年科学基金;
关键词
Forsyte广义正交多项式; 交错网格有限差分; 复杂非均匀介质; 地震波场; 数值模拟;
D O I
暂无
中图分类号
P631.4 [地震勘探];
学科分类号
0818 ; 081801 ; 081802 ;
摘要
地震波场模拟方法研究对于与波动现象有关的地震学问题的重要性是不言而喻的.就目前现有的各种正演算法来说,精度较高的算法(如有限元法、谱元法、高阶有限差分法等),其计算速度较慢;计算速度较快的算法(如低阶有限差分法、付氏伪谱法等)计算精度却比较低.为了兼顾地震波场模拟的精度与速度,本文推出了一种快速的、高精度地震波场模拟方法(基于Forsyte广义正交多项式的褶积微分算子法),该方法是以计算数学中的Forsyte广义正交多项式插值函数为基础,构建一个新的褶积微分算子,并将该算子引入到地震波动方程的一阶速度-应力方程的空间微分运算中去,采用时间交错网格有限差分算子替代普通的差分算子以匹配高精度的褶积微分算子,从而构造一种全新的地震波场数值模拟方法.该方法同时具有广义正交多项式方法的高精度和短算子低阶有限差分算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型中的波场数值模拟实验证实了该方法的可行性及优越性.
引用
收藏
页码:531 / 537
页数:7
相关论文
共 11 条
[1]   数值模拟在地球动力学中的研究进展 [J].
郑洪伟 ;
李廷栋 ;
高锐 ;
贺日政 .
地球物理学进展, 2006, (02) :360-369
[2]   地震波传播数值模拟 [J].
裴正林 ;
牟永光 .
地球物理学进展, 2004, (04) :933-941
[3]   有限元法与伪谱法混合求解弹性波动方程 [J].
马德堂 ;
朱光明 .
地球科学与环境学报, 2004, (01) :61-64
[4]   错格实数傅里叶变换微分算子及其在非均匀介质波动传播研究中的应用 [J].
赵志新 ;
徐纪人 ;
堀内茂木 .
地球物理学报, 2003, (02) :234-240
[5]   各向异性弹性波场的有限元数值模拟 [J].
张美根 ;
王妙月 ;
李小凡 ;
杨晓春 ;
王磊 .
地球物理学进展, 2002, (03) :384-389+413
[6]   一阶弹性波方程交错网格高阶差分解法 [J].
董良国 ;
马在田 ;
曹景忠 ;
王华忠 ;
耿建华 ;
雷兵 ;
许世勇 .
地球物理学报, 2000, (03) :411-419
[7]   各向异性介质中波动方程有限元法模拟及其稳定性 [J].
周辉 ;
徐世浙 ;
刘斌 ;
何樵登 .
地球物理学报, 1997, (06) :833-841
[8]   声波与弹性波场数值模拟中的褶积微分算子法 [J].
张中杰,滕言文,杨顶辉 .
地震学报, 1996, (01) :63-69
[9]   三维波动方程有限差分正演方法 [J].
金胜汶,陈必远,马在田 .
地球物理学报, 1994, (06) :804-810
[10]  
物探数据处理的数学方法[M]. 地质出版社 , 陈天与, 1981