地震记录数值模拟地褶积模型法.pdfVIP

本 科 生 实 验 报 告 实 验 课 程 数 值 模 型 模 拟 学 院 名 称 地 球 物 理 学 院 专 业 名 称 勘 测 技 术 与 工 程 学 生 姓 名 学 生 学 号 指 导 教 师 熊 高 君 实 验 地 点 5417 实 验 成 绩 2015 年 5 月 成都理工大学 《地震数值模拟》实验报告 实验时间 2015年 5 月 开课单位 地球物理学院 指导教师 熊高君 实验题目：地震记录数值模拟的褶积模型法 姓名 学号 班级 专业 勘测技术与工程（石油物探） 院（系） 地球物理学院 地球探测与信息技术系 内容理解 写作结构 单 程序设计 项 成 绩 模型设计 计算结果 结果分析 总成绩 实验报告 一、 实验题目： 地震记录数值模拟的褶积模型法 二、 实验目的： 掌握褶积模型基本理论、 实现方法与程序编制， 由褶积模型初步分析地震信 号的分辨率问题。 三、原理公式 1、褶积原理 地震勘探的震源往往是带宽很宽的脉冲，在地下传播、反射、绕射到测线， 传播经过中 高频衰减，能量被吸收。吸收过程可以看成滤波的过程，滤波可以 用褶积完成。在滤波中， 反射系数与震源强弱关联，吸收作用与子波关联。最 简单的地震记录数值模拟，可以看成反射系数与子波的褶积。通常， 反射系数是 脉冲，子波取雷克子波。 （1） 雷克子波： wave(t)=cos( 2 ft)* （2 ） 反射系数： （3 ） 褶积公式： 数值模拟地震记录 trace(t): trace(t) =rflct(t)*wave(t) ； 反射系数的参数由 z 变成了 t ，怎么实现？在简单水平层介质，分垂直和 非垂直入射两种实现，分别如图 1 和图 2 所示。 图 1 图 2 1） 垂直入射： t=2h/v ； 2 ）非垂直入射： t= 2 、褶积方法 （1） 离散化（数值化） 计算机数值模拟要求首先必须针对连续信号离散化处理。 反射系数在空间模 型中存在，不同深度反射系数不同， 是深度的函数。 子波是在时间记录上一延续 定时间的信号，是时间的概念。 在离散化时， 通过深度采样完成反射系数的离散 化，通过时间采样完成子波的离散化。如果记录是 Trace(t) ，则记录是时间的 函数，以时间采样离散化。 时间采样间距以 Δt 表示，深度采样间距以 Δz 表示。 在做多道的数值模拟时，还有横向 Δx 的概念，横向采样间隔以