基于Matlab实现的地震波场边界处理软件.doc

《基于Matlab实现的地震波场边界处理软件.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.基于Matlab实现的地震波场边界处理软件：嘉德学号：2015301130007院系：资源与环境科学学院摘要：用有限差分法模拟地震波场是研究地震波在地球介质中传播的有效方法。但我们在实验室进行波场数值模拟时有限差分网格是限制在人工边界里面，即引入了人工边界条件。本文采用Clayton_Engquist_Majda二阶吸收边界条件，通过MATLAB编程实现了这一算法。依靠MATLAB具有更加直观的、符合大众思维习惯的代码，为用户提供了友好、简洁的程序开发环境，方便同行们交流。利用Matlab本身所具有可视化功能以及像素识别功能，可以

2、将生成的动画电影进行识别，用于地震局实时分析有着深远意义。关键词：有限差分法，地震波场，吸收边界条件，MATLAB矢量帧,像素识别Abstract：ModelingseismicwavefieldwiththeFiniteDifferenceMethod(FDM)isaneffectivemethodtostudytheseismicwavepropagationintheearthmedium.Whenwemodelseismicwavefieldinthelaboratory,thefinitedifferencegridsar

3、erestrictedintheartificialboundary.Soitshouldintroducetheartificialboundaryconditions.ThispaperadoptsClayton_Engquist_MajdasecondabsorbingboundaryconditionsandrealizesthearithmeticwithMATLAB.TheMATLABcodesaredirectandaccordwithourthinkingcustom.Soitcanprovidethefriendl

4、yandsuccinctprogrammingenvironmentandiseasytocommunicatewithother.UsingthefunctionsofMatlabthatmakevisualizationcometrueandidentifythepixel,wecanidentifytheearthquakewavefield.Keywords:finitedifferencemethod,seismicwavefield,numericalmodeling,absorbingboundarycondition

5、s,MATLABWord资料.一、引言用有限差分法模拟地震波场是研究地震波在地球介质中传播的有效方法。但我们在实验室进行波场数值模拟时，只能在有限的空间进行，所以有限差分网格是限制在人工边界里面，即引入了人为的边界条件。这种人为边界条件的引入将对有限区域的波场值的计算带来严重影响，所以必须进行特殊的边界处理。边界条件处理的好坏直接影响地震正演模拟的最终效果。本文中我们采用Clayton_Engquist_Majda二阶吸收边界条件[2]。被称作是第四代计算机语言的MATLAB语言，利用其丰富的函数资源把编程工作者从繁琐的程序代码中解

6、放出来。MATLAB用更加直观的、符合大众思维习惯的代码，为用户提供了友好、简洁的程序开发环境。本文介绍运用MATLAB实现带有吸收边界条件的地震波场数值模拟方法和步骤，便于同行们交流，亦可用于本科地震理论的教学中，让学生们在程序演示中理解地震波的传播规律。二.、Clayton_Engquist_Majda二阶吸收边界条件我们给定二维标量声波波动方程（含震源）：（1）式中：是声波波场，是声波速度，是震源。对（1）式进行时间和空间2阶精度有限差分离散（见图1），整理后可得（2）式中，，为别为空间、时间离散步长，，，为震源函数。震源函数

7、：（3）Word资料.Clayton_Engquist_Majda二阶吸收边界条件的微分表达式可参见文献[2]，其左、右、上、下边界的差分格式分别为：三、基本算法步骤从图1可以看出，k+1时刻的波场值由k时刻和k-1时刻的波场值决定。所以在MATLAB里实现的基本算法步骤如下：(1)初始时刻的全波场值均为零，P(i,j,dt)=0（在MATLAB中初始从dt开始，不能从0开始）；(2)时刻2dt时，在炮点S(m,n)给出一个脉冲震源Src(t)（见式（3）），其它点波场P=0；(3)时刻t大于或等于3dt时，P(i,j,k+1)根据

8、式（2）计算，其它点波场P=0；(4)在波传播到四周边界时，左、右、上和下边界的波场值分别由式（4-1）、（4-2）、（4-3）和（4-4）计算出来。四、数值模拟由于是计算机模拟，为了能说明问题且便于计算，我们设地质模型边界为1，具体