Design and Implementation of High-performance Computing Platform Accounting Service System in Key Laboratory of Earthquake Prediction
-
摘要: 本文简要回顾了高性能计算在地震研究中的应用,介绍了中国地震局地震预测重点实验室高性能计算平台及其应用情况,并指出作业计费的必要性。其次,详细介绍了1种高性能计算平台计费服务系统的设计与实现。该系统基于TORQUE和LSF集群作业管理系统日志信息的分析处理,具备作业查询、作业统计、作业计费、作业预约、信息管理、留言管理等功能,有益于解决高性能计算平台计算收费和运维支出平衡的问题,提高了平台的透明度和利用率,促进信息共享和技术交流,并具有较好的扩展性和应用前景。最后,对高性能计算过渡到高性能云计算进行了展望。Abstract: This paper briefly reviews the application of high-performance computing in seismic research, introduces the high-performance computing platform and its application in the key laboratory of earthquake prediction of China Earthquake Administrator, and points out the necessity of cluster accounting. The design and implementation of a high-performance computing platform accounting service system is presented in detail. Based on job log analysis and processing of TORQUE and LSF cluster management system, the system provides with such functions as job query, job statistics, job accounting, job order function and online information sharing and exchange. The system is beneficial to solve the problem of calculating fee and platform cost balance on high-performance computing platform, improves the transparency and utilization of platform, promotes information sharing and technology exchange, and has better expansibility and application prospect. Finally, the development prospect of high-performance computing transition to high-performance cloud computing is briefly described.
-
Key words:
- Seismic research /
- High-performance computing /
- Job log /
- Accounting system
-
引言
地震预警是利用电磁波传播速度远大于地震动传播速度以及破坏性地震波(S波或面波)的传播速度小于破坏性较小的P波速度的原理,在地震发生后尽可能短的时间内,确定地震基本参数或直接估计预警目标区的地震动和破坏程度,在破坏性地震波到来前,对可能破坏区发出地震警报。地震预警系统可以快速地进行地震定位并根据震中距不同,提供几秒、十几秒、甚至几十秒的时间,将预警信息在破坏性地震波到达之前发送至用户,以帮助用户采取应急措施,达到减灾的目的(马强,2008;殷海涛等,2012;张红才等,2013)。目前,日本、墨西哥、罗马和中国台湾地区已经建立了地震预警系统并提供预警服务,美国、瑞士、意大利和中国大陆则已建立了试验性质的地震预警系统(Allen等,2009;赵兵,2011)。
天津是中国北方最大的港口城市和环渤海经济中心,是国家防震减灾的重点设防城市,具有发生强烈地震的背景和条件,历史上曾多次遭遇特别重大的地震灾害。1679年9月2日三河-平谷8级地震和1976年7月28日唐山7.8级地震都使天津市遭受了严重的破坏。其中,唐山7.8级地震使宁河和汉沽的地震烈度达到Ⅸ度,塘沽及市区大部分地区地震烈度达到Ⅷ度,从而使天津成为中国唯一遭受Ⅷ度地震灾害的特大型城市。因此,建设天津简易烈度计地震预警试验区,能够建立完善有效的地震预警机制,进一步提高天津地区的地震灾害防御能力,尽可能地避免未来强震造成的人员伤亡和经济损失。
1. 试验区建设概况
天津简易烈度计地震预警试验区的建设目标是在天津地区新建80个简易烈度计观测点,与天津行政区内具备实时传输能力的测震台站和强震动台站共同组建天津地震预警观测网络,其网络规模如图 1所示。观测网络东西孔径约127km,南北孔径约160km,平均台站间距约10km。结合天津地震观测台网现有布局特征,选择观测条件较好的中小学校内平房或二层楼房的一层作为简易烈度仪观测场址。仪器固定在观测室不高于30cm的承重墙上,采用锚固方式安装,如图 2所示。数据传输采用移动4G无线VPDN方式,VPDN为虚拟专用拨号网络,是建立在VPN基础上且隶属与于VPN的1种业务。根据天津市远程网络传输需求,烈度计端路由设备首先通过4G无线网络接入附近移动运营商的无线网,再通过移动运营商VPDN专用网接入地震行业网,实现节点端与地震系统内网的互联互通。移动4G无线VPDN组网流程如图 3所示。
2. 简易烈度计性能
综合考虑仪器的安装方式、性能指标以及京津冀地震烈度速报与预警示范区的整体规划,天津简易烈度计地震预警试验区的台站观测仪器选用GL-P2B型和MI3000型三分向型简易烈度计,均具备16位以上分辨率,水平向测量范围-2—2g,垂直向测量范围-1—3g,动态范围大于80dB。传感器采用微机械加速度传感器(MEMS),MEMS加速度传感器为微电子技术与微机械工程结合发展的1种惯性器件,与传统力加速度传感器相比,具有体积小、成本低等特点。将MEMS传感器应用于地震仪器中,可以较低的成本提高台网密度,并可实现烈度速报与地震预警功能(Chung等,2011;Wu,2013;D'Alessandro等,2014)。2种仪器的性能指标见表 1,仪器的幅频特性曲线见图 4。
表 1 简易烈度计性能指标Table 1. Performance index of simple intensity meter型号 数量 采集延迟 分辨率 动态范围 加速度平坦频率 线性误差 GL-P2B 40 <1s 28位 >80dB 0.1—20Hz <1% MI3000 40 <1s 16位 >80dB 0.1—20Hz <1% 3. 地震预警数据处理系统
首先,简易烈度计台站、强震台站和测震台站的数据实时汇集到地震数据汇集系统,再通过行业内网将数据实时传输至地震预警数据处理系统,数据传输流程如图 5所示。当地震发生时,地震预警数据处理系统利用多个台站收到的地震初至波信号,使用STA/LTA算法检查震相触发,使用AIC算法精确拾取震相到时,使用“着未着”算法快速分析多个台站震相到时数据,确定地震发生的时刻和位置,使用P波最大振幅和台站震中距估算地震震级,从而获得地震基本参数,并向地震波尚未到达的地区发布准确的地震预警信号。
2018年8月16日13时50分20.97秒,天津市蓟州区发生2.2级地震,距离震中较近的简易烈度计台站、强震实时台站和测震台站均监测到了P波信号。地震预警数据处理系统采用256字节的数据打包模式进行数据传输,地震发生后8s发出第一报,发震时刻为2018年8月16日13时50分20.7秒,震中地点为天津蓟州区,震级为2.0级(转换震级),震中地点距天津市中心约94km,横波预计到达时间为12s,触发台站39个(图 6)。由于此次地震的震级较小,尚需经历更大的地震对天津简易烈度计地震预警试验区的效能进行检验。
4. 简易烈度计数据记录
简易烈度计观测的数据存储在地震数据处理中心服务器上,采用定长512字节MiniSeed格式,对于任一时段的数据,可利用JOPENS6.0数据处理系统转换为SEED、ASCII等多种数据文件。在数据处理过程中,采用python语言编写相关的计算程序,主要使用了python中的Numpy、Matplotlib、Basemap、Obspy等扩展包,有效提高了系统的开发效率和数据处理的工作效率。
在2018年8月16日天津蓟州2.2级地震中,位于天津蓟州附近的简易烈度计台站A0055、A0056、A0058、A0063均获取了观测数据。根据《仪器地震烈度计算(征求意见稿)》1的要求,对0.1—10Hz频率范围内的观测数据进行带通滤波处理,然后三分向合成计算最大峰值加速度。计算得到距离震中4km的A0058台站的最大峰值加速度为7.4gal;距离震中10km的A0055台站的最大峰值加速度为6.46gal;距离震中10km的A0063台站的最大峰值加速度为6.78gal;距离震中15km的A0056台站的最大峰值加速度为6.27gal。
1 中国地震局,2015.仪器地震烈度计算暂行规程.
由于简易烈度计台站没有同址的测震或强震台站,故选择台站距离相邻、方位相差较小的强震台站BAJ与简易烈度计台站A0056进行波形数据对比。BAJ台距离震中12km,三分向合成的最大峰值加速度为6.89gal,使用的地震计为SLJ-100型加速度计,测量的频带范围为0.1—30Hz,采样率为100Hz。简易烈度计台站A0056使用GL-P2B型MEMES加速度计,测量频带范围为0.1—30Hz,采样率为100Hz。观测台站分布如图 7所示。
由于观测频带不同,采用4阶巴特沃斯(Buttorworth)带通滤波器对强震台BAJ与简易烈度计台A0056的0.1—20Hz观测数据进行滤波,结果如图 8所示。强震台BAJ三分向记录的波形峰值加速度(PGA)分别为3.99gal、3.88gal和1.43gal,简易烈度计台A0056三分向记录的波形峰值加速度(PGA)分别为1.23gal、1.28gal和0.78gal。2个台站的峰值加速度和波形数据时程均不相同,这是由于台站位置不同造成的,其中强震台BAJ距离震中12km,简易烈度计台A0056距离震中15km,二者直线距离相差3km,且相对震中的方位有差异,但二者NS向与EW向波形存在一定的相似性。
观测仪器的记录噪声水平在一定程度上能够直接反应仪器的记录能力,尤其是记录中小地震事件的能力。为此,截取了简易烈度计台站A0056无明显干扰的10分钟的数据记录,利用Welch算法计算加速度噪声功率谱(Peterson,1993),结果如图 9所示。由图可以看出,简易烈度计台A0056的加速度噪声功率谱处于高噪声模型(New High Noise Model)和低噪声模型(New Low Noise Model)之上,且远高于高噪声模型。由此可知,相对于测震台站对地脉动噪声的完整记录(在NHNM和NLNM之间),简易烈度计台站记录的噪声数据完全是仪器的自噪声,天然噪声则淹没在仪器的自噪声中无法分辨。由于简易烈度计的仪器自噪声已经超出地脉动水平,无法观测记录到有效地脉动,且仪器通常布设在人员比较密集的建筑物内,因此,简易烈度计台站的记录数据信噪比较低,地震数据初至震相会淹没在背景噪声中。
各类地震观测仪器的观测特点不同,测震台站主要服务于地球物理学的研究,旨在了解地球内部的构造和地震活动性,为此提供各种地震参数,包括地震基本参数、地震活动性、地震波传播、地震趋势等,地震观测的场地条件要求为基岩、深井且背景噪声应尽量小。强震观测主要为地震工程学研究和结构抗震提供资料,旨在测定地面和建筑物在地震作用下的运动过程,包括地震影响场、地面及结构的地震动特征,主要针对的是近场强震,仪器频带范围为-2—2g,通常要求背景噪声小于1cm/s2,城区或乡村均可布设观测仪器。而简易烈度计是以MEMS传感器为主的非专业强震仪,主要针对的是近场强震,要求背景噪声小于2cm/s2,记录质量较低,一般布设在自由地表或2层以下低矮结构建筑物内1。简易烈度计具有成本低、安装要求不高、维护简单等特点,适于高密度布设,以提高区域地震预警和烈度速报能力。
1 马强,2017.地震预警与烈度速报系统的试验示范.浙江省地震局内部报告.
5. 讨论与结论
本文对天津简易烈度计地震预警试验区的建设、台站分布情况、仪器性能指标、数据处理及数据记录情况进行了介绍,探讨了简易烈度计在预警台网中的应用,获得了以下初步认识:
(1)通过建设天津简易烈度计地震预警试验区,加密了天津强震台网的密度,初步实现了天津及邻区的地震预警功能。
(2)由测震台站、强震动实时传输台站及简易烈度计台站组成的天津地震预警观测网络,拓宽了各台网的数据应用范围,不仅可应用于地震应急快速响应,还可应用于地震工程和科学研究。
(3)简易烈度计台站造价低、建设周期短,易于大面积布设或规模化应用,可用较少的投资实现更好的应用(姚会琴等,2012),是现有测震台站和强震台站的有益补充。
(4)简易烈度计记录的噪声数据完全是仪器的自噪声,高于高噪声模型,但对于震中附近满足一定信噪比条件的台站,可获得较准确的初至震相到时信息,对于地震预警信息的准确产出具有意义。
(5)今后还需进一步积累地震数据,结合天津及邻区的地震分布特点优化地震预警参数模型,以期提升该地区的防震减灾能力以及公共服务水平。
-
表 1 计算资源价格及权重赋值
Table 1. The price and weight value of computing resoureces
资源指标 价格/元 权重/% CPU时 0.3 70 物理内存/GB 1.0 20 虚拟内存/GB 1.0 10 -
陈振娅, 杨英, 陈丽娟等, 2009.基于OpenPBS的高性能计算平台计费系统.山东科学, 22(5):98-102. http://d.old.wanfangdata.com.cn/Periodical/sdkx200905022 郭婷婷, 徐锡伟, 2011.有限元法在构造应力场与地震预测研究中的应用与发展.地震研究, 34(2):246-253. doi: 10.3969/j.issn.1000-0666.2011.02.023 何戈, 董小社, 钱德沛, 2003.网络计算环境下一种计费方法的研究与实现.计算机工程与应用, 39(17):185-187. doi: 10.3321/j.issn:1002-8331.2003.17.060 李惠欢, 杨敏, 吴汝明, 2016.基于TORQUE的高性能计算平台记账系统.计算机应用与软件, 33(8):126-130, 205. doi: 10.3969/j.issn.1000-386x.2016.08.028 李闽峰, 李圣强, 王斌等, 2009.地震预测研究所高性能集群计算系统.见: 中国地球物理2009.合肥: 中国地球物理学会. 李圣强, 李闽峰, 刘桂平等, 2012.高性能集群计算系统的构建.地震, 32(1):144-149. doi: 10.3969/j.issn.1000-3274.2012.01.014 廖凯宁, 郝永伟, 2015.地震数据处理和地震动力学仿真平台的优化研究.计算机工程与科学, 37(4):663-669. doi: 10.3969/j.issn.1007-130X.2015.04.007 林皎, 张武生, 徐伟平等, 2012.高性能计算平台预付费管理系统的研制与应用.实验技术与管理, 29(12):104-105, 128. doi: 10.3969/j.issn.1002-4956.2012.12.030 刘明良, 2008.基于校园网格系统计费方案研究.计算机系统应用, (3):60-63, 91. http://d.old.wanfangdata.com.cn/Periodical/jsjxtyy200803016 龙桂华, 赵宇波, 李小凡等, 2011.三维交错网格有限差分地震波模拟的GPU集群实现.地球物理学进展, 26(6):1938-1949. doi: 10.3969/j.issn.1004-2903.2011.06.007 陆新征, 卢啸, 许镇等, 2014.基于高性能计算的工程抗震与防灾: 从单体到城市.见: 徐正安, 任伟新, 丁克伟主编, 第六届结构工程新进展国际论坛文集.北京: 中国建筑工业出版社. 陆新征, 陈磊, 曾翔等, 2015.基于云计算的工程结构地震响应数值模拟.沈阳建筑大学学报(自然科学版), 31(5):769-777. http://www.cnki.com.cn/Article/CJFDTOTAL-SYJZ201505002.htm 吕弋培, 徐锐, 田文洪等, 2009.基于中国国家高性能计算与网格系统的抗震减灾应用平台.中国科技纵横, (12):472-476. doi: 10.3969/j.issn.1671-2064.2009.12.347 马新生, 谢志南, 王延伟等, 2010.高性能集群系统在地震波数值模拟中的应用.地震工程与工程振动, 30(6):43-48. http://d.old.wanfangdata.com.cn/Periodical/dzgcygczd201006007 牛铁, 朱鹏, 赵毅等, 2010.超级计算环境计算配额系统设计和实现.计算机应用, 30(S2):8-9, 39. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201003007085 沈瑜, 李娟, 常飚等, 2014.高性能计算机统一资源管理系统的设计与实现.计算技术与自动化, 33(1):83-90. http://d.old.wanfangdata.com.cn/Periodical/jsjsyzdh201401018 盛乐标, 周庆林, 2014.利用LSF API实现高性能计算机集群的机时统计.实验科学与技术, 12(2):50-52, 117. doi: 10.3969/j.issn.1672-4550.2014.02.017 王洪涛, 谢礼立, 2009.考虑楼板作用的钢筋混凝土框架有限元模型及并行计算效率.地震工程与工程振动, 29(1):63-69. http://d.old.wanfangdata.com.cn/Periodical/dzgcygczd200901008 王伟, 杨利, 刘建锋, 2008.一种高性能计算环境中的计费系统.计算机工程与科学, 30(1):148-150. doi: 10.3969/j.issn.1007-130X.2008.01.042 占昌宝, 罗川, 丁振坤等, 2016.高层建筑抗震性能预测仿真研究.计算机仿真, 33(8):397-402. doi: 10.3969/j.issn.1006-9348.2016.08.086 张军华, 臧胜涛, 单联瑜等, 2010.高性能计算的发展现状及趋势.石油地球物理勘探, 45(6):918-925. http://d.old.wanfangdata.com.cn/Periodical/zggckx200001016 赵毅, 朱鹏, 迟学斌等, 2007.浅析高性能计算应用的需求与发展.计算机研究与发展, 44(10):1640-1646. http://d.old.wanfangdata.com.cn/Periodical/jsjyjyfz200710002 -