《地震子波波形显示及一维地震合成记录》由会员分享,可在线阅读,更多相关《地震子波波形显示及一维地震合成记录(22页珍藏版)》请在人人文库网上搜索。
1、地震勘探原理实验一地震子波波形显示及一维地震合成记录姓名: 学号:专业:地球物理勘察技术与工程 级 一、实验目的1. 认识子波,对子波的波形有直观的认识。(名词:零相位子波,混合相位子波,最小相位子波;了解子波的分辨率与频宽的关系;)2. 利用褶积公式合成一维地震记录。二、实验步骤1. 雷克子波零相位子波(最小相位子波) n= m1/m2为最大波峰m1和最大波谷m2之比钟型子波 xw为初相m为时间域主波峰与次波峰之比w(t)=exp(-2*Fm2*t2*ln(n)*sin(T-2*pi*Fm*t) n=m1/m2 最大相位子波(最大相位子波请同学们自己查找相关文献完成,非必须完成)其中代表子波。
2、的中心频率, t =i*dt,dt为时间采样间隔,i为时间离散点序号; 这里可以为 = 10,25,40,100 Hz等,采样间隔dt=0.002秒,i为0256;2. 根据公式编程实现不同频率的零相位子波的波形显示;不同中心频率的零相位子波图= 25:= 100:3. 其地质模型为:震源检波器第一层介质r(100)=1.0第二层介质r(200)=-0.7第三层介质r(300)=0.5第四层介质r(400)=0.4第五层介质r(450)=0.6设计反射系数(n=512),n为地层深度,其中,为第一层介质深度;,为第二层介质深度;,为第三层介质深度;,为第四层介质深度;,为第五层介质深度;其它为。
3、0。地震波在介质中传播,当到达介质分界面时,发生反射和透射,反射波被检波器接受,生成地震记录。反射系数表示地震波在两层介质分界面的能量重新分配,如r(100)=1.0,表示地震波入射到分界面时,只有一种波,反射纵波(或反射横波)。反射系数不为1.0时,表示当地震波入射到分界面时,产生两种反射波。反射系数为正,表示反射波相位与入射波相位相差2;反射系数为负,表示反射波相位与入射波相位相差,存在半波损失。4. 应用褶积公式合成一维地震记录,并图形显示;应用褶积公式求f(n)的程序为:#include#include#define PI 3.#define FM 100void main() dou。
4、ble fac(double x,double y,double z,int m,int n);FILE *fp;int i,j,x;double W,dt=0.002,t,a256;double b512=0;double r512=0;r100=1.0;r200=-0.7;r300=0.5;r400=0.4;r450=0.6;fp=fopen(Date.txt,w+);printf(please input x:n);scanf(%d,&x);for(i=0;i0&i-j#include13KFFT.C#include#define PI 3.#define FM 100void main。
5、() double fac(double x,double y,double z,int m,int n);FILE *fp,*fpr,*fpre,*fpi,*fpamp,*fpha;int i,j,x;double W,dt=0.002,t,a256,pr512,pi512=0.0,fr512,fi512,amp512,pha512;double b512=0;double r512=0;r100=1.0;r200=-0.7;r300=0.5;r400=0.4;r450=0.6;fp=fopen(褶积结果.txt,w+);fpr=fopen(反射系数.csv,w+);fpre=fopen(实部.txt,w+);fpi=fopen(虚部.txt,w+);fpamp=fopen(振幅.csv,w+);fpha=fopen(相位.csv,w+);for(i=0;i0&i-j=256)sum+=xj*yi-j;zi=sum。