Application of Building Spatial Distribution Method Based on Remote Sensing Image in Earthquake Disaster Loss Pre-assessment
-
摘要: 房屋建筑数据是地震重点危险区预评估工作的基础,需基于获取到的房屋建筑信息开展人员伤亡、经济损失、救援物资需求等预评估工作。历年地震重点危险区预评估工作能够通过现场调查得到的房屋建筑信息占比极小,仅能进行抽样调查。因此,为批量完成危险区内全部房屋建筑损失估算,需基于遥感影像获取房屋建筑矢量数据,并建立数据库。为实现全国地震重点危险区预评估工作中大批量建筑物矢量化数据的快速获取,本文采用基于遥感影像的建筑物空间分布数据批量获取方法,得到地震重点危险区内建筑物空间矢量数据,结合现场抽样调查得到的建筑物属性信息,建立地震重点危险区建筑物空间分布数据库,进而为地震重点危险区灾害损失预评估工作提供数据基础。本文采用的方法可广泛应用于地震重点危险区房屋建筑信息获取工作中,可提高工作效率,降低工作成本,提升预评估工作的科学性和准确性。Abstract: Building data is the basis of pre-assessment of key earthquake risk areas. Based on the obtained building information, pre-assessment of casualties, economic losses and relief material needs can be carried out. In the pre-assessment of key earthquake risk areas over the years, the proportion of building information obtained through field investigation is very small, so it can only be sampled. Therefore, in order to estimate the loss of all buildings in the dangerous area in batches, the best way is to obtain the vector data of buildings based on remote sensing images and establish a database. In order to achieve the rapid acquisition of large quantities of building vector data in the pre-assessment of key seismic risk areas in China, this paper proposes a method of building spatial distribution data batch acquisition based on remote sensing images, and then obtains the building spatial vector data in the key seismic risk areas, combined with the building attribute information obtained from the field sampling survey, The spatial distribution database of buildings in key seismic risk areas is established to provide data basis for disaster loss pre assessment in key seismic risk areas. The method proposed in this paper can be widely used in the acquisition of building information in key seismic risk areas in the future, which can not only improve the work efficiency, but also reduce the cost, and improve the scientificity and accuracy of pre-assessment.
-
Key words:
- Key seismic risk area /
- Pre assessment /
- Remote sensing image /
- Building /
- Spatial distribution
-
引言
2017年4月12日2时25分,浙江省杭州市临安市(现为杭州市临安区,考虑地震发生时行政区划名称为临安市,本文沿用临安市名称)发生4.2级地震,震中(30.08°N,119.34°E)位于临安市河桥镇与潜川镇交界附近的九龙坑,震源深度16.2km。震中附近的潜川镇、河桥镇震感强烈,於潜镇、太阳镇、湍口镇等地也有较明显的震感。地震发生后,浙江省地震局启动了Ⅲ级响应,组织工作队赶赴地震现场,在对震区地质构造概况、震源机制解、实际震情等情况进行分析的基础上,开展了烈度调查,获取了震区震害分布基本资料,并绘制了地震烈度分布图。
1. 震区地形地貌及地质构造概况
震区在地貌上属于浙西低山丘陵区,地势自西北、西南向东部倾斜。震区在大地构造上属于扬子准地台的钱塘台褶带,新构造运动主要表现为大面积的间歇性升降运动,以整体抬升为主,震区发育第四系早期的坡洪积物、侏罗系凝灰岩、志留系、奥陶系粉砂岩、泥岩及燕山期花岗岩等地层。震中附近发育区域性北东走向的马金-乌镇断裂带(图 1)。马金-乌镇断裂带是浙江省内1条重要的区域性断层(浙江省地质矿产局,1989),控制前新生代地层分区,总体走向约40°,倾向南东,倾角约80°,震区一带断层带宽12km、由多条分支组成,在东山畈采石场观察点见其中一主要分支断层的出露剖面(图 2),剖面上断层切割寒武纪厚层灰色灰岩,剖面中见3个断面(f1、f2、f3),性质皆逆断。其中,f1为单一断面,f2与f3则构成宽约2m的黄褐色粗碎裂岩;f2、f3延伸到冲沟对面,断面及黄褐色粗碎裂岩仍出露清楚(图 3),且沿2个断面还形成1—2cm的灰绿色断层泥,断层泥ESR样品测试结果为(151±30)ka BP,表明断层中更新世有过活动迹象。
2. 震源机制
采用CAP方法(Zhao等,1994;Zhu等,1996;Wang,2011),计算获得了此次地震的震源机制解。其中,节面Ⅰ走向为298°,倾向87°,滑动角-12°;节面Ⅱ走向为29°,倾向78°、滑动角为-177°;主压应力轴(P轴)方位角平均为253°;主张应力轴(T轴)方位角平均为344°。震源机制解具体数值见表 1,拟合波形如图 4所示,图中红色曲线为理论地震波形,黑色曲线为观测地震波形。
表 1 临安4.2级地震震源机制解参数Table 1. The parameters of focal mechanism solution of Lin'an M 4.2 earthquake最佳双力偶解 节面Ⅰ/° 节面Ⅱ/° 走向 倾向 滑动角 走向 倾向 滑动角 298 87 -12 29 78 -177 应力轴 T轴/° B轴/° P轴/° 方位角 倾角 方位角 倾角 方位角 倾角 344 6 104 78 253 11 CAP方法求解震源机制解受震源深度影响较大。通过设置步长1km,搜索不同深度对应的震源机制解,拟合误差越小,深度越准确。最终,得到误差最小的16.2km为此次临安M 4.2地震最佳震源深度(图 5)。
震区现代构造应力场的特征以北东东向水平压应力、北北西向水平张应力为主,北东向断层容易发生右旋走滑错动,震源机制解显示本次地震受大范围构造应力的控制。
3. 现场调查与震害特征
3.1 现场调查
依据《地震现场工作(第3部分):调查规范》、《地震现场工作(第4部分):灾害直接损失评估》(中华人民共和国国家质量监督检验检疫总局等, 2012a, 2012b)及《中国地震烈度表》(GB/T 17742—2008)(中华人民共和国国家质量监督检验检疫总局等,2009)中的技术要求,主要根据震区人的感觉、器物反应、房屋震害程度等(夏坤等,2011)开展地震烈度调查与评定工作。根据当地政府反馈的震感情况,分别在潜川镇、河桥镇、於潜镇、太阳镇和湍口镇安排调查路线,共25个调查点,取得了震区大量的基础数据和第一手资料,调查范围基本上覆盖了本次地震震区。地震现场调查路线及调查点分布见图 6。
3.2 震区烈度分布情况
通过对25个调查点的震害评估与科学考察,同时收集了当地政府提供的人员伤亡、房屋受损排查情况以及山塘水库基础设施和地质灾害点排查情况,4个调查点的烈度评定为Ⅴ度,分别为潜川镇的五里亭、茶叶山、鱼潭和过水埠,主要分布在震中的东北和西南,Ⅴ度区面积约22.6km2;其余21个调查点烈度评定为Ⅳ度(图 7)。本次地震未造成人员伤亡和较大经济损失。
Ⅴ度区人和器物的反应特征主要为:①震区主要建筑为3层房屋,绝大多数调查点居民居住在2楼,反映震感强烈,多数人梦中惊醒,少数逃出户外,个别居民在室外过夜直至天明;②绝大多数被调查居民听到地声,且反映声音很响,如同打雷的轰轰声或放炮声,部分居民能明确指出声音由震源方向传来;③大多数被调查居民听到门窗、屋顶、屋架颤动作响,且门窗响声较大;④多数被调查居民感觉床和悬挂物有较大幅度的晃动。
Ⅴ度区房屋震害特征主要为:①个别土木、砖木结构房屋有老裂缝加宽、加长及墙皮脱落现象;②个别砌体结构房屋出现细微裂缝和老裂缝加宽、加长及灰土掉落现象(图 8)。
通过调查分析本次地震的房屋破坏情况,认为造成破坏的主要因素有以下方面:
(1)房屋质量
震区内产生破坏的个别房屋属于老旧房屋,建造时基本没有采取任何抗震措施,而且房龄长,一般多在20多年以上,有的甚至30、40年,这些房屋多为夯土(砖、石)墙木结构,或者某一面承重墙为空斗砖墙,由于房龄长、年久失修且普遍存在裂缝,有的墙体已倾斜,特别是房屋的屋顶、屋内木架构,梁木多有腐烂现象;另外,震区内部分建造于2000年以前的砌体房屋中,仅有少数有地梁或圈梁,多数无地梁、圈梁和构造柱,抗震性能较差,加之地基基础的不均匀,个别房屋墙体出现细微裂缝或抹灰缝。
(2)场地效应
震区位于山区,多数自然村分布于陡坡、边坡等地,导致部分调查点上的个别房屋在建造时为天然地基加垒石形成的不均匀基础(属于抗震不利地段),地震时容易产生沉降差异,抗震稳定性差,易出现轻微破坏。
4. 结语
通过对临安4.2级地震开展现场调查,得出如下认识:
(1)震源机制解显示长轴方向为北东向,与马金-乌镇断裂带走向基本一致;地震现场地震烈度调查结果表明,本次地震的影响烈度最高为Ⅴ度,面积22.6km2,烈度圈长轴方向为北东向,与地震震源机制解及北东向马金-乌镇断裂带方向一致。
(2)本次地震震中附近的潜川镇、河桥镇震感强烈,但未造成人员伤亡和较大经济损失,震中区少数房屋有轻微开裂和老裂缝加长、加宽现象。
(3)与浙江省内其它类似震级大小的地震相比,此次地震震源深度较深(16km),地面的地震破坏影响相对较小,因此虽然震感强烈,但震害相对较轻。
-
表 1 格网数据表结构
Table 1. Grid data table structure
字段名 中文名 字段类型 字段长度 是否可空 OBJECTID 要素ID Integer 14 否 GRIDID 网格ID VarChar 14 否 BUILDNUM 建筑物数量 Number 20 是 EARTHAREA 建筑物占地面积 Number 10 是 BUILDAREA 建筑物面积 Number 10 是 表 2 建筑物数据表结构
Table 2. Buildings data table structure
字段名 中文名 字段类型 字段长度 是否可空 OBJECTID 要素ID Integer 14 否 COUNTYID 区县ID VarChar 14 否 COUNTYNAME 区县名称 VarChar 40 是 CITYNAME 市名称 VarChar 40 是 PROVINCENAME 省名称 VarChar 40 是 BUILDID 建筑ID VarChar 14 是 GRIDID 网格ID VarChar 14 是 表 3 区县建筑数据表结构
Table 3. District and county building data table structure
字段名称 字段描述 字段类型 字段长度 OBJECTID 唯一标识符 Number 14 PROVINCE 省编码 nvarchar 14 XZQMC 区县行政区名称 nvarchar 40 XZQEMC 行政区名称拼音 nvarchar 40 FULLNAME 全名 nvarchar 40 CITY 市代码 nvarchar 10 CITYNAME 市名称 nvarchar 10 XZQDM 区县行政区代码 nvarchar 14 EARTHAREA 占地面积 Number 10 BUILDAREA 建筑面积 Number 10 BUILDNUM 建筑物数量 nvarchar 10 表 4 山西省临汾市尧都区和晋中市榆次区建筑物结构类型比例(单位:%)
Table 4. Proportion table of building structure types in Yaodu district of Linfen city and Yuci district of Jinzhong city, Shanxi province(Unit:%)
地区结构类型 砌体结构 钢混结构 钢结构 单层民宅 单层厂房 临汾市尧都区 70.67 5.48 0.22 21.08 2.54 晋中市榆次区 57.47 5.53 1.75 31.28 3.97 表 5 山西省临汾市尧都区和晋中市榆次区建筑物层数比例(单位:%)
Table 5. Proportion table of building floors in Yaodu district of Linfen city and Yuci district of Jinzhong city, Shanxi province(Unit:%)
地区层数 10层及以上 7~9层 3~6层 1~2层 临汾市尧都区 0.88 0.39 14.93 83.73 晋中市榆次区 1.42 0.35 16.37 81.86 表 6 山西省浑源县建筑物结构类型比例(单位:%)
Table 6. Proportion table of building structure types in Hunyuan county of Shanxi province(Unit:%)
城镇地区 农村地区 结构类型 占比 结构类型 占比 砖混结构 77.16 砖混结构 7.60 砖木结构 3.20 砖木结构 72.30 土木结构 0.49 土木结构 11.20 土石结构 0 土石结构 8.90 钢混结构 18.66 钢混结构 0 钢结构 0.49 钢结构 0 木结构 0 木结构 0 其他结构 0 其他结构 0 表 7 设定地震预评估结果
Table 7. Table of pre-evaluation results
震级 灾区面积/km2 灾区人口/万人 灾区GDP/亿元 预评估结果 死亡人数/人 受伤人数/人 需紧急安置人数/人 需紧急救援力量/人 需帐篷数/顶 5.5 1 096 36 17 0~20 10~80 14 600~16 600 2 000~2 200 2 500~2 700 6.0 4 723 88 36 5~60 50~200 39 800~41 800 5 200~5 400 6 700~6 900 6.5 14 956 261 160 15~100 600~2 000 97 400~99 400 14 600~16 600 15 400~17 400 -
[1] 邓宏宇, 孙柏涛, Dong W. M., 2013. 遥感技术在地震应急基础数据库建设中的应用. 地震工程与工程振动, 33(3): 81—87.Deng H. Y., Sun B. T., Dong W. M., 2013. The application of remote sensing technology in the construction of earthquake emergency foundational database. Journal of Earthquake Engineering and Engineering Vibration, 33(3): 81—87. (in Chinese) [2] 罗磊, 汪斌, 2018. 基于GIS的震后修复建筑物空间分布格网化方法研究. 地震工程学报, 40(6): 1362—1365. doi: 10.3969/j.issn.1000-0844.2018.06.1362Luo L., Wang B., 2018. GIS-based gridding method for spatial distribution of post-earthquake repaired buildings. China Earthquake Engineering Journal, 40(6): 1362—1365. (in Chinese) doi: 10.3969/j.issn.1000-0844.2018.06.1362 [3] 石玉成, 高晓明, 谭明等, 2013. 2013年岷县漳县6.6级地震灾害损失评估. 地震工程学报, 35(4): 717—723. doi: 10.3969/j.issn.1000-0844.2013.04.0717Shi Y. C., Gao X. M., Tan M., et al., 2013. Disaster loss assessment of the Minxian-Zhangxian MS6.6 Earthquake, 2013. China Earthquake Engineering Journal, 35(4): 717—723. (in Chinese) doi: 10.3969/j.issn.1000-0844.2013.04.0717 [4] 魏本勇, 聂高众, 苏桂武等, 2017. 基于公里格网的地震人员埋压风险评估——以云南昭通市为例. 见: 2017中国地球科学联合学术年会论文集. 北京: 中国地球物理学会, 1938—1940. [5] 熊俊楠, 韦方强, 苏鹏程等, 2013. 基于多源数据的四川省GDP公里格网化研究. 应用基础与工程科学学报, 21(2): 317—327. doi: 10.3969/j.issn.1005-0930.2013.02.013Xiong J. N., Wei F. Q., Su P. C., et al., 2013. Research on GDP Spatialization approach of Sichuan province supported by multi-source data. Journal of Basic Science and Engineering, 21(2): 317—327. (in Chinese) doi: 10.3969/j.issn.1005-0930.2013.02.013 [6] 徐国栋, 袁艺, 方伟华等, 2011. 玉树7.1级地震震后损失快速评估. 地震工程与工程振动, 31(1): 114—123.Xu G. D., Yuan Y., Fang W. H., et al., 2011. Fast loss assessment of M7.1 Yushu earthquake. Journal of Earthquake Engineering and Engineering Vibration, 31(1): 114—123. (in Chinese) [7] 杨海霞, 2015. 建筑物空间分布格网化方法研究. 北京: 中国地震局地震预测研究所.Yang H. X., 2015. Research on gridding method of buildings’ spatial distribution. Beijing: Institute of Earthquake Science, China Earthquake Administration. (in Chinese) [8] 杨海霞, 王晓青, 窦爱霞等, 2015. 基于RS和GIS的建筑物空间分布格网化方法研究. 地震, 35(3): 136—146. doi: 10.3969/j.issn.1000-3274.2015.03.014Yang H. X., Wang X. Q., Dou A. X., et al., 2015. Multi-source and Multi-factor gridding method of building distribution based on RS and GIS. Earthquake, 35(3): 136—146. (in Chinese) doi: 10.3969/j.issn.1000-3274.2015.03.014 [9] 周中红, 陈文凯, 何少林等, 2019. 基于人口公里格网的地震压埋人员分布评估方法的应用与评价——以甘肃岷县漳县6.6级地震为例. 地震研究, 42(2): 288—294. doi: 10.3969/j.issn.1000-0666.2019.02.018Zhou Z. H., Chen W. K., He S. L., et al., 2019. Application and evaluation on a evaluating method for distribution of earthquake buried personnel based on population kilometer grid: taking Minxian-Zhangxian M6.6 earthquake as an example. Journal of Seismological Research, 42(2): 288—294. (in Chinese) doi: 10.3969/j.issn.1000-0666.2019.02.018 [10] Sridharan H., Qiu F., 2013. A spatially disaggregated areal interpolation model using light detection and ranging-derived building volumes. Geographical Analysis, 45(3): 238—258. doi: 10.1111/gean.12010 -