F-K数字滤波器的设计和应用

将地震记录在二元正交基中表示为则其二维谱为然后将不需要的圆频率ω和圆波数K_x所对应的谱值充零,再用上述第一个公式作反傅氏变换,便完成了二维滤波。在程序设计中,根据公式A_p（K·2^p+j）=A_p-1(k·2^p-1+j)+A_p-1(K·2^p-1+j+2^m-1) A_p(K·2^p+j+2^p-1)=A_p-1(K·2^p-1+j)-A_p-1(K·2^p-1+j+2^m-1)]W_N^K·2p-1排列成快速傅氏变换子程序FFT,调用FFT对地震道分别进行纵向和横向快速傅氏变换,便得二维谱。程序具有时变和空变的功能。为了加强二维滤波的功能,还设计了F-K陷波器。陷波器的主要参数是斜率和频率,斜率参数既可保留正或负视速度波组,也可同时保留在正、负视速度域范围内的波组。频率取值应该比一维滤波的取值范围稍大一些。当要保留高速波组时,高频数值应取大一些;当要保留低速波组时,低频数值应取小一些。为了提高分辨率,频率采样间隔应小于7.8赫兹,波数采样间隔应小于0.0003125/米。