1、Seismic Unix基本操作介紹報告人：黃欣芮目錄3幫助工具4核心SU程序5繪圖6常見處理2關于SU7SU與SHELL語言18附錄幫助工具 一、前期準備工作： 1.設置軟件源 剛裝上ubuntu后，首先需要將網絡設置好，并選擇合適的軟件源。對于教育網來說，推薦用中國科大的軟件源。在百度上均可找到軟件源的代碼。 1)首先備份Ubuntu12.04源列表，終端中輸入一下命令： (備份下當前的源列表) 2)修改更新源 (打開Ubuntu12.04源列表文件) 注：軟件源的設置應該根據電腦上所裝ubuntu版本而定。這里以ubuntu12.04版本的安裝為準。 幫助工具 3)將中科大源的代碼粘貼進

2、去(可直接放在文件中的前面)： 4)通知ubuntu啟用新的更新源：終端中鍵入 至此，源的更新工作已經完畢。 注：軟件源的設置應該根據電腦上所裝ubuntu版本而定。這里以ubuntu12.04版本的安裝為準。 幫助工具 2.安裝相應軟件包 在終端中直接用 sudo apt-get install XXX-XXX命令安裝相應的軟件包。我主要安裝了如下的包： build-essential gcc g77 lesstif2 lesstif2-dev libx11-dev libglu1-mesa ligblu1-mesa-dev xlibmesa-glu freeglut3 freeglut3.

3、dev libxmu6 libxmu-dev libxi6 libxi-dev libxt-dev libglut3-dev gfortran 以安裝build-essential為例:在終端中鍵入sudo apt-get install build-essetial命令，若出現如下畫面，表示安裝正常啟動，稍后的工作就是等待安裝結束，進行下一個軟件包的安裝:后面的軟件包的安裝也依舊如此，只需在終端中輸入相應的安裝命令即可(見上面綠體字)。幫助工具 3.設置系統環境 1)創建su目錄： 2)在home目錄下，用vi或gedit打開.bashrc文件(終端下命令格式 gedit .bashrc):

4、 3)在文件最后加上以下三行(這里的yourname=huangxr)： export CWPROOT=/home/yourname/su #設置seismic unix的根目錄。 export PATH=$PATH:/home/yourname/su/bin #設置seismic unix的路徑。 export XAPPLRESDIR=/home/yourname/su/lib/X11/app-defaults #這是README_X中提示要設置的環境變量。 幫助工具 4)保存退出。在終端下運行： source .bashrc 完成環境設置生效。 echo $CWPROOT 檢查CWPROO

5、T設置是否正確。 echo $PATH 檢查CWPROOT設置是否正確。 如果環境變量設置沒問題，那么終端中的界面同下面的一樣： 至此，系統環境設置完畢。 幫助工具 二、安裝Seismic Unix 1.下載Seismic Unix 從/cwpcodes/下載當前的SU版本。我下載的是cwp_su_all_43R3.tgz (進入上面的網站，點擊下面界面紅色方框部分的內容即可下載)。 具體存放在是/home/yourname/su下： 幫助工具 2.解壓文件 在終端中進入/home/yourname/su目錄，然后鍵入tar -zxvf cwp_s

6、u_all_41_tar.Z 解壓文件包： 在/home/yourname/su/下生成/src目錄。SU軟件的相關文件都在其下。安裝說明Installation_Instructions和安裝配置文件Makefile.config，需要認真閱讀。幫助工具 *3.修改makefile.config文件(這是以前的SU版本需要進行的步驟) 進入/home/yourname/su/src 這個目錄，可以看到makefile.config 這個文件，要修改里面的三項,輸入命令：gedit makefile.config , 然后修改所打開的文件： ENDIANFLAG = -DCWP_LITTLE_

7、ENDIAN 句激活（即去掉#） include $(CWPROOT)/src/Rules/gnumake.rules 句激活，并屏蔽掉 include$(CWPROOT)/src/Rules/oldmake.rules 句,修改后 保存。 XDRFLAG = -DSUXDR 句激活（即去掉#） LARGE_FILE_FLAG = -D_FILE_OFFSET_BITS=64 -D_LARGEFILE_SOURCE 句激活(即去掉#，否則seismic unix處理數據有大小限制) OPTC = -O3 -Wall -Wno-long-long -ansi -pedantic -D_POSIX

8、_SOURCE 句激活(即去掉#，這一項要與cc編譯系統相配) CC = cc OPTF = -O -static -march=i686 -fno-automatic -fno-second-underscore 這一項要與fortran編譯系統相配 FC=g77, FC=f77,FC=ifort 屏蔽掉，激活FC=gfortran #現在的新的SU版本中的makefile.config該文件已經修改過了，因此不需要進行*3.修改makefile.config文件這一步。幫助工具 此時在進行正式安裝SU之前，建議先看一下/home/yourname/su/src目錄下的安裝說明Install

9、ation_Instructions文件，確保前面的所有工作都已經正確完成(因為新版本的SU的發布，可能會有一些變化，導致它在安裝上也會出現有一些變化，如果此時完全照著本文的安裝步驟進行可能會出錯)。幫助工具 二.正式安裝 1.編譯安裝 根據Installation_Instructions 提示，在/home/yourname/su/src目錄下分步進行安裝： make install (to install the basic set of codes) make xtinstall (to install the X-toolkit applications) make finstall

10、 (to install the Fortran codes) make mglinstall (to install the Mesa/ Open GL items) 先以make install為例，展示其安裝時的部分界面： 現在終端鍵入上述的命令: 。隨之安裝開始。 若安裝結束后出現下邊的畫面，則表示這一步的安裝成功： 其它幾步的安裝也同上，在此不必一一說明。 幫助工具 2.測試安裝是否成功 為了檢測安裝SU是否成功，在終端中鍵入下面的管程命令： suplane | suxwigb & 如果出現下圖，就表示安裝成功！ 目錄3幫助工具4核心SU程序5繪圖6常見處理27SU與SHELL語言1

11、SU安裝8附錄關于SU Seismic Unix地震數據處理系統（簡稱SU軟件包）由科羅拉多礦業學院(Colorado School of Mines,CSM)波場處理中心(the Center for Wave Phenomena,CWP)開發的處理程序。此軟件是在基于Unix的操作系統上創造一個地震處理環境(用C語言編寫)來擴展Unix操作系統來進行地震處理和研究。而且，該軟件的所有源碼免費提供給想用的人。 SU軟件包不僅限于進行地震處理。用SU軟件可以作與波場相關的一切大套處理，在一定意義上這是一個含有比“Seismic”(地震)一詞更廣義的軟件。SU目的是作為Unix操作系統的一個擴充

12、，因此它兼有Unix系統的很多特性，包括Unix的靈活性和擴展性。 就像Unix系統操作一樣，SU可以堪稱是一種語言(或者元語言) SU是什么關于SU SU軟件包目前沒有使用功能驅動的圖形用戶界面。它是Unix操作系統的一個延伸，一些SU的圖形界面只能是有限地訪問軟件包中的大部分功能。 SU不是商業地震軟件包的替代品，但是SU能夠作為使用的任何一種商業軟件包的重要助手。商業軟件包用于研究工作，SU常用作一個原型軟件包(a prototyping package)。還有如果需要寫進新的編碼，SU能夠作為新的軟件應用的開始的基礎。 SU并不限于地震應用。它可能用于勘探地球物理和通常的信號處理應用中

13、。它也一定能夠用于向學生介紹關于“與波有關”的信號處理(特別是傅立葉變換)教學中。這包括雷達、非地震聲波和圖像資料處理等。 再有，SU既不能說是一個3D軟件包，也不能說是一個2D軟件包，因為其中有很多的濾波和到操作任務與2D和3D中一樣。 SU不是什么關于SU 通過訪問CWP主頁：/cwpcodes/index.html可以下載最新版本的SU軟件包，也可以通過該網站了解更多的CWP/SU的內容。 因為SU包含有很多程序，剛開始接觸沒有那么容易。這里推薦讀者可以查閱書籍Seismic Unix地震數據處理系統-SU3.3處理操作系統指南。這本書好比一

14、本“字典”，可以很好地回答SU中一些不可避免的問題。 參考資料目錄34核心SU程序5繪圖6常見處理2關于SU7SU與SHELL語言1SU安裝8附錄幫助工具 有很多工具按不同的詳細程度對軟件包中的主程序、shell文本以及庫函數提供了內部文檔： SUHELP - 列出CWP/SU程序shell文本。 SUNAME - 自述文檔中各項的名字行以及源碼的位置。 The sefdoc - 是關于大多數可執行主程序和shell文本的一個內部文檔工具。通過在命令行鍵入不帶參數、沒有用管城（|）或Unix重定向符()進行輸入和輸出重定向的程序名和shell文本命令，就可以看到自述文檔內容。 SUDOC -

15、得到編碼的DOC列表。 SUFIND - 自述文檔里得到的信息。 SUKEYWORD - segy.h 文件中SU關鍵字說明。 Demos - SU中包含了一套demos，位于目錄$CWPROOT/src/demos下，主要是shell腳本。 $CWPROOT/src/su/main下，SU程序源碼。 $CWPROOT/src/su/examples下，相關的shell程序。 在此，只介紹一下sukeyword、The selfdoc、Demos。其它幫助工具請見附錄-幫助工具。幫助工具 使用方法：在Linux終端窗口命令行中輸入不帶參數或重定向到文件的程序命令。 用途：簡單說明程序的用途。

16、舉例：sustack # sustack的自述文檔段落就打印（顯示）到屏幕中 The selfdoc對該程序的用途進行了簡短說明 指出如何在命令行鍵入程序命令，關鍵字“stdin”和“stdout”分別表示輸入取自標準輸入和標準輸出。Unix這些術語的意思是用戶能夠從磁盤文件輸入和輸出參數，或者用Unix重定向輸入命令符“”，或管程“|”輸入和輸出參數。幫助工具 說明：查看SU道頭域中的關鍵字。 使用方法： sukeyword name # name - 道頭關鍵字名稱 舉例： 若要查看關鍵字“tracl”，在終端中鍵入: sukeyword tracl # tracl - 側線中道順序號 #

17、當從模擬程序或從野外數據中提取數據時，了解上面這個集合的內容是有好處的。SUKEYWORD幫助工具 說明：SU軟件包中有一套演示程序，它們是shell腳本，放在目錄$CWPROOT/src/demos中。它們可以用來驗證SU程序。 建議：如果是剛學SU軟件包，可以按以下的順序運行演示程序： Making_Data目錄中有用susynlv程序制作合成記錄炮點集和共偏移距剖面的基礎內容演示流程。 Filtering/Sufilter目錄中有一些實際數據處理中消除地滾波和初至演示程序。 Deconvolution目錄中演示了使用supef和其它工具簡單合成脈沖道集，來舉例說明去混響和脈沖反褶積處理。

18、演示程序包括使用loops系統檢驗濾波參數影響的命令。 Demos-演示程序幫助工具 Sorting_Traces指南是一個交互文本，加強了文獻中討論的一些UNIX和SU基本知識。 下一步是激活Selecting_Traces演示程序。再接著是NMO演示程序。除此之外，像SYNLV、Migaration等里面也包含重要的演示程序，也應該值得演示。 當更加熟習SU的邏輯知識后，就能更自如地按源碼運行任何其它演示程序，或者用戶按自己的目的修改shell文本來演示程序！ 目錄3幫助工具45繪圖6常見處理2關于SU7SU與SHELL語言1SU安裝8附錄核心SU程序 Seismic Unix的核心程序可

19、以處理多種多樣的任務，這些任務可以看作是關于數據研究和處理的。這些任務包括： 輸入/輸出。 數據格式轉換。 設置、查看和編輯道頭。 瀏覽查看SU數據。 開時窗顯示、抽道集和編輯數據。 一般操作。 變換和濾波操作。 對SU數據的地震處理操作。 在這一章節中主要介紹前四部分的內容。核心SU程序 輸入/輸出輸入輸入/ /輸出輸出1.SEGY格式和SU數據格式： 在20實際80年代早期，最通用的數據存儲格式是SEG-Y。這是美國勘探地球物理學家協會(Society of Exploration Geophysicists,SEG)的Y格式。這種格式現在仍然廣泛使用。 SEGY數據格式由三部分組成。第一

20、部分是3200字節的EBCDIC卡片頭，包括40個卡片對該磁帶進行描述的文本數據。第二部分是400個字節的二進制（卷）頭，描述了該磁帶卷上內容的相關信息。第三部是實際的地震道數據。每道含有240個字節的道頭，數據記錄在道頭后面。 SU數據格式基于SEGY格式的道部分。SEGY和SU道的主要不同在于SU格式的道數據是浮點型，使用運行SU程序的本機上的二進制浮點格式寫進。SU數據只由SEGY的道數據組成！SU格式中不保存EBCDIC和二進制卷頭。 用segyread程序可以把SEGY數據轉成SU程序能用的格式。核心SU程序 SEGY格式SU格式卷頭說明卷頭說明 (3200)卷頭信卷頭信息息(400

21、)道頭道頭1(240)數據數據1道頭道頭N(240)數據數據N3200個字節ebcid編碼區：40條記錄80字節=3200字節400字節binary二進制數區域卷頭文件240個字節道頭信息；每道數據核心SU程序 2.SEGYREAD - 將SEG-Y數據讀入SU 說明：程序segyread用來把數據從SEGY格式轉換成SU格式。 用法：segyread tape=*segy verbose=1 endian=0 conv=0 | segyclean*.su tape后面可以跟文件路徑，如 tape=/home/zc/*.sgy 參數：endian=0 or 1 #=0 表示“小頭”；=1 表示

22、“大頭” 注意事項： 1）當讀SEGY磁帶或數據文件時，需要使用機器的位序(byter-order)(即頭，endian)。SGI、SUN、IBM和所有Motorola芯片機使用的是“大頭”(big-endian)或高位(high-byte)IEEE格式，而基于Intel和Dec芯片的系統用的是“小頭”(little-endian)或低位(low-byte)； 2）在SEGY道頭里有可選的道頭字(181-240字節)。這些道頭字的使用沒有標準。SU也不例外。有些SU圖形程序使用的參數就存儲在這些道頭字節中。segyclean程序用于把這些可選的道頭字節信息清零，使得SU圖形程序不至于被這些信息

23、混淆。 #更多關于segyread的參數信息請鍵入自述文檔查看。核心SU程序 3.SEGYHDRS - 為segywrite創建SEGY ascii碼和二進制道頭 說明：segyhdrs命令創建“header”(道頭)和“binary”(二進制)文件。 用法：segyhdrs data.su 舉例：利用suplane命令制作一些測試數據，然后運行segyhdrs命令。 注意：“header”（道頭）文件是一個ascii碼文件。 4.SEGYWRITE - 寫SEGY磁帶或磁盤文件 說明：與segyread相聯系的命令segywrite程序，可以將SEGY格式按多種不同的方式把數據寫到磁帶或磁盤

24、文件。segywrite創建文件時需要有header部分和binary部分 。如果你沒有binary和header文件，你必須用上面提到的程序sgyhdrs來創建它們。最后通過segywrite命令來寫磁帶，總的命令是： segyhdrs *.su | segywrite tape=*.segy verbose=1 核心SU程序 輸入/輸出數據格式轉換數據格式轉換常常需要對其它系統的數據進行轉換，或需要輸入不同格式的數據。在SU里有很多工具和方法處理這些問題。 下面列出幾個相關的程序： A2B - 把ascii碼浮點型轉換換成二進制碼。 B2A - 把二進制浮點型轉換換成ascii碼。 FTN

25、STRIP - 把Fortran浮點型轉換成C浮點型。 FTNUNSTRIP- 把C浮點型轉換成Fortran浮點型。 H2B - 把8位16進制浮點型轉換成二進制。 RECAST - 改變數據類型(從一種數據類型轉到另一種)。 SWAPBYTES - 不同的數據類型(高、低位)字節(BYTE)交換。 SUSWAPBYTES - SU數據字節序從“大頭”(big endian)交換為“小頭”(little endian)，或者是反過來。 下面列舉幾個常用的程序。核心SU程序 1.A2B和B2A - ASCII到Binary、Binary到ASCII 說明：在所有的數據格式中，ASCII碼最具移

26、植形(也最耗空間)。任何系統均能輸入或輸出ASCII碼數據，還有文本編輯器支持ASCII碼。而二進制文件卻不能訪問。當需編輯二進制碼文件時， 需要用到A2B程序，將ASCII轉換成Binary。反之，需用到b2a，可查看或編輯ASCII碼文件。 用法：a2b *.bianry 相反的操作是： b2a *.ascii 注意：ASCII碼文件可以用通常的文本編輯器進行編輯。只要滿足總共個40行、每行80個字符的格式要求，可以在文件中放任何內容。因此上面的n1= number 中的number最好別太小(代表列數)。否則超過40行，ASCII碼文件不能正常顯示。 核心SU程序 2.FTNSTRIP

27、- Fortran數據輸入到C 說明：Fortran語言中二進制數據用“記錄頭”(BOR)和“記錄尾”(EOR)分隔符分開。而C程序語言創建的二進制數據沒有這些分隔符。要在C程序中使用Fortran數據需要去掉這些Fortran標識符。 用法：ftnstrip c.data 3.RECAST - 改變二進制數據類型 說明：C支持很多種數據類型，recast程序能對大量的數據類型集進行轉換處理。 舉例：把整型轉換成浮點型 用法：recast data.floats 可以改變的數據類型： float - 浮點型;double - 雙精度;int - (帶符號)整型;short - 短整型; cha

28、r - 字符;uchar - 無符號字符;long - 長整型;ulong - 無符號長整型.核心SU程序 輸入/輸出設置、查看和編輯道頭在SU數據(SEGY格式例子后面)這些參數的值儲存在道頭中。有很多程序可以訪問這些道頭，因此可以根據不同目的設置、瀏覽和修改這些道頭字。 SUADDHEAD - 把道頭加到空道(bare traces)上，設置tracl和ns道頭字。 SUSTRIP - 從道中去掉SEGY道頭。 SUPASTE - 已存在的數據粘貼上已存在的SEGY道頭。 SURANGE - 獲取非零道頭輸入值中最大值和最小值。 SUSHW - 設置一個或多個道頭字，利用道數、取余數(mo

29、d)或整除(integer divide)計算道頭值或從一個文件中輸入道頭字。 SUCHW - 利用一個或兩個道頭關鍵字改變道頭關鍵字。 SUGETHEW - 獲取SU數據中道頭關鍵字。 SUEDIT - 檢查segy磁盤文件并編輯道頭。 SUXEDIT - 檢查segy磁盤文件并編輯道頭核心SU程序 1.SUADDHEAD - 加SU(SEGY風格)道頭 說明：如果數據格式正確，也就是C浮點型排列中快速空間方向為每道采樣點數增加方向，就可以給這些數據加上道頭，以便其它SU程序訪問。 用法：suaddhead *.su 例1：加道頭到含有C二進制浮點型數據文件中,用surange檢查數據。 從

30、上圖可以看到,shot.dat文件已經成功從一般的地震數據格式轉換成SU數據格式的shot.su。相關的SU程序surange、suedit可以訪問該文件. 核心SU程序 例2：對于其它的類型(如整形)使用命令recast recast data.su 例3：如果數據首先是從Fortran轉換而來的整型，那么處理流程為 ftnstrip data.su核心SU程序 2.SUSTRIP - 將SU道頭與SU數據分離 說明：命令sustrip是命令suaddhead的逆操作。 用法：sustrip data.bin 上面的命令系列將去掉文件SU道頭，并把SU道頭存在文件data.headers 3

31、.SUPASTE - 將SU道頭粘貼到二進制數據 說明：對二進制數據完成處理操作后，可能要把道頭粘貼回去，這時用命令supaste。 用法：supaste data.su or suspate data.su heade=headers ns= #注意：在用supaste命令時，千萬不要把“ns= ”參數忘記，否則“粘貼”不成功。 核心SU程序 舉例：用之前由正演模擬轉換來的數據“shot.dat”上聯合使用sustrip和supaste命令。 核心SU程序 4.SURANGE - 獲得道頭字值的范圍 說明：對于道頭字一條有用的信息是看指定數據集中道頭字數值范圍，此時可用surange程序。

32、用法：susrange *.su 參數： traces=201 #共200道 ns=4000 #每道采樣點為4000 #注意：對于破壞的數據來說，大量的道頭中可能顯示出奇怪的數值。檢測這種問題是命令surange的主要用法之一。 核心SU程序 5.SUGETHW - 在SU數據中獲取道頭關鍵字的值 說明：前面的surange命令可以查看整個數據道頭字表示的最大和最小值。但是，我們嘗嘗要一道一道查看道頭字的數值，并按我們選擇的順序。命令sugethw就是這樣的一個工具。 用法：surange data.su key=keyword1,keyword2, | more 舉例： 核心SU程序 對于“

33、觀測系統設置”(geometry setting)，可以按照下面的命令操作： 可以看到，此時屏幕上沒有顯示道頭字，是因為 的緣故。 核心SU程序 6.SUSHW - 在SU數據中設定道頭字 說明：命令sushw可根據各種目的來設置地震道頭值。該命令可讓用戶設置一個或更多的道頭字。它常用來對一個特定道頭字設置一個值。 用法：sushw data.su 參數： key=cdp， 設置的道頭關鍵字 a=0， 第一道頭的值 b=0， 組內遞增量 c=0， 組增量 d=0， 道頭序號增減 j=ULONG_MAX,ULONG_MAX， 組中的元素的數量 核心SU程序 例1：地震數據“shot.su”中沒有

34、表示“dt”(采樣間隔)道頭字。通過sushw命令給“shot.su”加上道頭字“dt”。并鍵入surange檢查。 核心SU程序 對于a,b,c,d,j這些額外的可選參數可用來做更復雜的操作。道頭字的值用以下公式計算： 1） i=itr+d 2）val(key)=a+b*(i%j)+c*(i/j) 這里itr是道號(第一道是itr=0,而不是1)。上面百分號(%)表示取余函數(modulo)，/表示除法。 例2：1)設定道頭字dt=2000(按2ms采樣；2)設定5道中的道頭字sx=6400，第二個5道中sx=6300，后面是每5道遞減(-100；3)設定每5道數據集道頭字offset的值為

35、200、400、6400 步驟1）： 或者 核心SU程序 步驟2）： 或者 核心SU程序 步驟3）： 或者 核心SU程序 當然，可以值引用一次sushw命令就完成所有3中操作： 或者 核心SU程序 在前面的key=dt中，為什么大小為2毫秒的dt其值要用“a=2000”表示，而不是a=2或a=0.002。此時可以鍵入sukeyword查看。 通過上面補充信息，我們可以進一步了解sukeyword的作用。可見，熟習地運用sukeyword程序，對于準確地查找相應關鍵字的信息具有很大的幫助作用！ 核心SU程序 7.SUCHW - 在SU數據中改變(或計算)道頭關鍵字 說明：有些道頭字(如cdp)可

36、以用已有的道頭字計算而來，程序suchw就提供這種功能。 用法：suchw data.su 參數： key1=cdp， 輸出的關鍵字 key2=cdp， 輸入的關鍵字 key3=cdp， 輸入的關鍵字 a=0， 整體偏移一常量 b=1， 對第一個輸入道頭字進行乘除運算 c=0， 對第二個輸入道頭字進行乘除運算 d=1， 整體進行乘除運算 核心SU程序 這個程序使用兩個 道頭字的值(即key2和key3) ，利用下面的方程式來計算第三個道頭字(key1) ： val(key1)=(a+b*val(key2)+c*val(key3)/d 例1）道頭字cdp的值偏移一常量(比如-1) suchw o

37、utdata a=-1 例2） 增加一個常量(比如1000)道頭字(比如說ns) suchw key1=ns key2=ns a=100 outfile 例3） 設置道頭字gx，用sushw命令將偏移距(offset)和sx(shot point)道頭字的值相加，然后對sx和gx取平均來計算cdp道頭字值。核心SU程序 8.SUEDIT和SUXEDIT - 編輯SU數據中的道頭關鍵字 說明：suedit和suxedit可以檢查或改變某個道頭。 用法：suedit data.su（suxedit區別不大） 舉例： 目錄3幫助工具4核心SU程序56常見處理2關于SU7SU與SHELL語言1SU安裝

38、8附錄繪圖 Seismic Unix軟件包有小部分(普通C浮點型格式或SU格式)數據繪圖的圖形工具，可在X-Windows環境下在屏幕顯示或PostScript格式硬拷貝。 在SU中可能的繪圖類型有： contour plots(等值線繪圖) gray or colorscale image plots(灰色或彩色圖像繪圖) wiggle trace plots(地震道波形繪圖) line or symbol graph(側線或符號圖) movies(電影) 3D cube plots(PostScript only)(3D立體圖繪圖)(只限于PostScript格式) 這些程序自述文檔較長，

39、有大量的可選項用來選擇圖形的外觀和標注。然而，取視窗選擇數據的功能，應該用程序subset或suwind進行數據處理后才送到繪圖程序。 繪圖 輸入/輸出X-Windows繪圖程序1.繪制通用浮點型數據 在X-Windows環境下，用于瀏覽通用浮點型數據(沒有SU道頭)的程序： XCONTOUR - 用矢量繪圖程序對函數f(x1,x2)的X等值線繪圖。 XIMAGE - 相同采樣點函數f(x1,x2)的X圖像繪圖。 XWIGB - 用位圖對函數f(x1,x2)的X波形地震道繪圖。 XGRAPH - ni對坐標(x,y)的X GRAPHer圖形。 XMOVIE - 相同采樣點函數f(x1,x2)的

40、一幅(幀)或幾幅(幀)圖像。繪圖 舉例：用suplane生成二進制數據，數據大小為n1=64n2=32。然后從SU數據中分離處道頭。接著用上面的繪圖程序來查看數據(這里只列舉出了xcontour，其它的繪圖程序的寫法均可仿照xcontour，這里就不再逐一列出) 說明：符號&是一個Unix命令，告訴工作的shell在后臺運行該程序。當程序執行完后，可鍵入字幕q退出。繪圖 2.X-Windows中SU數據繪圖 為了SU格式數據繪圖，創建了很多程序。它們是： SUXCONTOUR - 用矢量繪圖程序對地震Unix到文件X等值線繪圖。 SUXIMAGE - SU數據X-Windows圖像繪圖。 SU

41、XWIGB - SU數據X-Windows位圖地震波形繪圖。 SUXGRAPH - SU數據X-Windows圖形繪圖。 SUXMOVIE - SU數據X電影繪圖。 SUXMAX - SU數據每道最大值、最小值和最大絕對值X-Windows圖。 舉例：用suplane命令生成的數據對上述各個程序進行測試 suplane | suxcontou title=contour & suplane | suximage title=image & 其它的SU繪圖程序均可仿照上面的suxcontour和suximage，這里不再逐一列出。 繪圖 輸入/輸出PostScript繪圖程序1.通常浮點型數據P

42、ostScript繪圖 可用與通常浮點型數據(數據沒有SU道頭)PostScript繪圖程序有： PSCONTOUR - 兩個空間函數f(x1,x2)PostScript等值線繪圖。 PSIMAGE - 等間隔函數f(x1,x2)PostScript圖像繪圖。 PSWIGB - 用位圖對函數f(x1,x2)的PostScript地震道繪圖。 PSGRAPH - ni對坐標(x,y)PostScript圖形。 PSMOVIE - 等間隔采樣函數f(x1,x2,x3)PostScript電影繪圖。 PSCUBE - 數據三維數據體PostScript圖像繪圖。 PSWIGP - 用多邊形對函數f(

43、x1,x2)的PostScript地震道繪圖。 舉例：用前面去掉道頭的data.bin(n1=64,n2=32)數據 pscontourdata.eps 在終端中鍵入gs data.eps可查看該圖片。繪圖 2.SU數據PostScript繪圖 這些程序有： SUPSCONTOUR - SU數據PostScript等值線繪圖。 SUPSIMAGE - SU數據PostScript圖像繪圖。 SUPSCUBE - SU數據PostScrip立體繪圖。 SUPSWIGB - SU數據PostScrip位圖地震波形繪圖。 SUPSGRAPH - SU數據PostScrip圖形繪圖。 SUPSWIGP

44、 - SU數據PostScrip多邊形填充波形繪圖。 SUPSMAX - SU數據每道最大值、最小值和最大絕對值PostScrip繪圖。 舉例：類同于前面的通常浮點型數據PostScript繪圖，這里不再敘述。 繪圖 3.另外的PostScript支持 在SU中還有幾種工具可用于PostScript操作，它們是： PSBBOX - 改變存在的PostScript文件的四周的框架(BoudingBOX)。 PSMERGE - 合并PostScript文件。 MERGE2 - 把多張PostScrip圖合并到一頁上。 MERGE4 - 把多張圖合并到一頁上MERGE4 figures onto o

45、ne page。 這些操作在這里就不詳細闡述。 可以鍵入下面的命令獲得更多的信息： more $CWPROOT/bin/merge2 或者下面的命令： more $CWPROOT/bin/merge4 另外，在終端中鍵入gs *.eps,可查看所有的“*.eps”類型的圖片。 繪圖 最后展示2張由有限差分正演模擬得來的地震記錄的結果。左邊一張是一般地震數據（無道頭）格式的shot.dat文件，由ximage得到的結果。右邊一張是轉換得來的加上道頭的SU數據格式的shot.su文件，由suximage的得到結果。 目錄3幫助工具4核心SU程序5繪圖62關于SU7SU與SHELL語言1SU安裝8附

46、錄常見處理 在第三章的“核心程序”中，主要介紹了一些道頭操作的命令。這一章中，將對一些常用的SU數據編輯操作、通用操作、地震模擬工具、地震處理工具做簡單介紹。 這些部分幾乎都有CWP自帶的DEMOS,參見$CWPROOT/src/demos目錄。對于一般情況，DEMOS的文件執行具有一定的順序性。常見處理 輸入/輸出編輯SU數據一旦數據讀入和道頭設置正確后，常常要對數據進行操作和編輯。這些操作任務常常不同，可以進行下面數據處理： Windowed(取時窗) Sorted(道分選) Truncated(道截取) Tapered(斜坡化處理) Zeroed(道充零) 作等間隔采樣 Concaten

47、ated(數據聯合) 在這里，只簡單介紹取視窗的操作suwind。其它處理的相關程序見附錄 編輯SU數據。常見處理 SUWIND - 用關鍵字選定一定時窗地震道 說明：很常見的操作是瀏覽和處理地震數據中的一小塊。suwind命令可以用于選定地震數據窗口大小。 參數： key=tracl # 設置窗口的關鍵道頭字 min=LONG_MIN # 數據關鍵道頭字最小值 max=LONG_MAX # 數據關鍵道頭字最大值 count=number # 計數參數 j=1 # 每個數據的第j道不選 s=0 # 以s為基數(如果(key-s)%j )= 0) reject=none # 按指定的關鍵字跳過的

48、地震道 accept=none # 用指定的關鍵字選擇地震道 常見處理 例1：用道數關鍵字對命令suplane產生的數據進行窗口大小的設置 常見處理 例2：每個2道選取數據，基數為1 常見處理 例3：用reject拒絕第3，8，9道。用accept接受第8道。 常見處理 輸入/輸出SU數據通用操作SU數據編輯以外的操作任務是進行下面的處理： gaining(增益) resampling(重采樣) unary operaions 一元操作(對單個文件的算術操作)。 binary operations 二元操作(對兩個文件的算術操作)。 進行這些處理的程序有： SUADDNOISE - 對地震道加

49、噪聲 SUGAIN - 按不同的增益方式顯示地震道 SUOP - 在segy數據做一元算術操作 SUOP2 - 對兩個數據做二元操作 在這里，不做詳細介紹。具體可參考其自述文檔“the selfdoc”。常見處理 輸入/輸出變換和濾波操作地震研究和地震處理的很多方面與數學變換方法有關。特別地，很多地震處理離不開數字Fourier變換。濾波是一個與此相關的主題，因為大多數濾波在頻率域進行，或者至少是可以用數學上進行頻率域操作來表示。1.1D Fourier變換 SUFFT - 從實數時間道變換到復數頻率道(正變換)。 SUIFFT - 從復數頻率道變換到實數時間道(逆變換)。 SUAMP - 從

50、(頻率，x)數據輸出振幅(amp)、相位(amp)、道實部和虛部。 SUFFT - 從實數時間道變換到復數頻率道(正變換)。 SUSPECFX - 地震數據(時間T到頻率F)Fourier頻譜分析。常見處理 例1：suplane | suxigb title=Original Data & #左圖 suplane | sufft | suifft | sushw key=d1,dt a=0,4000 | suxwigb & #右圖 #結果與輸入完全一樣，除了地震道結果中有更多的采樣點，因為變換是要補零(zero-padding) 常見處理 例2：瀏覽sufit程序輸出的振幅譜和相位譜、及實部和

51、虛部 suplane|sufft|suamp mode=amp|suxwigb title=amplitude& suplane|sufft|suamp mode=phase|suxwigb title=phases& suplane|sufft|suamp mode=real|suxwigb title=real& suplane|sufft|suamp mode=imag|suxwigb title=imaginary& 例3：SU數據格式可以存儲復數數據的實部和虛部。鍵入下面的命令查看： 例4：用suspecfx命令快速查看一個地震道或一塊地震數據的振幅譜 suplane|suspecf

52、x|suximage title=F-X Amplitude Spectrum& 常見處理 2.2D Fourier變換 SUSPECFK - 數據F-K Fourier頻譜分析。 SUSPECK1K2 - (x1,x2)數據2D(K1,K2)Fourier頻譜分析。 舉例： suplane|suspecfk|suximage title=F-K Amplitude Spectrum& suplane|suspeck1k2|suximage title=K1-K1 Amplitude Spectrum& 當然，還有很多其它變換方法，如： Hilbert變換 Gabor變換 Randon變換(T

53、an P濾波) 它們都包含在SU軟件包中。具體可以查閱目錄： $CWPROOT/src/demos/Time_Freq_Analysis $CWPROOT/src/demos/Filtering/Sugabor $CWPROOT/src/demos/Tau_P 常見處理 3.濾波操作 地震處理中很大一部分工作可以稱作是濾波。在SU軟件包中濾波應用程序可以進行簡單的濾波操作、和更復雜的反褶積和子波整型處理。 SUFILTER - 應用零相位正弦平方斜坡濾。 SUBFILT - Butterworth帶通濾波。 SUACOR - 自相關。 SUCONV,SUXCOR - 根據用戶提供的濾波器進行褶

54、積、相關r。 SUPEF - Wiener預測誤差濾波。 SUSHAPE - Wiener整型濾波。 SURESAMP - 時間域重采樣。 SUTVBAND - 時變帶通濾波(正弦平方濾波)。 這里不再仔細敘述，具體可以參考各程序的自述文檔The selfdoc和查閱目錄： $CWPROOT/src/demos/Deconvolution $CWPROOT/src/demos/Filtering 常見處理 輸入/輸出地震模擬工具地震勘探和研究的重要方面就是制作合成記錄。地震模擬包括兩個方面：第一部分是構造背景波動速度剖面(background wavespeed profile)，由均勻采樣系

55、列浮點數組成；第二部分是構造在波速剖面中傳播的合成地震信息。 1.背景速度剖面 SU中有幾個程可用于生成背景波速剖面，常常這些數據需要平滑： UNISAM - 由一對(x,y)左邊確定的函數y(x)均勻重采樣。 UNISAM2 - 2D函數f(x1,x2)均勻重采樣。 MAKEVEL - 制造一個速度函數v(x,y,z)。 UNIF2 - 產生層狀模型2D均勻采樣速度剖面。每層中速度位置是線性函數 SMOOTHINT2 - 非均勻采樣界面平滑，阻尼最小平方技術。 SMOOTH2 - 在用戶制定窗口內對均勻采樣2D數組數據平滑。 SMOOTH3D - 阻尼最小平方3D網格速度平滑。 這里不再詳細

56、列出所有相關的程序，僅以unif2為例。請查閱這些程序的自述文檔得到更多的信息。也可查看：$CWPROOT/src/demos/Velocity_Profiles目錄下的演示程序。 常見處理 舉例：UNIF2的第一種表達方法 常見處理 文件model里面的內容如下： 表示橫向坐標，單位m表示縱向坐標，單位m表示一個分界面常見處理 舉例：UNIF2的第二種表達方法 因為經unif2得到的vfile文件不是SU數據格式的，因此需經過suaddhead一處理，轉換為SU數據格式，才能使用suximage繪圖。常見處理 對于 ，這個文件里的內容如下： 參數： ninf=3 #表示分界面的數量 nx、n

57、z= #橫向和縱向的采樣數 dx、dz= #橫向和縱向的采樣間隔 npmax= #界面上的最大點數 v00=2000,3500,5000,2000 #每個界面的速度(包括表層) 可見，parfile.par文件里存放了背景波動速度剖面需要的參數，它可以精簡編輯界面。 常見處理 2.合成記錄生成器 在SU軟件包中有很多程序可用于生成地震記錄和與地震類似的記錄： SUPLANE - 產生可到3層(3 planes)共偏移距數據文件。 SUSPIKE - 制作小的脈沖數據。 SUIMP2D - 使用Born積分方程在3D空間直線上生成共炮點記錄。 SUIMP3D - 使用Born積分方程在3D空間點

58、上生成平面共炮點記錄(inplane shot records)。 SUFDMOD2 - 聲波方程2階有限差分模擬。 SUSYNCZ - 分段常速度V(Z)函數2.5D真振幅模擬(只限縱波)合成地震記錄。 SUSYNLV - 制作線性速度函數合成地震記錄。 SUSYNVXZ - V(X,Z)介質Kirchhoff模擬共偏移距合成地震記錄。 SUSYNLVCW - 制作線性速度函數模式轉換波合成記錄。 SUSYNVXZCS - V(X,Z)介質Kirchhoff模擬共炮點合成地震記錄。 只有程序sufdmod2、susynvxz和susnvxzcs要求輸入波速文件。而其它程序通過命令行鍵入模型參

59、數。請查閱目錄$CWPROOT/src/demos/Synthetic中演示程序獲取更多信息。這些程序的其它演示程序用于生成合成記錄數據。常見處理 輸入/輸出地震處理工具有一些操作是專門針對地震性質設計的，用于地震數據處理和將地震數據轉換成地層圖像。 stacking data 數據疊加。 picking data 數據拾取。 velocity analysis 速度分析。 normal moveout correction 正常時差校正。 dip moveout correction 傾角時差校正。 seicmic migration and related operation 地震偏移及相

60、關操作。 各種地震處理工具對應的程序請見附錄地震處理工具，這里不再詳細列出。 簡約的定義3幫助工具4核心SU程序5繪圖6常見處理2關于SU71SU安裝8附錄SU與SHELL語言 Shell程序能夠用來極大地擴充SU程序，而不必寫C語言程序編碼。例如，可查閱在$CWPROOT/src/su/examples目錄下的CvStack、FilterTest、FirstBreak和Velan程序。現以CvStack為例加以說明。這個shell程序的核心和關于速度和共深度點道(cdps)的兩個循環，產生速度數據集(velocity panel)，任何單個SU程序都沒有這個概念。SU與SHELL語言 She