计算机数值模拟实验报告,数值模拟实验报告..docx

数值模拟实验报告.

实验题目地震记录数值模拟的这几模型法二、实验目的掌握褶积模型基本理论、实现方法与程序编制，由褶积模型初步分析地震信号的分辨率问题三、实验原理1、褶积原理地震勘探的震源往往是带宽很宽的脉冲，在地下传播、反射、绕射到测线，传播经过中高频衰减，能量被吸收。吸收过程可以看成滤波的过程，滤波可以用褶积完成。在滤波中，反射系数与震源强弱关联，吸收作用与子波关联。最简单的地震记录数值模拟，可以看成反射系数与子波的褶积。通常，反射系数是脉冲，子波取雷克子波。雷克子波Wave(t) = (2)反射系数：? ?rflct(z)=(3)褶积公式：数值模拟地震记录trace(t):trace(t) =rflct(t)*wave(t)反射系数的参数由z 变成了t，怎么实现？在简单水平层介质，分垂直和非垂直入射两种实现，分别如图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(2)离散序列的褶积?四、实验内容1、垂直入射地震记录数值模拟的褶积模型；2、非垂直入射地震记录数值模拟的褶积模型；3、点绕射的地震记录数值模拟的褶积模型；五、方法路线 根据褶积模型的实验原理编写C++程序，完成对于垂直入射波的褶积。改变子波的长度与主频的大小，关注其对于实验结果的影响。通过增加一个地层来模拟地下两层界面的反射情况，通过改变界面的高度来说明其对于实验结果的影响，同时改变子波长度与主频。非垂直入射改变时间t来改变褶积结果显示地面情况。点绕射模型通过时间的改变，任意位置检波器的自激自搜时间来改变褶积结果，同时改变子波长度与主频来分析影响。六、实验结果：1.垂直入射：单层界面h=1000，fm=10，Nw=80v=2000v=3250v=4500v=5750v=7000h=1000，v=4000，fm=10Nw=80Nw=100Nw=120Nw=140Nw=160h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm-46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=46v=4000，fm=10，Nw=80h=800h=1000h=1200h=1400h=1600双层界面H12=20h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=46H12=40h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=46H12=60h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=46H12=80h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=462.非垂直入射1> h=1000，fm=10，Nw=80 v=3000v=4000v=5000v=6000v=70002> h=1000，v=4000，fm=10Nw=80Nw=100Nw=120Nw=140Nw=160h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=55h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=55v=4000，fm=10，Nw=80h=800h=1000h=1200h=1400h=16003.点绕射卷积模型图件Nw=80 fm=15相关程序：#include#include#includeintCnltn(float,float); intRflct(float,float,float); int Wave(float,float); #define Nx 128#define Nt