叠加地震记录的相移波动方程正演模拟素材.docxVIP

叠加地震记录的相移波动方程正演模拟素材.docx

  1. 1、本文档共25页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
本科生实验报告实验课程数值模型模拟学院名称地球物理学院专业名称勘测技术与工程(石油物探)学生姓名学生学号指导教师熊高君实验地点5417实验成绩2015年5月 成都理工大学《地震数值模拟》实验报告实验时间2015年 5月开课单位地球物理学院指导教师熊高君实验题目:叠加地震记录的相移波动方程正演模拟姓名学号班级专业勘测技术与工程(石油物探)院(系)地球物理学院 地球探测与信息技术系单项成绩内容理解写作结构程序设计模型设计计算结果结果分析总成绩实验报告实验题目叠加地震记录的相移波动方程正演模拟二、实验目的掌握各向同性介质任意构造、水平层状速度结构地质模型的相移波动方程正演模拟基本理论、实现方法与程序编制,由正演记录初步分析地震信号的分辨率。三、原理公式1、地震波传播的波动方程设(x,z)为空间坐标,t为时间,地震波传播速度为v(x,z),则二维介质中任意位置、任意时刻的地震波场为p(z,x,t):压缩波——纵波。则二维各向同性均匀介质中地震波传播遵循的声波方程:2、傅里叶变换的微分性质p(t)与其傅里叶变换的P()的关系:则有时间微分性质:ω为频率,,T为周期。同理有空间微分性质:k为波数,k=,为波长3、地震波传播的相移外推公式令速度v 不随x 变化,只随z 变化,则利用傅里叶变换微分性质(3)和(4)式,把波动方程(1)式变换到频率‐波数域,得:或:令:则(5)式的解为: 包括上行波和下行波两项,正演模拟取上行波:若和间隔为,速度v(z)为在此间隔内不随Z变的常数,(7)式实现波场从到的延拓,即:在深度开始向上延拓到,若延拓深度为零,即:,则对于任意深度到的延拓,可得正演模拟中地震波的传播方程(延拓公式4、初始条件和边界条件按照爆炸界面理论,反射界面震源在t=0 时刻同时起爆,此时刻的波场就是震源。根据不同情况,可直接使用反射系数脉冲或子波作震源。如果直接使用反射系数作震源脉冲,则初始条件可表示为:对时间t和空间x做二维傅立叶变换,则得频率-波数域的初始波场。边界条件: 其他参数都是在范围内定义的。5、边界处理(1)边界反射问题把实际无穷空间区域中求解波场的问题化为有穷区域求解时,左右两边使用零边界条件。物理上假设探区距与两个端点很远,在两个端点上收到的反射波很弱。但是,上述条件在实际中不能成立,造成零边界条件反而成为绝对阻止波通过的强反射面。在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。(2)边界强反射的处理镶边法、削波法、吸收边界都能有效消除边界强反射。削波法就是在波场延拓过程中,每延拓一次,在其两侧均匀衰减到零,从而消除边界强反射的影响。假设横向总长度为NX,以两边Lx道吸波为例,有以下吸波公式:6、数字化根据数字信号处理的采样定理,把连续的信号变为计算机能处理的数字信号,使相移法正演模拟得以实现。频域抽样定理:一个频谱受限信号,如果时间只占据的范围,若在频域以不大于频率间隔对信号的频谱采样,则抽样到的离散信号可以唯一表示原信号。时域抽样定理:一个时间受限信号,如果频谱只占据的范围,则信号可以用等间隔的抽样值唯一表示出来,而时间抽样间隔必须不大于,,。四、实验内容削波法相位移正演模拟(1)点绕射构造和水平层状速度模型(参数如图1所示)的正演数值模拟;1)削波的正演;2)无削波的正演;(2)计算中点和两个边界的信号位置,分析实验结果的正确性;(3)做同样模型的褶积模型数值模拟,对比分析两者的异同;图1 点绕射构造与水平速度模型(4)改变绕射点位置、速度,再做正演模拟;五、方法路线1、参数初始化;2、形成边界削波数据;3、波场初始化;3、Zmax层波场延拓到深度Zmax-1;5、Zi+1层波场延拓到深度Zi;6、重复5,从Iz=Nz-1开始,直到Iz=1,得测线上的频率—空间域波场;7、频率-空间域波场对频率做反傅里叶变换,得时间-空间波场;8、使用Fimage软件显示所得结果。六、实验结果1、结果显示及分析(点绕射结构)图2 削波(左)正演模拟结果与未削波正演模拟结果(右)由图2分析可知:削波前,由于在两个端点上收到的反射波很弱,造成零边界条件反而成为绝对阻止波通过的强反射面,在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。削波后,边界假反射的影响消除了,且曲线变得更光滑了。a、速度V1=5000、V2=5500,深度h=2000时,改变绕射点x的位置;图3x=Nx/4-1时,反射系数(左)与削波正演模拟结果(右)图4x=Nx/2-1时,反射系数(左)与削波正演模拟结果(右)由图3与图4对比分析可知:当速度,深度一定时,改变绕射点的位置,曲线随绕射点位置的改变而左右移动,且移动趋势与绕射点位置的移动一致。b、绕射点x=Nx/2-1,深度h=2000,速度V2=5500时,改变速度V1;V2=5500V1=4

文档评论(0)

希望之星 + 关注
实名认证
文档贡献者

我是一名原创力文库的爱好者!从事自由职业!

1亿VIP精品文档

相关文档