实时地震监测台网中Voronoi图生成算法的改进.pdfVIP

第37卷 第2期 地 震 研 究 VoL37．No．2 2014年4月 JOURNALOFSEISM0L0GICAL RESEARCH Apr．，2014 实时地震监测台网中Voronoi图生成算法的改进 水 赵 瑞，刘希强 (山东省地震局 ，山东 济南 250014) 摘要：地震预警是近年来发展起来的地震监测新技术和震灾防御新手段，实时地震监测台网中Voronoi图生成算 法研究可为地震早期预警系统提供有效的技术支持。在对Voronoi图逐点插入算法进行分析研究的基础上，结合 实时地震监测台网运行实际，提出了改进后的局部三角网格生长法 ，能够快速有效生成除中断台站之外正常运行 台站的Voronoi单元，同时有效约束凸壳的生成。与传统方法相比，该方法生成的Voronoi图不失精度 ，但运行时 间更短 。 关键词：实时地震监测；Voronoi图；Delaunay三角形 ；逐点插入算法 ；局部三角网格生长法 中图分类号：P315．6 文献标识码：A 文章编号：1000—0666(2014)02—0244—07 角网格互为对偶图，可通过Delaunay三角网格寻找 0 引言 与每个离散点相关的三角形，顺序连接每个三角形 的外接 圆心得到多边形，即为 V图 (孙继忠等， 建设地震监测台网的主要科学 目的之一就是一 2010)。经过众多专家学者的研究，Delaunay三角网 旦地震发生后，第一时间开展地震速报工作。当前 的生成算法 日趋成熟。目前主要有分割一归并法、 基于实时地震监测台网建立地震预警系统和地震应 逐点插入法和三角网格生长法及每种方法的改进算 急控制系统在 国际上 已形成潮流。日本、美 国、墨 法，其 中前两种 方法 应用 较 多 (Aurenhammer， 西哥、土耳其、意大利、瑞典、罗马尼亚、中国台 1991；Delaunay，1934；Okabeetal，1994)。每种算 湾等地震多发国家和地区发展了多套预警系统，如 法各有其优缺点，分割一归并算法适合离散点较多 日本铁路UrEDAS系统、墨西哥SAS和 SASO系统、 的情况，但该算法含有大量迭代运算，空间复杂度 美国ElarmS和PreSEIS系统、土耳其 PreSEIS和SO． 高，并且凸壳形成时，算法较复杂，时间复杂度也 SEWIN系统等，并有一些成功预警且收到减灾实效 较高 (周培德，2000)。逐点插入法算法简单，占用 的先例 (Gasparinietal，2007；Wuetal，2007；A1一 内存少，空间和时间复杂度不高，但没有很好地解 lenetal，2009；殷海涛等，2012)。在地震预警技 决凸壳问题。三角网格生长法算法简单，以基边为 术中，Voronoi图 (以下简称V图)在地震事件触 循环条件，产生凸壳时较困难，同时易产生边与边 发判断、地震定位等方面 已得到初步应用 (Ry— 相交，不利于 Delaunay三角网格优化。武晓波等 delek，Pujol，2004；Hofiuchietal，2005；Satriano (1999，2000)将逐点插入法与分割一归并法相结 etal，2008；王庆民等，2012)；同时在地理信息系 合，提高了计算速度，降低了算法复杂度，但同样 统、计算机图形学、气象学、虚拟实现、机械工程、 没有很好解决凸壳、初始三角网的生成问题。王强 机器人等领域也有着广泛的应用 (刘金义，刘爽， 等 (2010)在此合成算法的基础上，提出一种改进 2004；范会川等，2010)。 的Delaunay三角网生成算法，采用虚拟 四边形凸