云南中小地震矩震级测定.docVIP

云南中小地震矩震级测定 　　摘要：利用云南省测震台网的数字地震波形资料，反演了2008年1月至2012年9月675个中小地震的震源波谱参数，并计算了这些地震的地震矩M0和矩震级MW，最后利用回归分析方法得到了近震震级ML和MW之间的关系式。 关键词：中小地震；震源波谱参数；地震矩；矩震级；回归分析法 中图分类号：P31563文献标识码：A文章编号：1000-0666（2013）03-0306-07 0引言 震级是表征地震大小的量，是地震的基本参数之一，是地震预报和其它有关地震的研究中的一个重要参数。Richter（1935）在研究美国南加州的地震时引入了地方性震级ML，现称作里克特震级（简称里氏震级），尽管里氏震级完全是经验性的，物理意义不清楚，但使用方便，更重要的是它为震级的发展奠定了基础。 现在常用的震级主要有：近震震级ML、面波震级MS、体波震级mb和mB、持续时间震级MD、矩震级MW。与传统的震级标度相比，矩震级MW不仅物理含义较明确，而且不会出现震级饱和的现象，是目前量度地震大小的最理想的物理量（高景春等，2011）。它是一个均匀的震级标度，适用于震级范围很宽的统计，无论对大震、小震、微震，甚至极微震，还是是对浅震、深震，均可测量其地震矩，在一定震级范围内能与常规的面波震级转换（陈运泰，刘瑞丰，2004）。 本文利用云南省测震台网的近场数字地震波形资料，采用多台联合反演方法计算了2008年1月至2012年9月的ML30以上675个中小地震的震源波谱参数，得到了这些地震的地震矩M0和矩震级MW，并利用回归分析方法研究了近震震级ML和MW之间的关系。 1计算方法111原理从震源辐射出的地震波，经过地球介质传播到地震观测台，既携带着震源的信息，又携带着从震源到地震台站之间的传播介质的信息，因此地震仪记录到的地震波是一种综合信息，包含了地震震源效应、传播路径效应、台站场地响应及仪器响应。要测定地震矩参数，就必须首先扣除地震记录中传播路径效应、台站场地响应及仪器响应等，即先求出震源谱。 在任一台站观测的任一地震地面运动的傅里叶谱可以表示为Aij（f）=Ai0（f）·G（Rij）·Sj（f）·Ij（f）·e-πRijf1Q（f）β（1）式中，Aij（f）是在第j个台观测到的第i个地震的傅里叶振幅谱（观测谱）；Ai0（f）为第i个地震的震源振幅谱（震源谱）；Rij为震源距（第i个地震到第j个台站）；G（Rij）为几何衰减函数；Q（f）为频率依赖的品质因子；β为地震波速度（本文仅分析S波）；Sj（f）为第j个台站的场地响应项；Ij（f）为第j个台站的仪器响应项。 由式（1）可见，在频率域内对台站的观测位移谱Aij（f）进行几何衰减、非弹性衰减、场地响应及仪器响应校正后，才能得到震源谱Ai0（f）。笔者采用三段几何衰减模型来计算传播路径的影响（Atkinson，Mereu，1992；黄玉龙等，2003），采用Moya等（2000）关于多台、多地震联合反演的方法来计算台站的场地响应，仪器响应可通过仪器标定进行校正。 通过以上计算，对每个地震、每个台站记录的位移谱逐一校正后，就得到了每个地震多个台站的震源谱，求其平均值以消除方向性效应，得到每个地震的平均震源谱。这样对于某个地震，就可以得到第j个台站记录的该地震的震源位移谱[A0（f）]j，再对台站得到的震源谱求平均，可得到该地震的震源谱A0（f）。采用Brune（1970）的二次方衰减震源谱模型作为理论震源谱Atheo（f），利用遗传算法求解由A0（f）和Atheo（f）定义的残差为极小的震源谱参数Ω0和fc，其中，Ω0是零频振幅，fc是拐角频率。 1地震研究136卷第3期杨晶琼等：云南中小地震矩震级的测定由得到的震源谱参数Ω0及拐角频率fc，就可计算地震矩（Brune，1970）：M0=4πρβ3Ω012Rθφ（2）式中，ρ是密度，取29 g/cm3；β是S波速度，取35 km/s；Rθφ是辐射花样系数，由于没有每次地震的断层面解，所以令Rθφ为一常数，并取S波在震源球上的平均值为2/5。 矩震级MW由地震矩M0通过下式计算得到（Kanamori，1977；Hanks，Kanamori，1979）：MW=213lgM0-607 （3）式中，地震矩M0的单位为N·m。 12参数选取 （1）Q值 在计算台站观测震源谱时，必须首先扣除地震记录中的传播路径效应、台站场地响应及仪器响应等，其中最为关键的是介质品质因子Q。Q是地球介质的重要物理参数，描述由于介质的非完全弹性，地震波在介质里传播过程中的衰减，或能量耗损的特征，其定义为11Q=112π ΔE1E（4）式中，E为一定体积的地球介质在地震波一个周期T