用Wiener滤波方法提取台站接收函数  

Receiver Function Estimated by Wiener Filtering

作　　者：吴庆举[1] 田小波[1] 张乃铃[1] 李卫平[1] 曾融生[1] 

机构地区：[1]中国地震局地球物理研究所,北京100081

出　　处：《中国地震》2003年第1期41-47,共7页Earthquake Research in China

基　　金：国家自然科学基金资助项目 (4 9974 0 2 1);中国地震局地球物理研究所论著号 0 2AC10 36

摘　　要：本文提出了一种在时间域用Wiener滤波方法提取台站接收函数的方法 ,用远震P波波形的垂直分量为输入 ,接收函数作为滤波因子 ,远震P波波形的径向和切向分量作为期望输出 ,通过期望输出与实际输出的均方误差达极小 ,来提取接收函数。接收函数的计算可归结为Toeplitz方程的求解 ,可以采用Levinson递推算法。Toeplitz方程的非奇异性保证了Wiener滤波反褶积方法的稳定性。合成地震图与观测地震图的检验表明 ,用Wiener滤波方法测定台站接收函数是一种有效的时间域反褶积方法。Wiener filtering is provided to estimate receiver function in time-domain. With the vertical component of 3-component teleseismic P waveform as the input of Wiener filter, receiver function as the filter response, and radial and tangential components as the expected outputs, receiver function is estimated by minimizing the error between expected and actual outputs. Receiver function can be obtained by solving Toeplitz equation through Levinson algorithm. The non-singularity of Toeplitz equation keeps the stability of Wiener deconvolution. Both synthetic and observational seismogram checks show that Wiener deconvolution is an effective time-domain method to estimate receiver function from teleseismic P waveform.

关 键 词：Wiener滤波方法 接收函数 时间域 远震P波 波形 合成地震图 观测地震图 

分 类 号：P315[天文地球—地震学]