Res SU (15)道数据计算1

1. 显示第二道 2. 打印第二道数据 2.1数据存储到txt 3. 读取具体的su数据 3.1 获取某一道道数据 3.2 获取偶数道数据 3.3 获取奇数道数据 3.4 时窗数据（垂向）SUVLENGTH sutaper suramp 4. sukill、suzero、sunull和sumute–数据充零 4.1 sukill –道充零 4.2 sunull –产生空道数据 4.3 suzero –在一段时间窗内对数据充零 4.4 sumute –数据去除 4.4.1 使用默认的数据 4.4.2测试自定义数据

1. 显示第二道


suwind < data.su key=tracf min=2 max=4 | suximage title="Second Trace"

2. 打印第二道数据

下面这种打印的是乱码信息


suwind< data.su key=tracl min=2 max=2 | cat

用sustrip程序去掉道头，只保留数据块，然后用xxd程序以16进制格式显示每4个字节的数据，最后用sed程序只显示第一行的输出。


# 读取第一道
sustrip < data.su | xxd -c 4 | sed 1q
# 读取第二道
sustrip < data.su | xxd -c 4 | sed 2q | tail -1


00000000: 0000 0000  ....

2.1数据存储到txt

这样只能存储省略的数据


sustrip < data1.su | xxd -c 4 | sed 2q | tail -1 > data.txt

3. 读取具体的su数据


suascii <output14.su >output.txt
suascii < data.su >output.txt

a=2和b=2表示提取第二道数据。

第二行命令使用suascii命令将输出SU文件转换为文本文件output.txt

3.1 获取某一道道数据


suwind < data.su key=tracl accept=8 max=0 > output98.su

3.2 获取偶数道数据


suwind < data.su key=tracl j=2 > output98.su

3.3 获取奇数道数据


suwind < data.su key=tracl j=2 s=1 > output98.su

3.4 时窗数据（垂向）


suwind < data.su tmin=0.1 tmax=0.14 > output98.su


suwind < data.su itmin=1 itmax=10 > output98.su


suwind < data.su nt=10 > output98.su

Options for vertical windowing (time gating):

tmin = 0.0 min time to pass
tmax = (from header) max time to pass
itmin = 0 min time sample to pass
itmax = (from header) max time sample to pass
nt = itmax-itmin+1 number of time samples to pass

SUVLENGTH

将可变长度迹线调整为通用数字样品数量

SUVLENGTH - Adjust Variable Length Traces to a Common Number of Samples


suvlength ns=64 < data.su > data4.su


suxwigb < data.su title="Test of suvlength" &

sutaper

如线性在数据集的每一端逐渐减少 5 条迹线）


sutaper <data.su >output98.su ntr=5

suramp

和“suramp”以平滑迹线的开头和/或结尾（例如，从 0 上升到tmin=.05 秒，并从 tmax=1.15 秒斜坡下降到迹线末尾


suramp <data.su tmin=.05 tmax=1.15 >output98.su

显示


suxwigb < output98.su title="Test of suvlength" &

4. sukill、suzero、sunull和sumute–数据充零

对有噪声的道、位于数据边界的道进行充零是有益的。或者是产生一些空道，在绘图时最为联系数据集中间的分割标志。

4.1 sukill –道充零

/src/su/main/windowing_sorting_muting

对一批地震道充零键入：

sukill <stdin >stdout min=MIN_TRACE count=COUNT


sukill < data-tm-flt-dec-bal-cdp.su > outdata.su min=250 count=10

这里参数count是要充零的总道数，参数min是要充零的这些道中的最小道号。

上述数据显示，CMP范围为231~240


suwind < outdata.su key=cdp min=231 max=240 | suxwigb

4.2 sunull –产生空道数据

有时需要生成道值为零的地震数据块。下面的命令产生共有NTR道，时间采样数为NT的地震数据：

sunull nt=NT ntr=NTR <stdin >stdout min=MIN_TRACE count=COUNT


sunull nt=10 ntr=10 < data.su > output98.su min=5 count=10

上述数据显示，CMP范围为231~240


suxwigb < output98.su

4.3 suzero –在一段时间窗内对数据充零

/src/su/main/amplitudes

suzero itmin=MIN_TIME_SAMPLE itmax=MAX_TIME_SAMPLE <indata.su > outdata.su


suzero itmin=1000 itmax=1015 < data-tm-flt-dec-bal-cdp.su > outdata.su 

suzero itmin=1445 itmax=1480 < data-tm-flt-dec-bal-cdp.su > outdata.su

上述数据显示，CMP范围为231~240


suwind < outdata.su key=cdp min=231 max=240 | suxwigb

4.4 sumute –数据去除

/src/su/main/windowing_sorting_muting

要进行高精度的压制操作，可用命令sumute对SU数据进行处理。

4.4.1 使用默认的数据

下面用suplane生成数据，然后做去除处理与原始生成数据比较：


suplane | suxwigb &

suplane | sumute key=tracl xmute=1,10,12 tmute=.06,.1,.11 | suxwigb &

去除根据xmute=和tmute=参数确定的多边形曲线内每一个初至。

4.4.2测试自定义数据

sumute <indata.su >outdata.su key=KEYWORD xmute=x1,x2,x3,... tmute=t1,t2,t3,...


sumute < data-tm-flt-dec-bal-cdp.su key=cdp xmute=230,235,240 tmute=.2,.6,1.2 | suxwigb &

显示源数据，对比


suwind < data-tm-flt-dec-bal-cdp.su key=cdp | suxwigb
suwind < data-tm-flt-dec-bal-cdp.su key=cdp min=231 max=240 | suxwigb


sushw < data.su key=tracf a=200 > cmp100_time200.su


suchw < cmp100_time200.su | sed -e 's/.*=//' > amplitudes.txt