地震處理上機實習報告

1、地震資料處理上機實習姓名：李夢謹班級：061123指導老師：卞愛飛同組成員：鄭思源完成日期：2014.5.14第一部分 SU基礎實驗條件簡介：物探樓 315 實驗室提供安裝英文版 Linux64 位操作系統(tǒng)的微機工作站， 該系統(tǒng)自帶 GNUGCC 編譯套件。 Linux 系統(tǒng)提供 SU（ Seismic Unix） 上機。SU 為 CWP 小組開發(fā)的開源地震信號處理軟件包，具有通用常規(guī)地震資料預處理、速度分析、波場正演、偏移等數(shù)據(jù)處理和顯示功能， 支持 SU、 SEGY 等格式數(shù)據(jù)操作。 SU 提供 C、Fortran 語言編寫的核心算法程序源碼， 具有二次開發(fā)功能。 SU使用方法說明介紹：在

2、我的使用過程中經(jīng)常發(fā)生如果連續(xù)打開文件夾就會出現(xiàn)the filter cant be found.但是后來我一個一個打開之后就解決了這一問題。SU 運行環(huán)境檢查從命令控制臺加載 SU 運行環(huán)境步驟：（ 1） 從 Linux 登陸界面利用實習賬號登陸系統(tǒng)，進入系統(tǒng)后在桌面空白位置單機鼠標右鍵找到Open Terminal選項，打開新的命令控制臺窗口（類似 windows 系統(tǒng)下的 DOS 控制臺）；（ 2） 輸入命令export a=1；（ 3）輸入命令echo $a，檢查控制臺窗口輸出，如果輸出 1，繼續(xù)后續(xù)步驟 4，否則輸入命令bash后進行步驟 4；（ 4） 根據(jù)從系統(tǒng)管理員處獲得 SU

3、軟件安裝目錄情況，在控制臺窗口內(nèi)通過輸入命令echo $CWPROOT，如果輸出結(jié)果與 SU 軟件安裝目錄一致，可以進行后續(xù)實驗，否則輸入source /cwpsu/load_su.sh再進行后續(xù)試驗。命令控制臺下常用 Linux 命令pwd 獲取當前目錄who獲取當前用戶信息ps -a 獲取當前用戶進程信息ls la 獲取當前目錄下子目錄及文件列表chmod u+x filename向當前用戶開放名為 filename 的文件執(zhí)行權(quán)限cp filein fileout 文件復制more file 文件查看gvim filename 利用 gvim 軟件查看文本文檔kedit filename

4、 利用 kedit 軟件查看文本文檔常用 SU 命令（ 1） suhelp說明： 列出軟件包中每個可執(zhí)行程序(即主程序和 shell 文本)的總列表。用法： 直接鍵入 suhelp， 回車直至屏幕輸出結(jié)束（ 2） sudoc說明： 列出 SU 中各條目在線完全文檔用法： sudoc program ，回車直至屏幕輸出結(jié)束 program 為查詢模塊名稱（ 3） sufind說明： 對給定的字符串查找 SU 項目用法： sufind keyword ，回車直至屏幕輸出結(jié)束 Keyword 為查詢關(guān)鍵詞名稱（ 4） suname說明：列出軟件包中完整執(zhí)行程序功能的詳細說明。用法： 直接鍵入 sun

5、ame，回車直至屏幕輸出結(jié)束（ 5） selfdoc說明：列出軟件包具體模塊的使用說明。用法： 直接輸入執(zhí)行模塊名稱， 回車直至屏幕輸出結(jié)束。第二部分 濾波帶通濾波基本原理以及圖示：根據(jù)幅頻特性所表示的通過或阻止信號頻率范圍的不同，濾波器可分為低通濾波器（LPF）、高通濾波器（HPF）、帶通濾波器（BPF）、和帶阻濾波器（BEF）四種。圖分別為圖分別為四種濾波器的實際幅頻特性的示意圖。濾波器是對輸入信號的頻率具有選擇性的一個二端口網(wǎng)絡，它允許某些頻率（通常是某個頻率范圍）的信號通過，而其它頻率的信號幅值均要受到衰減或抑制。這些網(wǎng)絡可以由RLC 元件或RC 元件構(gòu)成的無源濾波器，也可由RC 元件

6、和有源器件構(gòu)成的有源濾波器。Su中帶通濾波模塊使用方法：f 為頻率控制點， amps 為對應頻率控制點的振幅值， &表示后臺運行， 表示方括號內(nèi)的參數(shù)有默認值，可選填。 對于標準的頻率域帶通濾波器， 4 個控制頻點及相應頻點振幅譜即可確定一個帶通濾波器的形態(tài)。實驗數(shù)據(jù)帶通濾波處理及分析：原始數(shù)據(jù)：通過低通濾波器，高頻成分被去除，例如面波等。高通濾波器去除低頻成分。帶通濾波器相當于兩個低通濾波器相減。第三部分 偏移FD,stolt，gazdag，pspi偏移方法的基本原理：疊前偏移: 使CSP道集記錄或COF道集記錄中的反射波歸位, 繞射波收斂疊后偏移: 基于水平疊加剖面，采用爆炸反射面

7、概念實現(xiàn)傾斜反射層歸位和繞射波收斂Stolt偏移法設 為 的二維傅里葉變換，對（1.2.1）式進行上述變換得到： 將（1.2.3）式代入上式有按上行波 取正號并對 微分得：對 做二維傅里葉逆變換得到： 就是要求取的偏移剖面。Gazdag相移法 對標量波動方程（1.2.1）相對x和t做二維傅里葉變換得到： 式中 。求解（1.2.9）式得出F-k域的向下外推公式 偏移成像公式是把上式變換回到空間-時間域，并取t=0時刻的波場值為成像值。即 FD,stolt，gazdag，pspi偏移SU模塊使用方法（ 1）按照 1.3 節(jié)步驟檢查 SU 運行環(huán)境；（ 2）在當前控制臺下運行cd $CWPROOT/

8、demos/filter命令進入 SU 濾波演示目錄;（ 3） 演示數(shù)據(jù)合成。 利用命令gvim README打開演示說明文檔查看演示步驟說明， 輸入命令./MakeDataLittle生成演示數(shù)據(jù)；（ 4） 原始數(shù)據(jù)質(zhì)量監(jiān)控。 輸入命令./XLook， 自動打開四個顯示窗口，利用 gvim 命令查看 XLook 內(nèi)容，分析不同模塊（ xmovie 、 sugain） 輸入顯示參數(shù)對剖面的影響， 分析完后選定顯示窗口按q退出顯示；（ 5）低通濾波。輸入命令./XFilter1 ， 自動打開四個顯示窗口， 利用 gvim 命令查看 XFilter1，分析不同模塊（ sugain、 suspecf

9、x、 sufilter） 功能及不同輸入?yún)?shù)對輸出剖面的影響，分析完后選定顯示窗口按q退出顯示；（ 6）高通濾波。輸入命令./XFilter2 ，自動打開四個顯示窗口，利用 gvim 命令查看 XFilter2，分析 sufilter 模塊濾波輸入?yún)?shù)對輸出剖面的影響，分析完后選定顯示窗口按q退出顯示；（ 7） 自己設計新的帶通/低通/高通濾波器， 并對濾波前后的數(shù)據(jù)加增益顯示，分析濾波效果（提示：改變 sufilter 輸入?yún)?shù) f 和 amps）；（ 8） 直達波切除。 輸入命令./XMute，自動打開四個顯示窗口， 分析直達波切除前以及使用模塊 sumute 及不同切除參數(shù)處理后剖面效果對比；（ 9） 使用不同增益對濾波、 切除直達波后的資料進行顯示。 輸入命令./XGain，自動打開四個顯示窗口，分析直達波切除后使用模塊 sugain 及不同增益參數(shù)處理后剖面效果對比。點脈沖處理分析：傾斜層狀介質(zhì)中偏移結(jié)果分析可以看出在改變了perc后perc越小那么他的偏移效果越不好，出現(xiàn)了很多噪音。第四部分 總結(jié) 通過這次實驗，我們初步掌握了SU的使用方法，了解了地震信號通過濾波器后的結(jié)果，形象的了解采樣頻率對地震信號的影響以及使用不同