SU(Seismic Unix)常用命令
(2014-08-08 22:53:26)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
$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
10.ximage-变密度显示二进制数据
#ximage n1=750 title=‘name’
11.xwigb-变波形+变面积显示二进制数据(用法与10相同)
#xwigb n1=750 title=‘name’
12.suaddnoise-对地震道加噪声
#suaddnoise 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
18. 画矢量图
supsimage
SU系统的命令可以和Unix系统的命令结合起来,因此可以进行su shell编程,设计特定的处理流程:
#vi plot
#!/bin/bash
#cd
/da
#suwind test_cmp.su key=cdp min=100 max=100
#suximage
#chmod +x plot
#./plot
另外,我们也可以修改或者扩充su的源代码,使其更加符合自己的处理要求。在此不详细介绍。

加载中…