SU(Seismic Unix)常用命令

SU(Seismic Unix)常用命令

SU（Seismic Unix）是由美国科罗拉多矿业学院开发的地震数据处理系统。虽然该系统占用资源很少，但是具有功能齐全、运算效率高、源代码公开等优点，在地球物理界应用广泛。但是也要注意到SU与诸如Promax、Omega、CGG等强大的处理软件相比，功能仍然非常有限，不过我们可以对SU进行修改或者加入自己的源代码，扩充它的功能。该系统没有可视化的软件界面，所有处理都是以命令的形式进行操作，因此需要对Linux系统有一定的了解。

学习SU的一个重要途径是利用强大的帮助工具，帮助命令有很多。 #name-得到该命令的相关用法及参数 e.g. #suximage

#sudoc-得到编码的DOC列表及命令和源程序的位置 e.g. #sudoc suximage #sufind-查找相关命令及用法 e.g. #sufind dmo

另外，SU系统有一套演示程序，放在$CWPROOT/src/demos目录下。进入这一目录，可以发现有很多文件夹，每个文件夹包含了SU各个模块的演示程序，进入这些文件夹，根据README文件的操作提示，一步一步执行里面的shell程序，就可以看到演示结果。

$cd $CWPROOT/src/demos $cd nmo $cat README $./MakeData

SU系统只能对SU格式的数据进行处理，因此对数据处理之前先要将该数据转化为SU数据。SU数据格式只含有SEGY格式的道头和数据记录，不含有EBCDIC卡片头和卷头。

1.Segyread-segy格式转换为su格式

#segyread tape=test.sgy endian=0 | segyclean > test.su su转segy

#segyhdrs < *.su | segywrite tape=*.sgy endian=0 2.Suedit-查看或编辑道头信息 #suedit

3.suaddhead-给二进制数据加道头 #suaddhead stdout ns=750

4.sushw-设置地震道头字

#sushw < stdin key=cdp,offset a=10,0 b=2,0 >stdout 5.sustrip-去掉道头 #sustrip stdout 6.susort-抽道集

#susort stdout cdp offset 7.suwind-截取数据

#suwind stdout key=cdp min=100 max=200 8.suximage-变密度显示su数据

#suximage

9.suxwigb-变波形+变面积显示su数据(用法与8相同) #suxwigb

#ximage n1=750 title=‘name’

11.xwigb-变波形+变面积显示二进制数据(用法与10相同) #xwigb n1=750 title=‘name’ stdout sn=20 13.sugazmig-Gazdag F-K域叠后偏移

#sugazmig stdout vfile=vel.bin dt=0.004 dx=15.0 14.sugethw 提取数据道头信息 #sugethw <*.su key=offset,cdp |more 15.Swapbytes 高低位格式转换 #swapbytes stdout 16.sugethw 获取道头字信息 显示道头字

sugethw

#sugethw < sudata key=key1,key2,... output=geom > hdrfile 17.velconv速度场转换

可以对深度域层速度、时间域层速度、时间域均方根速度之间进行转化 以深度域层速度转时间域均方根速度为例：

#velconv < mar_v.bin > mar_vrmst.bin intype=vintz outtype=vrmst nz=750 dz=4.0 nt=750 dt=0.004

SU系统的命令可以和Unix系统的命令结合起来，因此可以进行su shell编程，设计特定的处理流程：

#vi plot #!/bin/bash #cd /data

#suwind test_cmp.su key=cdp min=100 max=100 #suximage

另外，我们也可以修改或者扩充su的源代码，使其更加符合自己的处理要求。在此不详细介绍。