网内地震震中近实时动态定位方法的研究.pdfVIP

中国地球物理2011 ·483· ●l专题十五：信息技术与地球物理l◆ (15)信息技术与地球物理 网内地震震中近实时动态定位方法研究 NearReal-Time Location Methods EarthquakeEpicenter Study 刘希强’孙亚强赵冰王庆民赵大鹏赵瑞李铂 Liu Sun Zhao eta1． XiqiangYaqiangBingWangQingrnin 山东省地震局 济南250014 地震震中实时定位是“地震早期预警”(EEw)系统的基本要求。近年来EEW技术在国际上得到迅速 的发展(Alienet a1．，2009)，但我国的地震预警系统建设刚刚起步。在如何利用最少的台站直达波到时、最 大限度地提高地震定位的时效性和准确性问题上，科学家已经开展了一些研究工作(Satrianoetal。2008)。 基于此，本文综合应用触发台站震相和位置以及未触发台站位置信息，借助Voronoi剖分技术和地震波传 播规律探索依次应用1个台站到3个台站实时动态确定地震发生的区域、线段区间和震中的综合分析新方 法，并对地震定位的时效性和精度进行分析，以期推动我国地震预警技术向前发展。 方法：首先将测震台站两两相连形成Delaunay三角网，同时作连线的中垂线，中垂线相交形成若干个凸 多边形，即Voronoi图(称为矿图)。基于此，根据第f个台站触发记录可确定地震震中发生区；其次忽略地震 震源深度，可近似认为震中分布满足双曲线规律。根据台网中其它台站的初至波不是最先到达的事实，可将双 曲线上震中位置进一步约束在第一个触发台站所在的形区中一段双曲线上；最后根据三个台站记录的P波到时， 可得到3个由两两台站组合得到的地震震中双曲线分布，三条双曲线的交点围成的区域中心位置即是震中位置。 上述方法仅在理想的条件下适用，必须满足地震发生在网内、地震波传播介质性质均匀和距离实际地 震最近的台站首先触发这三个条件。但在实际地震实时动态定位过程中有许多问题具有不确定性，表现在 以下几个方面：一是地震观测台站偶有中断，如距离地震最近的台站在震时恰逢中断；二是矿区内或不同 的y区P波速度并不一致；三是地震震中的双曲线分布有时会出现三个以上交点；四是由于震相到时读取 的误差有时甚至不能形成双曲线；五是台站分布不合理，如三个台站呈近似线性分布等。 对此，本文采用以下改进思想进行地震震中测定：一是当双曲线间不相交时，则第一个台站所在的矿图内 双曲线间最近两点的中点为最大可能的震中位置；二是当双曲线相交只有一个交点，且距离该交点最近的三个 台中含有至少一个台站参与本研究震中定位时，则最大可能的震中位置为交点位置：三是当双曲线相交有一个 交点，但距离该交点最近的三个台中不含有一个台站参与本研究震中定位时，则最大可能的震中位置为第一个 台站所在的V图内双曲线间最近两点的中点；四是当双曲线相交有二个交点时，如果距离这两个交点最近的三 个台都一样且都参与本研究震中定位，那么最大可能的震中位置为距离第一个台站最近的交点：如果距离其中 一个交点(比如A)最近的三个台中含有至少一个台站参与本研究震中定位，而另一个交点不含有一个台站参 与本研究震中定位时，那么取A点为最大可能的震中位置：如果距离任一交点最近的三个台中都不含有一个 台站参与本研究震中定位时，那么最大可能的震中位置为第一个台站所在的V图内双曲线问最近两点的中点。 研究结果表明，对发生在山东虚拟测震台网网内的地震，基于有限台站和有限震相到时信息在震后平 均4秒钟左右即可确定地震发生的最小区域，在震后平均lO秒钟左右(不考虑数据传输延迟和数据处理 时间)可以给出比较准确的震中位置。425次定位结果的平均误差为6．3千米，误差小于20千米的地震频 次占总数的95．1％。因此定位方法可同时满足预警地震速报时效性和精度的双重要求。台站分布、速度模 型、到时读取精度等是影响震中测定精度的因素。为进一步提