su常用指令介绍.doc

SU开源代码 
一、home/xwd/SU/src/multiponent 关于多波 
/home/xwd/SU/src/su/main/multiponent目录下有关于多波得代码。
 1、 三分量旋转suhrot、c Z-N-E旋转变换到Z-R-T 
2、  线性变换程序 sullt、c R-T通过Li线性变换法变换到S1-S2 
3、 极化滤波 
4、 alford变化程序sualford、c R-T通过alford旋转到S1-S2、 
5、  极化分析 

二、/home/xwd/SU/src/Fortran  Fortran程序 
1  CMP道集NMO速度椭圆方位分析 由Vladimir Grechka所写代码。
 2 3D射线追踪 

三、/home/xwd/SU/src/par 代码挺多。 
1 计算HTI介质准确反射系数 refRealAziHti、c,由Ruger所写。
 2计算VTI介质反射系数 
RefRealVTI、c,由Ruger所写。 
3 二维射线追踪 Ray2D、c 
4各向异性参数与、弹性参数、刚度系数之间得转换 Thom2hti、c, thom2stiff、c 5 其它 
Regrid3、c,三维网格化  不同数据格式之间转换  Recast、c,此外还有很多其它得格式转换数据。 2D、3D最小平方平滑算法  
四、/home/xwd/SU/src/PVM 1、 偏移 2、 贝叶斯反演 3、 静校正 

五、/home/xwd/SU/src/refl 
1、 模拟层状介质中PSV与SH波反射率法 Sureflpsvsh、c  
六、/home/xwd/SU/src/su 
主要包括SU下得一套完整处理流程 1、 amplitude 
振幅处理,包括均衡、去非法值==。 
2、attribute-parameter-estimated 
参数属性计算,包括井资料得各向异性系数、瞬时振幅属性 
3.  correlation-convolution 相关与褶积计算。 
4.  data-conversion 
数据格式转换,均方根速度与层速度之间得计算
 5、 migration_inversion 偏移,速度反演 
5、 其它 
动校正、反褶积、去噪、滤波、道头操作 初至拾取、振幅拾取、静校正、 

七、/home/xwd/SU/src/synthetics-waveforms-testpatterns 波场模拟 

八、/home/xwd/SU/src/transform 
各种变换,包括傅里叶变换、拉冬变换、希尔伯特变化==。  
九、 /home/xwd/SU/src/velocity_analysis 速度分析。
SU学习笔记(一)
[这就是我学习SU时得笔记,其实主要内容就就是SU Manual得翻译。以后,我会陆续把我学习SU时得一些笔记拿出来,给大家分享,希望能对SU得初学者有所帮助。]
第一章 帮助工具
1、 suhelp 显示可执行得程序与Shell脚本。
2、 suname 列出SU中各项命令得名字与简短描述,以及编码得地址。
3、 sudoc 得到编码得DOC列表,列出SU中各条目得在线文档。
4、 sufind 在自述文档中得到信息,使用给定得字符串查找SU命令。
5、 Demo演示程序:
l SU软件包中有一套Shell脚本演示程序,放在目录$CWPROOT/src/demos下。$CWPROOT/src/demos/README文件就是这些演示程序得说明书。
l Making Data 演示程序显示使用susynlv程序制作合成记录炮集与共偏移距道集得基础内容。应当特别注意演示中好得标注风格。
l Filter/Sufilter 使用实际数据处理例子演示说明消除地滚波与初至。
l Deconvolution 使用supef与其它工具简单合成脉冲道集来实例说明去混响与脉冲反褶积处理。演示程序包括使用loops系统检验滤波参数影响得命令。
6、 sukeyword 列出SU道头中得关键字 sukeyword –o
7、 SU基本用法就是通过创建Shell脚本来实现相应得数据处理。su/examples目录下有很多这样得例子。
第二章 核心SU程序
一、从磁带上读写数据
下面得程序对于地球物理应用中特定得数据输入与输出任务就是有用得,对内部SU数据格式也一样。
BHEDTOPAR – 把二进制磁带HEaDer文件转换成PAR文件格式
DT1TOSU – 把Sensors & Software X、dtl GPR得地质雷达数据转成SU格式
SEGDREAD – 读取SEG-D磁带
SEGYCLEAN - zero out unassigned portion of header
SEGYREAD – 读SEG-Y磁带
SEGYHDRS – 为segywrite构造SEG-Y文件得ascii与二进制头文件SEGYWRITE –写SEG-Y磁带
SETBHED - 设置一个SEGY二进制磁带HEaDer文件得道头字
SUADDHEAD – 为裸道加上头文件并设置tracl与ns道头字
SUSTRIP – 从道中去掉SEGY头文件
SUPASTE – 为已存在得数据加上已存在得SEGY头文件
下面得程序可以用于一般得数据输入,输出与数据格式转换,在磁带读写中也可以使用。
A2B – 把ascii float转成二进制
B2A -把二进制float转成ascii
FTNSTRIP – 把Fortran得float转成C格式得float
H2B - convert 8 bit hexidecimal oats to binary
RECAST – 改变数据类型(从一种数据类型转到另一种)
TRANSP – 转置一个n1×n2个元素得矩阵
1. SEGY格式与SU数据格式
SEGY数据格式包括三个部分。第一个部分就是3200字节得EBCDIC卡片,包括40个卡片(等于每行包含80个字符得40行文本),用来磁带。 第二个部分时400个字节得二进制头文件,含有磁带卷内容得信息。SEGY格式得第三个部分就是真正得地震道数据。每道有240个字节得道头文件。接下来, 就是32位得IBM浮点型数据(在IBM Form GA 22-6821中定义)。注意,IBM格式与现代IBM PC上所用得IEEE格式就是不同得。
SU数据格式就是基于SEGY格式得道部分。SEGY道与SU道得主要不同在于SU格式得道数据就是浮点型,就是与您运行SU程序得计算机上得浮点格式一 致得。SU数据只含有SEGY得道部分！SU格式中不保存EBCDIC与二进制卷头,所以无法在任何SU程序中直接使用SEGY文件。
为了把SEGY数据转成SU程序所用得格式,需要使用segyread。
2. SEGYREAD 将SEGY数据读入SU
程序segyread用来把数据从SEGY格式转成SU格式。
当读取SEGY磁带或数据文件时,您需要知道您所使用得机器得byte-order(endian)。SGI,SUN,IBM RS6000使用得就是所谓“big-endian”或high-byte IEEE格式,而基于Intel与Dec芯片得系统用得就是“little-endian”或low-byte。
在big-endian机器上运行segyread得典型方式如下所示:
segyread tape=/dev/rmt0 verbose=1 endian=1 > data、su
更经常使用得就是如下得格式来为big-endian平台读入数据:
segyread tape=/dev/rmt0 verbose=1 endian=1 | segyclean >data、su
在SEGY道头里有可选得道头字(字节181-240)。这些道头字得使用没有标准,所以很多人按自己得需要来填写。SU也不例外。有几个SU图形程序使用得参数存储在这些道头字里。程序segyclean会把可选道头字里容易让SU图形程序产生误解得参数清零。
糟糕得就是有很多号称SEGY格式,但却不符合SEG得标准SEGY格式。最常见得情况就就是为了方便,道部分就是用IEEE格式。这种IEEE格式数据可用下面得命令来读取:
segyread tape=/dev/rmt0 verbose=1 endian=1 conv=0 | segyclean > data、su
这里conv=0就是告诉程序不进行IBM型到float浮点型得转换。
还有DOS SEGY格式,基本与前面相同,除了她得道与头都就是用little-endian格式写得。如果用big-endian机器来读得话用下面得命令:
segyread tape=/dev/rmt0 verbose=1 endian=0 conv=0 | segyclean > data、su
注意:endian=0就是设置交换字节(所有得字节,头与数据都就是交换格式)。在little-endian机器上,程序就是:
segyread tape=/dev/rmt0 verbose=1 endian=1 conv=0 | segyclean > data、su
endian=1会阻止交换字节。
不管哪种情况,如果我们得磁盘文件名为“”,那么应该使用“tape=”参数。
3. SEGYWRITE写SEGY格式得磁带或磁盘文件
与segyread相对应得命令就是segywrite。这个程序可以将SEGY格式按多种不同得方式把数据写到磁带或磁盘文件。该程序可用于把数据写成商业软件可以使用得形式。在学习如何使用segywrite命令之前,有几个需要得准备步骤必须要讨论一下。
4. SEGYHDRS – 为segywrite准备ascii与二进制头文件。
要写一个符合SEG数字磁带标准得SEGY格式,您需要提供ASCII与二进制得卷头文件,在SEGY磁带或文件里会变成EBCDIC与二进制得卷头文件。也就就是segywrite创建文件时需要有header部分与binary部分。
如果您没有binary与header文件,您必须用程序sgyhdrs(创建SEG Y文件)来创建它们。命令:
segyhdrs < data、su
会在当前工作目录下写header与binary文件。举例来说:利用suplane来做一个文件然后对它使用segyhdrs:
suplane > data、su
segyhdrs < data、su
您会瞧到binary与header文件出现在当前工作目录下。
这个命令带有选项,允许您设置binary header fields得值。这些字可以用如下得方式瞧到:
sukeyword jobid
文件header就是一个ASCII文件,可以用正常得文本编辑器来编辑。可以放任意内容,只有格式就是每行80个字符得40行。Segywrite会自动把segyhdrs产生得缺省头文件转成下面得格式:
C This tape was made at the
C
C Center for Wave Phenomena
C Colorado School of Mines
C Golden, CO, 80401
C
、、、
C
C
5. BHEDTOPAR,SETBHED – 编辑二进制头文件
要编辑二进制头文件,首先要转成ASCII格式。程序bhedtopar允许把binary文件写成“parfile”得格式:
bhedtopar < binary outpar=binary、par
可以编辑产生得ASCII码文件“binary、par”进行修改,然后通过setbhed来重新读入:
setbhed b par=binary、par
也可以单独设置头文件字。例如:
setbhed b par=binary、par lino=3
使用了binary、par得内容,但就是lino单独设为3。
最后,可以通过下面得命令来写磁带:
segywrite tape=/dev/rmtx verbose=1 < data、su
注意header文件与binary文件就是在当前目录下得。您也可以使用您自己得文件名。Segywrite得选项bfile=与hfile=用来输入您指定得文件名。
二、数据格式转换
经常需要把来自其它系统或不同格式得输入数据进行转换。在SU里有很多工具与方法处理这些问题。
下面得程序就就是用来处理这种转换得问题:
A2B – 把ascii浮点型转换成二进制
B2A - 把二进制转换成ascii浮点型
FTNSTRIP – 把Fortran floats转成C-style floats
FTNUNSTRIP – 把C-style floats转换成Fortran-style floats
H2B – 把8 bit hexidecimal floats to binary
RECAST – 改变数据类型(从一种数据类型转到另一种)
TRANSP – 转置一个n1×n2得矩阵
1. A2B与B2A - ASCII to Binary, Binary to ASCII
在所有得数据格式中,ASCII就是最常传送得(也就是最耗空间得)。不管您使用什么系统,都可能需要把ASCII转来或转去。而且,因为文本编辑器支持ASCII,因此经常可能要用文本编辑器来做数据输入或编辑。
这种格式大都就是多行格式,用空格或tab来隔开。要转换一个例如5行得数据到二进制,键入:
a2b < data、ascii n1=5 > data、binary
相反得操作就就是:
b2a < data、binary n1=5 > data、ascii
2. FTNSTRIP – 把Fortran数据转换成C语言格式
Fortran在地震数据处理中就是比较流行得语言,因此常常需要处理Fortran创建或处理过得数据。Fortan得二进制数据就是被beginning-of-record与end-of-record分隔符隔开得。而C程序创建得二进制数据没有这些分隔符。要在C程序中使用Fortran数据需要去掉这些Fortran标签,通过:
ftnstrip < fortdata > cdata
3. RECAST – 改变数据类型(从一种数据类型转到另一种)
float 浮点型, double 双精度, int (带符号)整型, char 字符, uchar 无符号字符, short 短整型, long 长整型, ulong 无符号长整型
例如,把整型转换成浮点型:
recast < data、ints in=int out=float > data、floats
4. TRANSP – 转置一个n1×n2个元素得矩阵
5. Farith - 对二进制数据作简单得算术运算
很多时候需要对文件做算术操作,或者在两个二进制数据文件之间。程序farith用来完成许多类似得任务、
Farith对单个文件得操作包括:
Scaling value(缩放数值)
Polarity reversal(极性反转)
Signum function(正负号函数)
absolute value (绝对值)
exponential(取指数)
logarithm(取对数)
square root (开平方根)
square(取平方)
inverse (punctuated), (取道数,带小数点)
inverse of square (punctuated), (取平方倒数,带小数点)
inverse of square root (punctuated) (取平方根倒数,带小数点)
二元操作(对两个文件得操作)包括
addition (加法)
subtraction (减法)
multiplication (乘法)
division (除法)
Cartesian product笛卡尔乘积:从给定得集合X与Y中构成得所有(x,y)元素对集合
地震操作()包括
slowness perturbation
sloth perturbation
使用farith得实例:
farith in=data、binary op=pinv out=data、out、bin
farith in=data1、binary in2=data2、binary op=add > data、out2、bin
三、道头操作
SU数据格式继承了SEGY得道头。如果您得数据不就是SEGY,而就是从其它格式转换得到得,就需要设置一些道头关键字,以使得数据与SU程序协调一致。
这一节涉及得问题(issues)包括:
SUADDHEAD – 在裸数据中写道头,设置道头字
SUSTRIP – 把SU格式数据中得道头切除并存放在文件中,形成裸数据
SUPASTE –把道头文件再粘贴回来,把裸数据添加道头文件生成ＳＵ格式得数据文件
SUKEYWORD – 文件“segy、h”中SU得关键字指南
SURANGE – 获得非零道头输入值得范围(最大值与最小值)
SUSHW – 设置一个或多个道头字,利用道数、取余数(mod)或整除(integer divide)计算道头字得值,或从一个文件中输入道头字得值
SUCHW – 利用一个或两个已存在得道头字来计算新得道头字
SUGETHW – 获取SU数据中道头字得值
SUEDIT – 检查Segy磁盘文件并编辑道头
SUXEDIT – 检查Segy磁盘文件并编辑道头
SUS – 把SU数据字节顺序从big endian 交换为 little endian,或者就是相反。
1. suaddhead – 给二进制数据加SU道头
如果我们得数据由二进制C浮点型文件组成(例如每道1024个采样点),那么下面命令系列将产生SU数据文件“data、su”:
suaddhead < data、bin ns=1024 > data、su
对于其它类型(如整型)使用命令recast:
recast < data、ints in=int out=float | suaddhead ns=1024 > data、su
如果数据首先就是从Fortran转换而来得整型,那么处理流程为:
ftnstrip < data、fortran | recast in=int out=float | suaddhead ns=1024 > data、su
2. sustrip – 将SU道头与SU数据分离
命令suaddhead得逆操作就就是sustrip、下面得命令行将去掉文件SU道头,并把SU道头存到“data、headers”中,形成裸数据data、bin:
sustrip < data、su head=data、headers > data、bin
3. supaste– 将SU道头与SU数据分离
对二进制数据完成处理后,我们可能需要把道头粘贴回去,这时用命令supaste。下面得命令行将把文件“data、headers”中得道头内容粘贴回数据中:
supaste < data、bin head=data、headers > data、su
4. sukeyword – 查瞧SU关键字
sukeyword -o
将显示SU关键字得列表,SU道头中确定有80多个关键字,大多时候只使用其中相对很小得一部分。
5. surange – 获得道头字值得范围
指定数据中道头字值得范围,对于道头字而言就是非常有用得信息。键入:
surange < data、su
将返回所有非零SU道头字值得范围。
请注意,对于损坏得数据来说,很多道头字中可能出现非常奇怪得值,检测这种问题也就是surange命令得主要用法之一。
6. sugethw – 获取SU道头字得值
surange命令可以查瞧整个数据道头字得范围(最大值与最小值)。但就是,我们往往需要按一定得顺序一道一道地查瞧道头字得值。命令sugethw就就是这样得一个工具。例如:
sugethw < data、su key=keyword1,keyword2,、、、 | more
下面给出一个具体得例子:
suplane | sugethw key=tracl,tracr,offset,dt,ns | more
sugethw对于道头字得排列顺序与个数没有要求,但至少要指定一个道头字。
如果需要把道头字输出成二进制得文件,可用下面得命令:
suplane | sugethw key=tracl,tracr,offset,dt,ns output=binary >
对于观测系统,您可以使用下面得命令输出到文件中去:
suplane | sugethw key=tracl,tracr,offset,dt,ns output=geom > hdrfile
7. sushw – 在SU数据中设定道头字得值
命令sushw可根据各种需要来设置道头字得值,该命令可以一次让用户设置一个或多个道头字得值。使用sushw对道头字指定一个固定得值,如我们给数据设定采样间隔:
sushw < data、su key=dt a=2000 > data、out、su
sushw指定得可选参数有:
key= 要修改得关键字
a= 第一道得值
b= 组内增加量
c= 组间增加量
d= 道偏移量
j= 组内元素个数
这些额外得可选参数可用来做更加复杂得操作。这样做就是非常重要得,因为道头字得值与道数据得位置常常有直接得关系。道头字得值用下面得公式计算:
i = itr + d
val(key) = a + b * (i % j) + c * (i / j)
这里itr就是道号(注意:第一道就是itr=0,而非1),%表示取余数,/表示除法。
例如,我们可以设定头五道得道头字sx=6400,第二个5道中sx=6300,依次类推,每5道递减100:
sushw < data、su key=sx a=6400 c=-100 j=5 > data、new、su
另一个例子,我们设置每5道得offset得值为200:200:1000,命令格式为:
sushw < data、su key=offset a=200 b=200 j=5 > data、out、su
我们可以只使用一个sushw命令就可以完成上面3个操作:
sushw < data、su key=dt,sx,offset a=2000,6400,200 b=0,0,200 c=0,-100,0 j=0,5,5 > newdata、su
下面就是一个实际得例子,tracl从1开始,每隔100道增加1;cdp从1开始,每道增加1,一直增加到100,然后重复从1开始;offset同tracl,只就是从0开始;sx同cdp;sy同offset;ns全部设为495;dt全部设为1000:
sushw < key=tracl,cdp,offset,sx,sy,ns,dt a=1,1,0,1,0,495,1000 b=0,1,0,1,0,0,0 c=1,0,1,0,1,0,0 j=100,100,100,100,100,0,0 >
8. suchw – 在SU数据中改变(或计算)道头字得值
有些道头字(如cdp)可以从已有得道头字计算而来,程序suchw就提供了这种功能。
参数有:key1 输出得关键字
key2 输入得关键字
key3 输入得关键字
a 偏移量
b key2关键字得倍数
c key3关键字得倍数
d overall scales
我们可以使用两个道头字得值(key2 与key3),利用下面得方程式计算第三个道头字得值(key1):
val(key1) = (a + b * val(key2) + c * val(key3)) / d
例如:
suchw<indata key1=gx,cdp key2=offset,gx key3=sx,sx b=1,1 c=1,1 d=1,2 >outdata
9. suedit与suxedit – 编辑SU数据中道头字得值
最后,您可能想检查或改变某个道头,suedit与suxedit命令就可以提供这个功能。允许交互浏览与编辑道头字。
例如:
suplane > data、su
suedit data、su
将得到下面得结果:
32 traces in input file
tracl=32 tracr=32 offset=400 ns=64 dt=4000
> <------- prompt for interactive use
suedit 与suxedit交互使用得命令可以通过在提示符后键入问号(?)显示。例如:
32 traces in input file
tracl=32 tracr=32 offset=400 ns=64 dt=4000
>？
n read in trace #n
<CR> step
+ next trace; step -> +1
- prev trace; step -> -1
dN adv N traces; step -> N
% percentiles
r ranks
p [n1 [n2]] tabplot
! key=val modify field
? print this file
q quit
>

该程序让用户将数据采样值按表格形式浏览到数据,或者浏览或改变单个道头字得值。
程序suxedit与suedit类似,但含有X-Windows图形绘图功能用于道数据绘图。
四、编辑SU数据
一旦SU格式得数据读入与道头设置正确后,常常要对数据进行操作与编辑。SU提供得数据处理与编辑得命令有:
SUWIND – 根据关键字开时窗
SUSORT – 基于segy道头关键字进行道分选
SURAMP – 从道起点到道终点进行线性斜坡化(Linearly taper)至零值
SUTAPER – 对一个数据窗进行斜坡化(taper)至零值
SUNULL – 产生空道(都为零值)
SUZERO – 在一时间窗内数据充零
SUKILL – 地震道充零
SUNUTE – 根据关键道头字确定得距曲线得距离,去除用户指定得多边形曲线上边或下边得值
SUVELENGTH – 调整变长度地震道到相同长度地震道
SUVCAT – 将一个数据附加到另一个数据(一道接一道进行)
1、 suwind – 用关键字选定一定时窗地震道
非常常见得操作就就是我们时常需要浏览或处理地震数据中得一小块数据。suwind命令使得我们可以根据大量得参数信息设置窗口来选定感兴趣得数据。
通过道头字来设定数据窗口得大小
suwind最简单得用法就就是,通过用户设定道头关键字得最小与最大值来选定数据:
key=tracl 设置窗口得道头字(see segy、h)
min=LONG_MIN 数据关键道头字得最小值
max=LONG_MAX 数据关键道头字得最大值
例如,用道数作关键字对命令suplane产生得数据进行窗口大小设置:
suplane | suwind key=tracl min=5 max=10 | sugethw key=tracl | more
对于一个大得数据体,应该使用计数(count)参数,而不就是最大值。如果直接设定最大值(max),命令suwind将查询全部数据后,再选择在 最小值与最大值之间得道集,这就是因为程序认为道得标志(trace labeling)可能多次出现。例如,可以比较下面得两个命令得结果:
suplane ntr=100000 | suwind key=tracl min=5 max=10 | sugethw tracl | more
suplane ntr=100000 | suwind key=tracl min=5 count=5 | sugethw tracl | more
更复杂得窗口大小选择如下(例如,使用分数,decimating data):
j=1 每个数据得第j道不选…
s=0 以s为基数(如果((key - s)%j) == 0)
下面得例子用suplane生成得数据,每两道抽取一道(即抽取道2、4、6…):
suplane | suwind key=tracl j=2 | sugethw key=tracl | more
或者每隔2道抽取1道,基数为1(即抽取道1、3、5…):
suplane | suwind key=tracl j=2 s=1 | sugethw key=tracl | more
也可以用suwind接受或拒绝某些地震道:
reject=none 按指定得关键字跳过得地震道
accept=none 用指定得关键字选择地震道
例如下面得命令所示,地震道中得第3、8、9道将不显示:
suplane | suwind key=tracl reject=3,8,9 | sugethw key=tracl | more
参数accept选项有点特殊,它表示接收这些道,即使这些道被拒绝过。例如:
suplane | suwind key=tracl reject=3,8,9 accept=8 | sugethw key=tracl | more
如果您只想接收列表中得道,这时需要设置参数“max=0”,下列所示只显示出第8道:
suplane | suwind key=tracl accept=8 max=0 | sugethw key=tracl | more
计数(count)参数覆盖接收(accept)参数,所以如果您想真正无条件接受道,就不能指定计数(count)参数。
选取时间窗口(Time gating)
有关窗口得第二个问题就就是时间窗得设定。垂直时间窗(time gating)得选项为:
tmin = 0、0 选取最小时间
tmax = (from header) 选取最大时间
itmin = 0 选取最小时间采样
itmax = (from header) 选取最大时间采样
nt = itmax-itmin+1 时间采样点数
设置得结果要么就是选取最小时间采样与选取最大时间采样(itmin与itmax),要么选取最小时间与选取最大时间(tmin与tmax),这样就将产生一个最接近采样得时窗。
2、 susort– 根据segy道头字分选道
susort命令利用Unix系统sort命令,通过道头关键字对地震道进行分选。
例如,根据两个关键字(cdp与offset)分选数据(数值按升序),可用下面得命令:
susort <indata、su >outdata、su cdp offset
如果offset关键字按降序,cdp关键字按升序进行分选,可用下面得命令:
susort <indata、su >outdata、su cdp -offset
3、 suramp与sutaper– 数据斜坡化(Tapering)
很多地震处理算法因为数据得突变边界而显示假象,在数据边界对振幅进行斜坡化(tapering),就就是压制这种假象得最容易得方法。在SU中,我们可以使用sutaper命令对数据边界进行斜坡化。
例如,对每5道数据从头到尾进行斜坡化:
sutaper <diskfile >stdout ntaper=5
suramp可以平滑从起始道与(或)末道。例如,斜坡化从0到tmin=0、05s,然后再向下从tmax=1、15s到末道:
suramp <disk tmax=1、15 >stdout
4、 sukill、suzero、sunull与sumute– 数据充零
对有噪声得道、位于数据边界得道进行充零就是有益得。或者就是产生一些空道,在绘图时最为联系数据集中间得分割标志。
sukill – 道充零
对一批地震道充零键入:
sukill <stdin >stdout min=MIN_TRACE count=COUNT
这里参数count就是要充零得总道数,参数min就是要充零得这些道中得最小道号。
sunull – 产生空道数据
有时需要生成道值为零得地震数据块。下面得命令产生共有NTR道,时间采样数为NT得地震数据:
sunull nt=NT ntr=NTR <stdin >stdout min=MIN_TRACE count=COUNT
suzero – 在一段时间窗内对数据充零
suzero itmin=MIN_TIME_SAMPLE itmax=MAX_TIME_SAMPLE <indata、su > outdata、su
sumute – 数据去除
要进行高精度得压制操作,可用命令sumute对SU数据进行处理。
sumute <indata、su >outdata、su key=KEYWORD xmute=x1,x2,x3,、、、 tmute=t1,t2,t3,、、、
下面用suplane生成数据,然后做去除处理与原始生成数据比较:
suplane | suxwigb &
suplane | sumute key=tracl xmute=1,10,12 tmute=、06,、1,、11 | suxwigb &
去除根据xmute=与tmute=参数确定得多边形曲线内每一个初至。
5、 suvcat与cat – 数据合并
有两种方法可以将一种数据附加到另一个上(合并)。第一种使用Unix命令cat,简单得将第二个文件得数据放到第一个文件中去。
cat data1、su data2、su > data3、su
另外,可能需要对地震道重编号:
cat data1、su data2、su | sushw key=tracl a=1 > data3、su
第二种就是将第二个数据集中得每一道垂直得附加到第一个数据集中每一道得末尾。这就需要使用suvcat命令:
suvcat data1、su data2、su > data3、su
在这个例子中,就不需要修改道头字。
6、 suvlength – 调整可变长度道到相同采样点数
有时数据中每道含有