su常用指令介紹

1、 SU開源代碼一、home/xwd/SU/src/multicomponent關(guān)于多波/home/xwd/SU/src/su/main/multicomponent目錄下有關(guān)于多波的代碼。1、三分量旋轉(zhuǎn)suhrot.cZ-N-E旋轉(zhuǎn)變換到Z-R-T2、 線性變換程序sullt.cR-T通過Li線性變換法變換到S1-S23、 極化濾波4、 alford變化程序sualford.cR-T通過alford旋轉(zhuǎn)到S1-S2.5、 極化分析二、/home/xwd/SU/src/FortranFortran程序1 CMP道集NMO速度橢圓方位分析由VladimirGrechka所寫代碼。23D射線追蹤三、

2、/home/xwd/SU/src/par代碼挺多。1計(jì)算HTI介質(zhì)準(zhǔn)確反射系數(shù)refRealAziHti.c，由Ruger所寫。2計(jì)算VTI介質(zhì)反射系數(shù)RefRealVTI.c，由Ruger所寫。3二維射線追蹤Ray2D.c4各向異性參數(shù)與、彈性參數(shù)、剛度系數(shù)之間的轉(zhuǎn)換Thom2hti.c,thom2stiff.c5其它Regrid3.c，三維網(wǎng)格化不同數(shù)據(jù)格式之間轉(zhuǎn)換 Recast.c，此外還有很多其它的格式轉(zhuǎn)換數(shù)據(jù)。2D、3D最小平方平滑算法四、/home/xwd/SU/src/PVM1.偏移2.貝葉斯反演3.靜校正五、/home/xwd/SU/src/refl1.模擬層狀介質(zhì)中PSV和S

3、H波反射率法Sureflpsvsh.c六、/home/xwd/SU/src/su主要包括SU下的一套完整處理流程1.amplitude振幅處理，包括均衡、去非法值=。2.attribute-parameter-estimated參數(shù)屬性計(jì)算，包括井資料的各向異性系數(shù)、瞬時(shí)振幅屬性3. correlation-convolution相關(guān)和褶積計(jì)算。4. data-conversion數(shù)據(jù)格式轉(zhuǎn)換，均方根速度與層速度之間的計(jì)算5.migration_inversion偏移，速度反演5.其它 動(dòng)校正、反褶積、去噪、濾波、道頭操作初至拾取、振幅拾取、靜校正、七、/home/xwd/SU/src/syn

4、thetics-waveforms-testpatterns波場模擬八、/home/xwd/SU/src/transform各種變換，包括傅里葉變換、拉冬變換、希爾伯特變化=。9、 /home/xwd/SU/src/velocity_analysis速度分析。SU學(xué)習(xí)筆記（一） 這是我學(xué)習(xí)SU時(shí)的筆記，其實(shí)主要內(nèi)容就是SU Manual的翻譯。以后，我會(huì)陸續(xù)把我學(xué)習(xí)SU時(shí)的一些筆記拿出來，給大家分享，希望能對SU的初學(xué)者有所幫助。 第一章 幫助工具1. suhelp 顯示可執(zhí)行的程序和Shell腳本。2. suname 列出SU中各項(xiàng)命令的名字和簡短描述，以及編碼的地址。3. sudoc 得到

5、編碼的DOC列表，列出SU中各條目的在線文檔。4. sufind 在自述文檔中得到信息，使用給定的字符串查找SU命令。5. Demo演示程序：l SU軟件包中有一套Shell腳本演示程序，放在目錄$CWPROOT/src/demos下。$CWPROOT/src/demos/README文件是這些演示程序的說明書。l Making Data 演示程序顯示使用susynlv程序制作合成記錄炮集和共偏移距道集的基礎(chǔ)內(nèi)容。應(yīng)當(dāng)特別注意演示中好的標(biāo)注風(fēng)格。l Filter/Sufilter 使用實(shí)際數(shù)據(jù)處理例子演示說明消除地滾波和初至。l Deconvolution 使用supef和其它工具簡單合成脈沖

6、道集來實(shí)例說明去混響和脈沖反褶積處理。演示程序包括使用loops系統(tǒng)檢驗(yàn)濾波參數(shù)影響的命令。6. sukeyword 列出SU道頭中的關(guān)鍵字 sukeyword o7. SU基本用法是通過創(chuàng)建Shell腳本來實(shí)現(xiàn)相應(yīng)的數(shù)據(jù)處理。su/examples目錄下有很多這樣的例子。 第二章 核心SU程序一、從磁帶上讀寫數(shù)據(jù)下面的程序?qū)τ诘厍蛭锢響?yīng)用中特定的數(shù)據(jù)輸入和輸出任務(wù)是有用的，對內(nèi)部SU數(shù)據(jù)格式也一樣。BHEDTOPAR 把二進(jìn)制磁帶HEaDer文件轉(zhuǎn)換成PAR文件格式DT1TOSU 把Sensors & Software X.dtl GPR的地質(zhì)雷達(dá)數(shù)據(jù)轉(zhuǎn)成SU格式SEGDREAD 讀取SEG

7、-D磁帶SEGYCLEAN - zero out unassigned portion of headerSEGYREAD 讀SEG-Y磁帶SEGYHDRS 為segywrite構(gòu)造SEG-Y文件的ascii和二進(jìn)制頭文件SEGYWRITE 寫SEG-Y磁帶SETBHED - 設(shè)置一個(gè)SEGY二進(jìn)制磁帶HEaDer文件的道頭字SUADDHEAD 為裸道加上頭文件并設(shè)置tracl和ns道頭字SUSTRIP 從道中去掉SEGY頭文件SUPASTE 為已存在的數(shù)據(jù)加上已存在的SEGY頭文件 下面的程序可以用于一般的數(shù)據(jù)輸入，輸出和數(shù)據(jù)格式轉(zhuǎn)換，在磁帶讀寫中也可以使用。A2B 把a(bǔ)scii float

8、轉(zhuǎn)成二進(jìn)制B2A -把二進(jìn)制float轉(zhuǎn)成asciiFTNSTRIP 把Fortran的float轉(zhuǎn)成C格式的floatH2B - convert 8 bit hexidecimal oats to binaryRECAST 改變數(shù)據(jù)類型(從一種數(shù)據(jù)類型轉(zhuǎn)到另一種)TRANSP 轉(zhuǎn)置一個(gè)n1n2個(gè)元素的矩陣 1 SEGY格式和SU數(shù)據(jù)格式SEGY數(shù)據(jù)格式包括三個(gè)部分。第一個(gè)部分是3200字節(jié)的EBCDIC卡片，包括40個(gè)卡片（等于每行包含80個(gè)字符的40行文本），用來磁帶。 第二個(gè)部分時(shí)400個(gè)字節(jié)的二進(jìn)制頭文件，含有磁帶卷內(nèi)容的信息。SEGY格式的第三個(gè)部分是真正的地震道數(shù)據(jù)。每道有240個(gè)

9、字節(jié)的道頭文件。接下來， 是32位的IBM浮點(diǎn)型數(shù)據(jù)（在IBM Form GA 22-6821中定義）。注意，IBM格式和現(xiàn)代IBM PC上所用的IEEE格式是不同的。SU數(shù)據(jù)格式是基于SEGY格式的道部分。SEGY道和SU道的主要不同在于SU格式的道數(shù)據(jù)是浮點(diǎn)型，是和你運(yùn)行SU程序的計(jì)算機(jī)上的浮點(diǎn)格式一 致的。SU數(shù)據(jù)只含有SEGY的道部分！SU格式中不保存EBCDIC和二進(jìn)制卷頭，所以無法在任何SU程序中直接使用SEGY文件。為了把SEGY數(shù)據(jù)轉(zhuǎn)成SU程序所用的格式，需要使用segyread。2 SEGYREAD 將SEGY數(shù)據(jù)讀入SU程序segyread用來把數(shù)據(jù)從SEGY格式轉(zhuǎn)成SU格

10、式。當(dāng)讀取SEGY磁帶或數(shù)據(jù)文件時(shí)，你需要知道你所使用的機(jī)器的byte-order（endian）。SGI，SUN，IBM RS6000使用的是所謂“big-endian”或high-byte IEEE格式，而基于Intel和Dec芯片的系統(tǒng)用的是“l(fā)ittle-endian”或low-byte。在big-endian機(jī)器上運(yùn)行segyread的典型方式如下所示：segyread tape=/dev/rmt0 verbose=1 endian=1 data.su更經(jīng)常使用的是如下的格式來為big-endian平臺讀入數(shù)據(jù)：segyread tape=/dev/rmt0 verbose=1 en

11、dian=1 | segyclean data.su在SEGY道頭里有可選的道頭字（字節(jié)181-240）。這些道頭字的使用沒有標(biāo)準(zhǔn)，所以很多人按自己的需要來填寫。SU也不例外。有幾個(gè)SU圖形程序使用的參數(shù)存儲(chǔ)在這些道頭字里。程序segyclean會(huì)把可選道頭字里容易讓SU圖形程序產(chǎn)生誤解的參數(shù)清零。糟糕的是有很多號稱SEGY格式，但卻不符合SEG的標(biāo)準(zhǔn)SEGY格式。最常見的情況就是為了方便，道部分是用IEEE格式。這種IEEE格式數(shù)據(jù)可用下面的命令來讀?。簊egyread tape=/dev/rmt0 verbose=1 endian=1 conv=0 | segyclean data.su這

12、里conv=0是告訴程序不進(jìn)行IBM型到float浮點(diǎn)型的轉(zhuǎn)換。還有DOS SEGY格式，基本與前面相同，除了他的道和頭都是用little-endian格式寫的。如果用big-endian機(jī)器來讀的話用下面的命令：segyread tape=/dev/rmt0 verbose=1 endian=0 conv=0 | segyclean data.su注意：endian=0是設(shè)置交換字節(jié)（所有的字節(jié)，頭和數(shù)據(jù)都是交換格式）。在little-endian機(jī)器上，程序是：segyread tape=/dev/rmt0 verbose=1 endian=1 conv=0 | segyclean dat

13、a.suendian=1會(huì)阻止交換字節(jié)。不管哪種情況，如果我們的磁盤文件名為“filename”，那么應(yīng)該使用“tape=filename”參數(shù)。3 SEGYWRITE寫SEGY格式的磁帶或磁盤文件與segyread相對應(yīng)的命令是segywrite。這個(gè)程序可以將SEGY格式按多種不同的方式把數(shù)據(jù)寫到磁帶或磁盤文件。該程序可用于把數(shù)據(jù)寫成商業(yè)軟件可以使用的形式。在學(xué)習(xí)如何使用segywrite命令之前，有幾個(gè)需要的準(zhǔn)備步驟必須要討論一下。4 SEGYHDRS 為segywrite準(zhǔn)備ascii和二進(jìn)制頭文件。要寫一個(gè)符合SEG數(shù)字磁帶標(biāo)準(zhǔn)的SEGY格式，你需要提供ASCII和二進(jìn)制的卷頭文件，

14、在SEGY磁帶或文件里會(huì)變成EBCDIC和二進(jìn)制的卷頭文件。也就是segywrite創(chuàng)建文件時(shí)需要有header部分和binary部分。如果你沒有binary和header文件，你必須用程序sgyhdrs（創(chuàng)建SEG Y文件）來創(chuàng)建它們。命令：segyhdrs data.susegyhdrs data.su你會(huì)看到binary和header文件出現(xiàn)在當(dāng)前工作目錄下。這個(gè)命令帶有選項(xiàng)，允許你設(shè)置binary header fields的值。這些字可以用如下的方式看到：sukeyword jobid文件header是一個(gè)ASCII文件，可以用正常的文本編輯器來編輯?？梢苑湃我鈨?nèi)容，只有格式是每行8

15、0個(gè)字符的40行。Segywrite會(huì)自動(dòng)把segyhdrs產(chǎn)生的缺省頭文件轉(zhuǎn)成下面的格式：C This tape was made at theCC Center for Wave PhenomenaC Colorado School of MinesC Golden, CO, 80401C.CC5 BHEDTOPAR，SETBHED 編輯二進(jìn)制頭文件要編輯二進(jìn)制頭文件，首先要轉(zhuǎn)成ASCII格式。程序bhedtopar允許把binary文件寫成“parfile”的格式：bhedtopar binary outpar=binary.par可以編輯產(chǎn)生的ASCII碼文件“binary.par”進(jìn)

16、行修改，然后通過setbhed來重新讀入：setbhed bfile=binay par=binary.par也可以單獨(dú)設(shè)置頭文件字。例如：setbhed bfile=binary par=binary.par lino=3使用了binary.par的內(nèi)容，但是lino單獨(dú)設(shè)為3。 最后，可以通過下面的命令來寫磁帶：segywrite tape=/dev/rmtx verbose=1 data.su注意header文件和binary文件是在當(dāng)前目錄下的。你也可以使用你自己的文件名。Segywrite的選項(xiàng)bfile=和hfile=用來輸入你指定的文件名。 二、數(shù)據(jù)格式轉(zhuǎn)換經(jīng)常需要把來自其它系統(tǒng)

17、或不同格式的輸入數(shù)據(jù)進(jìn)行轉(zhuǎn)換。在SU里有很多工具和方法處理這些問題。下面的程序就是用來處理這種轉(zhuǎn)換的問題:A2B 把a(bǔ)scii浮點(diǎn)型轉(zhuǎn)換成二進(jìn)制B2A - 把二進(jìn)制轉(zhuǎn)換成ascii浮點(diǎn)型FTNSTRIP 把Fortran floats轉(zhuǎn)成C-style floatsFTNUNSTRIP 把C-style floats轉(zhuǎn)換成Fortran-style floatsH2B 把8 bit hexidecimal floats to binaryRECAST 改變數(shù)據(jù)類型(從一種數(shù)據(jù)類型轉(zhuǎn)到另一種)TRANSP 轉(zhuǎn)置一個(gè)n1n2的矩陣1 A2B和B2A - ASCII to Binary, Binar

18、y to ASCII在所有的數(shù)據(jù)格式中，ASCII是最常傳送的（也是最耗空間的）。不管你使用什么系統(tǒng)，都可能需要把ASCII轉(zhuǎn)來或轉(zhuǎn)去。而且，因?yàn)槲谋揪庉嬈髦С諥SCII，因此經(jīng)?？赡芤梦谋揪庉嬈鱽碜鰯?shù)據(jù)輸入或編輯。這種格式大都是多行格式，用空格或tab來隔開。要轉(zhuǎn)換一個(gè)例如5行的數(shù)據(jù)到二進(jìn)制，鍵入：a2b data.binary相反的操作就是：b2a data.ascii 2 FTNSTRIP 把Fortran數(shù)據(jù)轉(zhuǎn)換成C語言格式Fortran在地震數(shù)據(jù)處理中是比較流行的語言，因此常常需要處理Fortran創(chuàng)建或處理過的數(shù)據(jù)。Fortan的二進(jìn)制數(shù)據(jù)是被beginning-of-recor

19、d和end-of-record分隔符隔開的。而C程序創(chuàng)建的二進(jìn)制數(shù)據(jù)沒有這些分隔符。要在C程序中使用Fortran數(shù)據(jù)需要去掉這些Fortran標(biāo)簽，通過：ftnstrip cdata 3 RECAST 改變數(shù)據(jù)類型(從一種數(shù)據(jù)類型轉(zhuǎn)到另一種)float 浮點(diǎn)型, double 雙精度, int （帶符號）整型, char 字符, uchar 無符號字符, short 短整型, long 長整型, ulong 無符號長整型 例如，把整型轉(zhuǎn)換成浮點(diǎn)型：recast data.floats4 TRANSP 轉(zhuǎn)置一個(gè)n1n2個(gè)元素的矩陣5 Farith - 對二進(jìn)制數(shù)據(jù)作簡單的算術(shù)運(yùn)算很多時(shí)候需要對

20、文件做算術(shù)操作，或者在兩個(gè)二進(jìn)制數(shù)據(jù)文件之間。程序farith用來完成許多類似的任務(wù).Farith對單個(gè)文件的操作包括：Scaling value（縮放數(shù)值）Polarity reversal（極性反轉(zhuǎn)）Signum function（正負(fù)號函數(shù)）absolute value （絕對值）exponential（取指數(shù)）logarithm（取對數(shù)）square root （開平方根）square（取平方）inverse (punctuated), （取道數(shù)，帶小數(shù)點(diǎn)）inverse of square (punctuated), （取平方倒數(shù)，帶小數(shù)點(diǎn)）inverse of square roo

21、t (punctuated) （取平方根倒數(shù)，帶小數(shù)點(diǎn)） 二元操作(對兩個(gè)文件的操作)包括addition （加法）subtraction （減法）multiplication （乘法）division （除法）Cartesian product笛卡爾乘積：從給定的集合X和Y中構(gòu)成的所有（x，y）元素對集合 地震操作（）包括slowness perturbationsloth perturbation 使用farith的實(shí)例：farith in=data.binary op=pinv out=data.out.binfarith in=data1.binary in2=data2.binary

22、 op=add data.out2.bin 三、道頭操作SU數(shù)據(jù)格式繼承了SEGY的道頭。如果你的數(shù)據(jù)不是SEGY，而是從其它格式轉(zhuǎn)換得到的，就需要設(shè)置一些道頭關(guān)鍵字，以使得數(shù)據(jù)與SU程序協(xié)調(diào)一致。這一節(jié)涉及的問題（issues）包括:SUADDHEAD 在裸數(shù)據(jù)中寫道頭，設(shè)置道頭字SUSTRIP 把SU格式數(shù)據(jù)中的道頭切除并存放在文件中，形成裸數(shù)據(jù)SUPASTE 把道頭文件再粘貼回來，把裸數(shù)據(jù)添加道頭文件生成格式的數(shù)據(jù)文件SUKEYWORD 文件“segy.h”中SU的關(guān)鍵字指南SURANGE 獲得非零道頭輸入值的范圍（最大值和最小值）SUSHW 設(shè)置一個(gè)或多個(gè)道頭字，利用道數(shù)、取余數(shù)(mo

23、d)或整除(integer divide)計(jì)算道頭字的值，或從一個(gè)文件中輸入道頭字的值SUCHW 利用一個(gè)或兩個(gè)已存在的道頭字來計(jì)算新的道頭字SUGETHW 獲取SU數(shù)據(jù)中道頭字的值SUEDIT 檢查Segy磁盤文件并編輯道頭SUXEDIT 檢查Segy磁盤文件并編輯道頭SUSWAPBYTES 把SU數(shù)據(jù)字節(jié)順序從big endian 交換為 little endian,或者是相反。1 suaddhead 給二進(jìn)制數(shù)據(jù)加SU道頭如果我們的數(shù)據(jù)由二進(jìn)制C浮點(diǎn)型文件組成（例如每道1024個(gè)采樣點(diǎn)），那么下面命令系列將產(chǎn)生SU數(shù)據(jù)文件“data.su”:suaddhead data.su對于其它類型

24、（如整型）使用命令recast:recast data.su如果數(shù)據(jù)首先是從Fortran轉(zhuǎn)換而來的整型，那么處理流程為：ftnstrip data.su 2 sustrip 將SU道頭與SU數(shù)據(jù)分離命令suaddhead的逆操作就是sustrip.下面的命令行將去掉文件SU道頭，并把SU道頭存到“data.headers”中，形成裸數(shù)據(jù)data.bin：sustrip data.bin3 supaste 將SU道頭與SU數(shù)據(jù)分離對二進(jìn)制數(shù)據(jù)完成處理后，我們可能需要把道頭粘貼回去，這時(shí)用命令supaste。下面的命令行將把文件“data.headers”中的道頭內(nèi)容粘貼回?cái)?shù)據(jù)中：supaste

25、 data.su4 sukeyword 查看SU關(guān)鍵字sukeyword -o將顯示SU關(guān)鍵字的列表，SU道頭中確定有80多個(gè)關(guān)鍵字，大多時(shí)候只使用其中相對很小的一部分。5 surange 獲得道頭字值的范圍指定數(shù)據(jù)中道頭字值得范圍，對于道頭字而言是非常有用的信息。鍵入：surange data.su將返回所有非零SU道頭字值得范圍。請注意，對于損壞的數(shù)據(jù)來說，很多道頭字中可能出現(xiàn)非常奇怪的值，檢測這種問題也是surange命令的主要用法之一。6 sugethw 獲取SU道頭字的值surange命令可以查看整個(gè)數(shù)據(jù)道頭字的范圍（最大值和最小值）。但是，我們往往需要按一定的順序一道一道地查看道頭

26、字的值。命令sugethw就是這樣的一個(gè)工具。例如：sugethw file.bin對于觀測系統(tǒng)，你可以使用下面的命令輸出到文件中去：suplane | sugethw key=tracl,tracr,offset,dt,ns output=geom hdrfile 7 sushw 在SU數(shù)據(jù)中設(shè)定道頭字的值命令sushw可根據(jù)各種需要來設(shè)置道頭字的值，該命令可以一次讓用戶設(shè)置一個(gè)或多個(gè)道頭字的值。使用sushw對道頭字指定一個(gè)固定的值，如我們給數(shù)據(jù)設(shè)定采樣間隔：sushw data.out.susushw指定的可選參數(shù)有：key= 要修改的關(guān)鍵字a= 第一道的值b= 組內(nèi)增加量c= 組間增加

27、量d= 道偏移量j= 組內(nèi)元素個(gè)數(shù)這些額外的可選參數(shù)可用來做更加復(fù)雜的操作。這樣做是非常重要的，因?yàn)榈李^字的值與道數(shù)據(jù)的位置常常有直接的關(guān)系。道頭字的值用下面的公式計(jì)算：i = itr + dval(key) = a + b * (i % j) + c * (i / j)這里itr是道號（注意：第一道是itr=0,而非1），%表示取余數(shù)，/表示除法。例如，我們可以設(shè)定頭五道的道頭字sx=6400,第二個(gè)5道中sx=6300,依次類推，每5道遞減100：sushw data.new.su另一個(gè)例子，我們設(shè)置每5道的offset的值為200：200：1000，命令格式為：sushw data.ou

28、t.su我們可以只使用一個(gè)sushw命令就可以完成上面3個(gè)操作：sushw newdata.su下面是一個(gè)實(shí)際的例子，tracl從1開始，每隔100道增加1；cdp從1開始，每道增加1，一直增加到100，然后重復(fù)從1開始；offset同tracl,只是從0開始；sx同cdp；sy同offset；ns全部設(shè)為495；dt全部設(shè)為1000：sushw filename_new.su 8 suchw 在SU數(shù)據(jù)中改變（或計(jì)算）道頭字的值有些道頭字（如cdp）可以從已有的道頭字計(jì)算而來，程序suchw就提供了這種功能。參數(shù)有：key1 輸出的關(guān)鍵字 key2 輸入的關(guān)鍵字 key3 輸入的關(guān)鍵字 a

29、偏移量 b key2關(guān)鍵字的倍數(shù)c key3關(guān)鍵字的倍數(shù)d overall scales我們可以使用兩個(gè)道頭字的值（key2 和key3），利用下面的方程式計(jì)算第三個(gè)道頭字的值(key1):val(key1) = (a + b * val(key2) + c * val(key3) / d例如：suchwoutdata9 suedit和suxedit 編輯SU數(shù)據(jù)中道頭字的值最后，你可能想檢查或改變某個(gè)道頭，suedit和suxedit命令就可以提供這個(gè)功能。允許交互瀏覽和編輯道頭字。例如：suplane data.susuedit data.su將得到下面的結(jié)果：32 traces in i

30、nput filetracl=32 tracr=32 offset=400 ns=64 dt=4000 ？ n read in trace #n step+ next trace; step - +1- prev trace; step - -1dN adv N traces; step - N% percentilesr ranksp n1 n2 tabplot! key=val modify field? print this fileq quit 該程序讓用戶將數(shù)據(jù)采樣值按表格形式瀏覽到數(shù)據(jù)，或者瀏覽或改變單個(gè)道頭字的值。程序suxedit和suedit類似，但含有X-Windows圖形

31、繪圖功能用于道數(shù)據(jù)繪圖。四、編輯SU數(shù)據(jù)一旦SU格式的數(shù)據(jù)讀入和道頭設(shè)置正確后，常常要對數(shù)據(jù)進(jìn)行操作和編輯。SU提供的數(shù)據(jù)處理和編輯的命令有：SUWIND 根據(jù)關(guān)鍵字開時(shí)窗SUSORT 基于segy道頭關(guān)鍵字進(jìn)行道分選SURAMP 從道起點(diǎn)到道終點(diǎn)進(jìn)行線性斜坡化(Linearly taper)至零值SUTAPER 對一個(gè)數(shù)據(jù)窗進(jìn)行斜坡化(taper)至零值SUNULL 產(chǎn)生空道（都為零值）SUZERO 在一時(shí)間窗內(nèi)數(shù)據(jù)充零SUKILL 地震道充零SUNUTE 根據(jù)關(guān)鍵道頭字確定的距曲線的距離，去除用戶指定的多邊形曲線上邊或下邊的值SUVELENGTH 調(diào)整變長度地震道到相同長度地震道SUVCA

32、T 將一個(gè)數(shù)據(jù)附加到另一個(gè)數(shù)據(jù)（一道接一道進(jìn)行）1. suwind 用關(guān)鍵字選定一定時(shí)窗地震道非常常見的操作就是我們時(shí)常需要瀏覽或處理地震數(shù)據(jù)中的一小塊數(shù)據(jù)。suwind命令使得我們可以根據(jù)大量的參數(shù)信息設(shè)置窗口來選定感興趣的數(shù)據(jù)。通過道頭字來設(shè)定數(shù)據(jù)窗口的大小suwind最簡單的用法就是，通過用戶設(shè)定道頭關(guān)鍵字的最小和最大值來選定數(shù)據(jù)：key=tracl 設(shè)置窗口的道頭字(see segy.h)min=LONG_MIN 數(shù)據(jù)關(guān)鍵道頭字的最小值max=LONG_MAX 數(shù)據(jù)關(guān)鍵道頭字的最大值例如，用道數(shù)作關(guān)鍵字對命令suplane產(chǎn)生的數(shù)據(jù)進(jìn)行窗口大小設(shè)置：suplane | suwind k

33、ey=tracl min=5 max=10 | sugethw key=tracl | more對于一個(gè)大的數(shù)據(jù)體，應(yīng)該使用計(jì)數(shù)（count）參數(shù)，而不是最大值。如果直接設(shè)定最大值（max），命令suwind將查詢?nèi)繑?shù)據(jù)后，再選擇在 最小值和最大值之間的道集，這是因?yàn)槌绦蛘J(rèn)為道的標(biāo)志（trace labeling）可能多次出現(xiàn)。例如，可以比較下面的兩個(gè)命令的結(jié)果：suplane ntr=100000 | suwind key=tracl min=5 max=10 | sugethw tracl | moresuplane ntr=100000 | suwind key=tracl min=5

34、 count=5 | sugethw tracl | more更復(fù)雜的窗口大小選擇如下(例如，使用分?jǐn)?shù)，decimating data)：j=1 每個(gè)數(shù)據(jù)的第j道不選s=0 以s為基數(shù)(如果(key - s)%j) = 0)下面的例子用suplane生成的數(shù)據(jù)，每兩道抽取一道（即抽取道2、4、6）：suplane | suwind key=tracl j=2 | sugethw key=tracl | more或者每隔2道抽取1道，基數(shù)為1（即抽取道1、3、5）：suplane | suwind key=tracl j=2 s=1 | sugethw key=tracl | more也可以用s

35、uwind接受或拒絕某些地震道：reject=none 按指定的關(guān)鍵字跳過的地震道accept=none 用指定的關(guān)鍵字選擇地震道例如下面的命令所示，地震道中的第3、8、9道將不顯示：suplane | suwind key=tracl reject=3,8,9 | sugethw key=tracl | more參數(shù)accept選項(xiàng)有點(diǎn)特殊，它表示接收這些道，即使這些道被拒絕過。例如：suplane | suwind key=tracl reject=3,8,9 accept=8 | sugethw key=tracl | more如果你只想接收列表中的道，這時(shí)需要設(shè)置參數(shù)“max=0”,下

36、列所示只顯示出第8道:suplane | suwind key=tracl accept=8 max=0 | sugethw key=tracl | more計(jì)數(shù)（count）參數(shù)覆蓋接收（accept）參數(shù)，所以如果你想真正無條件接受道，就不能指定計(jì)數(shù)(count)參數(shù)。選取時(shí)間窗口（Time gating）有關(guān)窗口的第二個(gè)問題就是時(shí)間窗的設(shè)定。垂直時(shí)間窗（time gating）的選項(xiàng)為：tmin = 0.0 選取最小時(shí)間tmax = (from header) 選取最大時(shí)間itmin = 0 選取最小時(shí)間采樣itmax = (from header) 選取最大時(shí)間采樣nt = itmax

37、-itmin+1 時(shí)間采樣點(diǎn)數(shù)設(shè)置的結(jié)果要么是選取最小時(shí)間采樣和選取最大時(shí)間采樣（itmin和itmax），要么選取最小時(shí)間和選取最大時(shí)間（tmin和tmax），這樣就將產(chǎn)生一個(gè)最接近采樣的時(shí)窗。2. susort 根據(jù)segy道頭字分選道susort命令利用Unix系統(tǒng)sort命令，通過道頭關(guān)鍵字對地震道進(jìn)行分選。例如，根據(jù)兩個(gè)關(guān)鍵字(cdp和offset)分選數(shù)據(jù)（數(shù)值按升序），可用下面的命令：susort outdata.su cdp offset如果offset關(guān)鍵字按降序，cdp關(guān)鍵字按升序進(jìn)行分選，可用下面的命令：susort outdata.su cdp -offset3. su

38、ramp和sutaper 數(shù)據(jù)斜坡化（Tapering）很多地震處理算法因?yàn)閿?shù)據(jù)的突變邊界而顯示假象，在數(shù)據(jù)邊界對振幅進(jìn)行斜坡化(tapering)，就是壓制這種假象的最容易的方法。在SU中，我們可以使用sutaper命令對數(shù)據(jù)邊界進(jìn)行斜坡化。例如，對每5道數(shù)據(jù)從頭到尾進(jìn)行斜坡化：sutaper stdout ntaper=5suramp可以平滑從起始道和(或)末道。例如，斜坡化從0到tmin=0.05s,然后再向下從tmax=1.15s到末道：suramp stdout 4. sukill、suzero、sunull和sumute 數(shù)據(jù)充零對有噪聲的道、位于數(shù)據(jù)邊界的道進(jìn)行充零是有益的?；蛘?/p>

39、是產(chǎn)生一些空道，在繪圖時(shí)最為聯(lián)系數(shù)據(jù)集中間的分割標(biāo)志。sukill 道充零對一批地震道充零鍵入：sukill stdout min=MIN_TRACE count=COUNT這里參數(shù)count是要充零的總道數(shù)，參數(shù)min是要充零的這些道中的最小道號。sunull 產(chǎn)生空道數(shù)據(jù)有時(shí)需要生成道值為零的地震數(shù)據(jù)塊。下面的命令產(chǎn)生共有NTR道，時(shí)間采樣數(shù)為NT的地震數(shù)據(jù)：sunull nt=NT ntr=NTR stdout min=MIN_TRACE count=COUNTsuzero 在一段時(shí)間窗內(nèi)對數(shù)據(jù)充零suzero itmin=MIN_TIME_SAMPLE itmax=MAX_TIME_S

40、AMPLE outdata.susumute 數(shù)據(jù)去除要進(jìn)行高精度的壓制操作，可用命令sumute對SU數(shù)據(jù)進(jìn)行處理。sumute outdata.su key=KEYWORD xmute=x1,x2,x3,. tmute=t1,t2,t3,.下面用suplane生成數(shù)據(jù)，然后做去除處理與原始生成數(shù)據(jù)比較：suplane | suxwigb &suplane | sumute key=tracl xmute=1,10,12 tmute=.06,.1,.11 | suxwigb &去除根據(jù)xmute=和tmute=參數(shù)確定的多邊形曲線內(nèi)每一個(gè)初至。5. suvcat和cat 數(shù)據(jù)合并有兩種方法可

41、以將一種數(shù)據(jù)附加到另一個(gè)上（合并）。第一種使用Unix命令cat，簡單的將第二個(gè)文件的數(shù)據(jù)放到第一個(gè)文件中去。cat data1.su data2.su data3.su另外，可能需要對地震道重編號：cat data1.su data2.su | sushw key=tracl a=1 data3.su第二種是將第二個(gè)數(shù)據(jù)集中的每一道垂直的附加到第一個(gè)數(shù)據(jù)集中每一道的末尾。這就需要使用suvcat命令：suvcat data1.su data2.su data3.su在這個(gè)例子中，就不需要修改道頭字。6. suvlength 調(diào)整可變長度道到相同采樣點(diǎn)數(shù)有時(shí)數(shù)據(jù)中每道含有不同的采樣點(diǎn)數(shù)。下面我

42、們對用suplane命令生成的數(shù)據(jù)，用命令suvlength處理的道相同的采樣點(diǎn)數(shù)：suplane nt=64 data1.susuplane nt=32 data2.sucat data1.su data2.su data3.su對于上面合并生成的數(shù)據(jù)data3.su，如果想用SU程序處理該數(shù)據(jù)文件將會(huì)失敗，因?yàn)榇蠖嗾fSU程序要求數(shù)據(jù)塊含有相同的采樣點(diǎn)數(shù)。使用命令suvlength可以解決該問題：suvlength ns=64 data4.susuxwigb data.susuxwigb data.su title=Ungained Data &sugain data.su scale=5.0 | suxwigb title=Scaled data &sugain data.su agc=1 wagc=.01 | suxwigb title=AGC=1 WAGC=.01 sec &sugain data.su agc=1 wagc=.2 | suxwigb title=AGC=1 WAGC=.1 sec &sugain data.su pbal=1 | suxwigb title=traces balanced by rms &sugain data.su qbal=1 | suxwigb title=t