关于SU（SeismicUnix）的使用——数字信号处理部分

一、绘制频谱曲线

首先我们需要一个数据，我们构造出一个模拟单炮数据data.su
用suwind将数据变成单道数据,这里我从单炮数据中抽取关键字tracl=60的数据： suwind<data.su key=tracl max=60 min= 60>singleline.su
对singleline.su作傅里叶变换（sufft），将振幅谱显示出来即可。命令行代码如下：

sufft<singleline.su |suamp mode=amp|suxgraph style=normal

注意：

仅仅使用sufft和suamp，而没有输出或显示是不可以的！会报错：’output can’t be tty‘，需要用suxgraph进行显示或者进行输出，就不会报错了！
通常将sufft与suamp组合起来使用，那么就可以将快速傅里叶变化后的结果中的mode=振幅（amp）、相位（phase）、实部和虚部(real/imag)显示或者输出。

二、使用sufilter进行滤波

*注：*在/su/src/demos/Filtering 中可以看到su自带的演示例子，但是由于demos里面需要OZ Yilmaz.H 和 .F 文件，但是在CWP的网页中无法打开对应链接进而下载数据。所以这可能就需要自己用su去构造一个模拟数据，然后分析demos里面的shell，对shell做一些修改并再应用了。

首先，我们sufilter的功能介绍，在命令行输入 sudoc sufilter，我们可以看到sufilter是应用零相位的滤波器，应用时必须在道头中有参数dt! 可选择的参数有 f 和 amps ，而且f和amps的参数是对应的，也就是说每一个频率对应一个滤波器的振幅值。
f和amps各有两个参数时，例如:在命令行输入：

sufilter< single.su f=30,40 amps=1,0 > tmp.su
//这行代码是说对 single.su 进行‘低通’滤波处理，
//滤波器的参数为小于30Hz的为幅值为1，大于40HZ的部分幅值为零，30-40Hz之间的频率即为滤波器的过渡带。

f 和amps各有四个参数时，就可随意构成带通，带阻，低通和高通的滤波器，例如一个30-40Hz的带通滤波器：

sufilter< single.su f=20, 30,40, 50 amps=0,1,1,0 > tmp.su

由此可见，sufilter的f=f1,f2,f3,…和amps=a1,a2参数表示的是:在频谱曲线上各个对应点(f1,a1), (f2,a2), (f3,a3), … 这些点的连线构成的滤波器，在这些点之外的曲线与第一个点和最后一个点的幅值相同。

最后注意：

sufilter 的f与amps 只是用一对参数是不允许的
对tmp.su可以使用上面一绘制频谱曲线的方法进行绘制频谱曲线，来查看效果
参数verbose待补充！参数verbose待补充！

三、绘制fk谱

OZ. Yilmaz 在他的书《Seismic data analysis》的二维傅里叶变换中说：
multichannel processing operation can be loosely defined as those that must operate on several data traces, simultaneously.

我们就用su/src/demos 中的例子来说明fk谱的问题，首先用suplane 构造一个地震剖面名字为suplane.su，如下图（左）所示，而右图即为fk域振幅谱。
suspecfk<suplane.su |suxwigb，这样就得到了右图

注意：

fk谱图中的斜率就表示视速度，当倾角越大时，视速度越小，对应的斜率也越小，这一点表现在上图（右）中就是：图中有四条短折线，但是对应三个斜率（其中最外侧两条线的斜率相同）。中间一根竖线斜率无穷大，对应的是左图中水平的同向轴（同向轴倾角为零）；中间那条同向轴斜率中等，我们将在下面将它去除；而倾角最大的同向轴，对应了fk图中最外侧的两个斜率相同的小折线，它产生两条斜率相同的折线，是因为产生了空间假频（spatial aliasing）
对于fk频谱图（右图），要翻转过来看，因为在图中y轴是向下为正的！所以左下角的一小段折线就表示了斜率最大处的空间假频。

四、进行fk滤波

4.1用sudipfilt进行倾角滤波

在上面的suplane.su的基础上，我们用sudipfilt进行倾角滤波，我们滤掉中间的那个同向轴。代码如下

slopes=-1.0,-0.5,0.5,1.0
amps=1,0,0,1
bias=0.0
sudipfilt=suplane.su dt=1 dx=1 slopes=$slopes amps=$amps bias=$bias >test.su
#注意slopes必须和amps对应，即这里个参数的个数必须相同，而且slopes必须单调递增
#对于sudipfilt，dt和dx的参数在道头中必须有！

上面的图可以看出，我们利用sudipfilt将中间的那一个地震道去除，而从fk谱图上面看，波数k为正的部分的中间的一个倾斜直线消失。