数值模拟实验一.docVIP

数值模型模拟实验报告 实验名称：地震记录数值模拟的褶积模型法实验 学院：地 球 物 理 学 院 学号：2010050603xx 教师：熊 高 军 姓名：Blackheart--Mike 日期：2010.6.13 实验一 一、 实验题目 地震记录数值模拟的褶积模型法 二、 实验目的 掌握褶积模型基本理论、实现方法与程序编制，由褶积模型初步分析地震信号的分辨率问题。 三、 原理公式 1、褶积原理 地震勘探的震源往往是带宽很宽的脉冲，在地下传播、反射、绕射到测线，传播经过中高频衰减，能量被吸收。吸收过程可以看成滤波的过程，滤波可以用褶积完成。在滤波中，反射系数与震源强弱关联，吸收作用与子波关联。最简单的地震记录数值模拟，可以看成反射系数与子波的褶积。通常，反射系数是脉冲，子波取雷克子波。 雷克子波 反射系数： （3） 褶积公式： 数值模拟地震记录trace(t): trace(t) =rflct(t)?wave(t) 反射系数的参数由z变成了t，怎么实现？在简单水平层介质，分垂直和非垂直入射两种实现，分别如图1和图2所示。 1）垂直入射： 2）非垂直入射： 2、褶积方法 （1） 离散化（数值化） 计算机数值模拟要求首先必须针对连续信号离散化处理。反射系数在空间模型中存在，不同深度反射系数不同，是深度的函数。子波是在时间记录上一延续定时间的信号，是时间的概念。在离散化时，通过深度采样完成反射系数的离散化，通过时间采样完成子波的离散化。如果记录是Trace(t)，则记录是时间的函数，以时间采样离散化。时间采样间距以?t表示，深度采样间距以?z表示。在做多道的数值模拟时，还有横向x的概念，横向采样间隔以?x表示。离散化的实现：t=It×?t；x=Ix×?x；z=Iz×?z或：It=t/?t; Ix=x/?x; Iz=z/?z 离散序列的褶积 四、 实验内容 1、 垂直入射地震记录数值模拟的褶积模型； 2、 非垂直入射地震记录数值模拟的褶积模型； 3、 点绕射的地震记录数值模拟的褶积模型； 五、 实验步骤 （一）实验一： 1、 根据垂直入射褶积模型理论算法，填充程序（附后）的下划线部分，使程序完整，调试程序，算出结果，用“Fimage”显示软件显示褶积结果； 2、 根据非零偏移距算法，编制非零偏移距褶积模型程序，算出结果，用“Fimage”显示软件显示褶积结果。（参考垂直入射褶积模型理论算法和程序，子波与反射层不变）； 3、 变换子波的长度Nw和主频fm，重复1和2； （1）Nw=80的情况下：1）fm=5，2）fm=10，3）fm=15，3）fm=20，3）fm=55； （2）Nw=160的情况下：1）fm=5，2）fm=10，3）fm=15，3）fm=20，3）fm=55； 六、 实验结果 结果显示 对比分析 对于垂直入射而言，当子波长度保持不变而频率逐渐升高时画出的图像都在0.5下方而且一次变得细起来，而对于子波长度不同频率相同时子波长度打的图像位于0.5上方，同样子波长度不变频率增加时图像变细 七、 讨论建议 1、 实验收获：获得了垂直入射地震记录数值模拟的褶积模型和非垂直入射地震记录数值模拟的褶积模型 2、 存在问题：开始的时候Nt取值较小让其范围超出了最大的容纳范围，最后经过改正把Nt改的比原来大一些才得出了正确的结果，开始非垂直绘图出来都是直线，经检查是程序出现咯一点点小问题，但却花了很长时间才找出来错误 3、 心得体会：在做程序的时候一定要细心，一个小小的错误都会耽误自己做程序的效率，而且写程序的时候一定要对其这样发现错误检查起来也方便得多. 实验二：两层反射界 如图3所示为两层反射界面的模型，对比不同间距H12情况下，不同长子波长度Nw和主频fm的褶积模型。 1、H12=20m时： （1）Nw=80的情况下： 1）fm=5，2）fm=10，3）fm=15，4）fm=20，5）fm=55； （2）Nw=160的情况下： 1）fm=5，2）fm=10，3）fm=15，4）fm=20，5）fm=55；2、H12=40m时： （1）Nw=80的情况下：1）fm=5，2）fm=10，3）fm=15，3）fm=20，3）fm=55； （2）Nw=160的情况下：1）fm=5，2）fm=10，3）fm=15，3）fm=20，3）fm=55； 3、H12=60m时： （1）Nw=80的情况下：1）fm=5，2）fm=10，3）fm=15，3）fm=20，