澜沧―耿马地区中小地震精确定位的研究.docVIP

澜沧―耿马地区中小地震精确定位的研究 　　摘要：概述了双差定位方法及其基本原理，并运用该方法对1979～2009年澜沧―耿马地区发生的4 279个中小地震进行精确定位，获得了4 124个地震基本参数，走时残差均方根的平均值由定位前的0.502 s降到了0.086 s。定位结果显示震中沿北北西向的构造呈明显的条带状分布，且更加集中于断层附近，在深度上改变了定位前的层状分布现象，震源深度的优势分布为3～15 km，表明该区的发震层较浅。虽然地表显示有部分段落不连续，但从深度剖面和震源机制等结果来看，该北北西向的断裂在地下可能是相通的。 　　 　　关键词：双差定位；中小地震；澜沧、耿马地震带；北北西向构造 　　 　　中图分类号：P315.6文献标识码：A文章编号：1000-0666（2011）04-0482-06 　　 　　 　　0引言 　　 　　1988年11月6日21时03分14.5秒和15分44.9秒，在云南省澜沧县境内（22°50′N，99°43′E）和耿马县境内（23°23′N，99°36′E）分别发生了7.6、7.2级强烈地震，两次地震的震中均位于北西向的耿马―澜沧―勐海强震带内。该区域地震活动密集，不同方向的断裂相互切割，活动断裂断续相连，且发生的5级以上地震在该强震带上呈条带状分布（姜葵，1993）。为了研究澜沧―耿马地震的发震构造及其孕震背景，进行地震活动性分析，需要了解震源之间精确的空间偏移情况。由于定位精度受台网分布形状、震相到时读数的精度和地下速度结构等因素的影响，使得常规定位方法典型的不确定性比事件本身的震源尺度大很多倍，因此限制了地震活动精细结构的研究（Waldhauser，Ellsworth，2000），为此需要用定位精度较高的定位方法对该地区发生的中小地震进行重新定位。 　　 　　双差定位法是一种相对定位方法，由于使用的是两两地震间的到时差，因此可以有效提高地震相对位置的精度。Waldhauser和Ellsworth（2000）利用该方法对1984～1998年发生在美国加州北海沃德断层上的地震重新进行定位，平面定位精度可达到几十米，深度可达到几百米。我国的许多地震工作者已将该方法广泛应用于国内的具体工作，并取得了一定的进展（杨智娴等，1999，2002，2003）。双差定位法基于：当两个地震震源之间的距离足够小时，震源区到同一个台站间的整个射线路径几乎相同，这两个地震事件的走时差就来自于事件之间高精度的空间偏移，这样就可以部分地消除速度模型的不确定性所带来的误差。因此双差定位法特别适合于地震发生密度高的地区，并能较好地显示地震发生区断层的特征。 　　 　　本文尝试运用双差地震定位法对1979年至2009年12月澜沧―耿马地区发生的4 279个中小地震进行重新定位，期望通过定位结果获得对发震构造的更多的认识。 　　 　　1方法简介 　　 　　双差地震定位方法是由Waldhauser和Ellsworth（2000）提出的一种相对定位方法。该法反演的是一组丛集地震中每个地震相对于丛集矩心的相对位置，有效地减少了由于对地壳结构了解不够精细而引起的误差，该法还能够应用在空间跨度比较大的地震事件群体，只要两个相邻事件的距离远小于震中距及速度不均匀体的线性尺度即可。 　　 　　在反演过程中，双差定位法基于传统的盖格法（Geiger，1912）的迭代求解技术，即通过对未知震源位置地震波走时方程的观测走时与理论计算走时差作泰勒展开，保留一阶项迭代求解，使得两者残差最小时获取震源位置和发震时刻。 　　 　　定义drijk是两个事件之间观测和理论计算走时差的残差： 　　drijk=(tik-tjk)obs-(tik-tjk)cal.（1） 　　式（1）表示第i个地震事件到第k个地震台站的地震波的走时tik与第j个事件到第k个地震台站的地震波走时tjk之差的观测值(tik-tjk)obs与理论计算值(tik-tjk)cal的残差。 　　 　　对地震到时t求震源参数的偏导数，联立双差残差就可以确定两个地震i和j的相对位置 　　??tik??mΔmi-??tjk??mΔmj=drijk.（2） 　　式中，Δmi=(Δxi,Δyi,Δzi,Δti)T是第i个地震震源参数：(xi,yi,zi,ti)T的改变量；Δmj是第j个地震震源参数的改变量。其中，x、y、z是震源东西、南北和垂直方向的坐标；t是发震时刻。结合一个台站的所有地震对建立式（2），再将所有台站的方程写成如下形式的方程组 　　WGm=Wd.（3） 　　式中，G是一个M×4N偏微分商矩阵；M是双差观测的数目；N是地震数；d是双差资料矢量；W是用以对方程加权的对角线矩阵，m为由待定的震源参数改变量构成的维数为4N