Structural Characteristics and Seismic Capacity Analysis of Rural Buildings in the Southeast of Qinghai Province
-
摘要: 我国农村民居的抗震能力具有显著的地域特点。本文基于现场调研结果总结了青海东南部黄河流域农村民居的结构特点,并对其抗震能力进行了评估。庄廓院是调研区内的典型农村民居,本文对调研区内庄廓院的木构架房屋进行研究,提出庄廓土墙围护木构架房屋的震害等级分类标准,根据实际震害资料和类比的方法确定了当地土墙围护木构架结构的震害矩阵和震害指数与峰值加速度的关系曲线。与云南鲁甸农村的土木结构房屋比较,青海庄廓院民居具有更高的抗震能力;与未设防的砖混结构比较,这种土墙围护的木构架结构在高烈度区造成人员伤亡的风险更低。Abstract: The anti-seismic capacity of rural buildings in China is of significant regional characteristics. Based on the investigation, this paper summarizes the structural characteristics and evaluates the seismic capacity of rural buildings in the Yellow River Basin in southeastern of Qinghai. Zhuangkuo folk house is a typical rural building in the research area. Through the investigation on the wooden frame structures in the rural buildings, we put forward to earthquake damage classification of the wooden frame structures with earth walls in Zhuangkuo folk house. According to the actual damage information, earthquake damage matrix of the wooden frame structures with earth walls in Zhuangkuo folk house by analogy and the relationship are determined between the earthquake damage index and the peak acceleration. Compared with the civil structure houses in the rural buildings of Ludian in Yunnan province, rural buildings have relative high anti-seismic capacity in Qinghai Province. The wooden frame structures with the earth wall are at a lower risk of casualties than the unfortified brick-concrete structures under the high-intensity condition.
-
引言
建筑物作为灾害的主要承灾体(陈振拓等,2012),其空间分布是震前开展地震灾害风险评估、震后进行灾害损失快速评估的基础数据(王龙等,2007;谷国梁等,2016;李皓等,2018)。基于精细可靠的建筑物空间分布数据开展评估,震前可给出抗震薄弱环节及相应的抗震措施、危旧房屋集中地段、高风险小区(孙柏涛等,2017;王志涛等,2019);震后可给出更可靠的灾情快速评估结果和可能被困人员的空间分布,为地震应急处置和现场应急救援提供可靠的信息支持,从而达到减轻地震灾害损失的目的(韩贞辉等,2013;丁文秀等,2015)。
目前,建筑物数据主要有两类:一类是建筑物统计数据,其广泛应用于震后快速评估中(王晓青等,2009;徐国栋等,2011),该数据与所属区域的空间特征脱节而无法反映区域内建筑物空间分布特点,仅反映统计单元内建筑物的总体情况,且其可靠性和精细化程度较低,数据更新速度慢,时效性较差;另一类是基于多源数据融合方式生成的建筑物公里格网数据(韩贞辉,2013;杨海霞等,2015;杨海霞,2015),该数据虽能较好地反映建筑物空间分布的宏观特征,但不能准确反映建筑物空间分布的细节特征,且未考虑不同地域建筑物空间分布的差异,无法反映区域间建筑物空间分布的异质性。然而建筑物空间分布具有明显的地域差异,采用同一建筑物空间分布格网化模型具有一定的局限性。
针对未考虑建筑物空间分布异质性问题,在充分考虑目标区域自然环境、社会经济条件等差异的基础上,开展建筑物特征一致性分区。根据分区结果,采用抽样的方法获取各区域建筑物样本数据。在此基础上,探究不同区域建筑物空间分布与各影响因子之间的关系,分区构建基于多因素的建筑物数据空间化模型,并根据模型生成建筑物格网数据,从而提高建筑物数据格网化精度。
1. 研究方法
1.1 总体技术思路
建筑物空间分布是受诸多因素影响的复杂非线性问题,不同地区的社会经济情况、自然地理条件等存在差异,建筑物空间分布特征也不相同(程晓亮等,2008;曾祥贵,2013)。因此,需将具有相似建筑物空间分布特征的区域划分为同一区域。在此基础上,分析不同区域建筑物空间分布与各影响因子之间的关系,分区构建建筑物空间分布格网化模型,反演建筑物空间分布情况。具体技术路线如图1所示。
首先,从自然地理和社会经济两方面出发选择分区指标,并以乡镇为基本单元进行统计,利用主成分分析方法提取主要特征,对其进行特征一致性分区。然后,在归纳总结目前基于多因素的建筑物空间分布格网化方法考虑影响因子的基础上,综合选择自然因子(如高程、坡度、坡向、地形起伏度及河流等)、社会经济因子(如道路、土地利用等),分析各影响因子对建筑物空间分布的影响,并对各因子进行分级设计。在此基础上,针对不同区域,从研究区中选取一定数量的格网作为建模区,采用目视解译的方式获取建模区的建筑物,统计分析不同建模区建筑面积密度与各影响因子之间的关系。采用相关系数法计算不同区域内各影响因子对建筑物空间分布的影响权重,并以归一化的相关系数作为该类因子的权重系数,对于同一类因子的不同子类,采用建筑面积密度作为权重的估计值。在此基础上,分区构建基于多因素的建筑物数据格网化模型。最后,随机选择若干格网作为检验样本,开展模型精度评估。
1.2 建筑物空间分布格网化方法
在借鉴已有研究的基础上,本文依据土地地表覆盖数据分类体系(Gong等,2013;Li等,2017),将研究区域划分为建设用地和非建设用地,其中建设用地指不透水层,包括城镇用地、工矿用地、商服用地、交通设施用地等,非建设用地包括耕地、森林、草地、灌木林、湿地、水域、裸地、冰川等。由于建筑物大多分布在建设用地内,少数或零星分布在非建设用地内(杨海霞等,2015),因此,针对这两类区域分别建立基于多因素的建筑物数据空间化模型。
1.2.1 单因子建筑物空间分布空间化权重
将研究区划分为规则格网,假定共有
$ {N}_{t} $ 个格网,如果格网t内存在影响建筑物空间分布的因子$ {F}_{i} $ ,其与建设用地和非建设用地的空间分布关系如图2所示。为探究因子$ {F}_{i} $ 对建筑物空间分布的影响程度,将其进一步划分为Q个子因子,在假设研究区内子因子$ {F}_{ij} $ 在第m类区域的建筑占地面积密度$ {D}_{mij} $ 不变的基础上,统计格网g内与第m类区域内第j个子因子$ {F}_{ij} $ 的占地面积为$ {A}_{mgij} $ ,则格网g内依据因子$ {F}_{i} $ 估计的建筑面积为:$$ {S}_{ mgi}=\sum _{j=1}^{Q}\left({D}_{mij}\cdot{A}_{mgij}\right) $$ (1) $ {\mathrm{D}}_{\mathrm{m}\mathrm{i}\mathrm{j}} $ 建筑占地面积密度常采用抽样统计的方式获取。假设抽样区域中,第m类区域中第i个因子第j个子因子$ {F}_{ij} $ 的占地面积为$ {A}_{mij} $ 、建筑面积为$ {S}_{mij} $ ,则建筑占地面积密度${{D}}_{{m}{i}{j}}$ 计算公式为:$$ {D}_{mij}={S}_{ mij}/{A}_{mij} $$ (2) 1.2.2 多因子建筑物分布空间化权重
建筑物空间分布是自然因子和社会经济条件共同作用的结果(江东等,2002),因此,基于单一因子的建筑物数据空间化模型可靠性及结果的准确度较低。本文选择多源数据融合的方法,在各类单因子格网建筑面积预测的基础上,根据各类因子之间的权重,通过权重比求和的形式综合同一格网内各类因子建筑面积预测结果。
假设影响建筑物空间分布的因子有M类,在m类区域中,第i类因子
$ {F}_{i} $ 的建筑面积分配权重为$ {W}_{mi} $ ,则格网g内建筑面积的计算公式为:$$ {S}_{ mg}=\sum _{i=1}^{N}\left({W}_{mi}\cdot{S}_{ mgi}\right) $$ (3) 目前,确定各类因子权重
$ {W}_{mi} $ 的方法主要有专家打分法、层次分析法、相关系数法等。其中专家打分法的主观性太强,与专家的偏好程度有很大关系。层次分析法依赖于专家经验,主观性相较于专家打分法有所降低。相关系数法是通过抽样确定的,目前在进行GIS空间分析时使用较多。2. 研究区及数据源
2.1 研究区概况
本文以雅安市为研究区开展试验。该区位于四川省中部,跨越四川盆地和青藏高原两大地形,地势呈北、西、南较高,东部和中部较低。全市辖2个市区和6个县,幅员面积为15 046 km2,其中山地占比94%,平原占比6%。
2.2 数据源与预处理
本文在分析研究区内建筑物空间分布特征及建筑物空间分布影响因子的基础上,模拟建筑物空间分布情况。因此,本文所需的基础数据主要包括行政区划数据、土地覆盖类型数据(FROM-GLC30)、DEM数据、道路交通数据、河流水系数据及人口统计数据。
由于建筑物空间分布具有显著的地域差异,若采用同一建筑物空间分布格网化模型具有一定的局限性,数据准确性难以保证。故本文在构建建筑物格网化模型时,首先进行建筑物空间分布特征一致性分区。建筑物空间分布受诸多因素影响,主要包括自然环境(如高程、坡度、坡向、水系等)及社会条件(如土地利用、道路、人口等),本文综合考虑上述方面进行建筑物一致性分区。
根据研究区建筑物空间分布特征,本文选择人口密度、路网密度、河网密度、平均高程、平均坡度、建设用地面积比例、耕地面积比例共7个指标进行建筑物特征一致性分区。首先,以乡镇为基本单元统计各指标的值,对其进行极差标准化处理,利用主成分分析方法进行降维。分析得到的第一主成分特征值和第二主成分特征值分别为3.743、1.483,且这两个主分量承载了74.68%的原始信息,其余主成分特征值均小于1,表明引入该主成分的影响程度不如1个基本变量。因此,选择第一主成分和第二主成分基本可表达建筑物分布特征。根据主成分分析得到的第一主成分和第二主成分计算公式,计算雅安市143个乡镇主第一主成分和第二主成分得分,再根据各主成分贡献率的比重作为权值,得到各乡镇综合得分,利用ArcGIS软件自然断点法将其划分为5个特征分区,如图3所示。
3. 模型构建与结果分析
3.1 建筑物空间分布影响因子分析
在归纳总结目前基于多因素的建筑物空间分布格网化方法考虑因素的基础上,结合研究区建筑物空间分布特征,综合选取土地利用、高程、坡度、坡向、河流距离、道路距离、地形起伏度7类因子,并对各类因子进行合理分类或分级,具体方案如表1所示。其中,高程等级划分原则为对建筑物分布密集段高程细分,对建筑物分布稀疏段高程粗分;坡度等级划分原则依据国际地理学联合会地貌调查与地貌制图委员制定的分类标准;坡向根据方位划分为阴坡、半阴坡、阳坡、半阳坡。河流距离和道路距离采用等间距划分等级;地形起伏度根据我国基本地貌类型划分等级。
表 1 影响因子子类分级Table 1. Classification of influencing factors土地利用 高程/m 坡度/(°) 坡向/(°) 河流距离/m 道路距离/m 地形起伏度/m 耕地 500~800 平原
(0~0.5)平缓坡
(−1)0~200 0~200 平原(<30) 森林 800~1 100 微斜坡(0.5~2) 向阳坡(135~225) 200~400 200~400 台地
(30~70)草地 1 100~1 400 缓斜坡(2~5) 向阳坡
(45~135,
225~315)400~600 400~600 丘陵
(70~200)灌木林 1 400~1 700 斜坡
(5~15)阴坡(0~45,
315~360)600~800 600~800 小起伏山地(200~500) 湿地 1 700~2 000 陡坡
(15~35)— 800~1 000 800~1 000 中起伏山地(500~1 000) 水体 2 000~2 400 峭坡
(35~55)— 1 000~1 200 1 000~1 200 大起伏山地(1 000~2 500) 不透水面 2 400~2 800 垂直壁
(55~90)— 1 200~1 400 1 200~1 400 极大起伏山地(>2 500) 裸地 2 800~3 200 — — 1 400~1 600 1 400~1 600 — 冰川 3 200~3 800 — — 1 600~1 800 1 600~1 800 — — 3 800~4 400 — — 1 800~2 000 1 800~2 000 — — >4 400 — — >2 000 >2 000 — 3.2 因子权重系数确定
以雅安市为研究区,将其划分为300 m×300 m的规则格网,基于建筑物特征一致性分区结果,从不同区域分别选取一定数量的格网作为建模区,如图4所示。根据天地图影像,采用目视解译的方式获取建模区内单体建筑物数据。
根据研究区各类因子数据(图5),针对不同区域,分别将抽样格网的影响因子与建设用地、非建设用地进行叠加分析,统计区域内各影响因子的子类或分级占地面积和建筑面积,根据式(2)确定各子类或分级的建筑占地面积密度。
针对各类影响因子权重的确定,本文选用相关系数法,根据同一特征分区内某一因子的子类或分级建筑面积密度与子类或分级占地面积比例之间的关系,确定相关系数,并以归一化的相关系数作为该类因子的权重系数,结果如表2所示。
表 2 雅安市不同特征分区抽样统计的各类因子权重Table 2. Weights of various factors in sampling statistics in different regions of Ya'an city分区 区域划分 权重 土地利用 高程 坡度 坡向 河流距离 道路距离 地形起伏度 一区 建设用地 — — 0.263 0.108 — 0.300 0.329 非建设用地 0.244 — 0.225 0.164 0.064 0.083 0.220 二区 建设用地 — — 0.230 0.231 0.144 0.159 0.236 非建设用地 0.163 — 0.178 0.154 0.187 0.141 0.177 三区 建设用地 — 0.216 0.213 0.156 — 0.217 0.198 非建设用地 0.227 0.177 0.224 — — 0.183 0.189 四区 建设用地 — 0.089 0.195 0.185 0.169 0.148 0.214 非建设用地 0.236 0.183 0.221 — — 0.175 0.185 五区 建设用地 — 0.176 0.183 0.169 0.156 0.147 0.170 非建设用地 0.156 0.168 0.192 0.180 — 0.148 0.157 3.3 建筑物空间分布格网化
在确定不同特征区各因子的不同子类或分级的建筑面积密度和各类因子权重的基础上,对雅安市内所有格网进行影响因子、建设用地和非建设用地统计,在此基础上,分别进行建设用地和非建设用地建筑物空间分布预测。考虑到1个格网内有可能出现预测面积小于1栋房屋的建筑面积,这与实际情况不符,需对这部分数据予以舍去。通过对雅安市建筑物开展野外调查,确定建筑用地单栋房屋建筑面积约为150 m2、非建筑用地单栋房屋建筑面积约为200 m2,因此,将建设用地预测面积<150 m2的格网建筑面积舍弃,将非建设用地预测面积<200 m2的格网建筑面积舍弃,即将格网的建筑面积设置为0。最后,将同一格网内的建设用地和非建设用地建筑面积相加,得到雅安市建筑面积预测结果,如图6所示。
3.4 精度验证
为验证建筑物空间分布格网化方法的预测精度,从研究区中随机选择若干格网作为验证样本,对验证格网内的建筑物进行解译,并统计格网内建筑面积,与模型得到的预测值进行对比,从而计算建筑面积相对误差,分级统计结果如表3所示。由表3可知,严重低估或严重高估的比例≤10%,较准确估计的比例接近75%,说明本文模型具有较高的精度。精度较差的格网主要分布在建筑分布特别密集或稀疏的地区,其原因可能在于建模时这部分样本抽取的比例较小,构建的模型对这类格网估计的结果不是很准确或未考虑因素之间的相关性,造成信息冗余。
表 3 相对误差分级统计Table 3. Statistics of relative error classification分级 数目 比例/% 严重低估,<−50% 794 4.7 一般低估,[−50%,−20%) 602 3.6 较准确估计,[−20%,20%] 12 568 74.7 一般高估,(20%,50%] 1 228 7.3 严重高估,>50% 1 626 9.7 4. 结论与讨论
本文考虑到建筑物空间分布的地域差异,根据影响建筑物空间分布的因素进行建筑物特征一致性分区。在分析建筑物空间分布影响因子的基础上,结合区域建筑物空间分布特征,综合选取土地利用、高程、坡度、坡向、河流距离、道路距离、地形起伏度7类因子,以归一化的相关系数作为各因子的权重系数,将建筑面积密度作为各因子子类或分级的权重,分区构建基于多因素的建筑物格网化模型,预测建筑物空间分布情况。以雅安市为研究区,得到了300 m格网的建筑物空间分布数据,并开展模型精度检验,结果表明该方法具有较好的精度。基于本文提出的方法获取的建筑物空间分布格网可为震前防御、震中救灾和震后重建提供可靠的数据支撑。
建筑物空间分布情况是考虑诸多影响因素的复杂问题,影响因子的选取是否全面、因子之间的相关性如何确定均会对建筑物空间分布产生影响。未来可深入分析影响建筑物空间分布的因素,进一步提高模型的精度。
-
表 1 调研区及地区破坏性地震情况一览表
Table 1. List of destructive earthquakes in the research area and the whole region
地点 时间 震级 震中烈度 青海托索湖西 1937.01.07 7.5 Ⅹ 青海阿兰克湖 1963.04.19 7.0 Ⅷ+ 青海共和 1990.04.26 6.9 Ⅸ 青海玛沁 1999.11.26 5.0 Ⅵ 青海玉树 2010.04.14 7.1 Ⅺ+ 青海门源 2016.01.21 6.4 Ⅷ 青海杂多 2016.10.17 6.2 Ⅶ 表 2 庄廓土墙围护木构架结构震害等级分类
Table 2. Earthquake damage classification of the wooden frame structures with earth walls
震害等级 震害描述 基本完好 木构架结构完好,或庄廓院墙轻微开裂,不加修理可继续使用。 轻微破坏 木构架结构出现轻微移位,庄廓院墙有明显裂缝,土墙体局部有少量崩塌,不需要修理或稍加修理即可继续使用。 中等破坏 木构架出现可见移位,庄廓院墙酥裂,裂缝严重,屋顶局部出现滑塌等,需要一般修理后可使用。 严重破坏 木构架移位较严重,木构架出现歪斜、部分断裂,庄廓院夯土院墙局部倒塌,房屋的围护土墙体酥裂、修复困难。 毁坏 木构架严重倾斜折断,房屋濒于崩溃或已倒毁,庄廓院夯土围墙倒塌、已无修复可能。 表 3 庄廓土墙围护木构架结构震害矩阵(%)
Table 3. Earthquake damage matrix of the wooden frame structures with earth walls
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 60 41 10 0 0 轻微破坏 32 30 22 8 0 中等破坏 8 22 42 38 20 严重破坏 0 7 20 42 50 毁坏 0 0 6 12 30 震害指数 0.096 0.197 0.412 0.582 0.73 表 4 未设防砖混结构的标准震害矩阵(%)
Table 4. Standard earthquake damage matrix of unfortified brick-concrete structures
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 77 51 13 4 0 轻微破坏 22 29 25 9 0 中等破坏 1 16 36 20 9 严重破坏 0 4 20 52 40 毁坏 0 0 6 15 51 震害指数 0.048 0.15 0.394 0.612 0.826 表 5 云南地区土木结构的震害矩阵(%)(1992—2008年模型)
Table 5. Earthquake damage matrix of civil structures in Yunnan province (the model in 1992—2008)
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 基本完好 68.96 38.59 6.91 0 中等破坏 29.8 53.61 65.41 40.9 毁坏 1.24 7.8 27.68 59.1 震害指数 0.132 0.292 0.538 0.755 表 6 地震动峰值加速度与地震烈度对照表
Table 6. Comparison of earthquake peak acceleration and earthquake intensity
地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ Ⅹ 地震动峰值加速度/g 0.04≤amax≤0.09 0.09≤amax≤0.19 0.19≤amax≤0.38 0.38≤amax≤0.75 amax≥0.75 平均加速度/g 0.063 0.125 0.25 0.5 1 -
崔文河, 王军, 2015.青海藏族庄廓民居及聚落研究—以巴麻堂村为例.华中建筑, (1):74-80. http://d.wanfangdata.com.cn/Periodical_hzjz201501015.aspx 葛学礼, 朱立新, 于文等, 2010.村镇木构架土坯围护墙房屋1/2缩尺模型振动台试验研究.土木建筑与环境工程, 32(增刊2):508-510, 515. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-CQDX201012001065.htm 郭光玲, 2016.汉中农村既有砌体房屋抗震加固技术研究.震灾防御技术, 11(1):111-116. doi: 10.11899/zzfy20160112 侯建盛, 李洋, 宋立军等, 2015.2014年云南景谷6.6级地震与云南鲁甸6.5级地震致灾因素分析.灾害学, 30(2):100-101, 143. https://www.wenkuxiazai.com/doc/628c9081cf84b9d529ea7a5b.html 李方一, 高晓路, 王英杰, 2012.中国农村地区住房结构的区域差异及其影响因素.地理科学进展, 30(12):1555-1563. http://www.cqvip.com/QK/96015A/201112/40264686.html 梁琦, 2005.青海少数民族民居与环境.西宁:青海人民出版社. 邱舒睿, 高惠瑛, 2015.青海省农居地震灾害易损性研究.震灾防御技术, 10(4):969-978. doi: 10.11899/zzfy20150415 孙柏涛, 陈洪富, 闫培雷等, 2014.南北地震带房屋建筑抗震能力分区特征研究—以四川省为例.土木工程学报, 47(S1):6-10. http://www.cnki.com.cn/Article/CJFDTotal-TMGC2014S1002.htm 王强, 王兰民, 袁中夏等, 2008.陇南地区农村民房遭受汶川8.0级地震震害分析及启示.世界地震工程, 24(3):20-24. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=sjdzgc200803004 肖梅玲, 叶燎原, 刘本玉等, 2012.云南农村土木结构房屋有限元抗震分析.地震研究, 35(2):276-281. http://d.old.wanfangdata.com.cn/Periodical/dzyj201202018 杨钦杰, 刘华贵, 高鹏飞等, 2016.广西农村民居抗震性能调查研究.震灾防御技术, 11(1):100-110. doi: 10.11899/zzfy20160111 中国地震局, 2009.GB/T 17742-2008中国地震烈度表.北京:中国标准出版社. 中国地震局, 2015.GB 18306-2015中国地震动参数区划图.北京:中国标准出版社. -