用多道反褶积方法测定台站接收函数
|
摘要
本文提出了一种在时间域用多道反褶积测定台站接收函数的方法,以提高台站接收函数的测量精度与分辨率.在单道反褶积的基础上,选取若干个质量较好的远震P波波形事件,构成多道信号,以垂直分量为输入,径向和切向分量作为期望输出,依据最小二乘,设计多道滤波器,提取接收函数.合成地震图与观测地震图的检验表明,多道滤波方法能有效地测定台站接收函数,特别是多道反褶积能够有效地恢复地壳上地幔间断面所产生的弱转换波震相.
Time-domain multi-channel deconvolution is put forward to estimate receiver function, to improve resolution. Based on the deconvolution results of individual teleseismic P waveforms, a number of events with good quality are selected to form multi-channel signals, the vertical components are regarded as inputs, and radial and transverse components are taken as expected outputs. Least-square error is used to design the multi-channel filter, to get the common filter factors, i.e., receiver function. Both synthetic and real data experiments show that multi-channel deconvolution is an effective approach to measure receiver function, especially, it can recover the weak converted phases from the upper mantle.
Time-domain multi-channel deconvolution is put forward to estimate receiver function, to improve resolution. Based on the deconvolution results of individual teleseismic P waveforms, a number of events with good quality are selected to form multi-channel signals, the vertical components are regarded as inputs, and radial and transverse components are taken as expected outputs. Least-square error is used to design the multi-channel filter, to get the common filter factors, i.e., receiver function. Both synthetic and real data experiments show that multi-channel deconvolution is an effective approach to measure receiver function, especially, it can recover the weak converted phases from the upper mantle.
引文
[1]Owens TJ,Zandt G,Taylor S R.Seismic evidencefor an ancient rift beneaththe Cumberland Plateau,Tennessee:Adetailed analysis of broadbandteleseismic P waveforms.J.Geophys.Res.,1987,89(B9):7783~7795
[2]刘启元,Rainer K,李顺成.接收函数复谱比的最大或然性估计及非线性反演.地球物理学报,1996,39(4):502~511Liu QY,Rainer K,Li S C.Maximal likelihood estimation and nonlinear inversion of the complex receiver function spectrumratio.Chinese J.Geophys(in Chinese),1996,39(4):502~511
[3]刘启元,Rainer K,李顺成.中国数字地震台网的接收函数及其非线性反演.地球物理学报,1997,40(3):356~368Liu QY,Rainer K,Li S C.The receiver functions at the stations of the Chinese digital seismic network(CDSN)and their nonlinear inversion.Chinese J.Geophys.(in Chinese),1997,40(3):356~368
[4]刘启元,李顺成,沈扬等.延怀盆地及其邻区地壳上地幔速度结构的宽频带地震台阵研究.地球物理学报,1997,40(6):763~773Liu QY,Li S C,Shen Y,et al.Broadband seismic array study of the crust and upper mantle velocity structure beneath Yanhuai basin and its neighbouring region.Chinese J.Geophys.(in Chinese),1997,40(6):763~773
[5]刘启元,陈九辉,李顺成等.新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨.地球物理学报,2000,43(5):356~365Liu QY,Chen J H,Li S C,et al.Passive seismic experiment in Xingjiang-Jiashi strong earthquake region and discussion on its seismic genesis.Chinese J.Geophys.(in Chinese),2000,43(5):356~365
[6]刘启元,Rainer Kind,陈九辉等.大别造山带壳幔界面的断错结构和壳内低速体.中国科学(D),2005,35(4):304~313Liu QY,Rainer K,Chen J H,et al.Dislocation structure of the crust-mantle boundary andlow-velocity body withinthe crust beneath the Dabie Shan collision orogen.Science in China(Series D),2005,48(7):875~885
[7]吴庆举,曾融生.用宽频带远震接收函数研究青藏高原的地壳结构.地球物理学报,1998,41(5):669~679Wu QJ,Zeng RS.The crustal structure of Qinghai-Xizang plateau inferredfrombroadbandteleseismic waveform.Chinese J.Geophys.(in Chinese),1998,41(5):669~679
[8]Zhu L P,Owens T J,Randall G E.Lateral variation in crustal structure of the Northern Tibetan Plateau inferred fromteleseismic receiver functions.Bull.Seismol.Soc.Am.,1995,85(6):1531~1540
[9]Levin V,Park J.Crustal anisotropy beneath the Ural mtns foredeep fromteleseismic receiver functions.Geophys.Res.Lett.,1997,24:1283~1286
[10]Zhu L P,Kanamori H.Moho depth variation in southern California fromteleseismic receiver functions.J.Geophys.Res.,2000,105:2969~2980
[11]Zhu L P.Crustal structure across the San Andreas Fault,SouthernCalifornia fromteleseismic converted waves.Earth and Planetary Sci.Lett.,2000,179:183~190
[12]吴庆举,曾融生,赵文津.喜马拉雅—青藏高原的上地幔倾斜构造与陆-陆碰撞过程.中国科学(D),2004,34(10):919~925Wu QJ,Zeng R S,Zhao WJ.The upper mantle structure of the Tibetan Plateau and its implication for the continent-continent collision.Science in China(Series D),2005,48(8):1158~1164
[13]Ai Yinshuang,Zheng Tianyu,Xu WW,et al.A complex660km discontinuity beneath northeast China.Earth and Planet.Sci.Lett.,2003,212:63~71
[14]Ai Yinshuang,Zheng Tianyu.The upper mantle discontinuity structure beneath eastern China.Geophys.Res.Lett.,2003,30(21):2089,doi:10.1029/2003GL017678
[15]Ai YS,Zhao D P,Gao X,et al.The crust and upper mantle discontinuity structure beneath Alaska inferred from receiver functions.Phys.Earth Planet.Inter.,2005,150:339~350
[16]Tian X B,Wu QJ,Zhang Z J,et al.Joint imaging by teleseismic converted and multiple waves andits applicationinthe INDEPTH-Ⅲpassive seismic array.Geophys.Res.Lett.,2005,32:L2131510.1029
[17]Gurrola H,Baker G E,Minster J B.Simultaneous time-domain deconvolution with application to the computation of receiver functions.Geophys.J.Int.,1995,120:537~543
[18]吴庆举,田小波,张乃铃等.用Wiener滤波方法提取台站接收函数.中国地震,2003,19(1):41~47Wu QJ,Tian X B,Zhang NL,et al.Receiver function estimated by Wiener filtering.Earthquake Research in China(in Chinese),2003,19(1):41~47
[19]吴庆举,田小波,张乃铃等.计算台站接收函数的最大熵谱反褶积方法.地震学报,2003,25(4):382~389Wu QJ,Tian X B,Zhang NL,et al.Receiver function estimated by Wiener filtering.Acta Seismologica Sinica,2003,19(1):41~47
[20]Park J,Levin V.Receiver functions from multiple-taper spectral correlation estimates.Bull.Seismol.Soc.Am.,2000,103(26):899~917
[21]Ligorria J P,Ammon C J.Iterative deconvolution and receiver-function estimation.Bull.Seismol.Soc.Am.,1999,89(1):395~400
[2]刘启元,Rainer K,李顺成.接收函数复谱比的最大或然性估计及非线性反演.地球物理学报,1996,39(4):502~511Liu QY,Rainer K,Li S C.Maximal likelihood estimation and nonlinear inversion of the complex receiver function spectrumratio.Chinese J.Geophys(in Chinese),1996,39(4):502~511
[3]刘启元,Rainer K,李顺成.中国数字地震台网的接收函数及其非线性反演.地球物理学报,1997,40(3):356~368Liu QY,Rainer K,Li S C.The receiver functions at the stations of the Chinese digital seismic network(CDSN)and their nonlinear inversion.Chinese J.Geophys.(in Chinese),1997,40(3):356~368
[4]刘启元,李顺成,沈扬等.延怀盆地及其邻区地壳上地幔速度结构的宽频带地震台阵研究.地球物理学报,1997,40(6):763~773Liu QY,Li S C,Shen Y,et al.Broadband seismic array study of the crust and upper mantle velocity structure beneath Yanhuai basin and its neighbouring region.Chinese J.Geophys.(in Chinese),1997,40(6):763~773
[5]刘启元,陈九辉,李顺成等.新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨.地球物理学报,2000,43(5):356~365Liu QY,Chen J H,Li S C,et al.Passive seismic experiment in Xingjiang-Jiashi strong earthquake region and discussion on its seismic genesis.Chinese J.Geophys.(in Chinese),2000,43(5):356~365
[6]刘启元,Rainer Kind,陈九辉等.大别造山带壳幔界面的断错结构和壳内低速体.中国科学(D),2005,35(4):304~313Liu QY,Rainer K,Chen J H,et al.Dislocation structure of the crust-mantle boundary andlow-velocity body withinthe crust beneath the Dabie Shan collision orogen.Science in China(Series D),2005,48(7):875~885
[7]吴庆举,曾融生.用宽频带远震接收函数研究青藏高原的地壳结构.地球物理学报,1998,41(5):669~679Wu QJ,Zeng RS.The crustal structure of Qinghai-Xizang plateau inferredfrombroadbandteleseismic waveform.Chinese J.Geophys.(in Chinese),1998,41(5):669~679
[8]Zhu L P,Owens T J,Randall G E.Lateral variation in crustal structure of the Northern Tibetan Plateau inferred fromteleseismic receiver functions.Bull.Seismol.Soc.Am.,1995,85(6):1531~1540
[9]Levin V,Park J.Crustal anisotropy beneath the Ural mtns foredeep fromteleseismic receiver functions.Geophys.Res.Lett.,1997,24:1283~1286
[10]Zhu L P,Kanamori H.Moho depth variation in southern California fromteleseismic receiver functions.J.Geophys.Res.,2000,105:2969~2980
[11]Zhu L P.Crustal structure across the San Andreas Fault,SouthernCalifornia fromteleseismic converted waves.Earth and Planetary Sci.Lett.,2000,179:183~190
[12]吴庆举,曾融生,赵文津.喜马拉雅—青藏高原的上地幔倾斜构造与陆-陆碰撞过程.中国科学(D),2004,34(10):919~925Wu QJ,Zeng R S,Zhao WJ.The upper mantle structure of the Tibetan Plateau and its implication for the continent-continent collision.Science in China(Series D),2005,48(8):1158~1164
[13]Ai Yinshuang,Zheng Tianyu,Xu WW,et al.A complex660km discontinuity beneath northeast China.Earth and Planet.Sci.Lett.,2003,212:63~71
[14]Ai Yinshuang,Zheng Tianyu.The upper mantle discontinuity structure beneath eastern China.Geophys.Res.Lett.,2003,30(21):2089,doi:10.1029/2003GL017678
[15]Ai YS,Zhao D P,Gao X,et al.The crust and upper mantle discontinuity structure beneath Alaska inferred from receiver functions.Phys.Earth Planet.Inter.,2005,150:339~350
[16]Tian X B,Wu QJ,Zhang Z J,et al.Joint imaging by teleseismic converted and multiple waves andits applicationinthe INDEPTH-Ⅲpassive seismic array.Geophys.Res.Lett.,2005,32:L2131510.1029
[17]Gurrola H,Baker G E,Minster J B.Simultaneous time-domain deconvolution with application to the computation of receiver functions.Geophys.J.Int.,1995,120:537~543
[18]吴庆举,田小波,张乃铃等.用Wiener滤波方法提取台站接收函数.中国地震,2003,19(1):41~47Wu QJ,Tian X B,Zhang NL,et al.Receiver function estimated by Wiener filtering.Earthquake Research in China(in Chinese),2003,19(1):41~47
[19]吴庆举,田小波,张乃铃等.计算台站接收函数的最大熵谱反褶积方法.地震学报,2003,25(4):382~389Wu QJ,Tian X B,Zhang NL,et al.Receiver function estimated by Wiener filtering.Acta Seismologica Sinica,2003,19(1):41~47
[20]Park J,Levin V.Receiver functions from multiple-taper spectral correlation estimates.Bull.Seismol.Soc.Am.,2000,103(26):899~917
[21]Ligorria J P,Ammon C J.Iterative deconvolution and receiver-function estimation.Bull.Seismol.Soc.Am.,1999,89(1):395~400
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心