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

一、绘制频谱曲线

1. 首先我们需要一个数据，我们构造出一个模拟单炮数据data.su

注意：

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

二、使用sufilter进行滤波

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

1. 首先，我们sufilter的功能介绍，在命令行输入 ， 我们可以看到sufilter是应用零相位的滤波器，应用时必须在道头中有参数dt! 可选择的参数有 f 和 amps ，而且f和amps的参数是对应的，也就是说每一个频率对应一个滤波器的振幅值。
2. f和amps各有两个参数时，例如:在命令行输入：
3. f 和amps各有四个参数时，就可随意构成带通，带阻，低通和高通的滤波器，例如一个30-40Hz的带通滤波器：

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

注意：

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

四、进行fk滤波

4.1用sudipfilt进行倾角滤波

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

文章知识点与官方知识档案匹配，可进一步学习相关知识CS入门技能树Linux入门初识Linux25045 人正在系统学习中