基于Voronoi图地震预警定位算法.docVIP

基于Voronoi图地震预警定位算法 　　摘 要：利用一种基于Voronoi图的震中定位方法，对日本K-NET和KiK-net的地震记录进行定位。结果表明，该算法不受速度结构模型的影响，只需要台站的波到时顺序，所需台站数目少、计算速度快。这些特性和地震预警系统需要尽快给出预警信息的特性相契合，因此适合于地震预警系统的震中定位。 　　关键词：地震预警；震中定位；Voronoi图；到时顺序；地震预警系统 　　中图分类号：TP312 文献标识码：A 文章编号：1672-7800（2016）003-0048-04 　　作者简介：王安（1990-），男，山东临沂人，中国地震局地震研究所硕士研究生，研究方向为地震预警算法、机器学习、图像处理。 　　0 引言 　　地震给人类带来了巨大灾难，人们一直期望能实现地震预报，然而迄今为止，尚没有一种可靠的短期地震预报方法，地震预报是一个世界公认的难题。在这种情况下，为了减少地震带来的损失，地震预警作为一种折衷的手段具有重大意义。地震预警系统利用的原理是电磁波的传播速度比地震波快[1]，以及破坏性较小的P波速度比破坏性大的S波速度快，将距离震中近的台站报警信息，通过电磁波传送给地震波还没有达到的区域，可为这些区域在地震到达之前预留出应急操作时间。因此，如何尽可能快速地给出报警信息是地震预警系统的关键问题。 　　针对地震预警系统中的震中定位问题，传统的地震定位方法需要综合所有台站信息，所需信息多，计算过程复杂，无法满足地震预警系统对时效性的要求。截至目前，针对地震预警提出的震中定位方法主要有以下几类：①根据地震动幅值和震中距的关系进行定位[2]；②利用台站的到时顺序进行定位[3]；③B-Δ方法，利用参数B和震中距的对应关系，通过对P波初始阶段求包络，得到参数B与震中距，从而进行定位[4]；④利用两个台站的到时差得到震中到两个台站的距离差，即震中到两个台站的距离差为常数，因此可以确定一条双曲线，震中则应在这条双曲线上[5]。同时利用未到台和已到台的信息对震中的约束，得到定位结果[6]。国内一些研究对以上方法进行了结合[7]，金星等[8]提出了单台、多台连续定位的方法。 　　Voronoi图定位算法是Rosenberge（2009）[9]最早提出的。该算法的唯一前提是时间与距离成正比，在地震预警网络中，该前提是可以满足的。因此其不受速度模型的影响，具有定位速度快、所需台站少等特点，适合地震预警网络中的震中定位。本文改进并初步实现了该算法，并利用日本K-NET和KiK-net的地震记录进行了验证。 　　1 Voronoi图定位算法 　　1.1 Voronoi图 　　以地震台站为例，假设有台站集S。对于集合的每一个台站Sx，都有一个关联区域，其中所有点到台站Sx的距离小于到集合S中其它台站的距离，把该关联区域称作Voronoi域。所有台站的Voronoi域用边界加以区分，组合在一起则得到了完整的Voronoi图（见图1）。 　　2 Voronoi图算法定位实例与结果分析 　　2.1 震例筛选 　　根据地震预警原理，地震预警网络的用途更多的是检测网内发生的地震事件，然后将该信息在地震达到目标之前传送给目标。如果事件发生在网外较远距离，那么对目标产生的破坏也是有限的。因此本文中选择的都是网内事件，详细地震目录见表1。 　　2.2 定位结果 　　根据以上筛选规则，笔者从日本K-NET和KiK-net的记录中选择了40个震例。对每一个事件，运用该地震定位算法，从第一个台站触发开始，每一个台站触发后给出新的定位结果，并与实际震中比较，得出每一次定位误差。从图4中可以看出，随着触发台站数目的增加，定位精度也随之提高，从地震发生到得到结果所需时间也随之增加，定位精度和所需时间是相矛盾的。因此，对地震预警系统而言，必须在定位精度和所需时间之间作一个折衷的选择。不过，这里还是对每个事件的前8个触发台站进行了定位以便分析。 　　结果如图4所示，定位结果与实际震中的偏差随触发台站数目的增加有了显著提高，首台触发时，到实际震中的平均距离接近11km，当4个台站触发时，平均距离略大于3km，当8个台站触发时，平均距离略大于1km。 　　图5显示了40个事件最终定位结果的误差分布，可以看出除一个事件的定位误差大于5km外，绝大部分事件的定位误差小于4km。图6中每个点代表一个事件，其横坐标为该事件用首台作为震中时与实际震中的距离，纵坐标为第8个台站触发后，用Voronoi图定位算法对事件进行定位得到的震中与实际震中的距离。对于一个事件，如果用Voronoi图定位算法得到的结果比用首台作为震中得到的结果好，即到震中的距离更近，则横坐标x大于纵坐标y，该点应在直线y=x下方。把40个事件和直线y=x绘制在