起伏地表下弹性波传播数值模拟方法对比研究

发布时间：2020-11-20 05:38

　　 随着地表速度复杂区域勘探力度的加大,地震勘探作业区域的速度结构越来越复杂。特别地,在我国西部探区和山前带区域,由于构造运动和风化作用的影响,地表构造特别复杂,低速异常普遍存在,由地表起伏引起的地震资料信噪比低、静校正问题复杂以及观测系统布置困难等问题正影响着复杂区域地震工作的开展。为了研究这一系列实际问题,必须对复杂地表下地震波的传播规律有清醒的认识。在地震波传播的研究中,地球物理界普遍认可地震波传播数值模拟,也就是通常所说的正演模拟。正演模拟的算法很多且各有优势,本次选择有限差分法和谱元法进行研究。有限差分法使用差分格式替代空间偏导项,该方法编程难度小、计算效率高以及对复杂地质模型适用的优点使其被广泛采用。而谱元法作为一种高阶有限元法,兼具有限元法剖分模型的贴体性和谱方法的计算效率,逐渐被更多人关注。本文先总结了复杂地表条件下地震波数值模拟方法,再从有限差分法和谱元法入手完成以下研究工作:1)通过应力镜像法实施自由边界条件后,使用交错网格差分完成了起伏地表下弹性波方程有限差分法数值模拟;2)在弹性波方程弱形式的基础上,使用Galerkin谱展开法、Lagrange插值基函数以及Gauss-Lobatto-Legendre积分方法,完成了起伏地表下弹性波方程谱元法数值模拟;3)对比有限差分法和谱元法的精度、计算效率和计算机内存占用情况,并使用MPI+CUDA的异构并行技术完成谱元法的加速。4)使用谱元法对地表起伏以及地表速度复杂模型进行弹性波传播数值模拟,并对波场特征做了详细分析。通过研究发现:1)谱元法的贴体网格划分能克服矩形网格剖分带来的角点数值散射,提高了数值模拟精度;2)有限差分法的计算时间和计算机内存占用均小于谱元法,但是并不能改善其由于矩形网格划分所带来的数值散射,使用MPI+CUDA的异构并行技术能极大的提高计算效率,使较大模型的正演变得现实;3)在起伏地表模型下弹性波的传播模式发生了很大的变化,散射波和多次波同时影响着记录,这也是起伏地表下地震资料信噪比低的原因。
【学位单位】：西南石油大学
【学位级别】：硕士
【学位年份】：2018
【中图分类】：P631.4
【部分图文】：

图２－３?ＰＭＬ吸收边界模型试算快照ｔ＝６６０ｍｓ：??

图２－５水平自由地表（地表过ｃｒＹＹ、ｃｒｚｚ方向）示意图??２０??

图２－６起伏自由地表（地表过（Ｊ？．、方向）示意图，黑色实线为为地表，实心方形代表切??应力，实心圆形代表正应力，垂直速度平行于切应力，水平速度垂直于正应力??
【参考文献】

相关期刊论文 前10条

1 刘有山;滕吉文;徐涛;刘少林;司芗;马学英;;三角网格谱元法地震波场数值模拟[J];地球物理学进展;2014年04期

2 刘有山;滕吉文;刘少林;徐涛;;稀疏存储的显式有限元三角网格地震波数值模拟及其PML吸收边界条件[J];地球物理学报;2013年09期

3 黄超;董良国;;可变网格与局部时间步长的交错网格高阶差分弹性波模拟[J];地球物理学报;2009年11期

4 黄超;董良国;;可变网格与局部时间步长的高阶差分地震波数值模拟[J];地球物理学报;2009年01期

5 杨旭明;裴正林;;复杂近地表弹性波波场特征研究[J];石油地球物理勘探;2007年06期

6 董良国;郭晓玲;吴晓丰;马在田;;起伏地表弹性波传播有限差分法数值模拟[J];天然气工业;2007年10期

7 王祥春;刘学伟;;起伏地表二维声波方程地震波场模拟与分析[J];石油地球物理勘探;2007年03期

8 王秀明;G.SERIANI;林伟军;;利用谱元法计算弹性波场的若干理论问题[J];中国科学(G辑:物理学 力学 天文学);2007年01期

9 董良国;复杂地表条件下地震波传播数值模拟[J];勘探地球物理进展;2005年03期

10 裴正林;任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟[J];石油地球物理勘探;2004年06期

本文编号：2891029