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

第38 卷 第1 期 煤田地质与勘探 Vol. 38 No.1 2010 年2 月 COAL GEOLOGY EXPLORATION Feb. 2010 文章编号: 1001-1986(2010)01-0062-04 Matlab 环境下瑞利波有限差分正演与曲线绘制 1 2 2 杨天春 , 朱自强 , 周 勇 (1. 湖南科技大学, 湖南省普通高等学校土木工程施工过程与质量安全控制重点实验室, 湖南 湘潭 411201; 2. 中南大学信息物理工程学院, 湖南 长沙 410083) : 以瑞利波频散方程为出发点的Abo-Zena 传递矩阵等方法，只能研究层状介质中瑞利波的传 播特性，对于非层状介质，传递矩阵无能为力。因此，不得不考虑有限元、有限差分等方法。根 据弹性动力学方程，采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟，并在 Matlab 环境下实现编程计算，再现了瑞利波在近地表的传播状态，实现了地震剖面曲线的正振幅 充填问题，从而获得了类似野外地震勘探的剖面记录。模拟结果表明，采用一阶差分格式所得到 的地震记录存在较严重的数值频散，在采用有限差分法研究瑞利波“之”字形成因时，须尽量压 制数值频散现象，才能获得更接近实际情况的频散曲线。 ：瑞利波；有限差分；交错网格；Matlab ：P631 ：A DOI: 10.3969/j.issn.1001-1986.2010.01.015 Finite-difference modeling of Rayleigh wave and drawing of seismograms in Matlab environment 1 2 2 YANG Tianchun , ZHU Ziqiang , ZHOU Yong (1. Hunan College Key Laboratory of the Construction Process, Quality and Safety of Civil Engineering, Hunan University of Science and Technology, Xiangtan 411201, China; 2. School of Info-Physics and Geomatics Engineering, Central South University, Changsha 410083, China) Abstract: The propagator matrix approach Abo-Zena based on the dispersion equation of Rayleigh wave could only be used for studying dispersion characteristics of layered media. For non-layered media, the propagator matrix approach was useless, and the finite-element method and the finite-difference method had to be considered. The authors applied the finit