地震數(shù)據(jù)處理課程設計(報告)

《地震數(shù)據(jù)處理課程設計(報告)》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在教育資源-天天文庫。

1、中國石油大學（華東）《地震資料數(shù)據(jù)處理方法》課程設計報告《地震資料數(shù)據(jù)處理》課程設計總結報告專業(yè)班級：姓名：學號：設計時間：指導老師：2011年5月30日37中國石油大學（華東）《地震資料數(shù)據(jù)處理方法》課程設計報告目錄一、設計內(nèi)容………………………………………………………………（1）褶積濾波………………………………………………（2）快變?yōu)V波………………………………………………（3）褶積濾波與快變?yōu)V波的比較…………………………（4）設計高通濾波因子……………………………………（5）頻譜分析……………………

2、…………………………（6）分析補零對振幅譜的影響……………………………（7）線性褶積與循環(huán)褶積…………………………………（8）最小平方反濾波………………………………………（9）零相位轉換……………………………………………（10）最小相位轉換…………………………………………（11）靜校正…………………………………………………二、附錄…………………………………………………………………………（1）附錄1：相關程序……………………………………（2）附錄2：相關圖件……………………………………37中國石油大學（華

3、東）《地震資料數(shù)據(jù)處理方法》課程設計報告【附錄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輸出濾波后的數(shù)據(jù)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中國石油大學（華東）《地震資料數(shù)據(jù)處理方法》課程設計報告WRITE(4,*)(H2(I),I=-50,50)CLOSE(4)CALLCON(X,H2,Y_BAND,100,101,200)CCCCCCCCCCCCCCCC輸出濾波后的數(shù)據(jù)CCCCCCCCCCCCCCCCCOPEN(5,FILE='Y_BAND.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(5,*)(Y_BAND(I),I=51,

8、150)CLOSE(5)ENDCCCCCCCCCCCCCCCCCCCCC褶積函數(shù)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.快變?yōu)V波CCCCCCCCCCCCCCC頻率濾波CC