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

(完整版)地震数据处理课程设计(报告) (完整版)地震数据处理课程设计(报告) PAGE PAGE64 (完整版)地震数据处理课程设计(报告) PAGE 中国石油大学〔华东〕?地震资料数据处理方法?课程设计报告 ?地震资料数据处理?课程设计 总结报告 专业班级： 姓 名： 学 号： 设计时间： 指导老师： 2021年5月30日 1 中国石油大学〔华东〕?地震资料数据处理方法?课程设计报告 目 录 一、内容????????????????????????〔1〕褶波?????????????????? 〔2〕快波?????????????????? 〔3〕褶波与快波的比?????????? 〔4〕高通波因子?????????????? 〔5〕分析?????????????????? 〔6〕分析零振幅的影响??????????? 〔7〕性褶与循褶????????????? 〔8〕最小平方反波??????????????? 〔9〕零相位????????????????? 〔10〕最小相位???????????????? 〔11〕静校正??????????????????? 二、附????????????????????????????〔1〕附1：相关程序?????????????? 〔2〕附2：相关件?????????????? 2 中国石油大学〔华东〕?地震资料数据处理方法?课程设计报告 【附录1：有关程序】 褶积滤波 褶积滤波CCCCCCCCCCCCCCCCC PROGRAMMAIN DIMENSIONX(100),H1(-50:50),H2(-50:50),Y_LOW(200),Y_BAND(200) PARAMETER(PI=3.141592654) H1是低通滤波因子，H2为带通滤波因子CCCCCCREALX,H1,H2,Y_LOW,Y_BAND REALdt,F,F1,F2 INTEGERIdt=0.002F=70.0F1=10.0F2=80.0 OPEN(1,FILE=INPUT1.DAT,FORM=FORMATTED,STATUS=UNKNOWN)READ(1,*)(X(I),I=1,100) 低通滤波器CCCCCCCCCCCCCCCCC DO10I=-50,50IF(I.EQ.0)THEN H1(I)=2*F*PI/PI ELSE H1(I)=SIN(2*PI*F*I*dt)/(PI*I*dt) ENDIF 10 CONTINUE 输出低通滤波因子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) 输出滤波后的数据CCCCCCCCCCCCCCCCOPEN(3,FILE=Y_LOW.DAT,FORM=FORMATTED,STATUS=UNKNOWN)WRITE(3,*)(Y_LOW(I),I=51,150) CLOSE(3) 带通滤波器CCCCCCCCCCCCCCCCCCCC DO20I=-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) ENDIF 20 CONTINUE 输出带通滤波因子CCCCCCCCCCCCCCCCCOPEN(4,FILE=H2_BAND.DAT,FORM=FORMATTED,STATUS=UNKNOWN) 3 中国石油大学〔华东〕?地震资料数据处理方法?课程设计报告 WRITE(4,*)(H2(I),I=-50,50) CLOSE(4) CALLCON(X,H2,Y_BAND,100,101,200) 输出滤波后的数据CCCCCCCCCCCCCCCCCOPEN(5,FILE=Y_BAND.DAT,FORM=FORMATTED,STATUS=UNKNOWN)WRITE(5,*)(Y_BAND(I),I=51,150) CLOSE(5) END 褶积函数CCCCCCCCCCCCCCCCCCCCSUBROUTINECON(A,B,C,I,J,K) DIMENSIONA(I),B(J),C(K)DO1K1=1,K C(K1)=0.0DO2I1=1,IDO2I2=1,JII=I1+I2-1 C(II)=C(II)+A(I1)*B(I2)*0.002 RETURNEND 快变滤波 频率滤波CCCCCCCCCCCCCCCCCCCC PROGRAMMAIN PARAMETER(PI=3.141592654) REALH_LOW(1:200),H_B