声波方程有限差分法数值模拟.pdfVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
声波方程有限差分法数值模拟

声波方程有限差分法数值模拟 学院:海洋地球科学学院 专业:地球信息科学与技术 年级:2012 级 姓名:王昊 指导教师:宋鹏 声波方程有限差分法数值模拟 2012 级 地球信息科学与技术 12040032037 王昊 摘要:本实验应用声波方程作为正演模拟的波动方程,将所提供震源函数离散 后绘图,并将给定两个二维速度-深度模型(一个小模型;一个大模型),绘出图 形。通过模拟地震波在介质中的传播,理解实际勘探中地震波在地层中的传播规 律,在模拟水平层状速度模型中,体会地震波在两种介质分界面的传播规律,并 能够从地震记录中识别出反射波,透射波,多次波,折射波和绕射波。并通过模 拟人工合成的地震记录,体会地震勘探基本原理和方法,验证地震波传播能量波 形变化趋势。 关键词:声波方程 数值模拟 地震波 有限差分 引言:地震波场模拟即地震正演,是指已知模型结构,通过物理或数值计算的 方法模拟该地质结构下的地震波的传播,最终合成地震记录,也可以认为其是野 外数据采集过程的室内再现。物理模拟花费昂贵,人们一般采用比较经济的数值 模拟技术。地震波场数值模拟是在给定数学模型(如弹性波方程,声波方程等)、 震源和地下几何界面、物性参数(岩层密度、速度等)情况下,研究弹性波或声 波的传播规律。地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶 变换法、有限元法和有限差分法等。相对于上述几种方法,有限差分法是一种更 为快速有效的方法。虽然其精度比不上有限元法,但因其具有计算速度快,占用 内存较小的优点,在地震学界受到广泛的重视与应用。本实验给定大小两个模型, 对于小模型,整个区域的速度值可设为常数,即只有一种介质,将震源点放在模 型中间,分别记录两个时刻的波前快照(即该时刻区域内所有网格点的波场值)。 第一时刻为地震波还未传播到边界上的某时刻,第二时刻为地震波已经传播到边 界上的某时刻,体会其人工边界反射;对于大模型,定义为水平层状速度模型(至 少两层);做两个实验,一是将震源点放在区域表层任一点,记录下某些时刻的 波前快照,体会地震波在两种介质的分界面上传播规律;二是合成一个地震记录, 即记录下与震源同一深度点的各点所有时刻的波场值,并指出记录上的同向轴分 别对应哪些波。 实验内容: 对于二维速度-深度模型,地下介质中地震波的传播规律可以近似地用声波方程描述: 2 2 2  u 2  u  u 2 v ( 2  2 )  S (t) (1) t x z u 其中v(x,z)是介质在点(x , z)处的纵波速度,本实验中v 为常数2500, 为描述速度位或者压力的波场, s(t) 为震源函数,本实验选用雷克子波(零相位子波) (2f / )2t2 s(t) e cos 2ft 上式中,t 为时间,fm 为中心频率,一般取为 20-40HZ,本实验取20,γ为控制频带宽度的参数, 一般取3-5,这里取4,在实际计算过程中,需把此震 源函数离散,参与波场计算。由图像可知子波能量 主要集中在前部,大致上t35 的时候S(t)都为零值, 故在绘制S(t)图像时k 取 100,来比较清晰得呈现雷 克子波。本实验不涉及吸收边界条件的使用。这样 声波方程数值模拟所需的(1)震源函数(2)地层 速度(波速)(3)边界条件均定义完全。 为求 (1)式的数值解,必须将此式离散化(包 括时间离散和空间离散,这里dt 取0.002,dh 取4, 即用有限差分来逼近导数,用差商代替微商。为此, 先把空间模型网格化,得到横向上 x=i*dh,纵向上

文档评论(0)

dajuhyy + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档