地震資料處理上機(jī)實(shí)習(xí)報(bào)告

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

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

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

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

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

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

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

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

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