接收函数复谱比的最大或然性估计及非线性反演

利用Ｓｈａｍｗａｙ和Ｄｅｒ的多道最大或然性反褶积原理，讨论了从单台三分量远震Ｐ波波形中分离接收函数径向与垂向分量复谱比的方法．根据Ｔａｒａｎｔｏｌａ的波形反演理论发展了接收函数复谱比的非线性反演方法．本文的反演方法除了拟合接收函数的复谱比之外，还需拟合时间域中接收函数垂向与径向分量的初至振幅比．合成波形数据反演的结果表明本文方法的反演结果不依赖初始模型．利用该方法和ＣＤＳＮ台网兰州台记录的宽频带远震Ｐ波波形数据，研究了该台站的接收函数及其随方位角的变化．接收函数复谱比非线性反演给出了该台站下方１４０ｋｍ深度岩石层的Ｓ波速度结构，得到兰州台下方地壳厚度为５４ｋｓ，岩石层厚度为９４ｋｍ．

MAXIMAL LIKELIHOOD ESTIMATION AND NONLINEAR INVERSION OF THE COMPLEX RECEIVER FUNCTION SPECTRUM RATIO

In this study, we describe a new method of estimating the complex spectrum ratio between the radial and vertical components of the receiver function from teleseismic P waveforms recorded by a three-component station.Our method is based on the multi-channel maximal likelihood deconvolution theory by Shumway and Der (1985).According to the theory by Tarantola(1987) we also present a nonlinear inversion algorithm of the complex spectrum ratio of the receiver function.Besides modeling the complex spectrum ratio,the onset amplitude ratio between the vertical and radial components of the receiver function in time domain is fitted also in our inversion procedore.The synthetic inversion test has demonstrated that the inverted model by our method is independent on the initial models.AS an example,the lithospheric receiver functions at the CDSN-LZH station have been obtained by using our method. Their variation with azimuths have been found. The newly developed inversion method has been used for investigating the fine lithosphere S velocity structure beneath the CDSN-LZH station.Our results show 54km thick crust and 94km thick lithosphere.