地震数据处理课程设计(报告)

《地震数据处理课程设计(报告)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、中国石油大学（华东）《地震资料数据处理方法》课程设计报告《地震资料数据处理》课程设计总结报告专业班级：姓名：学号：设计时间：指导老师：2011年5月30日37中国石油大学（华东）《地震资料数据处理方法》课程设计报告目录一、设计内容………………………………………………………………（1）褶积滤波………………………………………………（2）快变滤波………………………………………………（3）褶积滤波与快变滤波的比较…………………………（4）设计高通滤波因子……………………………………（5）频谱分析……………………

2、…………………………（6）分析补零对振幅谱的影响……………………………（7）线性褶积与循环褶积…………………………………（8）最小平方反滤波………………………………………（9）零相位转换……………………………………………（10）最小相位转换…………………………………………（11）静校正…………………………………………………二、附录…………………………………………………………………………（1）附录1：相关程序……………………………………（2）附录2：相关图件……………………………………37中国石油大学（华

3、东）《地震资料数据处理方法》课程设计报告【附录1：有关程序】1.褶积滤波CCCCCCCCCCCCCCCCC褶积滤波CCCCCCCCCCCCCCCCCPROGRAMMAINDIMENSIONX(100),H1(-50:50),H2(-50:50),Y_LOW(200),Y_BAND(200)PARAMETER(PI=3.141592654)CCCCCCCCH1是低通滤波因子，H2为带通滤波因子CCCCCCREALX,H1,H2,Y_LOW,Y_BANDREALdt,F,F1,F2INTEGERIdt=0.

4、002F=70.0F1=10.0F2=80.0OPEN(1,FILE='INPUT1.DAT',FORM='FORMATTED',STATUS='UNKNOWN')READ(1,*)(X(I),I=1,100)CCCCCCCCCCCCCCCCCC低通滤波器CCCCCCCCCCCCCCCCCDO10I=-50,50IF(I.EQ.0)THENH1(I)=2*F*PI/PIELSEH1(I)=SIN(2*PI*F*I*dt)/(PI*I*dt)ENDIF10CONTINUECCCCCCCCCCCCCCCC输

5、出低通滤波因子CCCCCCCCCCCCCCCCOPEN(2,FILE='H1_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(2,*)(H1(I),I=-50,50)CLOSE(2)CALLCON(X,H1,Y_LOW,100,101,200)CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCCOPEN(3,FILE='Y_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(3,*

6、)(Y_LOW(I),I=51,150)CLOSE(3)CCCCCCCCCCCCCCCCCC带通滤波器CCCCCCCCCCCCCCCCCCCCDO20I=-50,50IF(I.EQ.0)THENH2(I)=140ELSEH2(I)=SIN(2*PI*F2*I*dt)/(PI*I*dt)-SIN(2*PI*F1*I*dt)/(PI*I*dt)ENDIF20CONTINUECCCCCCCCCCCCCCC输出带通滤波因子CCCCCCCCCCCCCCCCCOPEN(4,FILE='H2_BAND.DAT',FO

7、RM='FORMATTED',STATUS='UNKNOWN')37中国石油大学（华东）《地震资料数据处理方法》课程设计报告WRITE(4,*)(H2(I),I=-50,50)CLOSE(4)CALLCON(X,H2,Y_BAND,100,101,200)CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCCCOPEN(5,FILE='Y_BAND.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(5,*)(Y_BAND(I),I=51,

8、150)CLOSE(5)ENDCCCCCCCCCCCCCCCCCCCCC褶积函数CCCCCCCCCCCCCCCCCCCCSUBROUTINECON(A,B,C,I,J,K)DIMENSIONA(I),B(J),C(K)DO1K1=1,K1C(K1)=0.0DO2I1=1,IDO2I2=1,JII=I1+I2-12C(II)=C(II)+A(I1)*B(I2)*0.002RETURNEND1.快变滤波CCCCCCCCCCCCCCC频率滤波CC