基于遥感影像的建筑抗震能力分析——以张家口万全区为例1

李 姜 张 合 刘志辉 许 冲

1）河北省地震局，石家庄 050021

2）河北省地震动力学重点实验室，河北三河 065201

3）中国地震局地壳应力研究所，北京 100085

引言

从近十几年国内外地震震害来看，房屋破坏是造成人员伤亡和经济损失的最主要原因。房屋建筑资料等信息一直是灾情快速评估的重要基础，尤其是房屋结构的合理分类，决定着评估结果的准确性（聂高众等，2002；帅向华等，2009；姜立新等，2012）。

近年来，随着遥感影像技术的迅速发展，高分卫星遥感分辨率提升至亚米级，为房屋基本信息的快速获取提供了途径。王晓青等（2004）和窦爱霞等（2012）通过遥感技术获取灾区震后遥感影像，并提出了遥感震害指数与定量分析模型，成功应用于汶川、玉树地震等的震害快速定量化评估中。郭建兴等（2020）提出了利用遥感影像结合居民地建筑物数据评估地震烈度的方法。马建等（2020）通过获取震后高分辨率无人机影像，对房屋进行结构分类和震害解译。李昌珑等（2019）利用三维图像模式识别各类建筑物的数量分布，并建立不同建筑物易损性模型。汶川、玉树地震等震害多次暴露出我国村镇房屋抗震能力较脆弱，快速、全面、准确掌握村镇房屋抗震能力一直是研究的重点方向（周强等，2016；李皓等，2018）。

目前，主要利用遥感影像进行震后房屋震害识别，较少用于房屋抗震能力的研究。因各地房屋结构复杂、形式多样，实际调查难度大，通过遥感技术快速获取房屋结构类型等基本信息，对分析房屋抗震能力具有现实意义。刘贾贾等(2019，2021)通过遥感影像对城镇和农村建筑物进行分类，并根据建筑物分类结果对震害风险进行初步划分。

本文通过无人机航拍和遥感影像获取房屋基本信息，通过人机交互方式提取村镇房屋结构类型，结合实地调查统计得到其他信息，建立当地村镇房屋抗震能力指数体系，分析村镇房屋抗震能力。利用遥感影像实现房屋基本信息的快速提取，提高房屋信息获取的便捷性，且具有较高的准确性和可靠性，对于地震灾害损失评估、风险普查、风险区划等工作具有重要意义。

1 研究区域

张家口市位于河北省西北部，京晋冀蒙交界地带，是“首都圈”地区地震重点监视防御区和晋冀蒙交界地震重点危险区，属于张家口−渤海地震构造带与汾渭地震构造带交汇地区，具备发生中等破坏性地震的条件和背景（李书进等，2010；张肖等，2019；茅远哲等，2019）。根据张家口市所处的地理位置、经济发展情况，部分地区新农村建设发展迅速，按抗震设防建造的房屋抗震能力存在差异。而部分地区发展较慢，存在无抗震能力的旧房，如土木、砖木、土窑结构房屋等（刘龙等，2019）。万全区位于张家口市西北部，地理位置险要，受地质构造控制，地形北高南低，东西向为西高东低，是川、丘、山并存的地区（图1），本文以该区域作为研究对象。

图1 万全区行政区划图Fig. 1 Administrative division map of Wanquan District

2 现场调查

调查数据作为万全区房屋抗震性能分析的重要基础数据，为后期房屋抗震能力分析奠定基础。基于随机均匀分布原则，对于研究区每个乡镇，选取不少于2 个调查点，调查范围确定了22 个村镇（图2），包括新河口村、代家房村、西湾村、武家庄村、新开口村、望虎台村、刘虎庄村、板山村、阳门堡村、旧羊屯村、吴家窑村、上营房村、西柳林村、洗马林镇、万全镇、沙家庄村、陈家沟村、高庙堡村、羊窖沟村、岸庄屯村、安家堡村、赵家梁村。采取拍取房屋照片、填写房屋信息、无人机测航调查点等方式开展调查，普查主要侧重于结构类型、抗震设防水平、构造措施等基本情况，详查是在普查的基础上，对调查区域内房屋进行随机抽样，更加详细地调查建造年代、层数、用途等信息（周强等，2016）。

图2 万全区地貌分区图与调查点分布Fig. 2 Geomorphic zoning map and distribution of survey points in Wanquan District

3 建筑物数据提取

目前，利用遥感影像提取房屋数据包括人工目视判读和计算机自动解译2 种方式。其中，人工目视判读精度较高，但费时费力，效率较低。计算机自动解译效率高，但易受周边背景的影响（赵春霞等，2004）。针对上述方式优缺点，本文采用eCognition 软件处理与人工目视判读提取无人机影像图的房屋信息（刘贾贾等，2019）。

本文以大疆无人机为测航平台，以Pix4D 软件为地理空间航片拼接工具，对万全区村镇进行抽样拍摄。所用无人机遥感影像具有高空间分辨率特点，可使地物几何结构和纹理信息显示更清晰。从拍摄影像图来看，研究区内房屋多为红瓦、青瓦房顶，且房顶多为尖顶，以屋脊为界，向两侧倾斜。无人机遥感影像中的房屋屋顶纹理清晰，形状规则，颜色与耕地、林地等有显著区别（胡勇等，2016）。图3所示为无人机拍摄的调查点影像图，主要技术流程如图4 所示。本文以赵家梁村为例进行详细介绍，由于无人机影像（像素分辨率约为0.12 m）与91 卫图助手遥感影像（像素分辨率为0.11 m）分辨率接近，其余村镇采用91 卫图助手遥感影像进行识别，识别步骤与无人机影像一致。

图3 无人机拍摄赵家梁村影像图Fig. 3 Image of zhaojialiang village taken by UAV

图4 无人机拍摄房屋侧面照片Fig. 4 Side photos of houses taken by UAV

影像处理流程如图5 所示，为提高正确率，采用eCognition 软件首先对遥感影像进行预处理，分别选择30、40、50 尺度进行分割。将3 种尺度的分割效果进行对比，发现分割尺度为40 时，建筑物能较好地同其他地物进行区分，提取建筑物较完整，未出现“粘连”现象。由此得出分割尺度为40，形状指数为0.4，紧密度为0.5。利用特征库的阈值分类法将建筑物与非建筑物进行区分，与调查结果对比修正后提取影像图房屋面积。

图5 影像处理流程Fig. 5 Flow chart of image processing technology

根据实地调查情况，张家口农村地区土木、砖木结构房屋普遍占地面积为40～60 m2，屋顶多为灰瓦小尖顶，颜色较暗沉，房屋进深多为3～4 m。砖混结构房屋占地面积为60～90 m2，屋顶多为灰瓦屋顶，颜色较鲜艳，房屋进深多为4～4.6 m。根据以上特征，采用样本邻近分类方法，将房屋类型识别为2 类，分别为土木、砖木结构房屋和砖混结构房屋。

通过人工目视校核和编辑进一步识别房屋类型，对照无人机拍摄的实地房屋侧面照片（图4）与无人机影像图（图3）进行对比匹配，最终给出房屋结构和层数提取结果。通过对图3 进行辨识，将房屋类型划分为4 类（图6），包括砌体结构房屋，工业厂房、公共空旷房屋，砖木结构房屋，土、木、石结构房屋。对于一般性居民地，可不考虑工业厂房和公共空旷房屋面积。表1 为计算得出的3 类房屋面积，赵家梁村砌体结构房屋面积为55 478.99 m2，砖木结构房屋面积为21 807.27 m2，土、木、石结构房屋面积为19 233.63 m2。

表1 22 个调查点房屋信息Table 1 House area summary o for 22 survey points

图6 赵家梁村房屋信息提取结果Fig. 6 Extraction results of house information in zhaojialiang Village

通过无人机遥感影像与91 卫图助手遥感影像对22 个调查点分别进行房屋结构与面积的提取，并结合实际调查和无人机拍摄房屋侧面照片，修正22 个调查点提取的房屋信息数据。因遥感影像无法判断砌体结构房屋是否设防，故砌体结构房屋中设防房屋与未设防房屋占比通过实地调查得到。

4 房屋抗震能力分析

获取房屋结构类型比例是分析房屋抗震能力的重要环节。通过建立遥感影像获取不同类型建筑物面积，得到研究区内建筑物结构比例，结合抽样调查得到的各影响因素类型和分布比例关系，得出不同时代建筑物面积所占比例及层数所占比例。利用样本建筑物房屋类型抗震能力指数，分析万全区建筑物综合房屋类型抗震能力指数。

4.1 影响建（构）筑物抗震能力的因素

无数震害实例表明，提高建（构）筑物抗震能力是减轻震害最有效的途径。因此，提高建（构）筑物抗震能力是防震减灾重要环节（张风华等，2004；刘莉，2009）。归纳震害经验，同时考虑建（构）筑物抗震能力影响因素和数据易获取性，将房屋抗震设防水平、结构类型、建造年代和层数作为抗震能力指标。

4.2 划分建（构）筑物抗震能力指数等级

根据《建（构）筑物地震破坏等级划分》（GB/T 24335−2009）（中华人民共和国国家质量监督检验检疫总局等，2009），将建（构）筑物地震破坏分为5个等级，即基本完好、轻微破坏、中等破坏、严重破坏、毁坏。将建（构）筑物抗震能力量化为抗震能力指数，其代表建（构）筑物抵御地震灾害的能力，取值范围为0～1，数值越大表明抗震能力越好，越小表明抗震能力越差（谢礼立，2006），基本完好、轻微破坏、中等破坏、严重破坏、毁坏对应的抗震能力指数分别为1.0、0.8、0.6、0.4、0.2。

4.3 建（构）筑物抗震能力指数的确定

《建筑抗震设计规范》（GB 50011−2010）（中华人民共和国住房和城乡建设部等，2010）规定，地震区建筑必须按规定设防，因此不同的基本烈度地区设防标准不同，震害矩阵也不同。本文依据《中国大陆地震灾害损失评估汇编》收集到的1989−2016 年华北片区成灾地震的震害记录和相关数据，参考河北省地震应急基础数据库震害矩阵，修改得到适用于本研究区的震害矩阵。根据《地震现场工作 第4 部分：灾害直接损失评估》（GB/T 18208.4−2011）（中华人民共和国国家质量监督检验检疫总局等，2011）将房屋类别划分为钢结构房屋，钢筋混凝土房屋，砌体房屋，砖木房屋，土、木、石结构房屋，工业厂房和公共空旷房屋。鉴于张家口地区农村房屋多为砌体房屋、砖木、土木结构房屋，可将房屋类别划分为A、B、C、D 类，其中A 类为土、木、石结构房屋，包括土墙木屋架的土坯房、砖柱土坯房、土坯窑洞、碎石砌筑房屋；B 类为砖木房屋，包括砖墙、木房架的多层砖木结构、砖木平房；C 类为未设防的砌体房屋，即未经设防的砖砌体房屋；D 类为设防的砌体房屋，即按照Ⅶ度设防的砖砌体房屋，选用的震害矩阵如表2～5 所示。

表2 A 类房屋震害矩阵（%）Table 2 Earthquake damage matrix of A type buildings（Unit:%）

表3 B 类房屋震害矩阵（%）Table 3 Earthquake damage matrix of B type buildings（Unit:%）

表4 C 类房屋震害矩阵（%）Table 4 Earthquake damage matrix of C type buildings（Unit:%）

表 5 D 类房屋震害矩阵（%）Table 5 Earthquake damage matrix of D type buildings（Unit:%）

因遥感影像无法区别设防与未设防的砌体房屋，所以通过实地调查对利用遥感影像提取到的砌体房屋进行分类，即分为C 类或D 类。一般来说，设防的砌体结构房屋所占比例高的地区建（构）筑物抗震能力强。不同地区建（构）筑物抗震能力在很大程度上取决于A、B、C、D 类房屋面积比例。因此，本文采用谢礼立（2006）确定的抗震能力指数公式计算不同地区建（构）筑物震害矩阵：

表6 建造年代修正参考值Table 6 Reference value of building age correction

表7 层数修正参考值Table 7 Corrected reference value of layers

基于上述分析，可以给出某地区建（构）筑物抗震能力指数计算公式为：

式中：a1、a2、a3、a4分别表示1979 年以前、1980−1989 年、1990−1999 年、2000−2019 年建筑物面积占比；b1、b2、b3分别表示平房、2～6 层、7 层以上建筑物面积占比。

建筑物抗震能力等级划分如表8 所示。

表8 建筑物抗震能力等级划分Table 8 Classification of seismic capacity of buildings

5 结果分析

本文选取不同类型的建筑物面积数据（表1），得到研究区内房屋结构类型比例，并结合实际调查中的详查资料，对研究区域内房屋建造年代、层数进行统计与分类，得到22 个调查点房屋建造年代、层数比例，进而计算房屋抗震能力指数。表9 所示为22 个调查点不同结构类型、建造年代和层数对应的抗震能力指数。根据《建筑工程抗震设防分类标准》（GB 50223−2008）（中华人民共和国住房和城乡建设部，2008），万全区村镇抗震设防水平应按Ⅶ度设防。

表9 22 个调查点房屋抗震能力指数Table 9 Seismic capacity index of 22 survey points

根据式（1）～（4）和修正系数得出22 个调查点建筑物抗震能力指数值，如图7 所示。

图7 22 个调查点建筑物抗震能力指数对比图Fig. 7 Comparison of seismic capacity index of buildings in 22 survey points

采用平均数算法得到万全区房屋抗震设防水平，不同结构类型、建造年代和层数对应的抗震能力指数如表10 所示。

表10 万全区房屋抗震能力指数Table 10 Seismic capacity index of buildings in Wanquan district

通过式（1）计算万全区建筑物震害矩阵，从而得到建（构）筑物平均抗震能力指数为：

万全区建（构）筑物抗震能力指数为：

可得万全区遭遇VI～Ⅹ度地震时，抗震能力指数为：[0.684 7，0.542 0，0.403 7，0.307 6，0.199 4]。

6 分析与讨论

本文利用无人机遥感技术获得高分辨率遥感影像图，通过eCognition 软件处理初步得到建筑物面积及分类，结合无人机拍摄的侧面照片，通过人工辅助进行识别与修正，并与实地调查相结合，使提取的村镇房屋信息数据更准确，为地震灾害风险评估和地震风险普查工作提供参考。本文在建筑物抗震能力指数计算中考虑了设防与未设防的砌体房屋情况，得到了较精准的建筑物抗震能力情况。但目前无法通过遥感影像识别砌体房屋是否设防，需通过实地调查确定。本研究村镇房屋结构较单一，需对进一步提高城市房屋识别正确率进行研究与探索。

由万全区村镇建筑物抗震能力指数计算结果可知，万全区房屋整体抗震能力较弱，以未设防的砌体结构和砖木、土木结构房屋为主。抗震能力指数为0.6～0.7，基本可抵御VI 度地震，当遭遇VI 度地震时，万全区村镇建筑物抗震能力为中等，除个别土木结构房屋集中的地区外，房屋不至于大面积损坏，也不会造成巨大损失。当遭遇Ⅶ度地震时，万全区村镇建筑物抗震能力较差。

通过对万全区房屋进行现场调查，发现交通较发达的地方砖木结构房屋较多，偏远难行的地方土木结构房屋较多。老旧的土窑洞、土木结构房屋由于建设年代较早，抗震性能普遍较差。砖木结构房屋普遍未采取抗震措施，木屋架多为硬山搁檩形式，与墙体无可靠拉结，整体性较差。同时，普遍存在门窗开洞过大、窗间墙过窄的现象，造成墙体抗剪能力减弱、房屋刚度不均匀，对房屋整体抗震不利。房屋多采用石砌基础，基础埋深为30～100 cm，部分基础砌筑不牢，影响房屋抗震性能。普遍存在附属设施安全隐患，其中出屋面烟囱、门楼、围墙、挑檐等震后均可能出现掉落、倒塌，进而导致伤人，应予以重视。

地震烈度为VI 度时，万全区房屋基本完好。地震烈度为Ⅶ度以上时，房屋会发生大面积损坏，由此看来，张家口地区农村房屋抗震性能仍不足，地震风险偏大，防震减灾形势严峻，应加大农村地区抗震减灾力度。