一种有效吸收边界条件MATLAB实现.pdf

一种有效吸收边界条件的MATLAB实现 陈敬国 中国地质大学（北京） 地球物理与信息技术学院 （100083） E-mail: chenjg_cugb@163.com 摘要：用有限差分法模拟地震波场是研究地震波在地球介质中传播的有效方法。但我们在实 验室进行波场数值模拟时有限差分网格是限制在人工边界里面，即引入了人工边界条件。本 文采用 Clayton_Engquist_Majda 二阶吸收边界条件，通过 MATLAB 编程实现了这一算法。依 靠 MATLAB 具有更加直观的、符合大众思维习惯的代码，为用户提供了友好、简洁的程序 开发环境，方便同行们交流。 关键词：有限差分法，地震波场，数值模拟，吸收边界条件，MATLAB 1. 引言 [1] 用有限差分法模拟地震波场是研究地震波在地球介质中传播的有效方法 。但我们在实 验室进行波场数值模拟时，只能在有限的空间进行，所以有限差分网格是限制在人工边界里 面，即引入了人为的边界条件。这种人为边界条件的引入将对有限区域内的波场值的计算带 来严重影响，所以必须进行特殊的边界处理。边界条件处理的好坏直接影响地震正演模拟的 最终效果。本文中我们采用Clayton_Engquist_Majda二阶吸收边界条件[2]。 被称作是第四代计算机语言的 MATLAB 语言，利用其丰富的函数资源把编程工作者从 繁琐的程序代码中解放出来。MATLAB 用更加直观的、符合大众思维习惯的代码，为用户提 供了友好、简洁的程序开发环境。本文介绍运用 MATLAB 实现带有吸收边界条件的地震波场 数值模拟方法和步骤，便于同行们交流，亦可用于本科地震理论的教学中，让学生们在程序 演示中理解地震波的传播规律。 2. Clayton_Engquist_Majda 二阶吸收边界条件 我们给定二维标量声波波动方程（含震源）： ∂2P ∂2P 1 ∂2P 2 + 2 +f (x ,z ) 2 2 （1） ∂x ∂z v ∂t v 式中：P 是声波波场P (x, z ,t ) ， 是声波速度v(x, z) ，f (x , z ) 是震源。 对（1）式进行时间和空间 2 阶精度有限差分离散（见图 1），整理后可得 k +1 k k −1 2 k k k P 2P =−P +A (P+ +P− +P, +1 i, j i, j i, j i 1, j i 1, j i j （2） k k +1 2 2 +P − −4P ) +v (Δt) Src(k ) i, j 1 i , j k vΔt 式中， P P (i=Δx, jΔz ,k Δt ) ， Δx , Δz , Δt 为别为空间、时间离散步长，