加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

SU(Seismic Unix)常用命令

(2014-08-08 22:53:26)
个人的学习使用笔记,从网上收集并应用,同时记下备忘。

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

EBCDIC卡片头和卷头。

1.Segyread-segy格式转换为su格式
segyread tape=vel_z6.25m_x12.5m_lw.segy endian=0|segyclean >vel_z6.25m_x12.5m_lw.su

2.sustrip-去掉道头
sustrip detached/vel.lw.bin
把道头和数据分开

http://einsteinliu1986.blog.163.com/blog/static/103289627200991510112220/

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

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 legend=1 threecolor=1 brgb=1.0,0.0,0.0 grgb=1.0,1.0,1.0 wrgb=0.0,0.0,1.0 label1='Depth(m)' label2='Tr        ace'  *.eps 

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

#vi plot

#!/bin/bash

#cd /data

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

#suximage

#chmod +x plot

#./plot

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

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有