工程与环境物探

工程与环境物探

实验报告及心得体会

2014~2015学年第二学期

教师：张玮

学生：

专业：地球物理学

班级：地物二班

学号：

实验一 地震信号的频谱分析 一 实验原理及要求

褶积合成地震反射记录：选择合适的编程语言或通用软件合成地震反射记录。结果用软件（Excel ）给出。 雷克子波公式为：

2

(f t)2

f (t)(12(f t))e

m m ππ-=-

真实的地震记录夹杂的干扰成分，通过频谱分析，做卷积，再通过反滤波得到的记录就清晰很多。雷克子波的宽度可以改变时间来实现，根据真实的记录要选择合适频率。

傅氏正反变换公式

二.实验程序

#include

#include

#include

#define v1 1000 #define v2 2000 #define v3 3500 #define B 0.31

#define pi 3.1415926 #define f 50 #define t 0.002

FILE *fp,*fp1,*fp2;

void cov(double Ref[],double Rec[],int N,int n_f,double x[]); void main() {

fp=fopen("fanshexishu.txt","w+"); fp1=fopen("zibo¨.txt","w+"); fp2=fopen("dizhenbo¨.txt","w+");

210210()()1()()N j nk N

n N j nk N k X k x n e x n X k e N ππ--=-=⎧=⎪⎪⎨⎪=⎪⎩

∑∑

int i;

double R[500];

double x[500];

double V[500],D[500];

double z[61];

double d1,d2,d3;

d1=B*sqrt(sqrt(v1));

d2=B*sqrt(sqrt(v2));

d3=B*sqrt(sqrt(v3));

for(i=0;i<10;i++)

{

V[i]=v1;

D[i]=d1;

}

for(i=10;i<100;i++)

{

V[i]=v2;

D[i]=d2;

}

for(i=100;i<250;i++)

{

V[i]=v3;

D[i]=d3;

}

for(i=0;i<499;i++)

{

R[i]=(V[i+1]*D[i+1]-V[i]*D[i])/(V[i+1]*D[i+1]+V[i]*D[i]);

R[499]=0;

}

for(i=0;i<500;i++)

{

fprintf(fp,"%lf\n",R[i]);

}

for(i=0;i<61;i++)

{

z[i]=(1-pi*f*(i-30)*t*pi*f*(i-30)*t)*exp(-(pi*f*(i-30)*t)*(pi*f*( i-30)*t));

fprintf(fp1,"%lf\n",z[i]);

}

cov(R,z,500,61,x);

for(i=0;i<500;i++)

{

fprintf(fp2,"%lf\n",x[i]);

}

fclose(fp);

fclose(fp1);

fclose(fp2);

}

void cov(double Ref[],double Rec[],int N,int n_f,double x[]) {

int i,j,m1,m2;

double sum;

for(i=0;i

{

m1=n_f/2-i;

m2=i-n_f/2;

sum=0;

if(m1>=0)

{

for(j=0;m1

sum+=Ref[j]*Rec[m1];

}

else if((N-m2)>=n_f)

{

for(j=0;j

sum+=Ref[m2]*Rec[j];

}

else

{

for(j=0;m2

sum+=Rec[j]*Ref[m2];

}

x[i]=sum;

}

}

三.实验结果

四 实验心得与体会

在本实验中，一个实际地震记录道就是由无数多个子波组成的复合振动。振动的幅值大小与界面的反射系数成正比。因此可近似认为，一个反射记录道是地层反射系数序列和地震子波的卷积的结果。但是如果岩层厚度薄到一定程度，雷克子波将会发生重叠现象，将不再能识别这一层。

实验二 工程地震正演

一 实验原理及要求

主要目的是对几种不同的波（直达波，反射波，折射波）进行了解，工程地震方法正演：使用计算机通用软件（Excel ），计算并绘制两层及两层以上的直达波、折射波、反射波时距曲线，并说明它们各自的特点和相互关系；

直达波:1

t V x

±

=； 单层反射波：21220

t V x t += 其中：t0是双层垂直时间1

0*2t V H

=；

单层折射波:)cos(H 2t 1

1

2i V V x *+=

其中:i 为临界角，H1为层厚度；

多层反射：22

20

t σ

V x t += 其中在当炮检距很小时，用均方根速度

代替上覆地层速度，即：

2/12)*(

∑∑=i

i

i

t

V t V σ

；

多层折射：∑

-=+=11n )

cos(*2V x t n k k

k k V H α；