地震地质 ›› 2004, Vol. 26 ›› Issue (3): 416-425.

• • 上一篇    下一篇

分离三分量远震接收函数的多道最大或然性反褶积方法

刘启元1, Rainer Kind2   

  1. 1. 中国地震局地质研究所, 地震动力学国家重点实验室, 北京, 100029;
    2. GeoForschung Zentrum Potsdam, Telegrafenberg, 14473 Potsdam, Germany
  • 收稿日期:2004-05-11 修回日期:2004-06-06 出版日期:2004-09-02 发布日期:2009-10-26
  • 作者简介:刘启元,男,1945年生,1969年毕业于北京大学地球物理系,1982年毕业于中国科技大学研究生院并获理学硕士学位,研究员,博士生导师,主要从事理论及应用地震学研究,电话:010-62009311,E-mail:qyliu@public.bta.net.cn.
  • 基金资助:
    国家自然科学基金(40074009)和德国科学基金会(DFG)联合资助。中国地震局地质研究所论著2004B0030。

MULTI-CHANNEL MAXIMAL LIKELIHOOD DECONVOLUTION METHOD FOR ISOLATING 3-COMPONENT TELESEISMIC RECEIVER FUNCTION

LIU Qi-yuan1, Rainer Kind2   

  1. 1. State Key Laboratory of Earthquake Dynamics, Institute of Geology, China Earthquake Administration, Beijing 100029, China;
    2. GeoForschungs Zentrum Potsdam, Telegrafenberg 14473 Potsdam, Germany
  • Received:2004-05-11 Revised:2004-06-06 Online:2004-09-02 Published:2009-10-26

摘要: 迄今为止,广泛应用的分离接收函数方法主要基于单个台站的观测数据,并依赖等效震源假定的有效性。当地壳内存在明显的间断面时,等效震源假定不适于宽频带地震波形数据,接收函数的估计将在不同程度上遭到破坏。解决这个问题的出路在于避免使用等效震源假定。依据多道最大或然性反褶积原理,给出了利用多道观测数据分离接收函数的方法。我们的方法避免了等效震源假定,这对于在复杂构造环境条件下改善接收函数的估计具有实际价值。合成地震图的数值检验结果验证了给出的方法。利用该方法,我们得到了跨越大别山造山带,从江西大冶附近的大箕铺(30°20′N,115°03′E)到河南兰考附近的崔林(34°40′N,114°49′E),总长度约500km的流动地震台阵剖面各台站的三分量接收函数估计。与基于单台数据的接收函数径向分量比较表明,两种方法给出的结果可能存在不可忽略的差别。

关键词: 接收函数, 最大或然性, 反褶积, 大别山

Abstract: Up to now,all of the methods of isolating the receiver function from teleseismic P waves are based on the data recorded at a single station. However,the estimation of the receiver function based on the so called source equalization scheme will not be valid for broadband data,when a sharp discontinuity exists within the crust. The outlet for solving this problem is to avoid using the source equalization assumption. In this study,we propose a method for isolating the 3 component receiver function from multi channel seismic data based on multi channel maximal likelihood deconvolution princple. The source equalization scheme is avoided in our method,and this is valuable for improving the estimation of the receiver function in case of complex crustal structure. The numerical test by synthetics proves our method to be effective. Using this new method,we obtain the 3 component receiver functions at stations along the passive seismic profile across the Dabieshan orogenic belt from Daqipu (30°20′N,115°03′E) nearby Daye,Jiangxi Province to Cuilin(34°40′N,114°49′E)nearby Lankao,Henan Province. In comparison with the radial components of the receiver function from the data of a single station,it is indicated that the results given by both methods have un negligible differences.

Key words: maximal likelihood, deconvolution, receiver function, Dabieshan

中图分类号: