Demo1
#! /bin/sh # 0. 设置参数，在屏幕显示命令 set -x # 1. 将输入的文件转重命名modeldata # mv $1 modeldata # 1. 加载叠后剖面 # suximage < modeldata legend=1 d2=tr.d2 d1=tr.d1 title="test_su" cmap=rgb1 & # 2. 在加载叠后剖面获取cdp、t时刻 cdp=208 tcdp=800 # 3. 获取动校正NMO后的数据 velpanel=/opt/seismic/data-tm-flt-dec-bal-cdp-nmo.su cdp_file="cmp_${cdp}_t_${tcdp}.su" suwind <$velpanel key=cdp min=$cdp max=$cdp >$cdp_file # 4. 将指定t时刻的数据归零 cdp_file_zero="cmp_${cdp}_t_${tcdp}_zero.su" suzero itmin=$tcdp itmax=$tcdp < $cdp_file > $cdp_file_zero # # 5. 反NMO计算出CMP道集数据 file_contents=$(cat velocity_model.txt) # echo "$file_contents" cdp_file_nnmo="cmp_${cdp}_nnmo.su" sunmo </opt/seismic/$cdp_file_zero invert=1 $file_contents >$cdp_file_nnmo
显示
suwind < outdata.su key=cdp min=231 max=240 | suxwigb