Matlab环境下瑞利波有限差分正演与曲线绘制

详细信息   

• 责任者：杨天春^c湖南科技大学湖南省普通高等学校土木工程施工过程与质量安全控制重点实验室^d1968^e博士^f副教授 ; 朱自强^c中南大学 ; 周勇
• 刊名：煤田地质与勘探
• 出版年：2010
• 卷期：38^b1总217
• 关键词：瑞雷波法20^a波动方程 ; 有限差分法 ; 正演问题 ; 数据处理 ; 数值模拟 ; 地震勘探
• 页码：62～65，70
• 索书号：P613.06/730-1
• 图表：^c参10插图
• 分类号：2100
• 文摘期号：2010/07

摘要

以瑞利波频散方程为出发点的Abo‐Zena传递矩阵等方法，只能研究层状介质中瑞利波的传播特性，对于非层状介质，传递矩阵无能为力。因此，不得不考虑有限元、有限差分等方法。根据弹性动力学方程，采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟，并在Matlab环境下实现编程计算，再现了瑞利波在近地表的传播状态，实现了地震剖面曲线的正振幅充填问题，从而获得了类似野外地震勘探的剖面记录。模拟结果表明，采用一阶差分格式所得到的地震记录存在较严重的数值频散，在采用有限差分法研究瑞利波“之”字形成因时，须尽量压制数值频散现象，才能获得更接近实际情况的频散曲线。