频率波数域matlab,频率-波数域方法的发展及其在台阵数据分析中的应用

0 引言

自1960年起，地震台阵作为一种新型的地震学研究工具，已经广泛应用于不同类型波场信号的检测及地球内部不同尺度结构的成像(

传统的F-K方法(

各种不同类型台阵的陆续建设，为F-K分析提供了大量的数据支持，使得该方法在理论和实践中都得到了极大的发展。人们对F-K方法的不断改进(

1 F-K分析方法

从原理上来说，F-K方法是基于频率波数域功率谱密度函数(F-K PSDF)的计算，将台阵观测信号从时间-空间域转换到频率-波数域，从而获取入射波场能量在不同慢度和方位角的分布。根据窗函数的不同，可分为聚束法(

1.1 F-K谱

$

E{P_b}({\rm{f,k}}) = \sum\limits_{n,m = 1}^N { {\phi _{nm}}} \exp \left[ { {\rm{i}}k\left( { { {\rm{X}}_m} - { {\rm{X}}_n}} \right)} \right]

$

(1a)

$

E{P_m}({\rm{f,k}}) = {[ {\sum\limits_{n,m = 1}^N {\phi _{nm}^{ - 1}} \exp \left[ { {\rm{i}}k\left( { { {\rm{X}}_n} - { {\rm{X}}_m}} \right)} \right]^{ - 1}} }

$

(1b)

此处，f表示频率，k为波数，N是台阵中所包含的地震台站数目，φnm是第n和第m个台站记录之间的互功率谱，Xn和Xm分别为第n和第m个台站的坐标。

上述两种方法最主要的不同是窗函数。聚束法在波数k0的窗函数可表示为

$

{W_b}(k,{k_0}) = {\rm{ }}\frac{1}{ { {N^2}}}\sum\limits_{n,m = 1}^N {\exp \left\{ {i\left( {k - {k_0}} \right)\left( { {X_n} - {X_m}} \right)} \right\}}

$

(2)

最大似然法的窗函数为

$

{W_m}(k,{k_0}) = {\rm{ }}\left| {\sum\limits_{j = 1}^N { {A_j}(f,{k_0})} } \right|{W_b}(k,{k_0})

$

(3)

其中

$

{A_j}(f,{k_0}) = \sum\limits_{n = 1}^N { { {({\phi _{jn}}\exp\{ i{k_0}({X_j} - {X_n})\} )}^{ - 1}}} {\rm{/}