fk一维变换滤波c语言程序

F-K滤波法压制多次波的基本原理为：将地震数据由时间-空间域（T-X）变换到频率-波数域(F-K)，根据F-K域内一次波和多次波能量分布特征在F-K域中进行多次波的滤波分离。此方法在考虑多次波频谱特性的同时，兼顾多次波信号在对应剖面的波长等特性，相比频率域切除，该方法通过增加波数域限制可以有效减少切除过程中有效波的损失。 Ryu（1982）首次利用F-K变换对多次波进行滤波处理，发现F-K滤波需要保证合适的压制带。F-K二维滤波不能简单地通过直线在频率波数域内进行一次波和多次波的分离，直线边界会引入Gibbs现象，为了有效选取切除函数，Zhou和Greenhalgh（1994）提出了非线性F-K滤波方法，利用波场外推之后的多次波来自动确定切除函数，但在近偏移距无法获取理想的数据分离效果。利用F-K滤波法压制多次波的最重要步骤是选取处于一次波和多次波速度之间的校正速度，校正速度选取不当则可能会影响一次波信息的完整性，导致原始数据信噪比的降低。