《地震资料数据处理》课程设计
总结报告
专业班级：
姓名：
学号：
设计时间：
指导老师：
2011年5月30日
目录
一、设计内容………………………………………………………………
（1）褶积滤波………………………………………………
（2）快变滤波………………………………………………
（3）褶积滤波与快变滤波的比较…………………………
（4）设计高通滤波因子……………………………………
（5）频谱分析………………………………………………
（6）分析补零对振幅谱的影响……………………………
（7）线性褶积与循环褶积…………………………………
（8）最小平方反滤波………………………………………
（9）零相位转换……………………………………………
（10）最小相位转换…………………………………………
（11）静校正…………………………………………………
二、附录…………………………………………………………………………
（1）附录1：相关程序……………………………………
（2）附录2：相关图件……………………………………
【附录1：有关程序】
1.褶积滤波
CCCCCCCCCCCCCCCCC 褶积滤波CCCCCCCCCCCCCCCCC
PROGRAM MAIN
DIMENSION X(100),H1(-50:50),H2(-50:50),Y_LOW(200),Y_BAND(200)
PARAMETER (PI=3.141592654)
CCCCCCCC H1是低通滤波因子，H2为带通滤波因子CCCCCC
REAL X,H1,H2,Y_LOW,Y_BAND
REAL dt,F,F1,F2
INTEGER I
dt=0.002
F=70.0
F1=10.0
F2=80.0
OPEN(1,FILE='INPUT1.DA T',FORM='FORMATTED',STATUS='UNKNOWN') READ(1,*)(X(I),I=1,100)
CCCCCCCCCCCCCCCCCC低通滤波器CCCCCCCCCCCCCCCCC DO 10 I=-50,50
IF (I.EQ.0)THEN
H1(I)=2*F*PI/PI
ELSE
H1(I)=SIN(2*PI*F*I*dt)/(PI*I*dt)
END IF
10 CONTINUE
CCCCCCCCCCCCCCCC输出低通滤波因子CCCCCCCCCCCCCCCC OPEN(2,FILE='H1_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(2,*)(H1(I),I=-50,50)
CLOSE(2)
CALL CON(X,H1,Y_LOW,100,101,200)
CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCC OPEN(3,FILE='Y_LOW.DA T',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(3,*)(Y_LOW(I),I=51,150)
CLOSE(3)
CCCCCCCCCCCCCCCCCC带通滤波器CCCCCCCCCCCCCCCCCCCC DO 20 I=-50,50
IF(I.EQ.0)THEN
H2(I)=140
ELSE
H2(I)=SIN(2*PI*F2*I*dt)/(PI*I*dt)-SIN(2*PI*F1*I*dt)/(PI*I*dt)
END IF
20 CONTINUE
CCCCCCCCCCCCCCC输出带通滤波因子CCCCCCCCCCCCCCCCC OPEN(4,FILE='H2_BAND.DAT',FORM='FORMA TTED',STATUS='UNKNOWN')
WRITE(4,*)(H2(I),I=-50,50)
CLOSE(4)
CALL CON(X,H2,Y_BAND,100,101,200)
CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCCC OPEN(5,FILE='Y_BAND.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(5,*)(Y_BAND(I), I=51,150)
CLOSE(5)
END
CCCCCCCCCCCCCCCCCCCCC褶积函数CCCCCCCCCCCCCCCCCCCC SUBROUTINE CON(A,B,C,I,J,K)
DIMENSION A(I),B(J),C(K)
DO 1 K1=1,K
1 C(K1)=0.0
DO 2 I1=1,I
DO 2 I2=1,J
II=I1+I2-1
2 C(II)=C(II)+A(I1)*B(I2)*0.002
RETURN
END
2.快变滤波
CCCCCCCCCCCCCCC频率滤波CCCCCCCCCCCCCCCCCCCC PROGRAM MAIN
PARAMETER (PI=3.141592654)
REAL H_LOW(1:200),H_BAND(1:200),Y_LOW(1:200),Y_BAND(1:200)
REAL X(1:200)
INTEGER I,NFFT,K,NP
COMPLEX C(1:200),TEMP(1:200)
REAL DT,DF,F,F1,F2
F=70.0
F1=10.0
F2=80.0
DT=0.002
OPEN(1,FILE='INPUT1.DA T',FORM='FORMATTED',STATUS='UNKNOWN') READ(1,*)(X(I), I=1,100)
NP=100
K=LOG(100.0)/LOG(2.0)
IF(2**K.LT.100)THEN
K=K+1
NFFT=2**K
END IF
DF=1.0/(NFFT*DT)
DO 10 I=101,128
X(I)=0.0
10 CONTINUE
DO 20 I=1,NFFT
20 C(I)=CMPLX(X(I),0.0)
CCCCCCCCCCC 向频率域转换CCCCCCCCCCCCCCCCC
CALL FFT(NFFT,C,1.0)
CCCCCCCCCCC 低通滤波因子的设计CCCCCCCCCCC
DO 30 I=1,NFFT/2
IF(I*DF.LE.F)THEN
H_LOW(I)=1.0
ELSE
H_LOW(I)=0.0
END IF
30 CONTINUE
CCCCCCCCC 使滤波器理想化（对称）CCCCCCCCCC
DO 40 I=NFFT/2+1,NFFT
F=I*DF
H_LOW(I)=H_LOW(NFFT-I+1)
40 CONTINUE
CCCCCCCCCCCCCCC 实现滤波CCCCCCCCCCCCCCCCC
DO 50 I=1 , NFFT
50 TEMP(I)=C(I)*H_LOW(I)
CCCCCCCCCCCCCCC 向时间域变换CCCCCCCCCCCCC
CALL FFT(NFFT,TEMP,-1.0)
DO 60 I=1 , NFFT
60 Y_LOW(I)=REAL(TEMP(I))
OPEN(2,FILE='LOWPASS.DAT',FORM='FORMA TTED',STATUS='UNKNOWN') WRITE(2,*)(Y_LOW(I),I=1,NFFT)
CLOSE(2)
CCCCCCCCCCC 带通滤波因子的设计CCCCCCCCCCC
DO 70 I=1,NFFT/2
IF((I*DF.GE.F1).AND.(I*DF.LE.F2))THEN
H_BAND(I)=1.0
ELSE
H_BAND(I)=0.0
END IF
70 CONTINUE
CCCCCCCCC 使滤波器理想化（对称）CCCCCCCCCC
DO 80 I=NFFT/2+1,NFFT
F=I*DF
H_LOW(I)=H_BAND(NFFT-I+1)
80 CONTINUE
CCCCCCCCCCCCCCC 实现滤波CCCCCCCCCCCCCCCCC
DO 90 I=1 , NFFT
90 TEMP(I)=C(I)*H_BAND(I)
CALL FFT(NFFT,TEMP,-1.0)
DO 100 I=1 , NFFT
100 Y_BAND(I)=REAL(TEMP(I))
OPEN(3,FILE='BANDPASS.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(3,*)(Y_BAND(I),I=1,NFFT)
CLOSE(3)
CLOSE(1)
END
CCCCCCCCCCCCCCCCC子程序CCCCCCCCCCCCCCCC
CCCCC LX 为输入数据的点数CCCCCCCCCCCCCCCCC
CCCCC CX 为复型数组CCCCCCCCCCCCCCCCCCCCCC
CCCCC SIGNI 为转换标志CCCCCCCCCCCCCCCCCCCC
SUBROUTINE FFT(LX,CX,SIGNI)
COMPLEX CX(LX),CARG,CEXP,CW,CTEMP
J=1
SC=1.0/LX
IF(SIGNI.EQ.1.0)SC=1.0
SIG=-SIGNI
DO 300 I=1,LX
IF(I.GT.J)GO TO 100
CTEMP=CX(J)*SC
CX(J)=CX(I)*SC
CX(I)=CTEMP
100 M=LX/2
200 IF(J.LE.M)GO TO 300
J=J-M
M=M/2
IF(M.GE.1)GO TO 200
300 J=J+M
L=1
400 ISTEP=2*L
DO 500 M=1,L
CARG=(0.0,1.0)*(3.14159265*SIG*(M-1))/L
CW=CEXP(CARG)
DO 500 I=M,LX,ISTEP
CTEMP=CW*CX(I+L)
CX(I+L)=CX(I)-CTEMP
500 CX(I)=CX(I)+CTEMP
L=ISTEP
IF(L.LT.LX)GO TO 400
RETURN
END
3.褶积滤波与快变滤波的比较
CCCCCCCCCCCCC 褶积滤波与递归滤波的比较CCCCCCCCCCCCCC CCCCCCCCCCCCCCCC 褶积滤波的设计CCCCCCCCCCCCCCCCCC PROGRAM MAIN
PARAMETER(DT=0.002)
C H_NZ为非零相位滤波因子,H_Z为零相位滤波因子 C
C Y_CNZ为非零相位滤波结果,Y_CZ为零相位滤波结果 C
REAL Y_CNZ(1:100),Y_CZ(1:200)
REAL H_Z(1:20),H_NZ(1:20)
REAL X(1:50)
INTEGER I,J,K,N
REAL A(0:100),B(0:100)
REAL Y_DNZ(1:100),Y_DZ(1:100)
REAL X3(1:100),X4(1:100),X1(1:100),X2(1:100)
REAL DF
DATA H_NZ/1.0,5.254,11.6,13.7,8.47,0.775,-3.325,-2.764,-0.364,
$ 1.099,0.97,0.15,-0.37,-0.344,-0.06,-0.126,0.122,0.025,-0.042,
$ -0.043/
OPEN(1,FILE='INPUT3.DAT',FORM='FORMATTED',STATUS='UNKNOWN') READ(1,*)(X(I),I=1,50)
CLOSE(1)
CALL CON(X,H_NZ,Y_CNZ,50,20,69)
OPEN(2,FILE='Y_CNZ.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(2,*)(Y_CNZ(i),I=1,69)
CLOSE(2)
CALL AUTCOR(20,20,H_NZ,H_HZ)
CALL con(X,H_CZ,Y_CZ,50,20,69)
OPEN(3,FILE='Y_CZ.DAT',FORM='FORMATTED',STATUS='UNKNOWN') write(3,*)(Y_CZ(i),I=1,69)
CLOSE(3)
CLOSE(2)
CCCCCCCCCCCCCCC 递归滤波的设计CCCCCCCCCCCCCCC
C X 存放INPUT3里数据的数组 C
C Y_DNZ 为非零相位递归滤波，Y_DZ 为零相位递归滤波C
C X1,2,X,3,X4 为运算过程中的过度变量 C
A(0)=1.0
A(1)=4.0
A(2)=6.0
A(3)=4.0
A(4)=1.0
b(0)=0.0
B(1)=-1.254
B(2)=0.987
B(3)=-0.341
B(4)=0.0524
CCCCCCCCCCCCCCC 对A 补零CCCCCCCCCCCCCCCCCC DO 40 i=5,49
40 A(I)=0.0
CCCCCCCCCCCCCCC 对B 补零CCCCCCCCCCCCCCCCCC DO 50 i=5,49
50 B(I)=0.0
CCCCCCC 从外部数据向X 导入数据CCCCCCCCCCCC
OPEN(1,FILE='INPUT3.DA T',FORM='FORMATTED',STATUS='UNKNOWN') READ(1,*)(X(I),I=1,50)
CCCCCCCCCCC 非零相位递归滤波的设计CCCCCCCCC
DO 10 I=0,49
X1(I)=0.0
X2(I)=0.0
DO 20 J=1,I
20 X1(I)=X1(I)+A(J)*X(I-J)
IF(I.EQ.0) THEN
X2(I)=0.0
ELSE
DO 30 K=1,I
30 X2(i)=X2(i)+B(k)*Y_DNZ(I-K)
END IF
Y_DNZ(I)=X1(I)-X2(I)
10 CONTINUE
OPEN(2,FILE='Y_DNZ.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(2,*)(Y_DNZ(I),I=1,49)
CLOSE(2)
CCCCCCCCCC 零相位递归滤波的设计CCCCCCCCCCCC
DO 60 I=49,0,-1
X3(I)=0.0
X4(I)=0.0
DO 70 J=0,49-I
70 X3(I)=X3(I)+A(J)*Y_DNZ(I+J)
IF(I.EQ.49) THEN
X4(I)=0.0
ELSE
DO 80 K=2,49-I
80 X4(I)=X4(I)+B(K)*Y_DZ(I+K)
END IF
Y_DZ(I)=X3(I)-X4(I)
60 CONTINUE
OPEN(3,FILE='Y_DZ.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(3,*)(Y_DZ(I),I=1,49)
CLOSE(3)
CLOSE(1)
END
CCCCCCCCCCCCCCCCCCC 褶积子程序CCCCCCCCCCCCCCCCCCC SUBROUTINE CON(A,B,C,I,J,K)
DIMENSION A(I),B(J),C(K)
DO 1 K1=1,K
1 C(K1)=0.0
DO 2 I1=1,I
DO 2 I2=1,J
II=I1+I2-1
2 C(II)=C(II)+A(I1)*B(I2)*0.004
RETURN
END
CCCCCCCCCCCCCCCCCC 自相关子程序CCCCCCCCCCCCCCCCCC SUBROUTINE AUTCOR(J,K,C,A)
DIMENSION C(J),A(K)
DO 7 N=1,K
A(N)=0.0
I=N
DO 6 IN=I,J
IL=IN-N+1
6 A(N)=A(N)+C(IN)*C(IL)
7 CONTINUE
RETURN
END
4.设计高通滤波因子
CCCCCCCCCCC 高通滤波器的设计CCCCCCCCCCC
PROGRAM MAIN
PARAMETER (N=101,DT=0.004,PI=3.141592654,F=30.0)
REAL H(150),H2_R(128),H2_I(128),H_S(128)
INTEGER K,NFFT
COMPLEX H2(128)
OPEN(1,FILE='H1_HIGHPASS.DAT',FORM='FORMATTED',STATUS='UNKNOWN') OPEN(2,FILE='H2_HIGHPASS.DA T',FORM='FORMATTED',STATUS='UNKNOWN') K=LOG(101.0)/LOG(2.0)
IF(2**K.LT.101)THEN
K=K+1
ENDIF
NFFT=2**K
DO 10,I=-50,50
IF(I.EQ.0)THEN
H(I+51)=1.0/DT-2*F
ELSE
H(I+51)=-SIN(2*PI*30.0*I*DT)/(PI*I*DT)
END IF
10 CONTINUE
DO 20,I=102,128
20 H(I)=0.0
DO 30,I=1,128
30 H2(I)=CMPLX(H(I),0.0)
CALL FFT(128,H2,1.0)
DO 40,I=1,128
H2_R(I)=REAL(H2(I))
H2_I(I)=AIMAG(H2(I))
H_S(I)=(H2_R(I)**2+H2_I(I)**2)
40 CONTINUE
WRITE(2,*)(H_S(I),I=1,128)
WRITE(1,*)(H(I),I=1,128)
CLOSE(1)
CLOSE(2)
END
CCCCCCCCCCCCC FFT 子程序CCCCCCCCCCCCCC SUBROUTINE FFT(LX,CX,SIGNI)
COMPLEX CX(LX),CARG,CEXP,CW,CTEMP
J=1
SC=1.0/LX
IF(SIGNI.EQ.1.0)SC=1.0
SIG=-SIGNI
DO 30 I=1,LX
IF(I.GT.J)GO TO 10
CTEMP=CX(J)*SC
CX(J)=CX(I)*SC
CX(I)=CTEMP
10 M=LX/2
20 IF(J.LE.M)GO TO 30
J=J-M
M=M/2
IF(M.GE.1)GO TO 20
30 J=J+M
L=1
40 ISTEP=2*L
DO 50 M=1,L
CARG=(0.0,1.0)*(3.14159265*SIG*(M-1))/L
CW=CEXP(CARG)
DO 50 I=M,LX,ISTEP
CTEMP=CW*CX(I+L)
CX(I+L)=CX(I)-CTEMP
50 CX(I)=CX(I)+CTEMP
L=ISTEP
IF(L.LT.LX)GO TO 40
RETURN
END
5.频谱分析
CCCCCCCCCCCCC 频谱分析CCCCCCCCCCCCCC
PROGRAM MAIN
PARAMETER (PI=3.141592654,DT=0.004)
REAL XSIN(200),X(200),WA VE(200)
REAL X_SIN_R(200),X_SIN_I(200)
REAL X_R(200),X_I(200)
REAL X_WA VE_R(200),X_WA VE_I(200)
REAL SINSPRT(200),XSPRT(200),WA VESPRT(200)
REAL DF
COMPLEX XSIN_C(200),X_C(200),X_WA VE_C(200)
CCCCCCCCCCCC 对SIN 的分析CCCCCCCCCCCC
DO 100 I=1,128
XSIN(I)=SIN(2*I*PI/64.0)
100 CONTINUE
CCCCCCCCC 是数据转换成复数形式CCCCCCCCCC
DO 200 I=1,128
200 XSIN_C(I)=CMPLX(XSIN(I),0.0)
CALL FFT(128,XSIN_C,1.0)
DO 300 I=1,128
300 X_SIN_R(I)=REAL(XSIN_C(I))
DO 400 I=1,128
400 X_SIN_I(I)=AIMAG(XSIN_C(I))
OPEN(1,FILE='XSINSPRT.DAT',FORM='FORMATTED',STATUS='UNKNOWN') DO 1 I=1,128
1 SINSPRT(I)=SQRT(X_SIN_R(I)**2+X_SIN_I(I)**2)
WRITE(1,*)(SINSPRT(I), I=1,128)
CLOSE(1)
CCCCCCCCCCCCCC 对脉冲信号的分析CCCCCCCCCCCCC
X(1)=1.0
DO 800 I=2,128
800 X(I)=0.0
DO 900 I=1,128
900 X_C(I)=CMPLX(X(I),0.0)
CALL FFT(128,X_C,1.0)
DO 1000 I=1,128
1000 X_R(I)=REAL(X_C(I))
DO 1100 I=1,128
1100 X_I(I)=AIMAG(X_C(I))
OPEN(5,FILE='XSPRT.DAT',FORM='FORMATTED',STATUS='UNKNOWN')
DO 2 I=1,128
2 XSPRT(I)=SQRT(X_R(I)**2+X_I(I)**2)
WRITE(5,*)(XSPRT(I), I=1,128)
CLOSE(5)
CCCCCCCCCCCC 对地震波W A VE 的分析CCCCCCCCCCC
OPEN(3,FILE='WA VE.DAT',FORM='FORMATTED',STATUS='UNKNOWN')
READ(3,*)(WA VE(I), I=1,128)
DO 500 I=1,128
500 X_WA VE_C(I)=CMPLX(WA VE(I),0.0)
CALL FFT(128,X_WA VE_C,1.0)
DO 600 I=1,128
600 X_WA VE_R(I)=REAL(X_WA VE_C(I))
DO 700 I=1,128
700 X_WA VE_I(I)=AIMAG(X_WA VE_C(I))
OPEN(4,FILE='WA VESPRT.DAT',FORM='FORMATTED',STATUS='UNKNOWN') DO 3 I=1,128
3 WA VESPRT(I)=SQRT(X_WA VE_R(I)**2+X_WA VE_I(I)**2)
WRITE(4,*)(WA VESPRT(I), I=1,128)
CLOSE(4)
CLOSE(3)
END
CCCCCCCCCCCCCCC FFT 子程序CCCCCCCCCCCCCCCCC
SUBROUTINE FFT(LX,CX,SIGNI)
COMPLEX CX(LX),CARG,CEXP,CW,CTEMP
J=1
SC=1.0/LX
IF(SIGNI.EQ.1.0)SC=1.0
SIG=-SIGNI
DO 30 I=1,LX
IF(I.GT.J)GO TO 10
CTEMP=CX(J)*SC
CX(J)=CX(I)*SC
CX(I)=CTEMP
10 M=LX/2
20 IF(J.LE.M)GO TO 30
J=J-M
M=M/2
IF(M.GE.1)GO TO 20
30 J=J+M
L=1
40 ISTEP=2*L
DO 50 M=1,L
CARG=(0.0,1.0)*(3.14159265*SIG*(M-1))/L
CW=CEXP(CARG)
DO 50 I=M,LX,ISTEP
CTEMP=CW*CX(I+L)
CX(I+L)=CX(I)-CTEMP
50 CX(I)=CX(I)+CTEMP
L=ISTEP
IF(L.LT.LX)GO TO 40
RETURN
END
6.分析补零对振幅谱的影响
CCCCCCCCCCC 分析补零对振幅谱的影响CCCCCCCCCCC
PROGRAM MAIN
CCCCCCCCCCC 不补零时的振幅普CCCCCCCCCCCCCCCCC CCCCCCCCCCC 对函数SIN 的分析CCCCCCCCCCCCCCCC
PARAMETER (PI=3.1415)
REAL XSIN(150),XWA VE(150)
REAL XWA VE_R(150),XWA VE_I(150)
REAL XSIN_Z(150),XWA VE_Z(150)
REAL XSIN_R(150),XSIN_I(150)
COMPLEX XSIN_60C(150),XWA VE_60C(150)
COMPLEX XSIN_64C(150),XWA VE_64C(150)
COMPLEX XSIN_128C(150),XWA VE_128C(150)
REAL XSIN_60S(150),XWA VE_60ZS(150)
REAL XSIN_64S(150),XWA VE_64ZS(150)
REAL XSIN_128S(150),XWA VE_128ZS(150)
OPEN(1,FILE='XSIN_60S.DA T',FORM='FORMATTED',STATUS='UNKNOWN') DO 1 I=1,60
1 XSIN(I)=SIN(2*PI*I/30.0)
DO 11 I=1,60
11 XSIN_60C(I)=CMPLX(XSIN(I),0.0)
CALL DFT(60,XSIN_60C,1.0)
DO 12 I=1,60
12 XSIN_R(I)=REAL(XSIN_60C(I))
DO 13 I=1,60
13 XSIN_I(I)=AIMAG(XSIN_60C(I))
DO 14 I=1,60
14 XSIN_60S(I)=SQRT(XSIN_R(I)**2+XSIN_I(I)**2)
WRITE(1,*)(XSIN_60S(I),I=1,60)
CLOSE(1)
CCCCCCCCCCCCC 对地震数据的分析CCCCCCCCCCCCCC
OPEN(2,FILE='WA VE.DAT',FORM='FORMATTED',STATUS='UNKNOWN')
READ(2,*)(XWA VE(I),I=1,60)
CLOSE(2)
DO 21 I=1,60
21 XWA VE_60C(I)=CMPLX(XWA VE(I),0.0)
CALL DFT(60,XWA VE_60C,1.0)
DO 22 I=1,60
22 XWA VE_R(I)=REAL(XWA VE_60C(I))
DO 23 I=1,60
23 XWA VE_I(I)=AIMAG(XWA VE_60C(I))
DO 24 I=1,60
24 XWA VE_60ZS(I)=SQRT(XWA VE_R(I)**2+XWA VE_I(I)**2)
OPEN(3,FILE='XWA VE_60ZS.DA T',FORM='FORMATTED',STA TUS='UNKNOWN') WRITE(3,*)(XWA VE_60ZS(I),I=1,60)
CLOSE(3)
CCCCCCCCCCCCC 对补零数据的分析CCCCCCCCCCCCCCCC CCCCCCCCCCCC 对函数SIN 的分析CCCCCCCCCCCCCCCC
OPEN(4,FILE='XSIN_64S.DA T',FORM='FORMATTED',STATUS='UNKNOWN')
DO 31 I=61,64
31 XSIN(I)=0.0
DO 32 I=1,64
32 XSIN_64C(I)=CMPLX(XSIN(I),0.0)
CALL DFT(64,XSIN_64C,1.0)
DO 33 I=1,64
33 XSIN_R(I)=REAL(XSIN_64C(I))
DO 34 I=1,64
34 XSIN_I(I)=AIMAG(XSIN_64C(I))
DO 35 I=1,64
35 XSIN_64S(I)=SQRT(XSIN_R(I)**2+XSIN_I(I)**2)
WRITE(4,*)(XSIN_64S(I),I=1,64)
CLOSE(4)
OPEN(5,FILE='XSIN_128S.DA T',FORM='FORMATTED',STATUS='UNKNOWN')
DO 41 I=65,128
41 XSIN(I)=0.0
DO 42 I=1,128
42 XSIN_128C(I)=CMPLX(XSIN(I),0.0)
CALL DFT(128,XSIN_128C,1.0)
DO 43 I=1,128
43 XSIN_R(I)=REAL(XSIN_128C(I))
DO 44 I=1,128
44 XSIN_I(I)=AIMAG(XSIN_128C(I))
DO 45 I=1,128
45 XSIN_128S(I)=SQRT(XSIN_R(I)**2+XSIN_I(I)**2)
WRITE(5,*)(XSIN_128S(I),I=1,128)
CLOSE(5)
CCCCCCCCCCCCC 对地震数据的分析CCCCCCCCCCCCCC
OPEN(6,FILE='XWA VE_64ZS.DAT',FORM='FORMATTED',STATUS='UNKNOWN') DO 51 I=61,64
51 XWA VE(I)=0.0
DO 52 I=1,64
52 XWA VE_64C(I)=CMPLX(XWA VE(I),0.0)
CALL DFT(64,XWA VE_64C,1.0)
DO 53 I=1,64
53 XWA VE_R(I)=REAL(XWA VE_64C(I))
DO 54 I=1,64
54 XWA VE_I(I)=AIMAG(XWA VE_64C(I))
DO 55 I=1,64
55 XWA VE_64ZS(I)=SQRT(XWA VE_R(I)**2+XWA VE_I(I)**2)
WRITE(6,*)(XWA VE_64ZS(I),I=1,64)
CLOSE(6)
OPEN(7,FILE='XWA VE_128ZS.DAT',FORM='FORMATTED',STATUS='UNKNOWN') DO 61 I=65,128
61 XWA VE(I)=0.0
DO 62 I=1,128
62 XWA VE_128C(I)=CMPLX(XWA VE(I),0.0)
CALL DFT(128,XWA VE_128C,1.0)
DO 63 I=1,128
63 XWA VE_R(I)=REAL(XWA VE_128C(I))
DO 64 I=1,128
64 XWA VE_I(I)=AIMAG(XWA VE_128C(I))
DO 65 I=1,128
65 XWA VE_128ZS(I)=SQRT(XWA VE_R(I)**2+XWA VE_I(I)**2)
WRITE(7,*)(XWA VE_128ZS(I),I=1,128)
CLOSE(7)
END
CCCCCCCCCCCCCCC DFT 子程序CCCCCCCCCCCCCCCC
SUBROUTINE DFT(N,C,S)
COMPLEX C(N),Y(512),W
IF(S.EQ.1.0)Z=1.0
IF(S.EQ.-1.0)Z=1.0/FLOAT(N)
DO 20 I=1,N
Y(I)=(0.0,0.0)
DO 20 J=1,N
A=COS(2.0*3.14159265*(I-1)*(J-1)/FLOAT(N))
B=SIN(2.0*3.14159265*(I-1)*(J-1)/FLOAT(N))*(-S)
W=CMPLX(A,B)
20 Y(I)=Y(I)+C(J)*W
DO 30 I=1,N
30 C(I)=Y(I)*Z
RETURN
END
7.线性褶积与循环褶积
CCCCCCCCC 线性滤波与循环褶积比较CCCCCCCCC
PROGRAM MAIN
CCCCCCCCCCCCC 线性滤波CCCCCCCCCCCCCCC
CCCCCCCCCCCC H是低通滤波因CCCCCCCCCCCC
PARAMETER (PI=3.141592654, DT=0.002,F=70.0)
REAL X(100),X2(101),H(-50:50),Y_LINEAR(200),Y_CIRCLE(200)
OPEN(1,FILE='INPUT1.DA T',FORM='FORMATTED',STATUS='UNKNOWN')
READ(1,*)(X(I),I=1,100)
CCCCCCCCCCCCCCC低通滤波器CCCCCCCCCCCCC
DO 10 I=-50,50
IF (I.EQ.0)THEN
H(I)=2*F*PI/PI
ELSE
H(I)=SIN(2*PI*F*I*DT)/(PI*I*DT)
END IF
10 CONTINUE
CALL CON(X,H,Y_LINEAR,100,101,200)
CCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCC
OPEN(2,FILE='Y_LINEAR.DA T',FORM='FORMATTED',STATUS='UNKNOWN')
WRITE(2,*)(Y_LINEAR(I),I=1,200)
CLOSE(2)
CCCCCCCCCCCCC 圆周滤波的设计CCCCCCCCCCCC
OPEN(3,FILE='Y_CIRCLE.DAT',FORM='FORMATTED',STATUS='UNKNOWN')
DO 20 I=1,100
20 X2(I)=X(I)
X2(100)=0.0
CALL CIRCON(X2,H,Y_CIRCLE,101)
WRITE(3,*)(Y_CIRCLE(I),I=1,101)
CLOSE(3)
CCCCCCC 线性卷积输出点等于圆周卷积CCCCCC
OPEN(4,FILE='Y_LINEARCUT.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(4,*)(Y_LINEAR(I),I=1,101)
CLOSE(4)
END
CCCCCCCCCCCCCC 线性卷积子程序CCCCCCCCCCCCCC
SUBROUTINE CON(A,B,C,I,J,K)
DIMENSION A(I),B(J),C(K)
DO 1 K1=1,K
1 C(K1)=0.0
DO 2 I1=1,I
DO 2 I2=1,J
II=I1+I2-1
2 C(II)=C(II)+A(I1)*B(I2)*0.004
RETURN
END
CCCCCCCCCCCCCCC 圆周卷积子程序CCCCCCCCCCCCCC
SUBROUTINE CIRCON(X,H,Y,N)
REAL X(N),H(N),Y(N)
DO I=1,N
Y(I)=0.0
ENDDO
DO I=1,N
DO J=1,N
IF(I-J.GT.0) Y(I)=Y(I)+X(J)*H(I-J)
ENDDO
ENDDO
END
8.最小平方反滤波
CCCCCCCCCC 两种反射系数序列的对比CCCCCCCCCC
PROGRAM MAIN
REAL B(12),R(200),X_CONV(211),RXX(211)
REAL A(211),A_CONV(270)
DATA B/17.5,12.5,7.0,0.0,-3.5,-7.0,-3.0,-1.0,0.0,1.4,3.5,2.0/
OPEN(1,FILE='INPUT8.DAT',FORM='FORMATTED',STATUS='UNKNOWN') READ(1,*)(R(I),I=1,200)
CCCCCCCCCCCCCC 求合成地震记录CCCCCCCCCCCCCC
CALL CON(R,B,X_CONV,200,12,211)
CCCCCCCCCCCC 求地震子波的自相关CCCCCCCCCCCC
CALL AUTCOR(211,211,X_CONV,RXX)
CCCCCCCCCCCCCCC 求反滤波因子CCCCCCCCCCCCCCC
CALL PEO(211,RXX,A)
OPEN(2,FILE='A.DAT',FORM='FORMA TTED',STATUS='UNKNOWN')
WRITE(2,*)(A(I),I=1,60)
CCCCCCC 通过反滤波因子求反射系数序列CCCCCCC
CALL CON(A,X_CONV,A_CONV,60,211,270)
OPEN(3,FILE='A_CONV.DA T',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(3,*)(A_CONV(I),I=1,270)
CLOSE(3)
CLOSE(2)
CLOSE(1)
END
CCCCCCCCCCCCCCC 褶积子程序CCCCCCCCCCCCCCCC SUBROUTINE CON(A,B,C,I,J,K)
DIMENSION A(I),B(J),C(K)
DO 1 K1=1,K
1 C(K1)=0.0
DO 2 I1=1,I
DO 2 I2=1,J
II=I1+I2-1
2 C(II)=C(II)+A(I1)*B(I2)*0.004
RETURN
END
CCCCCCCCCCCCCCC 自相关子程序CCCCCCCCCCCCCCC SUBROUTINE AUTCOR(J,K,C,A)
DIMENSION C(J),A(K)
DO 7 N=1,K
A(N)=0.0
I=N
DO 6 IN=I,J
IL=IN-N+1
6 A(N)=A(N)+C(IN)*C(IL)
7 CONTINUE
RETURN
END
CCCCCCCCCCCCCC 求反滤波因子CCCCCCCCCCCCCCCC SUBROUTINE PEO(LR,R,A)
C PEO IS THE AUXILIARY TOEPLITZ RECURSION
C IT GIVES THE PREDICTION-ERROR OPERA TOR A(I)
DIMENSION R(LR),A(LR)
V=R(1)
A(1)=1.0
IF(LR.EQ.1)RETURN
DO 6 L=2,LR
D=0.0
L3=L-1
DO 1 J=1,L3
K=L-J+1
1 D=D+A(J)*R(K)
C=D/V
IF(L.EQ.2)GOTO 4
L1=(L-2)/2
L2=L1+1
IF(L2.LT.2)GOTO 3
DO 2 J=2,L2
HOLD=A(J)
K=L-J+1
A(J)=A(J)-C*A(K)
2 A(K)=A(K)-C*HOLD
IF(2*L1.EQ.L-2)GOTO 4
3 LT3=L2+1
A(LT3)=A(LT3)-C*A(LT3)
4 A(L)=-C
6 V=V-C*D
RETURN
END
9.零相位转换
CCCCCCCCCCC 零相位的转换CCCCCCCCCCCCC
PROGRAM MAIN
REAL B(32),B1(32),BC_S(32),B_CONV(32)
REAL BC_R(32),BC_I(32)
DATA B/0,1,3,5,7,9,10,8,6,4,2,0,-1,0,1,2,3,1.5,0,-1,0,2,1,0,0/
COMPLEX B_C(32),BC_SZ(32)
DO 10,I=1,32
IF(I.GE.26)THEN
B_C(I)=(0.0,0.0)
ELSE
B_C(I)=CMPLX(B(I),0.0)
END IF
10 CONTINUE
CALL FFT(32,B_C,1.0)
DO 20,I=1,32
BC_R(I)=REAL(B_C(I))
BC_I(I)=AIMAG(B_C(I))
BC_S(I)=SQRT(BC_R(I)**2+BC_I(I)**2)
20 BC_SZ(I)=CMPLX(BC_S(I),0.0)
CALL FFT(32,BC_SZ,-1.0)
DO 30,I=1,32
30 B_CONV(I)=REAL(BC_SZ(I))
OPEN (1,FILE='B_CONV.DA T',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(1,*)(B_CONV(I),I=1,32)
CLOSE(1)
END
CCCCCCCCCCCCCC FFT子程序CCCCCCCCCCCCCC
SUBROUTINE FFT(LX,CX,SIGNI)
COMPLEX CX(LX),CARG,CEXP,CW,CTEMP
J=1
SC=1.0/LX
IF(SIGNI.EQ.1.0)SC=1.0
SIG=-SIGNI
DO 30 I=1,LX
IF(I.GT.J)GO TO 10
CTEMP=CX(J)*SC
CX(J)=CX(I)*SC
CX(I)=CTEMP
10 M=LX/2
20 IF(J.LE.M)GO TO 30
J=J-M
M=M/2
IF(M.GE.1)GO TO 20
30 J=J+M
L=1
40 ISTEP=2*L
DO 50 M=1,L
CARG=(0.0,1.0)*(3.14159265*SIG*(M-1))/L
CW=CEXP(CARG)
DO 50 I=M,LX,ISTEP
CTEMP=CW*CX(I+L)
CX(I+L)=CX(I)-CTEMP
50 CX(I)=CX(I)+CTEMP
L=ISTEP
IF(L.LT.LX)GO TO 40
RETURN
END
10.最小相位转换
CCCCCCCCCCCC 最小相位转换CCCCCCCCCCCC
PROGRAM MAIM
REAL B(25),B0(75)
DATA B/0,1,2,3,1.5,0,-1,0,1,3,5,7,9,10,8,6,4,2,0,-1,0,2,3,2,0/
REAL RBB(25),Y0(25),A0(25)
REAL SUMX,C(49),G(49)
OPEN(1,FILE='B0.DA T',FORM='FORMATTED',STA TUS='UNKNOWN') OPEN(2,FILE='B.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(2,*)(B(I),I=1,25)
CALL AUTCOR(25,25,B,RBB)
CALL PEO(25,RBB,Y0)
SUMX=0.0
DO 10,I=1,25
CALL CON(B,Y0,C,25,25,49)
10 SUMX=SUMX+(C(I))**2
B0(1)=1.0/SQRT(SUMX)
地震资料解释期末复习（王松版） 1地震资料解释——以地质理论和规律为指导，运用地震波传播理论和地震勘探方法原理，综合地质、测井、钻井和其它物探资料，对地震数据进行深入研究、综合分析的过程。 2地震子波（wavelet）：地震勘探过程中，爆炸产生的尖脉冲传播到一定距离时波形逐渐稳定。 3褶积模型的应用： 已知r(t)和w(t)，求s(t)：正演问题 已知w(t) 和s(t) ，求r(t) ：反演问题 已知s(t) 和r(t)，求w(t)：子波处理 4同相轴：指地震时间剖面上相同相位的连接线 5极性判断 6有效波的识别标志 1）强振幅： 叠后资料往往经提高信噪处理，反射波能量大于干扰波能量 2）波形相似性： 子波相同、同一界面反射波传播路径相近，传播过程影响因素相近，相邻地震道上的波形特征（主周期、相位数、振幅包络形状等）是相似的。 3）同相性： 同一个反射波的相同相位，在相邻地震道上的到达时间也是相近的，每道记录下来的振动图是相似的，形成一条平滑的、有一定长度的同相轴，也称相干性。 4）时差变化规律： 在共炮点道集上，直达波、折射波是直线，反射波、绕射波、多次波等为曲线。在动校正后的剖面上，原来直线的同相轴被校正成曲线，一次反射波成为直线，多次波、绕射波为曲线。 1、2用于识别波的出现； 3、4用于识别波的类型、特征及地层界面特征的判断。 7水平叠加剖面的特点 （1）在测线上同一点,根据钻井资料得到的地质剖面上的地层分界面,与时间剖面上的反射波同相轴在数量上、出现位置上,常常不是一一对应的。 （2）时间剖面的纵坐标是双程旅行时t0 ，而地质剖面或测井资料是以铅垂深度表示的,两者需经时深转换,其媒介就是地震波的传播速度,它通常随深度或空间而变化。 （3）反射波振幅、同相轴及波形本身包含了地下地层的构造和岩性信息,如振幅的强弱与地层结构、介质参数密切相关。但是反射波同相轴是与地下的分界面相对应,同相轴与界面两侧的地层、岩性有关。必须经过一些特殊处理（如声阻抗反演技术等）才能把反射波所包含的“界面”的信息转换成为与“层”有关的信息后,才能与地质和钻井资料进行直接地对比。 （4）地震剖面上的反射波是由多个地层分界面上振幅有大有小、极性有正有负、
安徽理工大学 一、名词解释（20分） 1、、地震资料数字处理:就是利用数字计算机对野外地震勘探所获得的原始资料进行加工、改进，以期得到高质量的、可靠的地震信息，为下一步资料解释提供可靠的依据和有关的地质信息。 2、数字滤波:用电子计算机整理地震勘探资料时，通过褶积的数学处理过程，在时间域内实现对地震信号的滤波作用，称为数字滤波。（对离散化后的信号进行的滤波，输入输出都是离散信号） 3、模拟信号：随时间连续变化的信号。 4、数字信号：模拟数据经量化后得到的离散的值。 5、尼奎斯特频率：使离散时间序列x(nΔt)能够确定时间函数x(t)所对应的两倍采样间隔的倒数，即f=1/2Δt. 6、采样定理： 7、吉卜斯现象：由于频率响应不连续，而时域滤波因子取有限长，造成频率特性曲线倾斜和波动的现象。 8、假频:抽样数据产生的频率上的混淆。某一频率的输入信号每个周期的抽样数少于两个时，在系统的的输出端就会被看作是另一频率信号的抽样。抽样频率的一半叫作褶叠频率或尼奎斯特频率fN；大于尼奎斯特频率的频率fN+Y，会被看作小于它的频率fN-Y。这两个频率fN+Y和fN-Y相互成为假频。 9、伪门：对连续的滤波因子h(t)用时间采样间隔Δt离散采样后得到h (nΔt)。如果再按h (nΔt)计算出与它相应的滤波器的频率特性，这时在频率特性图形上，除了有同原来的H (ω)对应的'门'外，还会周期性地重复出现许多门，这些门称为伪门。产生伪门的原因就是由于对h(t)离散采样造成的。 10、地震子波：由于大地滤波作用，使震源发出的尖脉冲经过地层后，变成一个具有一定时间延续的波形w（t）。 11、道平衡：指在不同的地震记录道间和同一地震记录道德不同层位中建立振幅平衡，前者称为道间均衡，后者称为道内均衡。 12、几何扩散校正：球面波在传播过程中，由于波前面不断扩大，使振幅随距离呈反比衰减，即Ar=A0/r，是一种几何原因造成的某处能量的减小，与介质无关，叫几何扩散，又叫球面扩散。为了消除球面扩散的影响，只需A0=Ar*r即可，此即为几何扩散校正， 13、反滤波(又称反褶积):为了从与干扰混杂的地震讯息中把有效波提取出来，则必须设法消除由于水层、地层等所形成的滤波作用，按照这种思路所提出的消除干扰的办法称为反滤波，即把有效波在传播过程中所经受的种种我们不希望的滤波作用消除掉。 14、校正不足或欠校正：如果动校正采用的速度高于正确速度，计算得到的动校正量偏小，动校正后的同相轴下拉。反之称为校正过量或过校正。 15、动校正:消除由于接受点偏离炮点所引起的时差的过程，又叫正常时差校正。 16、剩余时差：当采用一次波的正常时差公式进行动校正之后，除了一次反射波之外，其他类型的波仍存在一定量的时差，我们将这种进过动校正后残留的时差叫做剩余时差。
地震勘探在海洋石油勘探中的基本原理
————————————————————————————————作者: ————————————————————————————————日期:
本科生课外研学任务书及成绩评定表 题目＿_地震勘探在海洋石油勘探中的基本原理学生姓名＿___ 黄邦毅_________＿＿__＿__ 指导教师＿＿＿_ 严家斌___＿_＿＿＿__＿_ 学院_＿＿_ 地信院_＿＿_____________ 专业班级_＿＿地科0９0１__________＿____
地震勘探在海洋石油勘探中的基本原理 一、引言 国内外的勘探实践表明,没有物探技术的进步，就没有更多圈闭的发现,就没有钻探成功率的提高,也就更不会有油田和储产量的快速增长。宏观看，物探的作用在勘探阶段是客观的目标评价,在开发阶段是精细的油藏描述。因此，油气勘探开发离不开地震技术和地震技术的进步与发展。如果说勘探技术是石油工业的第一生产力，那么物探技术就是获得油气储量的第一直接生产力。 纵观近些年的勘探技术的具体运用，最常见的莫过于地震勘探，所谓地震勘探就是通过人工方法激发地震波,研究地震波在地层中传播的情形,以查明地下的地质构造,为寻找油气田或其它勘探目的服务的一种物探方法! 21世纪是海洋的世纪,海洋蕴藏着很多宝贵的资源，随着生产技术的日趋进步，世界各国(包括中国在内)目前都在积极寻求开发海洋资源，在海洋的勘探开发中离不开物探，而且运用最广泛也最有效的是地震勘探。 二、海洋地震勘探 在茫茫大海里寻找石油最有效的技术方法是地球物理方法,其中主要是地震勘探方法。近几十年来,随着电子计算机的广泛应用,海洋地震勘探的数据采集和装备得到了极大的改进，数据处理技术和解释方法也得到迅速的发展。在油气勘探中,利用地震资料不仅能确定地下的构造形态、断裂分布，而且能了解地层岩性、储层厚度、储层参数甚至能直接指示地下油气的存在。在油气开发中，地震资料同测井、岩芯资料以及其它地下地质资料相结合能对油藏进行描述和监测。地震技术远远超出了石油勘探领域，已向石油开发和生产领域渗透。 用于寻找海上石油的地震反射法，和陆地的地震反射法相比,在方法基本原理、资料处理和解释方法等方面基本上是一样的。其中， 测量原理 在这类方法中，地震波在介质中传播的物理模型如图1所示。从震源O激发出的弹性波投射到反射界面上产生反射波,其条件是:入射角α等于反射角β。能
中国地质大学（武汉）地空学院 地震实验报告 姓名：沈 班级：班 学号： 时间： 2015年05月 指导老师：张
一、实验目的 实验一： 1、浅层地震装备的基本组成； 2、认识GEODE96浅层地震仪的主要结构，并学会该类仪器的操作方法； 3、地震波认识。 实验二： 1、掌握浅层地震数据采集方法及注意事项 二、仪器介绍 1、仪器简介 全套美国GEOMETRICS公司生产的Geode96浅层地震仪（相当于四套独立的24道浅层地震仪）该仪器能满足折、反射地震勘探、井间勘探、面波调查等地震监测需要,应用Crystal公司的A/D转换器和高速过采样技术达到了24位地震仪的精度。频带从1.75Hz到20,000Hz,使得采样间隔可以从20毫秒到16微秒。采样到的数据叠加到32位的叠加器中，然后传回到主机的硬盘或其它介质上。内置预触发器，每道有16K的内存。用硬件相关器对震源信号进行实时相关运算。Geode包装坚固、防水、防震,有提手,重4.1公斤,用12V的外接电池可以连续工作10个小时。(如下图)
2、主要操作功能键及快捷键 注释： 1锁定与解锁；2清除界面；4检测噪声；7保存 3、操作步骤及注意事项 1、每个GEODE用数传线按规定串联，通过数传盒与笔记本电脑的USB口连 接。 2、每个GEODE接上12V电源。 3、开关接到与笔记本相连的第一个GEODE上。 4、传盒上的开关置于POWER UP处。 5、采集控制程序，并按工作需要设置好各项参数，然后进行正常数据采集工作。 6、出采集控制程序之前，应将数传盒上的开关置于POWER DOWN处。 7、卸下各连接线并清理整齐。 8、注意的是：在正常工作过程中，任何时候移动数传线与GEODE的连接头时，必须退出采集控制程序。另外Y型头上有红色标记的与GEODE的前12道相连接。而且采集控制软件运行的语言环境必须是英语(美国)。
地震解释基础 复习题 1.为什么并非每一个地质界面都对应一个反射同相轴？ 子波有一定的延续长度，若地层很薄，相邻分界面的信号可能会重叠到一起形成复合波，导致无法分辨界面。所以一个反射同相轴可能包含多个地质界面。 2.影响地震资料纵向分辨率的因素有哪些？提高分辨率的实质是什么？ 1）激发条件——激发宽频带子波——井深、药量、激发岩性、虚反射、激发组合 2）接收条件——检波器类型、地表岩性、检波器耦合、组合方式、仪器响应 3）近地表低降速带的影响 4）大地滤波作用、地层速度 实质：提高主频，拓宽频带 3.提高横向分辨率的方法是什么？为什么它能提高横向分辨率？ 偏移是提高地震勘探横向分辨率的根本方法 提高横向分辨率的核心是减小菲涅尔带的大小, 菲涅尔带的极限 : 要想减小菲涅尔带的大小就要减小h ，偏移将地表向下延拓到地下界面，使h=0，所以 菲涅尔带减小到极限L=λ/4,所以偏移能提高横向分辨率。 4.地震剖面的对比方法 1）掌握地质规律、统观全局 在对比之前，要收集和分析勘探区的各种资料。研究规律性的地质构造特征，用地质规律指导对比解释。了解地震资料采集和处理的方法及相关因素，以便准确识别和判断出剖面假象。 2）从主测线开始对比 在一个工区有多条地震剖面，应先从主测线开始对比工作，然后从主测线的反射层延伸到其他测线上去。（主测线：指垂直构造走向、横穿主要构造，并且信噪比高、反射同相轴连续性好的测线。它还应有一定的延伸长度，最好能经过钻探井位。） 3）重点对比标准层 对某条测线而言，可能有几个反射层，应重点对比目标层（或称为标准层，标准层：具有较强振幅、连续性较好、可在整个工区内追踪的目标反射层。它往往是主要的地层或岩性的分界面，与生油层或储集层有一定的关系，或本身就为生油层、储油层）。 4）相位对比 反射波的初至难以辨认，采用相位对比。若选振幅最强、连续性最好的某同相轴进行追()222042164h L O C h h h λλλλ??'==+-=+== ???
地震资料解释课程教学大纲 课程代码：74190110 课程中文名称：地震资料解释 课程英文名称：Seismic Interpretation 学分：2.0 周学时：1.5-1.0 面向对象： 预修要求：地层学、构造地质学、海洋沉积学、地球海洋物理学 一、课程介绍 （一）中文简介 《地震资料解释》是海洋科学专业的一门专业必修课，其总目标是结合地震资料解释实习课，使学生能够理解地震资料解释的基本原理和概念、掌握复杂地质条件下的层序地层、构造和地震相分析等地震资料解释的基本方法。 （二）英文简介 “Seismic Interpretation” is a compulsory course for the students majored in Marine Science. In combination with associated practice course, the students who attend this course would: (1) understand the fundamentals and basic concepts in interpreting the seismic data; (2) master basic skills and methodology to analyze the sequence stratigraphy, structure and seismic facies in the subsurface with complex geological conditions. 二、教学目标 （一）学习目标 通过本课程系统学习，要求学生全面掌握地震地质解释的地球物理基础和地震地质解释方法；学会应用地震资料进行地质解释的技能；了解地震资料地质解释的现状及发展方向。 （二）可测量结果 （1）掌握沉积层序的概念、沉积层序的边界类型、层序划分的原则和方法，能在地震剖面
Vista 5.5的基本使用方法 数据输入 地震分析窗口 一维频谱 二维频波谱 观测系统 工作流 一、数据输入 1.1 把数据文件加入Project 首先选择File/New Project，新建一个Project，按住不放，出现按钮组合，可以选择不同类型 的数据集，选择，向Project中增加一个新的2-D数据集，按住不放，出现按钮组合， 可以选择加入不同类型的地震数据，选择，选择一个SEG-Y数据，即可将该数据文件加入新建的数据集。 1.2 命令流中数据的输入 双击进入如下界面 1.2.1 Input Data List 数据输入列表，选择已加入到Project的数据集，下面的文本框中会显示选择的数据的基本信息。 1.2.2 Data Order 选择输入数据的排列方式，对不同的处理步骤可以选择不同的数据排列方式 Sort Order a. NO SORT ORDER 输入数据原始排列方式 b. SHOT_POINT_NO 输入数据按炮点排列方式 c. FIELD_STATION_NUMBER d. CMP_NO 输入数据按共中心点排列方式 e. FIELD_STATION_NUMBER 1.2.3 Data Input Control 数据输入控制 右键-->Data Input Control a. Data Input 进入Flow Input Command（见上） b. Data Sort List 查看数据排列方式的种类 c. Data/header Selection 输入数据的选择，可以控制输入数据的道数和CMP道集 查看所有已经选择的数据 如果没有定义任何可选的数据信息，则如下图所示： 可以选择一种选择方式，单击并设置选择信息。定义有可选的数据信息后，在查看，则如下图所示，会显示选择的信息。 选择共炮点集 单击后，会弹出如下界面：
三维地震勘探技术及其应用 [摘要] 本文应用三维地震勘探技术对某矿南三采区进行探测，探测区内解释断层71条，其中可靠断层61条，较可靠断层10条，31个无煤带。为煤矿安全生产提供了科学依据,节约了生产成本的投入。 [关键词] 三维地震采区 [abstract] this paper introduces the application of three dimensional seismic exploration method on the south third mining area of a certain coal mine. 71 faults were showed in this exploration area, in which there are 61 reliable faults, 10 relatively reliable faults and 31 areas without any coal. those information provides scientific foundation for the production safty of the coal mine and saves the cost. [key words] three dimensional seismic mining area 0.引言 随着煤炭地震勘探技术的提高，尤其是九十年代以来三维地震勘探在煤炭系统的应用与推广，三维地震勘探技术在煤矿采区进行小构造勘探成为现实，给煤矿建设和生产带来了巨大的效益。 近年来，随着我国煤炭资源勘查理论和技术的不断发展，已形成了中国煤炭地质综合勘查理论与技术新体系，其中三维地震勘探技术是五大关键技术之一。[1]
三维地震解释课程设计 一、摘要，关键词---------------------------------------------------1 二、引言----------------------------------------------------------------2 三、设计部分----------------------------------------------------------4 1、三维地震解释的方法与原理---------------------------------------------4 2、三维地震解释的工作流程------------------------------------------------6 3、North sea 三维地震资料解释--------------------------------7 四、结束语------------------------------------------------------------14 五、参考文献---------------------------------------------------------15
摘要 利用垂直时间剖面、水平切片等多视角、全方位的三维可视化手段对采区高分辨率三维地震数据进行解释，为精细地质构造的解释提供了有效地手段。 Summary Using vertical time profile、horizontal slice multi-angel、omnibearing 3d visual device to interpret high-resolution 3d seismic data in mining area provides effective device to interpret elaborate geological structure. 关键词 三维地震勘探垂直时间剖面顺层切片水平切片
地震勘探实验报告记录
————————————————————————————————作者：————————————————————————————————日期：
中国地质大学（武汉）地空学院 地震实验报告 姓名：沈 班级：班 学号： 时间： 2015年05月 指导老师：张
一、实验目的 实验一： 1、浅层地震装备的基本组成； 2、认识GEODE96浅层地震仪的主要结构，并学会该类仪器的操作方法； 3、地震波认识。 实验二： 1、掌握浅层地震数据采集方法及注意事项 二、仪器介绍 1、仪器简介 全套美国GEOMETRICS公司生产的Geode96浅层地震仪（相当于四套独立的24道浅层地震仪）该仪器能满足折、反射地震勘探、井间勘探、面波调查等地震监测需要,应用Crystal公司的A/D转换器和高速过采样技术达到了24位地震仪的精度。频带从1.75Hz到20,000Hz,使得采样间隔可以从20毫秒到16微秒。采样到的数据叠加到32位的叠加器中，然后传回到主机的硬盘或其它介质上。内置预触发器，每道有16K的内存。用硬件相关器对震源信号进行实时相关运算。Geode包装坚固、防水、防震,有提手,重4.1公斤,用12V的外接电池可以连续工作10个小时。(如下图)
2、主要操作功能键及快捷键 注释： 1锁定与解锁；2清除界面；4检测噪声；7保存 3、操作步骤及注意事项 1、每个GEODE用数传线按规定串联，通过数传盒与笔记本电脑的USB 口连接。 2、每个GEODE接上12V电源。 3、开关接到与笔记本相连的第一个GEODE上。 4、传盒上的开关置于POWER UP处。 5、采集控制程序，并按工作需要设置好各项参数，然后进行正常数据采集工作。 6、出采集控制程序之前，应将数传盒上的开关置于POWER DOWN处。 7、卸下各连接线并清理整齐。 8、注意的是：在正常工作过程中，任何时候移动数传线与GEODE的连接头时，必须退出采集控制程序。另外Y型头上有红色标记的与GEODE的前12道相连接。而且采集控制软件运行的语言环境必须是英语(美国)。
本科生实验报告 实验课程基于 Vista 系统的地震资料处理学院名称地球物理学院 专业名称勘查技术与工程（石油物探）学生姓名 学生学号 指导教师唐湘蓉 实验地点5417 实验成绩 2015年3月- 2015年5月
基于 Vista 系统的地震资料处理 一、实验目的及要求 1）认知熟悉地震资料处理软件系统--vista软件的基本功能，了解其并熟练掌握vista软件运行的基本操作； 2）了解并掌握地震数据处理的基本流程，掌握地震数据处理的流程和基本方法，选择合适的处理参数以提高地震数据处理的精度； 3）对比地震资料处理与解释的理论与实际资料处理的结果，深入理解理论，并在理论指导下提高处理解释的水平、提高资料处理的质量； 4）提高综合分析问题的能力与编写实验报告或生产报告的能力。 二、实验内容 总流程 图1 总流程图 1）加载数据 打开Vista软件后选择加入2D的SEG-Y格式的原始地震数据，本实验
所用数据为给定的SHOT-20。加载后的原始地震数据如图2： 图2 原始地震数据显示 2）道均衡 各个道由于炮检距的不同，导致的反射波的振幅的变化，因为在共反射点叠加中，要求每一个叠加道的振幅都应该相等，每一道对叠加所做的贡献是等价的，无特殊情况，一般就以记录图中间的振幅为基准，使近激发点的地震道振幅减少，增加远离激发点的地震道记录的振幅。道均衡流程模块如图3，道均衡结果如图4： 图3 道均衡流程模块
3）建立观测系统 图5 观测系统显示4）初至拾取 初至拾取结果显示如图6：
图6 初至拾取结果显示 5）初至切除 地震记录上的初至波包括直达波和浅层折射波，它们能量强且有一定延续时间，对紧接而来的浅层反射波有干涉和破坏作用。另外，动校正后会引起波形畸变，浅层尤其厉害。对这些强能量初至波和动校正畸变引起的处理办法是“切除”，即将这些波的采样值全部变为零值（充零）。初至切除流程模块如图7，初至切 除结果如图8： 图7 初至切除流程模块
地震资料解释报告 序言 勘查技术与工程卓越班的实践性很强，加强实践教学可以提高学生的动手能力和处理实际问题、分析解决实际问题的能力、使之能更好的适应毕业后实际工作，是一个非常重要的教学环节，也是进一步提高教学质量的重要途径之一。 我们的地震资料解释实践共分两步完成，第一是在学校手工地震资料构造解释课程设计，第二是在东营对news软件的学习。此次实习是在完成了《地震勘探原理》和《地震资料解释》的基础上完成的实习，通过此次实习的机会我们得以理论联系实际并用实践以检验所学理论，各项安排有条不紊的展开。 在每一步的实习过程中都有老师的带领，手工地震资料构造解释课程设计由杨国权老师负责，news软件的学习由张繁昌老师负责。实习过程中注意理论和实际的结合，在老师的带领及同学的相互帮助下，我们顺利的完成了实践所要求的所有内容。
目录 一、实习目的及意义 (4) 二、实习内容 (4) 三、地震资料构造解释 (5) 四、News学习 (7) 五、结论与建议 (26)
一、实习目的及意义 通过课程的学习，对解释软件系统、数据的地质地球物理解释过程等有基本的认识和掌握，通过实习熟悉了勘探方法的整个工作原理和处理解释流程以及实习报告编写等过程。 了解到了反射波的追踪对比、地震资料的地质解释、构造图的绘制、以及研究成果的提交等过程。培养实际技能及对分析和解决实习问题的能力；掌握仪器的工作原理，并学会操作和使用；掌握各方法的基本数据分析和处理技能。 对本专业所从事工作的性质、手段、方法以及新技术、新方法有一个全面的了解，培养学生的实际操作和计算技能以及综合分析问题的独立工作能力，巩固已学过的专业知识，为下一步进入专业课程和毕业论文阶段以及今后走上本专业的工作岗位打下基础。 二、实习内容 地震自资料的构造解释内容主要有工区的地质情况总结、地震资料解释流程、对地震构造解释的分析、体会和建议等。News 的实习内容主要在理论学习好的基础上，学会利用软件完成地震资料解释的整个过程，并得出理论成果。 三、地震资料的构造解释 构造解释是以水平叠加时间剖面为主要资料，利用由地震资料提供的反射波旅行时间、速度等信息，查明地下地层的构造形态、埋藏深度、解除关系等，通过构造解释成果，即使提供钻井井位。 构造解释的三大环节：
《浅层折射波勘探》实验报告
《浅层折射波勘探》实验成绩评定表班级姓名学号
一、实验名称：浅层折射波勘探 二、实验目的 加深对地震勘探基本概念的理解，巩固已学的理论知识，了解数字地震仪的使用和仪器工作参数的选择；了解地震勘探人工震源激发，检波器的安置条件；地震折射波法野外资料的采集技术及方法，并进行资料的整理与解释；了解地震勘探野外工作施工的过程以及组织管理工作。 三、实验原理 1、折射波法基本原理 以水平界面的两层介质进行简要的说明，假设地下深度为h ，有一个水平的速度分界面R ，上、下两层的速度分别为V 1和V 2，且V 2>V 1。 如图1所示。从激发点O 至地面某一接收点D 的距离为X ，折射波旅行的路程为OK 、KE 、ED 之和，则它的旅行时t 为： 图1 水平两层介质折射波时距曲线 1 21V ED V KE V OK t ++= 式1 为了简便起见，先作如下证明：从O ，D 两点分别作界面R 的垂线，则OA ＝DG ＝h ，再自A 、G 分别作OK ，ED 的垂线，几何上不难证明∠BAK ＝∠EGF ＝i ，因
已知2 1 sin V V i = ，所以： 2 1 V V EG EF AK BK == 式2 即 21V AK V BK = 和 2 1V EG V EF = 式3 上式说明，波以速度V 1旅行BK （或EF ）路程与以速度V 2旅行AK （或EC ）路程所需的时间是相等的。将式3的关系和式1作等效置换，并经变换后可得： 2 121222122cos 2V V V V h V x V i h V x t -+=+= 式4 这就是水平两层介质的折射波时距曲线方程。它表示时距曲线是一条直线，若令x ＝0，则可得时距曲线的截距时间t 0（时距曲线延长与t 轴相交处的时间值） 2 12122102cos 2V V V V h V i h t -== 式5 式5表示出界面深度h 和截距时间t0之间的关系，当已知V 1和V 2时，可以求出界面的深度h 。 2、折射波分层解释的t 0法 折射波t 0解释法是常用的地震折射波解释方法，它是针相遇时距曲线观测系统采集发展起来的解释方法。 t 0法解释的主要原理与方法如下： t 0法又称为t 0差数时距曲线法，是解释折射波相遇时距曲线最常用的方法之一。当折射界面的曲率半径比其埋深大得很多的情况下，t 0法通常能取得很好的效果，且具有简便快速的优点。 如图2所示，设有折射波相遇的时距曲线S 1和S 2，两者的激发点分别是O 1 和O 2，
吉林大学实验教学大纲 教学单位名称：吉林大学地球探测科学与技术学院 课程名称：地震勘探数据处理与解释 课程代码：08262026 课程类别：专业课 课程性质：必修课 学时/学分：32/2（其中实验8学时） 面向专业：勘查技术与工程 一．实验课程的教学任务、要求和教学目的 《地震数据处理与解释》课程是应用地球物理系列课程中的一个重要方向，是地球物理勘探中的重要方法之一，与地震勘探原理一起构成了地震勘探研究方向的一个完整体系。是勘查技术与方法专业中应用地球物理方向本科生的一门重要选修课。 本实验课是与理论课紧密联系在一起的。通过实验课的教学，使学生加深对理论理解和将理论知识应用于实践的能力，熟悉基本的数据处理流程，并进行实际的地震资料处理。本实验课实际上是地震勘探数据处理与解释课程的重要组成部分。 二．学生应掌握的实验技术及基本技能 1、掌握常用地震数据处理系统的基本操作方法 2、了解常用地震记录的数据格式及剖面显示方式； 3、掌握动、静校正及水平叠加处理的方法； 4、掌握地震信号的频谱分析和一维、二维滤波； 5、掌握预测反褶积处理技术； 6、了解速度分析的方法和步骤； 7、了解地震波场偏移处理的目的和方法； 8、掌握合成地震记录的制作和分析方法； 9、掌握波动方程地震记录的正演模拟； 10、能编写简单的地震数据处理程序。 三．实验项目内容、学时分配和每组人数
四．实验教材或指导书或主要参考资料 教材采用《应用地球物理教程—地震勘探》。另外可参考以下文献： 1.《地震资料分析—地震资料处理、反演和解释》，渥.伊尔马滋 2.CWP/SU：Seismic Un*x用户手册 五．考核要求、考核方式及成绩评定标准 实验成绩可通过写实验报告，或总结性考核而定，占学生学期总成绩的20%～30%。 六．制定人、审核人、日期 制定人：王德利 审核人：潘保芝 审核日期：2009年9
地震勘探原理课程设计 报告精选文档 TTMS system office room 【TTMS16H-TTMS2A-TTMS8Q8-
地震勘探原理课程设计报告班级： 学号： 姓名： 指导老师： 日期：
目录 Contents 前言 (1) 一、工区概况 (1) 1、工区位置 (1) 2、勘探概况及石油地质特征 (2) 3、T1层位地震地质层位特征 (2) 4、钻井深度及地震层位的相应关系 (2) 5、地震剖面资料描述 (2) 二、完成工作量 (3) 三、成果（资料）解释 (4) 1、层位标定 (4) 2、地震反射时间剖面对比解释 (4) 3、断层识别解释 (5) 4、上数据 (5) 5、断层平面组合 (5)
6、勾绘T0等值线 (6) 7、空间校正，将等T0图转换为真深度图 (6) 8、解释两张图并作报告 (6) 四、成果分析 (7) 五、体会和建议 (8)
前言 作为一门专业基础课程，地震勘探原理在资源勘查工程专业中有着不可或缺的重要地位。对地震勘探原理较好的掌握将使我们在实际工作中能运用地震勘探方法进行矿产资源勘查,工程地质勘查, 地质灾害调查等方面的工作，为进一步深造及研究工作奠定基础。通过学习地震勘探原理, 初步学会如何运用所学的基础理论知识解决专业中的问题, 提高分析问题, 解决实际问题的能力, 训练逻辑思维能力和科学思维方法, 渗透学科前沿问题,懂得所学的基本理论的意义及价值。地震勘探原理课程设计则是将理论知识运用与实际，通过对地震课设的学习，我们将掌握以下内容： 1、地震剖面的对比解释； 2、绘制等t0构造图，包括断点组合，等值线的勾绘等； 3、绘制真深度构造图的一种方法，即将等t0构造图转换为真深度构造图； 4、地震成果的地质分析； 5、编写解释文字报告。
地震数据整体流程 不同软件的地震数据处理方式不同，但是所有软件的处理流程基本是固定不变的，最多也是在处理过程中处理顺序的不同。整体流程如下： 1 数据输入（又称为数据IO） 数据输入是将野外磁带数据转换成处理系统格式，加载到磁盘上，主要指解编或格式转换。 解编：将多路编排方式记录的数据（时序）变为道序记录方式，并对数据进行增益恢复等处理的过程。如果野外采集数据是道序数据，则只需进行格式转换，即转成处理系统可接受的格式。 注：早期的时序数据格式为记录时先记录第一道第一个采样点、第二道第一个采样点、……、第一道第二个采样点、第二道第二个采样点、……直至结束。现在的道序记录格式为记录时直接记录第一道所有数据、第二道所有数据、……直至结束，只是在每一道数据前加上道头
数据。将时序数据变为道序数据只需要对矩阵进行转置即可。 2 置道头 2.1 观测系统定义 目的为模拟野外，定义一个相对坐标系，将野外的激发点、接收点的实际位置放到这个相对的坐标系中。即将SPS文件转换为GE-Lib文件，包括1）物理点间距2）总共有多少个物理点3）炮点位置4）每炮第一道位置5）排列图形。 2.2 置道头 观测系统定义完成后，处理软件中置道头模块，可以根据定义的观测系统，计算出各个需要的道头字的值并放入地震数据的道头中。当道头置入了内容后，我们任取一道都可以从道头中了解到这一道属于哪一炮、哪一道？CMP号是多少？炮间距是多少？炮点静校正量、检波点静校正量是多少？等等。 后续处理的各个模块都是从道头中获取信息，进行相应的处理，如抽CMP道集，只要将数据道头中CMP号相同的道排在一起就可以了。因此道头如果有错误，后续工作也是错误的。 GOEAST软件有128个道头，1个道头占4个字节，关键的为2（炮号）、4（CMP号）、17（道号）、18（物理点号）、19（线号）、20（炮检距）等。 2.3 观测系统检查 利用置完道头的数据，绘制炮、检波点位置图、线性动校正图。 3 静校正（野外静校正） 静校正为利用测得的表层参数或利用地震数据计算静校正量，对地震道进行时间校正，以消除地形、风化层等表层因素变化时对地震波旅行时的影响。 静校正是实现共中心点叠加的一项最主要的基础工作。直接影响叠加效果，决定叠加剖面的信噪比和垂向分辨率，同时影响叠加速度分析的质量。 静校正方法： 1）高程静校正 2）微测井静校正-利用微测井得到的表层厚度、速度信息，计算静校正量 3）初至折射波法 4）微测井（模型法）低频+初至折射波法高频 4 叠前噪音压制 干扰波严重影响叠加剖面效果。在叠前对各种干扰进行去除，为后续资料处理打好基础。 常见干扰有：面波、折射波、直达波、多次波、50Hz工业电干扰及高能随机干扰等多种情况。不同干扰波有不同特点和产生原因，根据干扰波和一次反射波性质（如频率、相位、视速度等）上的不同，把干扰和有效波分离，从而达到干扰波的去除，提高地震资料叠加效
地震勘探:通过人工方法激发地震波,研究地震波在地层中传播的情况,以查明地下的地质构造,力寻找油气田或其他勘探目的服务的一种物探方法. 水平叠加:将不同接收点收到的来自地下同一反射点的不同激发点的信号,经动校正后叠加起来,这种方法可以提高信噪比,改善地震记录的质量,特别是压制一种规则干扰波效果最好 波形曲线:选定一个时刻t1,我们用纵坐标表示各质点离开平衡位置的距离,就得到一条曲线,这条曲线就叫做波在t1时刻沿x方向的波形曲线. 动校正:在水平界面情况下,从观测到的波的旅行时中减去正常时差Δt1得到x/2处的t0时间,这一过程叫动校正或正常时差校正. 多次覆盖:对被追踪的界面进行多次观测. 剖面闭合:是检查对比质量,连接层位,保证解工作正确进行的有效办法,他包括测线交点闭合,测线网的闭合,时间闭合等. 几何地震学:地震波的运动学是研究地震波,波前的空间位置与传播时间的关系，他与几何光学相似,也是引用波前,射线等几何图形来描述波的运动过程和规律,因此又叫几何地震学. 水平分辨率:指沿水平方向能分辨多大的地质体,其值为根号下0.5λh. 时距曲线:从地震源出发,传播主观测点的时间t与观测中点相对于激发点的距离x之间的关系 剩余时差:把某个波按水平界面一次反射波作动校正后的反射波时间与共中心点处的时间tom之差. 绕射波:地震波在传播过程中,如遇到一些岩性的突变点,这些突变点就会成为新震源,再次发出球面波,想四周传播,这就叫绕射波. 三维地震:就是在一个观测面上进行观测,对所得资料进行三维偏移叠加处理,以获得地下地质体构造在三维空间的特征. 水平切片:就是用一个水平面去切三维数据体得出某一时刻tk各道的信息,更便于了解地下构造形态个查明某些特殊地质现象. 同相轴:一串套合很好的波峰或波谷. 相位:一个完整波形的第i个波峰或波谷. 纵波:传播方向与质点振动方向一致的波. 转换波:当一入射波入射到反射界面时,会产生与其类型相同的反射波或透射波,也会产生类型不同的,与其类型不同的称为转换波. 反射定律:入射波与反射波分居法线两侧,反射角等于入射角,条件为:上下界面波阻抗存在差异,入射波与反射波类型相同. 地震子波：震源产生的信号传播一段时间后，波形趋于稳定，我们称这时的地震波为地震子波。 爆炸时产生的尖脉冲，在爆炸点附近的介质中以冲击波的形式传播，当传播到一的距离后，波形逐渐稳定，我们称这时的地震波为地震子波。 正常时差的定义第一种定义:界面水平情况下,对界面上某点以炮检距x进行观测得到的反射波旅行时同以零炮检距(自激自收)进行观测得到的反射波旅行时之差,这纯粹是因为炮检距不为零引起的时差. 第二种定义:在水平界面情况下,各观测点相对于爆炸点纯粹是由于炮检距不同而引起的反射波旅行时间差. 1.简述地震勘探原理 地震勘探根据岩石的弹性差别进行工作的,波遇到障碍物会发生反射和透射,折射.通过测反射波和透射波的性质,可以确定障碍物的距离.地震勘探是人工激发地震波.通过在地面布置测线,接收反射波,然后进行一些处理,从而来反映地下构造情况,为寻找油气和其他勘探目的的服务,生产工作包括三个环节:1野外数据采集2室内数据处理3地震资料解释,与其他方法
名词解释： 1.褶积模型：地震记录的褶积模型是当今地震勘探中三大环节的主要理论基础之一，其应用十分广泛，主要表现在三大方面：正演、反演和子波处理。层状介质的一次反射波通常用线性褶积模型表示,即：式中:w(t)为系统子波;r(t)为反射系数函数,符号“*”表示褶积运算。 2.分辨率：分辨能力是指区分两个靠近物体的能力。度量分辨能力强弱的两种表示：一是距离表示，分辨的垂向距离或横向范围越小，则分辨能力越强；二是时间表示，在地震时间剖面上，相邻地层时间间隔dt 越小，则分辨能力越强。时间间隔dt 的倒数为分辨率。垂向分辨率是指沿地层垂直方向所能分辨的最薄地层厚度。横向分辨率是指横向上所能分辨的最小地质体宽度。 3.薄层解释原理：Dt 《地震资料数据处理》课程设计 总结报告 专业班级： 姓名： 学号： 设计时间： 指导老师： 2011年5月30日 目录 一、设计内容……………………………………………………………… （1）褶积滤波……………………………………………… （2）快变滤波……………………………………………… （3）褶积滤波与快变滤波的比较………………………… （4）设计高通滤波因子…………………………………… （5）频谱分析……………………………………………… （6）分析补零对振幅谱的影响…………………………… （7）线性褶积与循环褶积………………………………… （8）最小平方反滤波……………………………………… （9）零相位转换…………………………………………… （10）最小相位转换………………………………………… （11）静校正………………………………………………… 二、附录………………………………………………………………………… （1）附录1：相关程序…………………………………… （2）附录2：相关图件…………………………………… 【附录1：有关程序】 1.褶积滤波 CCCCCCCCCCCCCCCCC 褶积滤波CCCCCCCCCCCCCCCCC PROGRAM MAIN DIMENSION X(100),H1(-50:50),H2(-50:50),Y_LOW(200),Y_BAND(200) PARAMETER (PI=3.141592654) CCCCCCCC H1是低通滤波因子，H2为带通滤波因子CCCCCC REAL X,H1,H2,Y_LOW,Y_BAND REAL dt,F,F1,F2 INTEGER I dt=0.002 F=70.0 F1=10.0 F2=80.0 OPEN(1,FILE='INPUT1.DA T',FORM='FORMATTED',STATUS='UNKNOWN') READ(1,*)(X(I),I=1,100) CCCCCCCCCCCCCCCCCC低通滤波器CCCCCCCCCCCCCCCCC DO 10 I=-50,50 IF (I.EQ.0)THEN H1(I)=2*F*PI/PI ELSE H1(I)=SIN(2*PI*F*I*dt)/(PI*I*dt) END IF 10 CONTINUE CCCCCCCCCCCCCCCC输出低通滤波因子CCCCCCCCCCCCCCCC OPEN(2,FILE='H1_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(2,*)(H1(I),I=-50,50) CLOSE(2) CALL CON(X,H1,Y_LOW,100,101,200) CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCC OPEN(3,FILE='Y_LOW.DA T',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(3,*)(Y_LOW(I),I=51,150) CLOSE(3) CCCCCCCCCCCCCCCCCC带通滤波器CCCCCCCCCCCCCCCCCCCC DO 20 I=-50,50 IF(I.EQ.0)THEN H2(I)=140 ELSE H2(I)=SIN(2*PI*F2*I*dt)/(PI*I*dt)-SIN(2*PI*F1*I*dt)/(PI*I*dt) END IF 20 CONTINUE CCCCCCCCCCCCCCC输出带通滤波因子CCCCCCCCCCCCCCCCC OPEN(4,FILE='H2_BAND.DAT',FORM='FORMA TTED',STATUS='UNKNOWN')地震数据处理课程设计(报告)
相关主题
文本预览