Earthquake Disaster Risk Assessment in Shenzhen and Related Discussion
-
摘要: 以第一次全国自然灾害综合风险普查的总体要求为框架,结合深圳市特点,以深圳市(深汕特别合作区除外)单体调查房屋为评估数据基础,选择适合深圳市房屋特点的地震易损性矩阵,在开展单体房屋震害评估的基础上,进行人员死亡和经济损失评估。深圳市地震灾害风险评估考虑4个超越概率水平,即50年63%(常遇地震)、50年10%(设防地震)、50年2%(罕遇地震)、100年1%(极罕遇地震)的地震作用影响。文中讨论了100年1%(极罕遇地震)地震作用下深圳市建筑物破坏所导致的人员死亡风险和建筑物直接经济损失风险,并给出风险治理的参考建议。Abstract: This study utilizes the overall requirements of the first national comprehensive natural disaster risk survey as its framework and incorporates the specific characteristics of Shenzhen City. Using individual surveyed houses in Shenzhen City (excluding the Shenzhen-Shantou Special Cooperation Zone) as the basis for assessment data, it presents a seismic vulnerability matrix tailored to Shenzhen's housing characteristics. The study assesses personnel mortality and economic losses based on seismic damage evaluations of individual houses. The seismic disaster risk assessment in Shenzhen considers four exceedance probability levels: 63% in 50 years (common earthquakes), 10% in 50 years (design earthquakes), 2% in 50 years (rare earthquakes), and 1% in 100 years (extremely rare earthquakes). It discusses the risks of personnel mortality and direct economic losses to buildings from the 1% in 100 years (extremely rare earthquake) event and provides recommendations for risk management.
-
Key words:
- Disaster risks /
- Earthquake disaster risk assessment /
- Shenzhen
-
引言
福建省地震预警项目简易烈度计台网系统用于评估地震对乡镇的影响,每个乡镇布置一个测站,按平均台间距10km左右设置,单台定位精度平均小于5km,最大不超过8km,运行维护率≥90%。目前福建省地震预警系统的总体规划分两期实施:一期在地震活动性较强或者地震影响较大的地区,选择人口较多和经济发达区域进行设备安装;二期在建设地区,如闽西与闽北农村山区等地稀疏布点。地震台网数据处理系统主要采用Earthworm进行搭建,拟在全省布设900台简易烈度计。该数据处理系统由C语言编写,由不同模块组建而成,每个模块实现不同的功能。该系统自从2015年6月正式运行以来,能够正常接收各地震台站的实时波形数据,从而实现地震台网中心地震速报、数据产出、数据管理等各项功能(王喆,2006)。
本文介绍了由美国地质勘探局开发的Earthworm系统的设计目标、思想和原则,重点介绍了Earthworm框架结构与各模块功能的搭建与配置、Earthworm系统的架构特点、主要工作原理和数据处理流程,以及其应用于福建省地震预警项目的相关情况(吴永权等,2010)。
1. Earthworm简介
Earthworm项目始于1993年,主要目的是为了解决美国地震区域台网出现的问题。当时,区域内的地震台网存在的主要问题有:观测设备陈旧,自动处理系统老化,维修费用逐年增长;仪器性能落后,地震研究的发展需要从精密传感器中获取研究数据;社会公众需要新的、显而易见的实时产品了解详情;最后,经费不足导致大多数台网不再支持区域系统的发展工作。
为了更好解决服务器老旧与新地震观测仪器相结合的问题,Earthworm项目应运而生。
1.1 Earthworm的特点
Earthworm系统是以避免早期地震处理系统的设计缺陷,保证地震预警的顺利实施为目标而建设的。系统具有以下5个方面的优点:
(1)模块化:在硬件以及软件方面,系统把所执行的每一个函数封装成一个模块,该模块能独立于其他模块运行。模块化使得系统内的任何一个关键性模块能与系统内的其他模块保持独立。因此,新的试验性需求模块可以在不破坏现有系统的基础上添加并保证良好的性能。
(2)系统独立性:由于不同的模块可独立运行在不同的计算机硬件和操作系统中,并作为一个系统运转,因此该系统可以无中断、在任意操作系统的计算机之间移植。在实际运用中,只使用各种计算机系统标准化的部分,而不使用任何无法规避的系统特殊功能。
(3)可扩展性:该系统提供较好的性价比以适应各种类型层次的网络需求。同时,用户可定制个性的需求,编写新的模块运行。
(4)连通性:该系统能够实现与其它自动实时系统、互动分析系统以及各种发布方式之间快速可靠地通信。其目标是在各层次提供自动和交互式的接口,使系统的配置范围可以从完全单机独立操作到分布式系统单个节点的运行。
(5)鲁棒性:在地震危机时刻,由于输入数据和电源可能中断、系统负荷急剧增加,只有高鲁棒性的地震报警系统才能向媒体和应急机构提供实时准确信息。该系统针对地震处理数据的错误检测及恢复、恢复时间、故障弱化和负荷控制等问题具有极好的可靠性(马强,2008)。
1.2 Earthworm框架的特点
Earthworm系统中每个模块都执行特定的任务,如数据采集,震相拾取等。模块之间通过广播和监听(如波形数据包和相位捡拾等)各种消息来通信。消息传递类似于无线电通讯:它由用于传递消息的“传输媒介”和一套标准模块程序组成。这些模块程序类似于运行在传输媒介中的多频率双向无线电设备。在传输媒介上模块可以使用这些标准程序广播并监听。
2. Earthworm地震数据处理系统的架构与功能模块
福建省简易烈度台网设备一期规模为300台左右,通过中国电信的MPLS VPN组网方式将数据传输至厦门地震勘测研究中心与福建省地震局。每个观测台站提供2M SDH链路带宽供烈度计上报数据使用,省局与厦门地震勘测研究中心分别采用一条10M光纤汇聚专线接入中国电信VPN平台,用于收集各观测台站数据。对于部分未铺设SDH链路的区域,采用3G信号进行辅助接入Earthworm的地震预警软件框架系统。
基于Earthworm地震数据处理系统的架构采用积木式组合的方式进行搭建(图 1),主要包括:
第1部分,数据接收,通过MPLS VPN组网方式将各台网观测的烈度计产生的数据在数据中心进行汇聚,主要采用Paler2ew模块进行功能搭建;
第2部分,地震波形展示,主要采用wave_server模块进行功能搭建;
第3部分,地震数据处理,主要采用pick_eew、tcpd、ShakeMap的模块进行功能搭建,pick_eew模块负责自动震相拾取功能,tcpd模块负责地震定位与震级计算功能,ShakeMap模块负责烈度图绘制功能;
第4部分,地震数据存储,主要采用Tbuf2mseed、Ewmseed Archiver的模块进行功能搭建,Tbuf2mseed模块负责进行数据格式转换功能,Ewmseed Archiver模块负责进行数据文件存档的功能。
3. Earthworm的安装步骤与配置
3.1 Earthworm源码包的获取
可以通过网址http://folkworm.ceri.memphis.edu/ew-doc/#overview下载最新的Earthworm安装源码包进行安装。目前福建省烈度计系统采用的最新版已更新为v7.9版,版本之间的主要区别在于系统模块功能的差异(Earthworm官方网址,2015)。
3.2 Earthworm源码包目录结构说明
Earthworm源码包目录结构主要有bin,environment,ewdoc,include,include_app,params,src,lib等,下面分别对该源码目录功能进行解释:
(1)bin目录为系统可执行程序存放位置;
(2)environment目录为系统环境变量存放位置;
(3)include目录为系统文件存放位置;
(4)include_app目录为系统程序应用存放位置;
(5)params目录为各程序参数配置文件存放位置;
(6)src目录为软件源码存放位置,可单独重新编译;
(7)lib目录为程序执行所依赖库文件存放位置。
3.3 Earthworm的安装与调试
(1)安装Earthworm前的环境准备
在安装Earthworm前,系统需要先安装make、gcc、java这3个软件,用于对Earthworm的程序包进行编译,3个软件的安装均可以通过yum源进行安装,命令如下:
yum install –y make & & yum install –y gcc & & yum install–y java
(2)创建安装目录结构
当前系统采用的是CentOS 7,64位操作系统进行搭建,我们在根目录下创建一级目录xmzk,然后再创建二级目录earthworm,earthworm目录下分别创建run目录与earthworm_7.9源码包,然后在run目录下分别创建params目录与logs目录,目录结构如下所示:
/xmzk
/earthworm
/run
/params(各程序参数配置位置)
/logs(各程序日志文件存放位置)
/data(各台站接收数据存放位置)
/earthworm_v7.9(源码包)
(3)拷贝全局参数与环境变量至安装目录
将/xmzk/earthworm/earthworm_v7.9/environment目录中的earthworm_global.d、earthworm. d、ew_linux.bash这3个文件拷贝至/xmzk/earthwom/run/params中。
(4)修改安装环境配置文件
安装Earthworm需要变更第二步拷贝至/xmzk/earthwom/run/params中的ew_linux.bash这个环境配置文件中几处环境变量的值,以适应我们自己的安装环境与配置需求,主要有以下6处需要变更:
① export EW_HOME=/home/xmzk/earthworm #定义Earthworm安装位置的环境变量
② export EW_VERSION=earthworm_7.9 #定义Earthworm版本的环境变量
③ export EW_RUN_DIR=/home/xmzk/earthworm/run #定义Earthworm运行位置的环境变量
④ export EW_PARAMS=/home/xmzk/earthworm/run/params #定义Earthworm各模块参数档位置的环境变量
⑤ export EW_LOG=/home/xmzk/earthworm/run/logs #定义Earthworm日志文件存放位置的环境变量
⑥ export EW_DATA_DIR=/home/xmzk/earthworm/run/data #定义Earthworm波形数据文件存放位置的环境变量
(5)拷贝启动文件
拷贝Earthworm源码包里params文件夹中的“startstop_unix.d”启动程序至/home/xmzk/ earthworm/run/params中,并根据自己的环境在文件中增加与删除功能模块。
(6)Earthworm运行与开机自启动
打开终端命令行,键入“startstop”命令开启Earthworm,如果命令行中提示各模块为alive说明earthworm运行正常,如有些模块出现zoombi的字样说明该模块未正常工作,需单独检查该模块,逐一排错,对其重新加载或编译。
4. 模块功能搭建与配置档案设置
4.1 单一模块功能在系统中的添加安装步骤
可以通过以下几个步骤添加任一模块功能:
第一步:编译该模块的源码,编译模块功能命令:make–f makefile name(模块名称),执行完可生成该模块的可执行程序;
第二步:将编译完成后的可执行程序拷贝至/xmzk/earthworm/earthworm_v7.9/bin目录下;
第三步:修改earthworm.d文件,在earthworm.d文档中新增写入该模块名称id与该模块数据类型,使得Earthworm系统可识别出该模块;
第四步:复制该模块的参数配置文件拷入/home/xmzk/earthworm/run/params中;
第五步:修改startstop.d文件,增加新的模块,在startstop.d文档中写入该模块的执行文件名称与配置文件名称。startstop.d参数档的作用主要是定义ring的个数与定义需要开启多少模块,需要开启的模块功能均需写入该文件中;
第六步:执行recon命令,重载startstop_unix.d文件,即可将该模块加载到现在运行的Earthworm环境中。
4.2 数据接收模块的关键参数配置与说明
该模块主要通过修改/home/xmzk/earthworm/run/params/palert2ew.d的参数文件进行配置,其主要参数设置说明如下:
4.3 数据流服务模块的关键参数配置与说明
该模块主要通过修改/home/xmzk/earthworm/run/params/Wave_serverV.d的参数文件进行配置,其主要参数设置说明如下:
4.4 数据处理—P波拾取模块的关键参数配置与说明
该模块主要通过修改/home/xmzk/earthworm/run/params/pick_eew.d的参数文件进行配置,其主要参数设置说明如下:
4.5 数据处理—地震报告生成的关键配置与说明
该模块主要通过修改/home/xmzk/earthworm/run/params/tcpd.d的参数文件进行配置,其主要参数设置说明如下:
4.6 数据处理—烈度图生成模块的关键配置与说明
该模块主要通过修改/home/xmzk/earthworm/run/params/shakemap.d的参数文件进行配置,其主要参数设置说明如下:
4.7 数据存储模块的关键配置与说明
该模块的主要通过修改/home/xmzk/earthworm/run/params/ewmseed Archiver.d的参数文件进行配置,其主要参数设置说明如下:
5. Earthworm系统常用命令总结
在Earthworm应用的过程中,对一些常用命令进行了总结,如表 1(Earthworm官方网址,2010)。
表 1 Earthworm常用命令Table 1. Commonly used commands of Earthworm序号 命令 作用 示例 1 startstop 在计算机上启动和停止所有Earthworm模块,这个模块是Earthworm系统的核心 2 restart 手动重新启动单个模块 3 recon 允许向正在运行的Earthworm添加新模块或RING 4 pau 完全关闭Earthworm和所有模块/RING pau startstop_nt.d 5 stopmodule 给定一个Earthworm模块进程ID,stopmodule停止它,startstop将其标记为“Stop”,以防止statmgr重新启动它 用法:Stopmodule < process_id > < 可选:startstop配置文件 > 6 quit 关闭earthworm系统 7 restart 重启系统中任何一支模块 restart module id 8 status 输出到屏幕EARTHWORM RING和模块的状态。 status < optional:startstop config file > 9 sac2tb 将sac文档转成tank文档 sac2tb [-n max-samples] infile > > outfile 10 ms2tb 将miniseed文档转成tank文档 ms2tb [-n max-samples] infile > > outfile 11 remux_tbu 将地震观测数据资料按时间进行排序 remux_tbuf < demuxed_file > < outputfile > 12 tankcut 将地震观测数据进行剪切 tankcut -s StartTime [-e EndTime|-d Duration] intank outtank all times for -s and -e options must be in YYYYMMDDHHMMSS format 13 dumpwave 将tank二进制文档转成ASCII文档 dumpwave < filename > < optional:pinno > 14 Findwave 查看数据池中是否已接收到地震波形数据 Findwave WAVE_RING 80 1 out.txt w 15 Sniffwave 显示服务器每秒收到封包值 Sniffwave WAVE_RING L001 HLZ TW–y 16 getmenu 检查waveserverV是否存活,有哪些资料 Getmenu 127.0.0.1:16001 6. 系统运行情况
6.1 波形数据接收展示
Earthworm地震数据处理系统可以通过swarm实时波形查看工具实时查看接收到的地震波(图 2),下载地址:https://volcanoes.usgs.gov/software/swarm/download.php。
6.2 Earthworm系统的硬件资源使用
福建省地震局Earthworm系统硬件配置采用一台HP DL580 Gen9服务器,4颗Inter至强E7-4820v3型号CPU,内存128GB,从系统截图(图 3)观察Earthworm各模块的资源利用率相对较低。
6.3 地震共享池数据结构
Earthworm数据录入规则是一秒钟一个封包,采样数据分别录入到各自的RING中,均可通过命令sniffwave查看获取的数据结构,通过sniffwave命令查得RING中每秒钟的数据结构如图 4所示,主要包括台站名称、数据封包时间、封包序列号、包长度、波长电频信号等信息。
6.4 地震报告产出流程
Earthworm系统在地震到来时获取4个台站数据后就可以触发计算,但还需验算,通过pick_eew进行P波拾取,再通过tcpd模块产出地震报告,主要工作流程如下:
第一步:进行autopicking工作,若autopicking的pa及pv大于执行模块参数文件设定值,则将数据放进PICK_RING里。执行模块的参数文件中提到的其中一个档案就是设定picking条件的档案;
第二步:tcpd模块会先以PICK_RING的数据进行grouping。采用Geiger method,如果RMS太大,则剔除RMS最大的测站,若小于4个测站则停止计算,若是grouping的触发测站达到6站则继续执行计算;
第三步:计算震源坐标(x,y,z)以及发震时间t。利用(x,y,z)计算到各个测站的理论走时;
第四步:发震时间加上理论走时,就是理论到时。如理论到时与autopicking的残差residual太大,则回到第三步,重新计算新的震源(x',y',z'),执行后续的步骤。如果residual一直无法收敛,则会将residual最大的测站剔除,再重新计算,若测站数小于4站,则停止计算。
通过该系统接收数据并进行分析,得到地震报告的产出如图 5。
6.5 地震图件生成样式
通过系统中产生的数据得到的PGA地震加速度峰值分布图如图 6所示。
7. 结束语
Earthworm系统引进了搭积木式模块化架构理念,将地震数据处理功能需求整合为6个软件模块,即数据接收模块、实时波形查看、自动震相拾取、定位与震级计算、数据格式转换、数据文件存档,各模块各司其职又相互关联。各模块之间分工明确,结构简单,具有优秀系统架构的稳定性、跨平台性、安全性、可扩展性、可定制化等特点。Earthworm系统运行稳定,功能强大,深受用户好评。其模块化设计,将庞大的系统功能分布运行于多台服务器,降低软件系统对服务器的性能要求,且架构配置简单,适合在各种不同规模的地震台网进行部署,用户可以根据自己的需求采用不同方法组建系统。
-
表 1 深圳市(不含深汕)单体房屋结构类型归并对应表
Table 1. Reclassification for single building structure types in Shenzhen (Excluding SSCZ)
调查结构类型 归并后结构类型 钢结构 钢结构 钢筋混凝土结构 钢筋混凝土结构 砌体结构 砌体结构 木结构 砖(土、石)木、木、简易结构 其他 子类钢混+钢结构 钢结构 子类砖木 砖(土、石)木、木、简易结构 子类土、石木结构 砖(土、石)木、木、简易结构 子类砌体+钢架 钢结构 子类简易结构 砖(土、石)木、木、简易结构 子类混合结构 8层及以上的非住宅钢混+钢结构 钢筋混凝土结构 7层及以下的住宅砌体+其他结构 砌体结构 非住宅钢混+钢结构 钢筋混凝土结构 子类组合结构 钢混+钢结构 钢结构 砌体+其他结构 砌体结构 子类其他住宅 1~2层 砖(土、石)木、木、简易结构 3~8层 砌体结构 8层以上 钢筋混凝土结构 子类其他非住宅 1~2层 砖(土、石)木、木、简易结构 3~8层 砌体结构 8层以上 钢筋混凝土结构 子类无 8层及以上(非住宅) 钢筋混凝土结构 7层及以下(住宅和非住宅) 砌体结构 表 2 深圳市(不含深汕)单体房屋栋数分区统计表
Table 2. Statistical results of single building structure types in each district of Shenzhen (Excluding SSCZ)
区划名称 房屋数量/栋 钢筋混凝土结构 钢结构 砌体结构 砖(土、石)木、木、简易结构 合计 大鹏新区 5023 1332 18968 166 25489 盐田区 3228 96 1661 7 4992 坪山区 31092 1469 8925 14 41500 光明区 26756 1500 19650 41 47947 罗湖区 15243 628 1634 23 17528 福田区 25606 876 4060 44 30586 南山区 23622 1675 5944 148 31389 龙华区 50675 2508 17477 76 70736 龙岗区 113887 5372 35962 4918 160139 宝安区 102287 3894 45550 77 151808 合计/栋 397419 19350 159831 5514 582114 表 3 深圳市(不含深汕)各区人口表
Table 3. Population in each district of Shenzhen (Excluding SSCZ)
区划名称 人数/人 面积/km2 密度/(人·km-2) 大鹏新区 156236 295.31 529 盐田区 214225 77.89 2750 坪山区 551333 163.83 3365 光明区 1095289 155.22 7056 罗湖区 1143801 78.98 14482 福田区 1553225 73.50 21134 南山区 1795826 178.09 10084 龙华区 2528872 175.74 14390 龙岗区 3979037 387.92 10257 宝安区 4476554 384.09 11655 合计 17494398 1970.57 8778 表 4 基岩地震动衰减关系模型系数
Table 4. Coefficients of attenuation relation model of bedrock ground motion
反应谱周期T/s 方向 A1 B1 A2 B2 C D E 标准差S PGA 长轴 2.024 0.673 3.565 0.435 2.329 2.088 0.399 0.245 短轴 1.204 0.664 2.789 0.420 2.016 0.944 0.447 0.245 0.20 长轴 2.558 0.643 3.680 0.470 2.309 2.088 0.399 0.261 短轴 1.779 0.628 2.918 0.454 1.999 0.994 0.447 0.261 1.00 长轴 0.226 0.895 2.409 0.559 2.157 2.088 0.399 0.300 短轴 −0.599 0.895 1.644 0.550 1.873 0.994 0.447 0.300 2.00 长轴 −0.666 0.936 1.247 0.641 2.047 2.088 0.399 0.342 短轴 −1.449 0.934 0.516 0.632 1.779 0.994 0.447 0.342 6.00 长轴 −1.432 0.859 −1.432 0.859 1.857 2.088 0.399 0.333 短轴 −2.041 0.841 −2.041 0.841 1.617 0.994 0.447 0.333 表 5 钢结构房屋易损性矩阵
Table 5. Vulnerability matrix of steel structure house
烈度 基本完好 轻微破坏 中等破坏 严重破坏 毁坏 Ⅵ 100% 0 0 0 0 Ⅶ 96.72% 3.03% 0.25% 0 0 Ⅷ 50.76% 45.96% 3.03% 0.25% 0 表 6 钢筋混凝土结构房屋易损性矩阵
Table 6. Vulnerability matrix of reinforced concrete house
烈度 基本完好 轻微破坏 中等破坏 严重破坏 毁坏 Ⅵ 99.75% 0.25% 0 0 0 Ⅶ 88.18% 11.82% 0 0 0 Ⅷ 49.00% 38.10% 12.07% 0.83% 0 表 7 砌体结构房屋易损性矩阵
Table 7. Vulnerability matrix of masonry house
烈度 基本完好 轻微破坏 中等破坏 严重破坏 毁坏 Ⅵ 76.46% 20.19% 2.95% 0.40% 0 Ⅶ 48.48% 31.36% 15.32% 4.46% 0.38% Ⅷ 16.01% 31.85% 35.23% 14.61% 2.30% 表 8 砖(土、石)木、木、简易结构房屋易损性矩阵
Table 8. Vulnerability matrix of brick (raw-soil, stone)with timber, timber, and simple structure houses
烈度 基本完好 轻微破坏 中等破坏 严重破坏 毁坏 Ⅵ 45% 36% 13% 6% 0 Ⅶ 16% 28% 30% 16% 10% Ⅷ 7% 12% 29% 27% 25% 表 9 不同分区分类的修正系数
Table 9. Correction factors for different partition classification
分区 一类地区 二类地区 三类地区 四类地区 五类地区 六类及以上地区 修正系数 0.6 0.76 0.85 0.9 0.95 1 表 10 人员死亡地震灾害风险等级分级指标
Table 10. Earthquake disaster risk grading index by fatality
风险等级 分级指标/人 Ⅰ级 死亡人数≥300 Ⅱ级 300>死亡人数≥150 Ⅲ级 150>死亡人数≥50 Ⅳ级 50>死亡人数≥10 Ⅴ级 死亡人数<10 表 11 建筑物直接经济损失地震灾害风险等级分级指标
Table 11. Earthquake disaster risk grading index by direct economic loss of buildings
风险等级 分级指标 Ⅰ级 (直接经济损失/区域内上年度GDP)≥75% Ⅱ级 75%>(直接经济损失/区域内上年度GDP)≥45% Ⅲ级 45%>(直接经济损失/区域内上年度GDP)≥25% Ⅳ级 25%>(直接经济损失/区域内上年度GDP)≥15% Ⅴ级 (直接经济损失/区域内上年度GDP)<15% 表 12 峰值加速度区间与烈度对照表
Table 12. Relationship between seismic intensity and peak ground acceleration
烈度 峰值加速度/Gal Ⅵ以下 PGA<50 Ⅵ 50≤PGA<100 Ⅶ 100≤PGA<200 Ⅷ 200≤PGA<400 Ⅸ PGA≥400 表 13 100年超越概率1%地震作用下不同破坏状态建筑面积评估统计结果
Table 13. Evaluation and statistical results of building area under different damage states with the 1% probability of exceedance in 100 years in Shenzhen
烈度 建筑面积/m2 基本完好 轻微破坏 中等破坏 严重破坏 毁坏 VIII 48556.40 ×10438477.72 ×10412708.26 ×1041255.31 ×10485.75×104 VII 14226.85 ×1042045.74 ×104117.42×104 34.20×104 3.35×104 总计 62783.25 ×10440523.46 ×10412825.68 ×1041289.50 ×10489.10×104 -
《2016−2025中国大陆地震危险区与地震灾害损失预测研究》项目组,2020.2016−2025年中国大陆地震危险区与地震灾害损失预测研究. 北京:中国地图出版社,245−261.《Forecasting Researches on Earthquake Risk Regions and Disaster Loss of Chinese Mainland During 2016 to 2025》Project Team, 2020. Forecasting researches on earthquake risk regions and disaster loss of Chinese mainland during 2016 to 2025. Beijing: China Cartographic Publishing House, 245−261. (in Chinese) 曹彦波,李永强,2017. 云南地震应急关键技术与信息服务. 昆明:云南科技出版社,79−80. 陈伟光,赵红梅,李富光等,2001. 深圳市断裂构造的活动性及其对地质环境的影响. 热带地理,21(1):45−50,60. doi: 10.3969/j.issn.1001-5221.2001.01.010Chen W. G., Zhao H. M., Li F. G., et al., 2001. Fault activities and their influence upon geologic environment in Shenzhen city. Tropical Geography, 21(1): 45−50,60. (in Chinese) doi: 10.3969/j.issn.1001-5221.2001.01.010 傅征祥,李革平,1993. 地震生命损失研究. 北京:地震出版社. 胡少卿,孙柏涛,王东明等,2017. 经验震害矩阵的完善方法研究. 地震工程与工程振动,27(6):46−50.Hu S. Q., Sun B. T., Wang D. M., et al., 2017. Approach to making empirical earthquake damage Matrix. Journal of Earthquake Engineering and Engineering Vibration, 27(6): 46−50. (in Chinese) 姜慧,郭恩栋,林旭川等,2022. 城市群地震灾害风险评估方法的一点探索−以粤港澳大湾区为例. 地震学报,44(5):868−880. doi: 10.11939/jass.20220096Jiang H., Guo E. D., Lin X. C., et al., 2022. A new exploration of the risk assessment method of earthquake disasters in urban agglomerations: taking the Guangdong-Hong Kong-Macao greater bay area as an example. Acta Seismologica Sinica, 44(5): 868−880. (in Chinese) doi: 10.11939/jass.20220096 马小平,林旭川,朱瑞等,2023. 城市建筑物情景构建及地震风险评估−−以玛曲县为例. 地震工程与工程振动,43(3):46−55.Ma X. P., Lin X. C., Zhu R., et al., 2023. Scenario construction of urban buildings and earthquake risk assessment: an example of Maqu County. Earthquake Engineering and Engineering Vibration, 43(3): 46−55. (in Chinese) 唐丽华,李山有,宋立军,2016. 地震灾害风险评估方法的对比分析−−以乌鲁木齐市为例. 地震工程学报,38(5):838−845. doi: 10.3969/j.issn.1000-0844.2016.05.0838Tang L. H., Li S. Y., Song L. J., 2016. Comparative analysis of earthquake risk assessment methods: a case study on Urumqi City. China Earthquake Engineering Journal, 38(5): 838−845. (in Chinese) doi: 10.3969/j.issn.1000-0844.2016.05.0838 王东明,高永武,2019. 城市建筑群概率地震灾害风险评估研究. 工程力学,36(7):165−173.Wang D. M., Gao Y. W., 2019. Study on the probabilistic seismic disaster risk assessment of urban building complex. Engineering Mechanics, 36(7): 165−173. (in Chinese) 徐俊,余成华,汤德刚等,2012. 深圳市活断层探测与地震危险性评价. 城市勘测,(1):161−166. doi: 10.3969/j.issn.1672-8262.2012.01.057Xu J., Yu C. H., Tang D. G., et al., 2012. Active fault exploration and seismic hazard assessment in Shenzhen city. Urban Geotechnical Investigation & Surveying, (1): 161−166. (in Chinese) doi: 10.3969/j.issn.1672-8262.2012.01.057 姚新强,孙柏涛,陈宇坤等,2016. 基于震害预测的动态震害矩阵方法研究. 地震工程学报,38(2):318−322. doi: 10.3969/j.issn.1000-0844.2016.02.0318Yao X. Q., Sun B. T., Chen Y. K., et al., 2016. Study on method of dynamic earthquake damage matrix based on seismic damage prediction. China Earthquake Engineering Journal, 38(2): 318−322. (in Chinese) doi: 10.3969/j.issn.1000-0844.2016.02.0318 尹之潜,杨淑文,2004. 地震损失分析与设防标准. 北京:地震出版社. 张桂欣,孙柏涛,陈相兆,2017. 分区分类的生命线工程地震直接经济损失研究. 地震,37(4):69−79. doi: 10.3969/j.issn.1000-3274.2017.04.007Zhang G. X., Sun B. T., Chen X. Z., 2017. Earthquake direct economic loss of lifeline engineering based on sub-area classification. Earthquake, 37(4): 69−79. (in Chinese) doi: 10.3969/j.issn.1000-3274.2017.04.007 张令心,刘琛,刘洁平,2015. 砖砌体房屋三水准抗震能力判别方法. 土木工程学报,48(12):34−40.Zhang L. X., Liu C., Liu J. P., 2015. Three level seismic capacity identification method of brick buildings. China Civil Engineering Journal, 48(12): 34−40. (in Chinese) 周光全,谭文红,施伟华等,2007. 云南地区房屋建筑的震害矩阵. 中国地震,23(2):115−123. doi: 10.3969/j.issn.1001-4683.2007.02.001Zhou G. Q., Tan W. H., Shi W. H., et al., 2007. Seismic hazard matrix of house construction in Yunnan. Earthquake Research in China, 23(2): 115−123. (in Chinese) doi: 10.3969/j.issn.1001-4683.2007.02.001 Sun B. T., Zhang G. X., Chen X. Z., 2018. The distribution of seismic capacity of buildings in Mainland of China. In: 16th European Conference on Earthquake Engineering.Greece:ECEE. -