非对称地震矩张量时间域反演:理论与方法
摘要
本文在对称地震矩张量反演的基础上,进一步研究了非对称地震矩张量时间域反演的理论与方法,结果表明:非对称地震矩张量反演与对称地震矩张量反演类似,只需将对称地震矩张量反演方法略加改动,即增加3个待解参数,便可实现非对称地震矩张量反演.为了判断非对称地震矩张量反演相对于对称地震矩张量反演是否存在过度拟合,运用了AIC准则(赤池信息准则).为了定量地描述地震矩张量之间的差异,引入了地震矩张量的矢量表示法.通过分析格林函数与地震矩张量各分量之间的关系,得出:在非对称地震矩张量反演时,若仅用垂直向地动位移数据,将无法区分M_(xy)与Myx这两个分量,需要同时运用垂直向与水平向地动位移数据进行联合反演才能区分M_(xy)与Myx;若采用不同的速度结构模型或不同的格林函数计算方法,则需重新评估地震矩张量各分量的分辨度问题.为检验非对称地震矩张量反演方法的可行性,利用合成地震图进行了一系列数值试验.数值试验结果表明,在非对称地震矩张量反演中,有必要引入S波进行P波与S波联合反演以提高反演的准确性和判定断层面的能力.
In this paper the theory and method of time-domain inversion for the asymmetric seismic moment tensor are studied based on the symmetric seismic moment tensor inversion.It is shown that the asymmetric seismic moment tensor inversion is analogous to the symmetric one,and that only a minor modification on the method of symmetric seismic moment tensor inversion,i.e.,increasing three unknown parameters is required to realize the asymmetric seismic moment tensor inversion.The Akaike Information Criterion(AIC)is introduced to judge whether the asymmetric seismic moment tensor inversion is over fitting as compared with the symmetric seismic moment tensor inversion.In analyzing the asymmetric seismic moment tensor the vector representation of the moment tensor is introduced to describe the difference between seismic moment tensors quantitatively.By analyzing the relation between Green's function and each of moment tensor's components,it is found that,in asymmetric seismic moment tensor inversion,the components M_(xy)and Myxcould not be distinguished if only vertical displacement data were used in inversion.The vertical and horizontal displacement data should be used jointly in inversion to distinguish M_(xy)from Myx.It is also found that if a different velocity structure model or a different method for calculating the Green's function was adopted for the inversion,the resolutions between any seismic moment tensor components need to be re-analyzed.A series of numerical tests were performed to test the feasibility of the present inversion method by using synthetic data.It is shown that the present inversion method is feasible and that it is necessary to use jointly both the P and S wave data to invert for the asymmetric seismic moment tensor.
引文
李旭.1993.用地震波波形资料反演1990年青海共和地震的震源过程[D].北京:中国地震局地球物理研究所:1-59.Li X.1993.Inversion of Seismic Body Wave Data for the Source Process of the Gonghe,Qinghai,China Earthquake[D].Beijing:Institute of Geophysics,China Earthquake Administration:1-59(in Chinese).
李旭,陈运泰,王培德.1994.水平层状介质中理论地震图计算程序的使用说明[R].北京:中国地震局地球物理研究所:1-20.Li X,Chen Y T,Wang P D.1994.User’s Manual of Computer Program of Theoretical Seismogram in Horizontally Layered Media[R].Beijing:Institute of Geophysics,China Earthquake Administration:1-20(in Chinese). 李旭,陈运泰.1996a.合成地震图的广义反射透射系数矩阵方法[J].地震地磁观测与研究,17(3):1-20.Li X,Chen Y T.1996a.The generalized reflection-transmission coefficient matrix method for synthetic seismograms[J].Seismological and Geomagnetic Observation and Research,17(3):1-20(in Chinese). 李旭,陈运泰.1996b.用长周期地震体波波形资料反演1990年青海共和地震的震源过程[J].地震学报,18(3):279-286.Li X,Chen Y T.1996b.Inversion of long period body wave data for the source process of the Gonghe,Qinghai,Chinaearthquake[J].Acta Seismologica Sinica,9(3):4-12. 刘超,张勇,许力生,陈运泰.2008.一种矩张量反演新方法及其对2008年汶川MS8.0地震序列的应用[J].地震学报,30(4):329-339.Liu C,Zhang Y,Xu L S,Chen Y T.2008.A new technique for moment tensor inversion with applications to the 2008Wenchuan MS8.0earthquake sequence[J].Acta Seismologica Sinica,21(4):333-343. 刘超.2011.断层厚度的地震效应和非对称矩张量[D].北京:中国地震局地球物理研究所:1-110.Liu C.2011.Seismic Effect of Fault Thickness and Asymmetric Seismic Moment Tensor[D].Beijing:Institute of Geo-physics,China Earthquake Administration:1-110(in Chinese). 刘超,陈运泰.2014.断层厚度的地震效应和非对称地震矩张量[J].地球物理学报,57(2):509-517.Liu C,Chen Y T.2014.Seismic effect of fault thickness and asymmetric seismic moment tensor[J].Chinese Journal ofGeophysics,57(2):509-517(in Chinese). Akaike H.1974.A new look at the statistical model identification[J].IEEE Trans Automat Contr,19(6):716-723. Aki K,Richards P G.1980.Quantitative Seismology:Theory and Methods[M].San Francisco:W H Freeman:1-932. Backus G,Mulcahy M.1976a.Moment tensors and other phenomenological descriptions of seismic sources:I.Continuous displacements[J].Geophys J R A S,46(2):341-361. Backus G,Mulcahy M.1976b.Moment tensors and other phenomenological descriptions of seismic sources:II.Discontinuous displacements[J].Geophys J R A S,47(2):301-329. Burridge R,Knopoff L.1964.Body force equivalents for seismic dislocations[J].Bull Seismol Soc Am,54(6A):1875-1888. Dziewonski A M,Chou T A,Woodhouse J H.1981.Determination of earthquake source parameters from waveform data for studies of global and regional seismicity[J].J Geophys Res,86(B4):2825-2852. Ekstr9m G.1989.A very broad band inversion method for the recovery of earthquake source parameters[J].Tectonophysics,166(1/2/3):73-100. Fitch T J,McCowan D W,Shields M W.1980.Estimation of the seismic moment tensor from teleseismic body wave data with applications to intraplate and mantle earthquakes[J].J Geophys Res,85(B7):3817-3828. Gilbert F,Dziewonski A M.1975.An application of normal mode theory to the retrieval of structural parameters and source mechanisms from seismic spectra[J].Phil Trans Roy Soc London,278(1280):187-269. Heaton T H.1990.Evidence for and implications of self-healing pulses of slip in earthquake rupture[J].Phys Earth Planet Interi,64(1):1-20. Hudson J A,Pearce R G,Rogers R M.1989.Source type plot for inversion of the moment tensor[J].J Geophys Res,94(B1):765-774. Kagan Y Y,Knopoff L.1985.The two-point correlation function of the seismic moment tensor[J].Geophys J R A S,83(3):637-656. Kagan Y Y.1991.3-D rotation of double-couple earthquake sources[J].Geophys J Int,106(3):709-716. Kennett B L N,Kerry N J.1979.Seismic waves in a stratified half space[J].Geophys J R A S,57(3):557-583. Kennett B L N.1980.Seismic waves in a stratified half space:II.Theoretical seismograms[J].Geophys J R A S,61(1):1-10. Kennett B L N.1983.Seismic Wave Propagation in Stratified Media[M].Cambridge:Cambridge University Press:1-339. Knopoff L,Chen Y T.2009.Single-couple component of far-field radiation from dynamical fractures[J].Bull Seismol Soc Am,99(2B):1091-1102. McCowan D W.1976.Moment tensor representation of surface wave sources[J].Geophys J R A S,44(3):595-599. Nielsen S,Madariaga R.2003.On the self-healing fracture mode[J].Bull Seismol Soc Am,93(6):2375-2388. Takei Y,Kumazawa M.1994.Why have the single force and torque been excluded from seismic source models?[J].Geophys J Int,118(1):20-30. Takei Y,Kumazawa M.1995.Phenomenological representation and kinematics of general seismic sources including the seismic vector modes[J].Geophys J Int,121(3):641-662. Willemann R J.1993.Cluster analysis of seismic moment tensor orientations[J].Geophys J Int,115(3):617-634. |