一种用于地震预警的三维实时定位方法.doc

一种用于地震预警的三维实时定位方法 韩雪君 马延路 刘瑞丰 中国北京 1000 中国地震 摘要 研究使用波前追踪算法，计算中国几个地区的三维地壳模型的走时网格，利用中国地震台网资料，对区域内发生的地震进行定位分析。结果显示，在一定的台网密度条件下，这种三维实时定位方法能在震后数秒给出震源位置，可满足地震预警的要求。 关键词 地震预警；实时定位；“着未着”定位方法；三维走时网格 0 引言 地震预警是指一个地震发生后，在具有强烈破坏性的S波到达之前发出警报信息（Kanamori，2005）。目前世界上多个国家或地区都已经建立了预警系统，比如日本的实时地震信息系统（REIS），能在第一个台站触发的34s内确定震源参数，震中定位误差小于20km（Nakamura et al，2009）；美国南加州的ElarmS预警系统在台间距小于20km的台网中采用单台开始定位，35台触发后定位结果稳定（Allen et al，2009）。我国也计划在“十二五”期间建成530s的预警系统。 地震发生后使用P波的信息迅速判定震中对增加预警时间非常重要。另外，预警时间也与震中距有关，对于震中附近区域，预警时间可能仅有数秒钟，但对于比较远的地区，可达到1以上。大地震的震源破裂时间通常比较长，其积累的位错能量往往并非在断层错动的初期就立刻辐射出来，利用初始破裂时辐射出的P波，来尽可能快地定位出震中，估计震级并发出预警，在目前地震预测尚未解决的情况下，能够最大限度地降低地震带来的人员伤亡和财产损失。因此，一个有效的地震预警系统首先必须能既快且准地确定震源位置。 地震定位是地震预警的第一步。常规的定位算法通常需要等到四、五个台站触发后才开始启动，在台网所含台站不够密集时往往无法满足预警的需求。为了提高定位速度，从第一个台站触发之后就开始定位的方法的研究和软件开发变得十分必要。Horiuchi等提出了“着未着”定位方法用来快速确定震源位置，这一方法的特点是同时利用了已触发台站到时和未触发台站位置信息。Satriano等通过使用已触发台与未触发台之间到时差的约束不等式，在定位早期已触发台站数不足的情况下，增加了对震中位置的约束，实现了从单台定位到多台定位的平稳过渡。 目前通过体波、面波层析成像和反射勘探剖面结果的研究和积累，已经得到地球结构的模型。随着计算机技术的不断发展，将这些三维模型用于地震定位也成为可实现的目标。本文以crust2.0模型（Bassin et al2000）为例，构造中国不同区域的三维模型，应用FMTOMO软件（Marthijn de Kool et al，2006；Rawlinson et al，2006）中的地震波前追踪程序计算了空间三维区域走时网格，并使用中国测震台网资料，用本文实现的定位方法，对若干地震震例进行了定位和分析，给出震源位置的实时概率空间分布。 1 方法和资料 1.1 实时定位方法概述 台网中有N个正常工作的台站：S1……SN，在tnow时刻S1、……、SK（KN），SK+1、SK+2……SN为未触发台站，对于已触发台站之间的到时差有以下约束： （1） 上式中，ttm(n)为震源到Sm(n)台站的走时， tm(n)为Sm(n)台站P波初至观测到时。这一等式要求震源到两个台站的走时差等于这两个台站初至的到时差。常规双曲线定位算法（Pujol，2004）就是根据该方程在三维空间里决定的K（K-1）/2个双曲面来限定震源的可能位置（见图1，图中三角形为已触发台站，三角形为未触发台站，为P波波前）。 图1 常规等时差（双曲线）算法确定震源区域 Figure 1 Earthquake location using normal equal differential time (EDT) formulation method 对于未触发台站Sl与已触发台站Sn间的走时差满足以下不等式，即： （2） 其中，ttl为震源到未触发台站Sl的走时，ttn为震源到已触发台站Sn的走时， tn为Sn台站初至观测到时。该不等式在三维空间里确定了一个以双曲面为界的区域。由每个已触发台站和未触发台站进行组合后所决定的K（N-K）个区域的交集就可以限定出震源位置（图2）。有了这一条件，使得我们在第一个台站触发后，即可开始定位。此时由于tnow=tn，不等式右边为0，震源的可能区域实际上就是第一个触发台站所在的Voronoi单元（如图3阴影区所示）。所谓的Voronoi单元（Voronoi，1908），是由一组由连接两相邻点直线的垂直平分线所组成的多边形构成。N个台站会将区域分割成为N个Voronoi单元，其中每个Voronoi单元中的点到该单元中台站的距离小于到任