场景一:数据拼接
- 从多个二进制炮记录各抽一道拼在一起最后画频谱
suaddhead < file1.dat ns=100 | suwind dt=0.001 key=tracl min=1 max=1 >f1.su suaddhead < file2.dat ns=100 | suwind dt=0.001 key=tracl min=2 max=2 >f2.su suaddhead < file3.dat ns=100 | suwind dt=0.001 key=tracl min=3 max=3 >f3.su cat f1.su f2.su f3.su | sushw key=tracl a=1 >all.su suxwigb < all.su style=normal
场景二:地震记录重采样
问题描述
现有地震记录640道(无道头),每道256个采样点,采样间隔5ms.现在需要重采样成1280个采样点
命令
suaddhead < input.dat ns=256 | sushw key=dt a=5000 >output1.su suresamp <output1.su nt=1280 dt=0.001 >output2.su sustrip < output2.su >output2.dat
场景三:地震记录抽稀
# 从in.su每隔20道抽出一道存到out.su suwind < in.su key=tracl j=20 >out.su
场景四:保存矢量图
img=$1 psimage < ${img}.dat n1=240 \ label1="Depth(km)" label2="Trace" labelfont=Times-Roman labelsize=24 \ f1=0 d1=0.055 f2=0 d2=1 perc=99 \ legend=0 lstyle=vertright lbeg=-0.1 lend=0.1 ldnum=0.05 lgrid=dot \ height=6.0 width=10.0 \ verbose=1 \ lwidth=0.2 lheight=6 \ n1tic=2 f1num=0 d1num=2 n2tic=2 f2num=0 d2num=100 \ axeswidth=2 labelsize=20 \ >${img}.eps
场景五:使用SU做波数域滤波
# 滤出低波数 suaddhead < pp_lap_pow.dat n1=234 n2=663 | sukfilter k=0.,0.3,0.7,1 amps=1,1,0,0 verbose=1 |sushw key=dt a=1000 | sustrip> pp_lap_pow_lk.dat # 滤出中波数 suaddhead < pp_lap_pow.dat n1=234 n2=663 | sukfilter k=0.,0.3,0.7,1 amps=0,1,1,0 verbose=1 |sushw key=dt a=1000 | sustrip> pp_lap_pow_mk.dat # 滤出高波数 suaddhead < pp_lap_pow.dat n1=234 n2=663 | sukfilter k=0.,0.3,0.7,1 amps=0,0,1,1 verbose=1 |sushw key=dt a=1000 | sustrip> pp_lap_pow_hk.dat
场景六:使用SU画波数谱
suaddhead < pp.dat n1=234 n2=663 |sushw key=dt a=1000| suspeck1k2 | suximage cmap=rgb1 legend=1 perc=99