Analysis of Seismic Capacity of Buildings Based on Remote Sensing ImageTaking Zhangjiakou Wanquan District as An Example
-
摘要: 据震害统计，房屋抗震能力是影响灾害的主要因素，抗震能力一般由房屋抗震设防水平、结构类型、建造年代和房屋层数等因素决定，通过对房屋抗震能力的综合评定，以便采取更有针对性的地震灾害对策。本文通过遥感影像实现房屋基本信息的快速提取，提高房屋结构类型获取的便捷性，具有较高的准确度和可靠度，并利用抽样调查得到张家口万全区房屋建造年代和层数分布特点，结合当地抗震设防水平，建立房屋抗震能力指数指标体系，阐述房屋抗震能力现状，为地震灾害损失评估、风险普查、风险区划等工作提供参考。Abstract: According to the statistics of earthquake damage, the anti-seismic capacity of buildings is the main factor affecting the disaster. The anti-seismic capacity is generally determined by the level of anti-seismic fortification, the type of structure, the construction age and the number of floors of buildings. Through the comprehensive evaluation of the anti-seismic capacity of buildings, we can take more targeted earthquake disaster countermeasures. In this paper, the basic information of buildings is extracted quickly by remote sensing image, which improves the convenience of obtaining the type of building structure, and has high accuracy and reliability. The distribution characteristics of building construction age and story number are obtained by sampling survey in Wanquan District of Zhangjiakou. Combined with the local seismic fortification level, the index system of building seismic capacity is established, and the current situation of building seismic capacity is expounded, It is of great significance for earthquake disaster loss assessment, risk census and risk zoning.
-
引言
从近十几年国内外地震震害来看，房屋破坏是造成人员伤亡和经济损失的最主要原因。房屋建筑资料等信息一直是灾情快速评估的重要基础，尤其是房屋结构的合理分类，决定着评估结果的准确性（聂高众等，2002；帅向华等，2009；姜立新等，2012）。
近年来，随着遥感影像技术的迅速发展，高分卫星遥感分辨率提升至亚米级，为房屋基本信息的快速获取提供了途径。王晓青等（2004）和窦爱霞等（2012）通过遥感技术获取灾区震后遥感影像，并提出了遥感震害指数与定量分析模型，成功应用于汶川、玉树地震等的震害快速定量化评估中。郭建兴等（2020）提出了利用遥感影像结合居民地建筑物数据评估地震烈度的方法。马建等（2020）通过获取震后高分辨率无人机影像，对房屋进行结构分类和震害解译。李昌珑等（2019）利用三维图像模式识别各类建筑物的数量分布，并建立不同建筑物易损性模型。汶川、玉树地震等震害多次暴露出我国村镇房屋抗震能力较脆弱，快速、全面、准确掌握村镇房屋抗震能力一直是研究的重点方向（周强等，2016；李皓等，2018）。
目前，主要利用遥感影像进行震后房屋震害识别，较少用于房屋抗震能力的研究。因各地房屋结构复杂、形式多样，实际调查难度大，通过遥感技术快速获取房屋结构类型等基本信息，对分析房屋抗震能力具有现实意义。刘贾贾等(2019，2021)通过遥感影像对城镇和农村建筑物进行分类，并根据建筑物分类结果对震害风险进行初步划分。
本文通过无人机航拍和遥感影像获取房屋基本信息，通过人机交互方式提取村镇房屋结构类型，结合实地调查统计得到其他信息，建立当地村镇房屋抗震能力指数体系，分析村镇房屋抗震能力。利用遥感影像实现房屋基本信息的快速提取，提高房屋信息获取的便捷性，且具有较高的准确性和可靠性，对于地震灾害损失评估、风险普查、风险区划等工作具有重要意义。
1. 研究区域
张家口市位于河北省西北部，京晋冀蒙交界地带，是“首都圈”地区地震重点监视防御区和晋冀蒙交界地震重点危险区，属于张家口—渤海地震构造带与汾渭地震构造带交汇地区，具备发生中等破坏性地震的条件和背景（李书进等，2010；张肖等，2019；茅远哲等，2019）。
根据张家口市所处的地理位置、经济发展情况，部分地区新农村建设发展迅速，按抗震设防建造的房屋抗震能力存在差异。而部分地区发展较慢，存在无抗震能力的旧房，如土木、砖木、土窑结构房屋等（刘龙等，2019）。万全区位于张家口市西北部，地理位置险要，受地质构造控制，地形北高南低，东西向为西高东低，是川、丘、山并存的地区（图1），本文以该区域作为研究对象。
2. 现场调查
调查数据作为万全区房屋抗震性能分析的重要基础数据，为后期房屋抗震能力分析奠定基础。基于随机均匀分布原则，对于研究区每个乡镇，选取不少于2个调查点，调查范围确定了22个村镇（图2），包括新河口村、代家房村、西湾村、武家庄村、新开口村、望虎台村、刘虎庄村、板山村、阳门堡村、旧羊屯村、吴家窑村、上营房村、西柳林村、洗马林镇、万全镇、沙家庄村、陈家沟村、高庙堡村、羊窖沟村、岸庄屯村、安家堡村、赵家梁村。采取拍取房屋照片、填写房屋信息、无人机测航调查点等方式开展调查，普查主要侧重于结构类型、抗震设防水平、构造措施等基本情况，详查是在普查的基础上，对调查区域内房屋进行随机抽样，更加详细地调查建造年代、层数、用途等信息（周强等，2016）。
图 2 万全区地貌分区图与调查点分布Figure 2. Geomorphic zoning map and distribution of survey points in Wanquan District3. 建筑物数据提取
目前，利用遥感影像提取房屋数据包括人工目视判读和计算机自动解译2种方式。其中，人工目视判读精度较高，但费时费力，效率较低。计算机自动解译效率高，但易受周边背景的影响（赵春霞等，2004）。针对上述方式优缺点，本文采用eCognition软件处理与人工目视判读提取无人机影像图的房屋信息（刘贾贾等，2019）。
本文以大疆无人机为测航平台，以Pix4D软件为地理空间航片拼接工具，对万全区村镇进行抽样拍摄。所用无人机遥感影像具有高空间分辨率特点，可使地物几何结构和纹理信息显示更清晰。从拍摄影像图来看，研究区内房屋多为红瓦、青瓦房顶，且房顶多为尖顶，以屋脊为界，向两侧倾斜。无人机遥感影像中的房屋屋顶纹理清晰，形状规则，颜色与耕地、林地等有显著区别（胡勇等，2016）。图3所示为无人机拍摄的调查点影像图，主要技术流程如图4所示。本文以赵家梁村为例进行详细介绍，由于无人机影像（像素分辨率约为0.12 m）与91卫图助手遥感影像（像素分辨率为0.11 m）分辨率接近，其余村镇采用91卫图助手遥感影像进行识别，识别步骤与无人机影像一致。
影像处理流程如图5所示，为提高正确率，采用eCognition软件首先对遥感影像进行预处理，分别选择30、40、50尺度进行分割。将3种尺度的分割效果进行对比，发现分割尺度为40时，建筑物能较好地同其他地物进行区分，提取建筑物较完整，未出现“粘连”现象。由此得出分割尺度为40，形状指数为0.4，紧密度为0.5。利用特征库的阈值分类法将建筑物与非建筑物进行区分，与调查结果对比修正后提取影像图房屋面积。
根据实地调查情况，张家口农村地区土木、砖木结构房屋普遍占地面积为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。
图 6 赵家梁村房屋信息提取结果Figure 6. Extraction results of house information in zhaojialiang Village通过无人机遥感影像与91卫图助手遥感影像对22个调查点分别进行房屋结构与面积的提取，并结合实际调查和无人机拍摄房屋侧面照片，修正22个调查点提取的房屋信息数据。因遥感影像无法判断砌体结构房屋是否设防，故砌体结构房屋中设防房屋与未设防房屋占比通过实地调查得到。
表 1 22个调查点房屋信息Table 1. House area summary o for 22 survey points序号 调查点 砌体结构房屋
面积/m2砖木结构
房屋面积/m²土、木、石结构房屋
面积/m2序号 调查点 砌体结构
房屋面积/m2砖木结构
房屋面积/m2土、木、石结构房屋
面积/m21 新河口村 17 181 40 492 14 168 12 上营房村 49 995 16 078 7 481 2 代家房村 7 335 15 503 3 693 13 西柳林村 2 044 5 492 821 3 西湾村 1 727 3 121 709 14 洗马林镇 263 798 98 653 54 821 4 武家庄村 48 671 111 185 31 716 15 万全镇 36 337 13 895 3 686 5 新开口村 70 999 44 171 13 381 16 沙家庄村 69 798 21 931 8 306 6 望虎台村 7 107 7 934 3 042 17 陈家沟村 20 917 43 403 15 590 7 刘虎庄村 15 150 39 950 13 182 18 高庙堡村 212 719 84 412 50 936 8 板山村 19 045 8 172 9 802 19 羊窖沟村 48 847 44 772 9 339 9 阳门堡村 124 920 130 124 42 232 20 岸庄屯村 84 151 109 636 21 021 10 旧羊屯村 124 071 127 114 43 988 21 安家堡村 277 443 92 840 50 190 11 吴家窑村 40 342 13 275 5 845 22 赵家梁村 55 479 21 807 19 234 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:%）破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 32.00 16.50 7.00 2.50 0 轻微破坏 26.50 18.50 12.00 8.50 1.50 中等破坏 22.50 20.00 16.50 14.00 7.50 严重破坏 16.50 26.00 27.00 25.00 18.50 毁坏 2.50 19.00 37.50 50.00 72.50 表 3 B类房屋震害矩阵（%）Table 3. Earthquake damage matrix of B type buildings（Unit:%）破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 47.35 26.85 11.53 7.45 1.98 轻微破坏 27.09 21.01 15.94 10.36 4.48 中等破坏 15.83 21.88 22.43 17.29 11.45 严重破坏 7.84 20.84 29.90 25.97 17.28 毁坏 1.89 9.42 20.20 38.93 64.81 表 4 C类房屋震害矩阵（%）Table 4. Earthquake damage matrix of C type buildings（Unit:%）破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 56.11 42.49 25.59 13.47 2.85 轻微破坏 27.76 23.75 21.43 16.76 7.30 中等破坏 10.89 16.65 21.58 23.16 17.40 严重破坏 4.29 12.38 20.36 22.15 24.86 毁坏 0.96 4.73 11.04 24.46 47.59 表 5 D类房屋震害矩阵（%）Table 5. Earthquake damage matrix of D type buildings（Unit:%）破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 69.29 67.24 53.28 33.65 10.96 轻微破坏 24.91 21.71 23.33 23.41 16.14 中等破坏 4.62 7.56 14.40 22.22 25.53 严重破坏 1.06 2.80 6.67 13.67 26.27 毁坏 0.12 0.69 2.32 7.05 21.10 因遥感影像无法区别设防与未设防的砌体房屋，所以通过实地调查对利用遥感影像提取到的砌体房屋进行分类，即分为C类或D类。一般来说，设防的砌体结构房屋所占比例高的地区建（构）筑物抗震能力强。不同地区建（构）筑物抗震能力在很大程度上取决于A、B、C、D类房屋面积比例。因此，本文采用谢礼立（2006）确定的抗震能力指数公式计算不同地区建（构）筑物震害矩阵：
$$ P\left[ {{D_i}/J,I} \right] = \sum\limits_{s = 1}^4 {{\omega _s}} {P_s}\left[ {{D_i}/J,I} \right] $$ (1) 式中：
$P\left[ {{{{D_i}} / {J,I}}} \right]$ 表示设防烈度为J的建筑物遭遇I烈度地震时的震害矩阵；${P_s}\left[ {{{{D_i}} / {J,I}}} \right]$ 表示s类结构设防烈度为J的建筑物遭遇I烈度地震时的震害矩阵；${\omega _s}$ 表示s类结构建筑物面积比例。不同地区建（构）筑物平均抗震能力计算公式如下：
$$ IL1(J,I) = {\boldsymbol{K}}\times P\left( {{D_i}/J,I} \right) $$ (2) 式中：
$IL1\left({J,I} \right)$ 表示城镇设防烈度为J的建（构）筑物遭受I烈度地震时的平均抗震能力指数；K表示抗震能力等级矩阵，为$ \{1,0.8,0.6,0.4,0.2\}$ 。通过收集的单体样本数据库考虑建造年代和层数对抗震能力的影响。考虑某年代对抗震能力指数的影响时，将单体样本数据库中去除这个年代后得到震害矩阵，比较原震害矩阵抗震能力指数和去除该年代后得出的震害矩阵抗震能力指数的差别，得出该年代对抗震能力的修正值，公式如下：
$$ {K_j} = \frac{{{F_{\rm{V}}}({{ {\rm{Newly}} }})}}{{{F_{\rm{V}}}({{ {\rm{Original}} }})}} $$ (3) 式中：
${{F_{\rm{V}}}({{ {\rm{Newly}} }})}$ 表示在单体样本数据库中去除j因素（年代）后得出的抗震能力指数；${{F_{\rm{V}}}({{ {\rm{Original}} }})}$ 表示根据原有震害矩阵得出的抗震能力指数；Kj表示第j个因素对抗震能力的影响值。根据式（3），得到建造年代和层数对抗震能力的修正值，如表6、7所示。
表 6 建造年代修正参考值Table 6. Reference value of building age correction建造年代 1979年以前 1980—1989年 1990—1999年 2000—2019年 修正值 0.700 0.800 0.934 1.000 表 7 层数修正参考值Table 7. Corrected reference value of layers层数 平房 2～6层 7层以上 修正值 0.822 0.863 1.000 基于上述分析，可以给出某地区建（构）筑物抗震能力指数计算公式为：
$$ IL = IL1 \times (0.7 \times a1\% + 0.8 \times a2\% + 0.934 \times a3\% + 1 \times a4\% )\\ \times (0.822 \times b1\% + 0.863 \times b2\% + 1.0 \times b3\% ) $$ (4) 式中：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建筑物抗震能力分级 优良 中等 差 范围 $ 0.8 \leqslant IL < 1$ $ 0.6 \leqslant IL < 0.8$ $ IL < 0.6$ 5. 结果分析
本文选取不同类型的建筑物面积数据（表1），得到研究区内房屋结构类型比例，并结合实际调查中的详查资料，对研究区域内房屋建造年代、层数进行统计与分类，得到22个调查点房屋建造年代、层数比例，进而计算房屋抗震能力指数。表9所示为22个调查点不同结构类型、建造年代和层数对应的抗震能力指数。根据《建筑工程抗震设防分类标准》（GB 50223—2008）（中华人民共和国住房和城乡建设部，2008），万全区村镇抗震设防水平应按Ⅶ度设防。
表 9 22个调查点房屋抗震能力指数Table 9. Seismic capacity index of 22 survey points调查点 新河口村 代家房村 西湾村 武家庄村 新开口村 望虎台村 刘虎庄村 板山村 阳门堡村 旧羊屯村 吴家窑村 结构类型 A类房屋 0.1851 0.3616 0.1348 0.1590 0.7604 0.1640 0.1855 0.2562 0.3438 0.1460 0.0940 B类房屋 0.5416 0.5368 0.6266 0.5574 0.1934 0.4278 0.5622 0.2136 0.5814 0.4219 0.2135 C类房屋 0.2632 0.1016 0.2386 0.2694 0.0462 0.2543 0.2523 0.5302 0.0515 0.3496 0.5876 D类房屋 0.0101 0 0 0.0142 0 0.1539 0 0 0.0233 0.0825 0.1049 建造年代 1979年以前 0.1234 0.1163 0.1047 0.1211 0.1011 0.1798 0.1221 0.1246 0.3563 0.1452 0.1254 1980—1989年 0.3164 0.3397 0.3451 0.3385 0.1384 0.3269 0.3199 0.3144 0.2364 0.3749 0.3154 1990—1999年 0.3423 0.3543 0.2456 0.3521 0.5916 0.3164 0.3384 0.3287 0.2963 0.3465 0.3336 2000—2019年 0.2179 0.1897 0.3046 0.1885 0.1689 0.1769 0.2196 0.2323 0.111 0.1334 0.2256 层数 平房 0.9659 0.9975 0.9834 0.7489 0.9746 0.8461 0.9644 0.9784 0.8831 0.8313 0.7290 3层以下楼房 0.0341 0.0025 0.0166 0.1746 0.0254 0.1539 0.0356 0.0216 0.1169 0.1569 0.1846 4层以上楼房 0 0 0 0.0765 0 0 0 0 0 0 0 调查点 上营房村 西柳林村 洗马林镇 万全镇 沙家庄村 陈家沟村 高庙堡村 羊窖沟村 岸庄屯村 安家堡村 赵家梁村 结构类型 A类房屋 0.0952 0.0969 0.1287 0.0674 0.0809 0.1908 0.1447 0.0894 0.0956 0.1165 0.1834 B类房屋 0.2046 0.6482 0.2316 0.2541 0.2136 0.5312 0.2398 0.4286 0.4986 0.2155 0.2080 C类房屋 0.5966 0.2549 0.5236 0.5461 0.6189 0.2780 0.5123 0.4131 0.3264 0.5784 0.6086 D类房屋 0.1036 0 0.1161 0.1324 0.0865 0 0.1032 0.0689 0.0794 0.0896 0 建造年代 1979年以前 0.1033 0.1661 0.093 0.0454 0.0765 0.1215 0.0821 0.0965 0.1569 0.1023 0.1651 1980—1989年 0.2985 0.3216 0.3469 0.2947 0.3621 0.3946 0.3416 0.3765 0.3263 0.2989 0.3547 1990—1999年 0.3045 0.3611 0.3136 0.3399 0.3611 0.3108 0.3656 0.3248 0.3625 0.3516 0.2812 2000—2019年 0.2937 0.1512 0.2465 0.3200 0.2003 0.1731 0.2107 0.2022 0.1543 0.2472 0.1990 层数 平房 0.8558 0.9347 0.6725 0.6214 0.8764 0.9036 0.6652 0.8996 0.8088 0.7102 0.9658 3层以下楼房 0.1442 0.0653 0.2036 0.2158 0.1236 0.0964 0.1879 0.1004 0.1123 0.1695 0.0342 4层以上楼房 0 0 0.1239 0.1628 0 0 0.1469 0 0.0789 0.1203 0 根据式（1）～（4）和修正系数得出22个调查点建筑物抗震能力指数值，如图7所示。
图 7 22个调查点建筑物抗震能力指数对比图Figure 7. Comparison of seismic capacity index of buildings in 22 survey points采用平均数算法得到万全区房屋抗震设防水平，不同结构类型、建造年代和层数对应的抗震能力指数如表10所示。
表 10 万全区房屋抗震能力指数Table 10. Seismic capacity index of buildings in Wanquan district结构类型 A类其他房屋 0.135 8 B类旧式房屋 0.380 7 C类砌体房屋 0.411 9 D类砌体房屋 0.062 2 建造年代 1979年以前 0.118 6 1980—1989年 0.343 1 1990—1999年 0.332 8 2000—2019年 0.239 4 层数 平房 0.827 5 3层以下楼房 0.123 7 4层以上楼房 0.048 8 通过式（1）计算万全区建筑物震害矩阵，从而得到建（构）筑物平均抗震能力指数为：
$$\begin{split} & IL 1(J,I) = K\times P\left( {{D_i}/J,I} \right)=\\ &[0.902 \;6,0.714 \;4,0.532 \;2,0.405 \;5,0.262 \;9] \end{split} $$ (5) 万全区建（构）筑物抗震能力指数为：
$$ \begin{split} & I L=I L 1 \times(0.7 \times a 1 \%+0.8 \times a 2 \%+0.934 \times a 3 \%+1 \times a 4 \%) \times\\ & (0.822 \times b 1 \%+0.863 \times b 2 \%+1.0 \times b 3 \%)=\\ & [0.684\;7,0.542\;0,0.403\;7,0.307\;6,0.199\;4] \end{split} $$ (6) 可得万全区遭遇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度时，万全区房屋基本完好。地震烈度为Ⅶ度以上时，房屋会发生大面积损坏，由此看来，张家口地区农村房屋抗震性能仍不足，地震风险偏大，防震减灾形势严峻，应加大农村地区抗震减灾力度。
-
图 2 万全区地貌分区图与调查点分布
Figure 2. Geomorphic zoning map and distribution of survey points in Wanquan District
图 6 赵家梁村房屋信息提取结果
Figure 6. Extraction results of house information in zhaojialiang Village
图 7 22个调查点建筑物抗震能力指数对比图
Figure 7. Comparison of seismic capacity index of buildings in 22 survey points
表 1 22个调查点房屋信息
Table 1. House area summary o for 22 survey points
序号 调查点 砌体结构房屋
面积/m2砖木结构
房屋面积/m²土、木、石结构房屋
面积/m2序号 调查点 砌体结构
房屋面积/m2砖木结构
房屋面积/m2土、木、石结构房屋
面积/m21 新河口村 17 181 40 492 14 168 12 上营房村 49 995 16 078 7 481 2 代家房村 7 335 15 503 3 693 13 西柳林村 2 044 5 492 821 3 西湾村 1 727 3 121 709 14 洗马林镇 263 798 98 653 54 821 4 武家庄村 48 671 111 185 31 716 15 万全镇 36 337 13 895 3 686 5 新开口村 70 999 44 171 13 381 16 沙家庄村 69 798 21 931 8 306 6 望虎台村 7 107 7 934 3 042 17 陈家沟村 20 917 43 403 15 590 7 刘虎庄村 15 150 39 950 13 182 18 高庙堡村 212 719 84 412 50 936 8 板山村 19 045 8 172 9 802 19 羊窖沟村 48 847 44 772 9 339 9 阳门堡村 124 920 130 124 42 232 20 岸庄屯村 84 151 109 636 21 021 10 旧羊屯村 124 071 127 114 43 988 21 安家堡村 277 443 92 840 50 190 11 吴家窑村 40 342 13 275 5 845 22 赵家梁村 55 479 21 807 19 234 
下载: 导出CSV
表 2 A类房屋震害矩阵（%）
Table 2. Earthquake damage matrix of A type buildings（Unit:%）
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 32.00 16.50 7.00 2.50 0 轻微破坏 26.50 18.50 12.00 8.50 1.50 中等破坏 22.50 20.00 16.50 14.00 7.50 严重破坏 16.50 26.00 27.00 25.00 18.50 毁坏 2.50 19.00 37.50 50.00 72.50
下载: 导出CSV
表 3 B类房屋震害矩阵（%）
Table 3. Earthquake damage matrix of B type buildings（Unit:%）
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 47.35 26.85 11.53 7.45 1.98 轻微破坏 27.09 21.01 15.94 10.36 4.48 中等破坏 15.83 21.88 22.43 17.29 11.45 严重破坏 7.84 20.84 29.90 25.97 17.28 毁坏 1.89 9.42 20.20 38.93 64.81
下载: 导出CSV
表 4 C类房屋震害矩阵（%）
Table 4. Earthquake damage matrix of C type buildings（Unit:%）
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 56.11 42.49 25.59 13.47 2.85 轻微破坏 27.76 23.75 21.43 16.76 7.30 中等破坏 10.89 16.65 21.58 23.16 17.40 严重破坏 4.29 12.38 20.36 22.15 24.86 毁坏 0.96 4.73 11.04 24.46 47.59
下载: 导出CSV
表 5 D类房屋震害矩阵（%）
Table 5. Earthquake damage matrix of D type buildings（Unit:%）
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 69.29 67.24 53.28 33.65 10.96 轻微破坏 24.91 21.71 23.33 23.41 16.14 中等破坏 4.62 7.56 14.40 22.22 25.53 严重破坏 1.06 2.80 6.67 13.67 26.27 毁坏 0.12 0.69 2.32 7.05 21.10
下载: 导出CSV
表 6 建造年代修正参考值
Table 6. Reference value of building age correction
建造年代 1979年以前 1980—1989年 1990—1999年 2000—2019年 修正值 0.700 0.800 0.934 1.000
下载: 导出CSV
表 8 建筑物抗震能力等级划分
Table 8. Classification of seismic capacity of buildings
建筑物抗震能力分级 优良 中等 差 范围 $ 0.8 \leqslant IL < 1$ $ 0.6 \leqslant IL < 0.8$ $ IL < 0.6$
下载: 导出CSV
表 9 22个调查点房屋抗震能力指数
Table 9. Seismic capacity index of 22 survey points
调查点 新河口村 代家房村 西湾村 武家庄村 新开口村 望虎台村 刘虎庄村 板山村 阳门堡村 旧羊屯村 吴家窑村 结构类型 A类房屋 0.1851 0.3616 0.1348 0.1590 0.7604 0.1640 0.1855 0.2562 0.3438 0.1460 0.0940 B类房屋 0.5416 0.5368 0.6266 0.5574 0.1934 0.4278 0.5622 0.2136 0.5814 0.4219 0.2135 C类房屋 0.2632 0.1016 0.2386 0.2694 0.0462 0.2543 0.2523 0.5302 0.0515 0.3496 0.5876 D类房屋 0.0101 0 0 0.0142 0 0.1539 0 0 0.0233 0.0825 0.1049 建造年代 1979年以前 0.1234 0.1163 0.1047 0.1211 0.1011 0.1798 0.1221 0.1246 0.3563 0.1452 0.1254 1980—1989年 0.3164 0.3397 0.3451 0.3385 0.1384 0.3269 0.3199 0.3144 0.2364 0.3749 0.3154 1990—1999年 0.3423 0.3543 0.2456 0.3521 0.5916 0.3164 0.3384 0.3287 0.2963 0.3465 0.3336 2000—2019年 0.2179 0.1897 0.3046 0.1885 0.1689 0.1769 0.2196 0.2323 0.111 0.1334 0.2256 层数 平房 0.9659 0.9975 0.9834 0.7489 0.9746 0.8461 0.9644 0.9784 0.8831 0.8313 0.7290 3层以下楼房 0.0341 0.0025 0.0166 0.1746 0.0254 0.1539 0.0356 0.0216 0.1169 0.1569 0.1846 4层以上楼房 0 0 0 0.0765 0 0 0 0 0 0 0 调查点 上营房村 西柳林村 洗马林镇 万全镇 沙家庄村 陈家沟村 高庙堡村 羊窖沟村 岸庄屯村 安家堡村 赵家梁村 结构类型 A类房屋 0.0952 0.0969 0.1287 0.0674 0.0809 0.1908 0.1447 0.0894 0.0956 0.1165 0.1834 B类房屋 0.2046 0.6482 0.2316 0.2541 0.2136 0.5312 0.2398 0.4286 0.4986 0.2155 0.2080 C类房屋 0.5966 0.2549 0.5236 0.5461 0.6189 0.2780 0.5123 0.4131 0.3264 0.5784 0.6086 D类房屋 0.1036 0 0.1161 0.1324 0.0865 0 0.1032 0.0689 0.0794 0.0896 0 建造年代 1979年以前 0.1033 0.1661 0.093 0.0454 0.0765 0.1215 0.0821 0.0965 0.1569 0.1023 0.1651 1980—1989年 0.2985 0.3216 0.3469 0.2947 0.3621 0.3946 0.3416 0.3765 0.3263 0.2989 0.3547 1990—1999年 0.3045 0.3611 0.3136 0.3399 0.3611 0.3108 0.3656 0.3248 0.3625 0.3516 0.2812 2000—2019年 0.2937 0.1512 0.2465 0.3200 0.2003 0.1731 0.2107 0.2022 0.1543 0.2472 0.1990 层数 平房 0.8558 0.9347 0.6725 0.6214 0.8764 0.9036 0.6652 0.8996 0.8088 0.7102 0.9658 3层以下楼房 0.1442 0.0653 0.2036 0.2158 0.1236 0.0964 0.1879 0.1004 0.1123 0.1695 0.0342 4层以上楼房 0 0 0.1239 0.1628 0 0 0.1469 0 0.0789 0.1203 0
下载: 导出CSV
表 10 万全区房屋抗震能力指数
Table 10. Seismic capacity index of buildings in Wanquan district
结构类型 A类其他房屋 0.135 8 B类旧式房屋 0.380 7 C类砌体房屋 0.411 9 D类砌体房屋 0.062 2 建造年代 1979年以前 0.118 6 1980—1989年 0.343 1 1990—1999年 0.332 8 2000—2019年 0.239 4 层数 平房 0.827 5 3层以下楼房 0.123 7 4层以上楼房 0.048 8
下载: 导出CSV
-
[1] 窦爱霞, 王晓青, 丁香等, 2012. 遥感震害快速定量评估方法及其在玉树地震中的应用. 灾害学, 27(3): 75—80. doi: 10.3969/j.issn.1000-811X.2012.03.016Dou A. X., Wang X. Q., Ding X., et al., 2012. Quantitative methods of rapid earthquake damage assessment using remote sensing and its application in Yushu earthquake. Journal of Catastrophology, 27(3): 75—80. (in Chinese) doi: 10.3969/j.issn.1000-811X.2012.03.016 [2] 郭建兴, 张宇翔, 姬建中等, 2020. 利用居民地建筑物数据和高分遥感影像评估地震烈度的方法初探. 地震地质, 42(4): 968—980. doi: 10.3969/j.issn.0253-4967.2020.04.013Guo J. X., Zhang Y. X., Ji J. Z., et al., 2020. A preliminary study on the method of seismic intensity assessment based on residential building data and high resolution remote sensing images. Seismology and Geology, 42(4): 968—980. (in Chinese) doi: 10.3969/j.issn.0253-4967.2020.04.013 [3] 胡勇, 张孝成, 马泽忠等, 2016. 无人机遥感影像中农村房屋信息快速提取. 国土资源遥感, 28(3): 96—101.Hu Y., Zhang X. C., Ma Z. Z., et al., 2016. Rural residential area extraction from UAV remote sensing imagery. Remote Sensing for Land Resources, 28(3): 96—101. (in Chinese) [4] 姜立新, 帅向华, 聂高众等, 2012. 地震应急指挥协同技术平台设计研究. 震灾防御技术, 7(3): 294—302. doi: 10.3969/j.issn.1673-5722.2012.03.008Jiang L. X., Shuai X. H., Nie G. Z., et al., 2012. Study on the design of earthquake emergency command collaboration technology platform. Technology for Earthquake Disaster Prevention, 7(3): 294—302. (in Chinese) doi: 10.3969/j.issn.1673-5722.2012.03.008 [5] 李昌珑, 李宗超, 吕红山等, 2019. 基于三维图像模式识别的西藏东南部地震灾害损失风险评估. 地球物理学报, 62(1): 393—410. doi: 10.6038/cjg2019M0360Li C. L., Li Z. C., Lv H. S., et al., 2019. Seismic disaster loss risk assessment for southeastern Tibet based on 3D image pattern recognition. Chinese Journal of Geophysics, 62(1): 393—410. (in Chinese) doi: 10.6038/cjg2019M0360 [6] 李皓, 张合, 吕国军, 2018. 基于遥感影像的建筑数据构建研究. 震灾防御技术, 13(1): 168—176. doi: 10.11899/zzfy20180115Li H., Zhang H., Lv G. J., 2018. Construction of building data based on remote sensing images. Technology for Earthquake Disaster Prevention, 13(1): 168—176. (in Chinese) doi: 10.11899/zzfy20180115 [7] 李书进, 毛羚, 陶礼龙等, 2010. 湖北农村民居抗震性能调查与分析. 震灾防御技术, 5(1): 116—124. doi: 10.3969/j.issn.1673-5722.2010.01.014Li S. J., Mao L., Tao L. L., et al., 2010. Seismic performance investigation and analysis of rural buildings in Hubei province. Technology for Earthquake Disaster Prevention, 5(1): 116—124. (in Chinese) doi: 10.3969/j.issn.1673-5722.2010.01.014 [8] 刘贾贾, 刘志辉, 刘龙等, 2019. 基于遥感影像的农村建筑物分类. 华北地震科学, 37(4): 65—72. doi: 10.3969/j.issn.1003-1375.2019.04.012Liu J. J., Liu Z. H., Liu L., et al., 2019. Classification of rural buildings in Zhangjiakou area based on remote sensing images. North China Earthquake Sciences, 37(4): 65—72. (in Chinese) doi: 10.3969/j.issn.1003-1375.2019.04.012 [9] 刘贾贾, 刘志辉, 刘龙等, 2021. 基于遥感影像的城镇建筑物分类. 测绘与空间地理信息, 44(1): 130—133. doi: 10.3969/j.issn.1672-5867.2021.01.036Liu J. J., Liu Z. H., Liu L., et al., 2021. Classification of town buildings based on remote sensing images. Geomatics & Spatial Information Technology, 44(1): 130—133. (in Chinese) doi: 10.3969/j.issn.1672-5867.2021.01.036 [10] 刘莉, 2009. 城市防震减灾能力标定及可接受风险研究. 哈尔滨: 中国地震局工程力学研究所.Liu L., 2009. Calibration capability of the urban seismic prevention and disaster mitigation & research on acceptable risk level. Harbin: Institute of Engineering Mechanics, CEA. (in Chinese) [11] 刘龙, 刘志辉, 刘晓丹等, 2019. 张北地区农村房屋抗震性能分析. 华北地震科学, 37(2): 82—88. doi: 10.3969/j.issn.1003-1375.2019.02.013Liu L., Liu Z. H., Liu X. D., et al., 2019. Analysis of the seismic performance of rural houses in the areas of Zhangbei. North China Earthquake Sciences, 37(2): 82—88. (in Chinese) doi: 10.3969/j.issn.1003-1375.2019.02.013 [12] 马建, 常想德, 黄帅堂等, 2020. 无人机摄影技术在精河地震房屋震害定量评估中的应用. 震灾防御技术, 15(1): 208—215. doi: 10.11899/zzfy20200121Ma J., Chang X. D., Huang S. T., et al., 2020. Application of UAV photography technology in quantitative assessment of building damage of the Jinghe earthquake. Technology for Earthquake Disaster Prevention, 15(1): 208—215. (in Chinese) doi: 10.11899/zzfy20200121 [13] 茅远哲, 曹筠, 高晨等, 2019. 京西北地区地应变观测与小震震源机制解一致性研究. 中国地震, 35(4): 709—717. doi: 10.3969/j.issn.1001-4683.2019.04.011Mao Y. Z., Cao J., Gao C., et al., 2019. Consistency of ground strain observation and source mechanism solution of small earthquakes in northwest Beijing. Earthquake Research in China, 35(4): 709—717. (in Chinese) doi: 10.3969/j.issn.1001-4683.2019.04.011 [14] 聂高众, 陈建英, 李志强等, 2002. 地震应急基础数据库建设. 地震, 22(3): 105—112. doi: 10.3969/j.issn.1000-3274.2002.03.017Nie G. Z., Chen J. Y., Li Z. Q., et al., 2002. The construction of basic database for earthquake emergency response. Earthquake, 22(3): 105—112. (in Chinese) doi: 10.3969/j.issn.1000-3274.2002.03.017 [15] 帅向华, 姜立新, 刘钦等, 2009. 地震应急指挥技术系统设计与实现. 测绘通报, (7): 38—41, 54.Shuai X. H., Jiang L. X., Liu Q., et al., 2009. Earthquake emergency command system: design and implementation. Bulletin of Surveying and Mapping, (7): 38—41, 54. (in Chinese) [16] 王晓青, 丁香, 2004. 基于GIS的地震现场灾害损失评估系统. 自然灾害学报, 13(1): 118—125. doi: 10.3969/j.issn.1004-4574.2004.01.019Wang X. Q., Ding X., 2004. Estimation system of disastrous losses in seismic site based on GIS. Journal of Natural Disasters, 13(1): 118—125. (in Chinese) doi: 10.3969/j.issn.1004-4574.2004.01.019 [17] 谢礼立, 2006. 城市防震减灾能力的定义及评估方法. 地震工程与工程振动, 26(3): 1—10. doi: 10.3969/j.issn.1000-1301.2006.03.001Xie L. L., 2006. A method for evaluating cities' ability of reducing earthquake disasters. Earthquake Engineering and Engineering Vibration, 26(3): 1—10. (in Chinese) doi: 10.3969/j.issn.1000-1301.2006.03.001 [18] 张风华, 谢礼立, 范立础, 2004. 城市防震减灾能力评估研究. 地震学报, 26(3): 318—329. doi: 10.3321/j.issn:0253-3782.2004.03.011Zhang F. H., Xie L. L., Fan L. C., 2004. Study on Evaluation of cities' ability reducing earthquake disasters. Acta Seismologica Sinica, 26(3): 318—329. (in Chinese) doi: 10.3321/j.issn:0253-3782.2004.03.011 [19] 张肖, 齐玉妍, 云萌等, 2019. 基于分量钻孔应变资料分析测区构造应变变化. 地震, 39(4): 76—83.Zhang X., Qi Y. Y., Yun M., et al., 2019. Analysis on changes of the tectonic strain in survey area based on multi-component borehole straingauge. Earthquake, 39(4): 76—83. (in Chinese) [20] 赵春霞, 钱乐祥, 2004. 遥感影像监督分类与非监督分类的比较. 河南大学学报(自然科学版), 34(3): 90—93.Zhao C. X., Qian L. X., 2004. Comparative study of supervised and unsupervised classification in remote sensing image. Journal of Henan University (Natural Science), 34(3): 90—93. (in Chinese) [21] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会, 2009. GB/T 24335—2009 建(构)筑物地震破坏等级划分. 北京: 中国标准出版社.General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, Standardization Administration, 2009. GB/T 24335—2009 Classification of earthquake damage to buildings and special structures. Beijing: Standards Press of China. (in Chinese) [22] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会, 2011. GB/T 18208.4—2011 地震现场工作第4部分: 灾害直接损失评估. 北京: 中国标准出版社.General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, Standardization Administration, 2011. GB/T 18208.4—2011 Post-earthquake field works—Part 4: assessment of direct loss. Beijing: Standards Press of China. (in Chinese) [23] 中华人民共和国住房和城乡建设部, 2008. GB 50223—2008 建筑工程抗震设防分类标准. 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People’s Republic of China, 2008. GB 50223—2008 Standard for classification of seismic protection of building constructions. Beijing: China Architecture & Building Press. (in Chinese) [24] 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 2010. GB 50011—2010 建筑抗震设计规范(附条文说明)(2016年版). 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People’s Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China. 2010. GB 50011—2010 Code for seismic design of buildings. Beijing: China Architecture & Building Press. (in Chinese) [25] 周强, 邵峰, 孙柏涛, 2016. 江西村镇房屋抗震能力调查与分析. 地震工程与工程振动, 36(6): 188—197.Zhou Q., Shao F., Sun B. T., 2016. Investigation and analysis of seismic capacity of rural buildings in Jiangxi. Earthquake Engineering and Engineering Dynamics, 36(6): 188—197. (in Chinese) 期刊类型引用(9)
1. 施唯，王东明. 全国房屋设施抗震设防信息化管理平台设计与实现. 震灾防御技术. 2024(02): 229-240 . 
本站查看2. 管庆伟，任金铜，樊祖洪. 基于RS和经验估计的区域房屋抗震能力评估——以纳雍县勺窝乡为例. 贵州工程应用技术学院学报. 2023(03): 139-146 . 百度学术3. 张树君，朱耿青，郑超，郑韵，林岩钊，王辉山. 遥感技术在房屋抗震设防能力初判中的应用探析. 地理空间信息. 2023(08): 118-121 . 百度学术4. 赵杰，王跃杰，程紫燕，王飞剑，郑树平. 基于遥感影像和经验估计初判山西省建筑物抗震能力. 山西地震. 2023(03): 18-22 . 百度学术5. 郑韵，姜立新，张颖，朱耿青，张树君，黄玲珠. 基于遥感影像的仙游县建筑物抗震水平初判. 地震地磁观测与研究. 2023(04): 48-56 . 百度学术6. 杨钦杰，李蕾，聂冠军，赵修敏. 广西南部农村群体建筑抗震能力评价. 世界地震工程. 2023(04): 116-123 . 百度学术7. 李姜，李铁纯，张兴，郭彦宇. 农村人口流动对地震灾害人员伤亡评估的影响研究-以张家口为例. 震灾防御技术. 2022(03): 557-568 . 本站查看8. 黄文涛，熊永良，张莹，赵真，赵雪慧，周琦. 川南地区民居建筑地震易损性研究. 华北地震科学. 2022(04): 40-49 . 百度学术9. 陈波，王芳，肖本夫. “情景-应对”型理论体系的发展及其在地震灾害应急管理中的应用探讨. 震灾防御技术. 2021(04): 605-616 . 本站查看其他类型引用(1)
-
下载:
点击查看大图
计量
- 文章访问数: 232
- HTML全文浏览量: 78
- PDF下载量: 19
- 被引次数: 10
