原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
云南大学学报 自然科学版
一种弱数值频散的四阶
方法及二维声波场模拟
张朝元 宋国杰
大理学院数学与计算机学院云南 大理 清华大学 数学科学系北京
清华大学 计算机系北京
摘要针对二维声波方程利用近似解析离散化方法对空间高阶偏导数进行八阶离散并采用四阶
方法对时间导数进行四阶离散得到了八阶方法将该方法应用于双层介质模型和三层介
质模型中进行波场数值模拟同时与八阶方法和八阶方法进行了比较结果表明八阶方
法具有弱数值频散和高计算模拟效果等优势
关键词声波方程近似解析离散化四阶方法波场模拟数值频散
中图分类号 文献标志码文章编号
为了在认识地球内部结构和对油气资源进行勘探研究时更好地掌握地球内部信息我们需要研究声
波在地球介质中的正演传播过程 发展好的正演数值模拟方法一直是地球物理学和计算数学中的重要研
究课题之一当前求解波动方程的数值模拟方法主要有有限元方法 伪谱法 反射率法 有限差
分法 和射线追踪法 等这些方法的优缺点在等 研究中已有详细比较有限差分方法是最常用
的数值计算方法但是有限差分方法存在严重的数值频散 进而使人们难以识别真波和假波而影响
了波场的观测和后续反演成像中地质结构的识别等基于 引入的近似解析离散化算子 构造的
数值模拟方法有着很好的低数值频散特点目前发展有基于泰勒公式进行时间离散的
等数值方法为了进一步提高波场模拟精度本文针对二维声波方程
在时间离散上采用简化的二级四阶方法 在空间离散上采用基于泰勒公式推导
出来的具有八阶精度的近似解析离散算子 构造出一种时间精度为四阶 空间精度为八阶的数
值方法称为八阶方法针对二维声波进行了波场数值模拟同时与八阶
方法和八阶位移 应力交错网格 方法进行了比较 结果表明 本文建立的八阶 方法具有
有效压制数值频散和提高模拟计算效率等优点
八阶方法
四阶 方法设常微分方程为
这里 为需求解的二维函数 为已知函数 则求解方程 的简化的二级四阶 格
式 为
收稿日期
基金项目国家自然科学基金 云南省教育厅科学研究基金重点项 目 云南省教育厅科学研究基金
文档评论（0）