基于 Python 和天地图的地震应急快速反应系统研究?以江苏地区为例

余其鹏　董建平

摘要为提升地震应急能力，方便震后快速获取地震影响场、地震应急专题图件和快速评估报告，为地震应急决策和救援提供专业技术支持。通过 QT 开发可视化窗体，利用天地图 API 调用天地图作为底图，天地图逆地理服务 API 获取基础地理信息，最后叠加地震专题数据。用户输入震中经纬度坐标等信息后，Arcpy库自动批量创建震中地区地震专题图件，并利用 Python-docx 库自动产出快速评估报告。实际应用结果表明，Python 脚本封装后系统程序只有一个单文件，即拷即用，部署简单，专题图批量生产、自动化产出震情简报等，应急效率显著提升，地震应急信息化服务水平显著提高。

关键词地震；地震专题图；地震应急；Python

中图分类号： P315.9文献标识码： A文章编号：2096-7780（2023）07-0296-07

doi：10.19987/j.dzkxjz.2022-139

Research on earthquake emergency rapid response system based on Python and Tianditu： A case study of Jiangsu

Yu Qipeng，Dong Jianping

（Zhenjiang Seismic Station， Jiangsu Zhenjiang 212000， China）

AbstractIn order to enhance earthquake emergency response capability，and to conveniently and quickly obtain seismic influence field，earthquake emergency thematic maps，and to rapidly assess post-earthquake reports，it is able to provide professional technical support for earthquake emergency decision-making and rescue. We develop visualization form through QT，use the Tianditu map API to call Tianditu map as the base map，and use Tianditu map geocoding API to obtain the basic geographic information，and finally overlay the seismic thematic data. After users input information such as epicenter latitude and longitude coordinates，arcpy library automatically creates a batch of earthquake thematic maps for the epicenter area，and a rapid assessment report is automatically output by using the Python-docx library. The actual application results show that the system program has only a single file after the Python script encapsulated，whichis ready to be copied and used，making it easy to deploy. It has features such as mass production of thematic maps and automatedoutputof earthquake briefings. Thesefunctions haveled toasignificant increaseinemergency response efficiency. The level of seismic emergency information service has been significantly advanced.

Keywordsearthquake; earthquake thematic map; earthquake emergency; Python

引言

在我國，地震灾害被列为群灾之首，地震灾害死亡人数居各类自然灾害死亡人数之首，比率高达54%，据统计我国的地震死亡人数约占全球地震死亡人数的一半[1]。破坏性地震发生后，快速产出的震中地区地震专题图以及震情简报能够直观的反应出震中情况，为地震应急决策和救援提供第一手的参考资料和强有力支持。所以地震相关专题图的制作以及震情简报都是各级地震工作部门在震后必须要面临的一项重要的工作任务，也有着较高要求，因此高效准确的制作地震影响场、震中地区历史地震和构造图等专题图件，以及快速形成一份专业震情简报以及后续的震情详细情况报告具有非常重要的现实意义。

从国内外研究来看，关于地震应急专题图以及震情评估报告的制作，许多学者通过手动、自动、本地网络等多种方法进行过探索。马宝君等[2]利用 VB6.0 开发系统，建立独立的数据库和独立的地震应急图件产出平台获取专业图件并批量输出。樊琪和刘洪光[3]通过建立预存储的地图切片数据拼接出需要的应急专题图地理底图，以达到快速批量出图的目的。韩贞辉等[4]利用 Python 语言编写脚本程序调用预制的地震应急指挥专题图模板进行出图。程显洲等[5]利用 ArcGIS 模型构建器快速调取数据和处理数据。使用Arcpy自动出图脚本进行自动批量制图，从而加快不同要求下地震应急专题图的制图速度，提高专题图的出图效率。利用天地图作为基础数据地图是研究者们常用的方法，孙哲等[6]通过研究“天地图—陕西”的切片方式和服务方式，在 ArcGIS 中实现了天地图瓦片数据与矢量数据相结合并基于 Python 语言和ArcPy包的出图技术。Zhou 等[7]基于ArcPy进行渔业海洋图集自动化制作，有效提高了绘制效率。Buckley和 Watkins[8]基于 GIS 数据库驱动制图，总结了具体的自动化地图制作工作流程。

针对前述研究，地震应急快速反应系统的关键问题有：一是基础数据的时效性，即能够及时调整更新系统内的基础数据，使其始终保持最新的数据；二是个性化定制性，即针对不同的需求，可以快速方便的调整产品的产出。本文利用 Python 语言，通过 QT 开发可视化窗体，结合天地图 API 调用天地图作为底图，天地图逆地理服务 API 获取基础地理信息，保证了基础数据的时效性，并叠加本地地震专题数据。用户输入震中经纬度坐标等信息后，Arcpy库自动批量创建震中地区地震专题图件，并利用 Python- docx 库自动产出快速评估报告。本文通过天地图 API 获取基础地理信息，摆脱保密数据限制，同时将 Python 脚本封装后系统程序只有一个单文件，无需安装，即拷即用，具有部署简单，专题图批量生产、自动化产出震情简报等优点，让软件产出更加丰富和可调整，让应急效率显著提升，地震应急信息化服务水平得到显著提高。

1 研究区概况与数据

1.1 研究区概况

江苏位于长江三角洲地区，总面积10.72万 km2，下辖13个地级行政区。江苏及其周边地区位于华北、华南两大地震区的交接部位， 多年来一直被列为我国地震重点监视区[9]。江苏地区在大地构造上跨越中朝断块、扬子断块和秦岭—大别褶皱带3大一级单元构造[10]。该地区内有郯庐断裂、茅山断裂带等多条断裂带，2012—2022年10年来，江苏地区大约发生 ML ≥2地震106次，其分布如图1所示。

1.2 数据来源与处理

本文所用到的数据主要是基础地理信息数据和地震专题数据，基础地理数据是利用天地图 API 调用天地图作为底图，天地图逆地理服务 API 获取基础地理信息，如行政区、交通、POI（兴趣点）等。另外一个是地震专题数据，包括历史地震和断层数据。本研究的地震目录数据和断层数据均源于中国地震台网中心国家地震科学数据中心（http：//data.earthquake.cn），地震目录数据记录了研究时间段内所有 ML ≥2.0的地震，数据列出了发震时间、经纬度、震源深度、震级、参考位置等详细记录，为xls存储格式，断层数据为1/400万中国活动构造空间数据库，存储格式为shapfile矢量格式。为了方便研究和存储，研究建立了 Geodatabase 空间数据库[1 ， 11]，将地震目录数据转换为坐标点矢量数据，然后建立数据集并导入空间数据库，同时将断层数据建立数据集也导入空间数据库。

2 研究框架

本文利用 Python 编程语言，使用天地图 API 调用天地图作为底图，天地图逆地理服务 API 获取基础地理信息，并叠加本地地震专题数据。使用Arcpy库自动批量创建震中地区地震专题图件，并利用 Python-docx 库自动产出快速评估报告。通过 QT 开发可视化窗体，QT 是 Python 开发窗体的工具，具有良好的兼容性，支持 Windows、Linux、MacOS X 三大操作系统。用户只需要输入震中经纬度坐标等信息，程序即可自动运行。最后利用Pyinsatller对程序进行封装成单文件文本，方便使用。其技术路线图如图2所示。系统界面如图3所示。

2.1 专题图制作

地震发生后，快速产出地震专题图件是地震应急的重要部分，为了规范破坏性地震应急专题图件的制作，中国地震台网中心发布了《破坏性地震应急专题图件产出与产出标准》。本研究专题图部分主要包括地震影响场快速评估图、震中位置图、历史地震分布图、震中地区断层分布图、震中与主要城市距离图以及震区人口密度图等图件（图4）。

2.1.1 地震影响场生成

自动生成影响场的具体计算步骤为：

（1）通过聂高众等[12]在2018年提出的震中烈度与震级、震源深度的关系式计算出震中烈度大小，公式如下：

式中，Ie是震中地区的烈度，M 是震级，H 是震源深度。

（2）确定影响场的方向，本研究通过查看距离震中最近的地震活动断层的走向来确定烈度圈的方向。

（3）通过中国东部地区烈度衰减模型[13]，计算出烈度圈的长轴和短轴，公式如下：

式中，Ra 为长轴，Rb 为短轴，I 为烈度值，M 为震级，标准差取值为0.517。

（4）通过公式（2），可以计算出烈度圈的长轴和短轴，具体公式如下：

通过上述方法输入震级和震源深度，以等效圆的方式展现。其生成的烈度图如图5所示。

2.1.2 地震专题图生成

地震发生后的情况比较复杂，本研究按照《破坏性地震应急专题图件产出与产出标准》进行出图，利用Arcpy库将不同需求目錄下的mxd模板文件转化成可打印的 jpg 文件，快速出图。除了震后第一时间产出震中分布图、震区历史地震图、烈度分布图和行政区划图等外，还尽可能满足多种用途，并不断补充各类专题地图，更好地服务于地震应急（图6）。

2.2 快速评估报告产出

本文按照《江苏省防震减灾条例》以及本地区的地震应急预案，制作快速评估报告模板，内容包括震情基本信息、距主要城市距离、震中地区历史地震和构造情况、影响区域以及相应建议等，系统利用 Python-docx 库将计算的数据和前期生成的专题图插入模版自动生成快速评估报告（图7）。为领导和专家的地震救援决策提供参考和服务，以支撑应急救援的协调指挥。

目前该系统已经在镇江市防震减灾工作中开始应用，为全市、全省多次地震应急演练以及小震的处置工作提供辅助决策支持，快速产出专题图件和快速评估报告，软件也在多次演练反馈中得到了不断改进和优化。

3 结论和展望

本文利用 Python 语言，通过 QT 开发可视化窗体，采用天地图 API 调用天地图作为底图，天地图逆地理服务 API 获取基础地理信息，叠加地震专题数据，生成数据库。用户输入震中经纬度坐标等信息后，系统自动批量创建地震专题图件，并自动产出快速评估报告。实际成果已经应用于镇江市多次地震应急演练中，效果显著，专题图和快速评估报告均能在3 min 内自动生成，且程序只有一个单文件，即拷即用，部署简单，让应急效率显著提升，地震应急信息化服务水平得到显著提高。

地方地震应急系统的建设通常是长期补充优化的过程。虽然本文中基础数据使用天地图，在更新上比较方便，但系统仍有需要完善和改进的地方，如调用天地图 API 依赖网络，地震评估模型还不完善，与其他部门数据共享不尽完善等。下一步将在开发单机版系统，基础数据高效更新，以及不同行业的数据共享方面进行深入探索。

参考文献

［1］余其鹏，程晓东.基于 Geodatabase 面向对象的地震应急数据库的设计与实现?以镇江市为例[J].防灾减灾学报，2018，34（1）：66-71

Yu Q P，Cheng X D. Design and realization of earthquake emergency database based on Geodatabase：A case study of Zhenjiang City[J]. Seismological Research of Northeast China，2018，34（1）：66-71

［2］马宝君，李太岩，赵建通，等.地震应急专题图快速产出系统开发?以黑龙江省为例[J].地震地磁观测与研究，2021，42（5）：172-179

Ma B J，Li T Y，Zhao J T，et al. Development of rapid output system for thematic maps of earthquake emergency：Take Heilongjiang Province as an example[J]. Seismological and Geomagnetic Observation and Research，2021，42（5）：172-179

［3］樊琪，刘洪光.地震应急专题图离线自动出图研究与应用?以海南省为例[J].地震科学进展，2021，51（7）：304-308

Fan Q，Liu H G. Research and application of off-line automatic mapping of earthquake emergency thematic map：A casestudy of Hainan Province[J]. Progress in Earthquake Sciences，2021，51（7）：304-308

［4］韩贞辉，何重阳，谢恒义，等.基于 Python 的地震应急专题图快速出图技术研究[J].地震科学进展，2020，50（4）：34-37

Han Z H，He C Y，Xie H Y，et al. Research on rapid mapping technology of seismic emergency thematic map of Henan Province based on Python[J]. Progress in Earthquake Sciences，2020，50（4）：34-37

［5］程显洲，张勇，王小明，等.地震应急專题图批量快速出图模型研究[J].测绘与空间地理信息，2018，41（8）：135-138

Cheng X Z，Zhang Y，Wang X M，et al. Study on batch rapid map model of earthquake emergency thematic map[J]. Geomatics & Spatial Information Technology，2018，41（8）：135-138

［6］孙哲，赵曦，贾宁.基于“天地图”的地震专题图自动出图技术研究与实现?以陕西省为例[J].震灾防御技术，2016，11（2）：396-402

Sun Z，Zhao X，Jia N. Research and implementation of automatic mapping for seismic thematic map based on TianDiTu（ Map World）of Shaanxi Province[J]. Technology for Earthquake Disaster Prevention，2016，11（2）：396-402

［7］Zhou W F，Li Y X，Hou J，et al. Automated cartography of fisheries oceanographic atlas using ArcPy based on global time series griddata of marine environment[J]. IOP Conference Series：Earth and Environmental Science，2019，234（1）：012023

［8］Buckley A，Watkins D. Automated Map Production Workflows[C]//24th International Cartographic Conference，2007

［9］王斌，梁雪萍，周健.江苏及其周边地区断裂活动性与地震关系的分析[J].高原地震，2008，20（1）：38-43

Wang B，Liang X P，Zhou J. Analysis on relationship between fault activity and earthquakes in Jiangsu Province and its adjacent areas[J]. Plateau Earthquake Research，2008，20（1）：38-43

［10］黄耘，孙业君，杨云，等.江苏及邻区震群活动特征分析[J].中国地震，2011，27（1）：72-82

Huang Y，Sun Y J，Yang Y，et al. The characteristics of earthquake swarms in Jiangsu Province and its adjacent area[J]. Earthquake Research in China，2011，27（1）：72-82

［11］余其鹏，董建平.基于 GIS 的地震时空统计分析?以长三角地区为例[J].防灾减灾学报，2022，38（2）：66-72

YuQ P， Dong J P. Spatialand temporalstatisticalanalysisof earthquakes basedonGIS： Acasestudyof Yangtze river delta region[J]. Seismological Research of Northeast China，2022，38（2）：66-72

［12］聂高众，徐敬海.基于震源深度的极震区烈度评估模型[J].地震地质，2018，40（3）：611-621

Nie G Z，Xu J H. Seismic intensity evaluation model in meizoseismal area based on focal depth[J]. Seismology and Geology，2018，40（3）：611-621

［13］汪素云，俞言祥，高阿甲，等.中国分区地震动衰减关系的确定[J].中国地震，2000，16（2）：5-12

Wang S Y，Yu Y X，Gao A J，et al. Development of attenuation relations for ground motion in China[J]. Earthquake Research in China，2000，16（2）：5-12