根据高密度地震监测台网生成的实测烈度值，开发1套基于Arcpy的计算程序，利用ArcGIS平台提供的Geostatistical Analyst模块进行克里金插值，得到烈度等值线分布图。在程序功能实现基础上，选用震例进行检验。研究结果显示，对于研究区域内仪器烈度在Ⅱ度以上的地震，基于本程序绘制的等值线图能够直观反映震后烈度分布，可为震后快速评估震害影响提供参考依据。
Based on the measured intensity values generated by the high-density earthquake monitoring network, a set of Arcpy-based computational programs was developed. The Geostatistical Analyst module provided by the ArcGIS platform was used for kriging interpolation to obtain the intensity contour distribution map. After functional implementation of the program, selected earthquake cases were tested. The research results show that for earthquakes with instrument intensity above magnitudeⅡ in the study area, the contour map drawn based on this program could intuitively reflect the post-earthquake intensity distribution, providing a reference for rapid hazard assessment after the earthquake.
2023,44（6）: 42-47 收稿日期：2022-11-15
DOI：10.3969/j.issn.1003-3246.2023.06.007
基金项目：北京市地震局科技项目（项目编号：BJWC-2022004，BJWC-2023005）
作者简介：梁芳（1973—），女，理学硕士，副研究员，主要从事地震烈度速报与预警工作。E-mail：liangfang@bjseis.cn
*通讯作者：白立新(1969-),男,理学博士,副研究员,主要从事强震动观测研究工作。E-mail:bailixin@bjseis.cn
参考文献：
白立新,成云辉,张杰, 等.河北永清M 4.3地震北京烈度仪台网记录分析[J].震灾防御技术,2019,14(1):210-219.
胡杨.非均匀分布强震动台站的仪器地震烈度等震线生成方法[D].哈尔滨:中国地震局工程力学研究所,2017.
金星,张红才,韦永祥,等.基于地震监测台网资料近实时插值计算震动图的初步研究[J].防灾减灾学报,2010,26(1):1-11.
李海涛,邵泽东.空间插值分析算法综述[J].计算机系统应用,2019,28(7):1-8.
檀斌,张洁,马犇,等.基于Arcpy的地震烈度自动生成技术的研究[J].电脑编程技巧与维护,2022,(8):63-66.
王士成,金星,张红才,等.台网密度对地震烈度速报的影响研究[J].地震工程与工程振动,2017,37(6):162-168.
夏苏琼,李乃强.基于Python与Arcpy的电子地图自动化制图研究[J].测绘与空间地理信息,2022,44(9):221-224.
薛蕾,王遹其,付萍.基于Python的地震烈度速报软件的设计与应用[J].科技与创新,2019,(19):146-148.
张红才.基于地震监测台网资料的震动图及震动烈度研究[D].哈尔滨:中国地震局工程力学研究所,2008.
周崴,徐东炯,沈丽娟.基于Python和Arcpy的专题图批量制作方法[J].测绘与空间地理信息,2020,43(10):17-20.
埃里克·马瑟斯[美]. Python编程:从入门到实践[M]. 袁国忠,译.北京:人民邮电出版社,2020.
Pimpler E[美]. 基于ArcGIS的Python编程秘笈[M].牟乃夏,张灵先,张恒才,译.北京:人民邮电出版社,2017.
Zandbergen P A[美].面向ArcGIS的Python脚本编程[M]. 李明巨,刘昱君,陶旸,等,译.北京:人民邮电出版社,2014.
