资源描述：
《1982年卢龙ms6．1地震余震序列震源机制的矩张量反演》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、1982年卢龙Ms6．1地震余震序列震源机制的矩张量反演第28卷第2期2006年3月(111～118)地震ACTASEISM0L0G1CASINICAVo】.28,NO.2Mar.,2006文章编号:0253—3782(2006)02—011卜O81982年卢龙Ms6.1地震余震序列震源机制的矩张量反演李文军王培德¨陈棋福1)中国北京100081中国地震局地球物理研究所2)中国北京100036中国地震局地震预测研究所摘要用矩张量反演的方法,通过对1982年1O月19日河北卢龙Ms6.1地震的余震序列
2、的震源机制求解,认为卢龙盆地的主压应力轴方位角为N74.E,近东西向;而盆地的北边变为N43.E;在盆地中心有些部位主压应力轴转为北西方向,与盆地周边的主压应力轴方位几乎垂直.这说明在卢龙地震中构造的不同部位地壳应力方位很不一致,反映出卢龙地区构造活动的复杂性.一方面,由于卢龙位处中国大陆东部,受到来自东部的日本海盆地在太平洋板块驱动力作用下的挤压;另一方面,卢龙又夹持在燕山地块与华北板块之间,因而受块体边界的限制也是很明显的.关键词卢龙盆地矩张量反演余震序列主压应力轴构造应力场中图分类号:P315
3、.33文献标识码:A引言卢龙位于河北省唐山市东北部约70km,1982年1O月19日在卢龙县城东北方向约4km处发生了6.1地震.据中国地震台网测定,这次地震的发震时刻为北京时间1O月19日2O时45分,震中位于39.57N,l19.O4E.按照中国大陆活动块体划分方法(张培震等,2003),卢龙处于燕山块体与华北平原块体的结合部位,在渤海一张家口一带存在一个活动强烈的WNW向的地震构造带(张家口一渤海地震带,以下简称张渤带).该地震构造带由多条WNW向延伸的活动断裂分布组成,1976年发生7.8级
4、强烈地震的唐山就位于这一地震带的东部.唐山地震,其后发生的滦县7.1级,宁河6.9级地震,以及本次卢龙地震.在空问上联系起来呈NE向展布,构成了唐山地震活动分布区.卢龙地区位于该活动区的北部端点,再往北地震很少发生.因此,卢龙实际上是处于唐山断裂带与张渤带的东部交叉部位.了解卢龙地区的构造应力分布,有利于理解河北东部地区的地震成因和地震构造.卢龙地震发生后,中国地震局地球物理研究所组织了对余震的观测(王培德,1987).使用5台DCS～302地震仪,围绕震中进行了近场观测,观测时间为1O月21—29
5、日.地震仪架设的地点为卢龙县城(XCH),康时各庄(KSG),下寨(XZH),白庄子(BZZ)和相公庄(XGZ),台站分布见图1.DCS-302为数字化地震数据记录仪,字长为12位,加上4个放*国家自然科学基金(40234038)资助.中国地震局地球物理研究所论着06AC10052004—12—08收到初稿,2005—09—20决定采用修改稿.十E-mail:liwjI@seis.ac.cn112地震28卷图1卢龙地震余震震中和台站分布F.～F5为根据震中分布确定的断层大倍率1,5,25和100的自
6、动增益转换,总的动态范围为112dB.采样速率为每信道100点/秒.与SSA一302伺服加速摆配合记录地面运动加速度,频率范围为O～3OHz.记录系统优良的记录特性和高质量的数字记录,为实行地震矩张量反演方法求解震源机制提供了基本条件.卢龙地震余震观测中所记录的地震多为3级以下的中小地震,使用地震矩张量反演方法可以对每次地震求解,从而得到了沿不同断层面,不同位置的地震震源机制和应力分布.这是其它方法,如P波初动方法,SH/P振幅比的方法很难做到的.从本文讨论的卢龙Ms6.1地震的余震序列,我们将看到
7、不同位置之间的震源机制的差别.研究余震之间的震源机制的差别,对于较深入地了解一个地区的活动构造具有积极的作用.在进行本研究之前,我们已经用Hypoinverse程序对该地震进行了定位(李文军等,2005),结果见图1.根据该定位的结果和所用的模型参数,我们用矩张量反演的方法计算各个余震事件的震源机制.1矩张量反演方法弹性动力学方程解用广义Green函数可以表示为(Madariaga,1983)U(r,￡)一G,(r,t;0,O)*M(0,￡)(1)反演所用的Green函数G.(,,t;0,O),可用
8、水平层状地壳介质模型计算得到(Kennett,Kerry,1979;Kennett,1980;周云好,2002).将上式的时间和褶积离散化,通过求解线性方程组,可反演出地震矩张量(0,￡).对时间积分(0,￡),即可求出随时间变化的地震矩张量(0,￡)的6个独立分量.由于时间域中的计算耗时长,占用内存大,我们把它转化为频率域中的计算.在频率域中,用体波反演地震矩张量的问题可转化为求解线性方程组UN×l((cJ)一GN×6((cJ)?M6×l((cJ)(2)式中,为离散
