基于单体建筑数据的震后人员埋压分布快速估计方法初探

曹颖; 安立强; 高武平; 张安东; 杨菲; 王志胜

doi:10.11899/zzfy20240105

基于单体建筑数据的震后人员埋压分布快速估计方法初探

doi: 10.11899/zzfy20240105

天津市地震局, 天津 300201

基金项目: 地震应急青年重点任务（CEAEDEM202201）

详细信息

作者简介:
曹颖，女，生于1990年。硕士，工程师。主要从事地震灾害风险评估工作。E-mail：wscaoying@163.com

通讯作者:
安立强，男，生于1984年。博士，高级工程师。主要从事地震灾害损失评估及地震安全韧性工作。E-mail：anliqiang@sina.com

12 https://data.earthquake.cn/
23 高武平，姚新强，张文朋等，2020. 宝坻建（构）筑物抗震性能调查与评价—天津市地震灾害风险调查与评估试点工作技术报告.
计量
- 文章访问数: 103
- HTML全文浏览量: 72
- PDF下载量: 49
- 被引次数: 1
出版历程
- 收稿日期: 2023-12-14
- 刊出日期: 2024-03-31

A Preliminary Study on a Rapid Estimation of People Trapped in the Single Building after Earthquake

Tianjin Earthquake Agency, Tianjin 300201, China

摘要

摘要: 最大限度降低人员伤亡，是震后开展救援工作的首要目标。当地震发生后，快速精确地估计震后人员埋压分布，为救援行动提供方向指导，是震后需要紧急开展的工作。本文在总结前人工作的基础上，提出了一种利用高精度数据开展震后人员埋压分布估计的方法。首先根据实验目的进行数据采集和整理；然后将单体建筑物按照结构类型和使用类型进行划分，在此基础上进行精细尺度人口空间化；最后依据典型建筑物易损性曲线和设定的PGA计算出建筑物破坏情况，进而实现对人员埋压位置的快速估计。本文以天津宝坻城区为研究区，选取了3种典型的震后埋压人员分布评估方法。在设定的2种地震情景下，计算结果显示人员埋压分布在数量上存在差异，但人员埋压分布区域和相对数量高值区一致，均为砖木结构，与之相对的抗震能力较强且建筑年龄较新的框架剪力墙住宅建筑发生震后人员埋压的可能性不大，本研究可为地震应急救援工作提供决策参考。
- 震后人员埋压 /
- 快速评估 /
- 单体建筑数据 /
- 精细化
Abstract: Minimizing casualties is the primary goal of post-earthquake relief efforts. When an earthquake occurs, it is urgent to quickly and accurately estimate the location distribution of buried person after the earthquake and to provide guidance for rescue operations. Based on the summary of previous work, this paper proposes a method for estimating the location distribution of victims buried after earthquake with high resolution data. First, data collection and sorting is carried out according to the purpose of the experiment. Then, the individual buildings are divided into structural types and usage types, and fine-scale population spatialization is carried out. Finally, the building damage situation is calculated according to the typical building vulnerability curve and the setting PGA, and then the buried position is quickly estimated. In this paper, taking the urban area of Baodi in Tianjin as the case study, three typical methods for assessing the distribution of buried people were selected. In earthquake scenarios settings, the calculation results showed that there are differences in the distribution of buried personnel in terms of quantity. However, the distribution area and the high-value region of the relative number of buried people were consistent, which is mainly in the brick-timber houses. In comparison, the relatively younger residential buildings with relatively stronger seismic resistance have a minimal likelihood of buried person after earthquake. Our work may provide decision-making reference for earthquake emergency rescue work.
- Earthquake buried person /
- Rapid assessment /
- Individual building data /
- Refinement
注释:

1) 12 https://data.earthquake.cn/

2) 23 高武平，姚新强，张文朋等，2020. 宝坻建（构）筑物抗震性能调查与评价—天津市地震灾害风险调查与评估试点工作技术报告.

HTML全文

引言

结构的抗倒塌能力一直都是地震灾害的研究重点。以往的震害经验表明，没有经过适当抗震设计的建筑结构在地震过程中可能遭遇更加严重的灾害。同时，近年来的地震调查表明，农村在地震中遭受了严重的破坏和生命、财产损失。如2008年汶川M_S 8.0地震，造成了建筑结构大面积倒塌，其中倒塌的楼房约有90%分布在农村和乡镇地区（孙景江等，2008；霍林生等，2009；李钢等，2009）；2014年鲁甸M_S 6.5地震也在农村、乡镇地区造成了重大损失（周铁钢等，2014）；类似的情况再次发生在2017年九寨沟M_S 7.0地震中，共造成73671间房屋不同程度受损，76栋农村地区房屋完全倒塌。中国有将近50.3%的人口生活在经济相对落后的农村地区，建筑结构缺少抗震设防措施，现存的大量年久失修和设计简陋的建筑结构（同济大学土木工程防灾国家重点实验室，2008），存在极大的安全隐患。

随着抗震经验的积累和抗震意识的加强，近年来农村地区开始重视建筑结构的抗震能力，意识到自行修建的多层框架结构缺乏科学的规划和设计，隐藏着巨大的潜在地震灾害危险，居民也纷纷采取一定的结构维护措施对结构进行维护和加固。但是，这些维护和加固措施受到当地经济和技术条件的限制，且缺乏专业的指导和科学的设计评估。针对这一现象，本文基于Ⅷ度设防区，设计了某典型的农村地区3层2跨非正规设计框架结构。利用增量动力分析（Incremental Dynamic Analysis，IDA）对结构进行非线性时程分析，并利用结构易损性方程评估结构在地震作用下的抗倒塌能力。通过对比分析地震作用下不同的支撑维护位置对非正规设计框架结构抗倒塌能力和抗震性能的影响，为农村地区框架结构的简单支撑维护提供一定参考。

1. 数值模型说明

针对Ⅷ度设防区（常遇地震PGA＝0.1g，设防地震PGA＝0.2g，罕遇地震PGA＝0.4g），设计了某城镇和乡村中常见的3层2跨非正规设计钢筋混凝土框架结构模型。首层层高4.5m，其余层高3m，总高10.5m。平面尺寸12m×6m，框架跨度6m，进深6m，如图 1所示。混凝土强度等级为C35，受力钢筋均采用HRB400。梁截面尺寸均为200mm×500mm，柱截面尺寸300mm×300mm，低于规范规定标准（中华人民共和国住房和城乡建设部等，2008）。

图 1 结构配筋示意（a）和柱网布置（b）

Figure 1. Reinforcement configurations of beams and columns (a) and layout of columns (b)

下载: 全尺寸图片幻灯片

采用有限元软件OpenSEES对结构进行非线性分析。OpenSEES中提供的非线性梁柱单元可较为精确地模拟梁柱构件和剪力墙构件在大变形下的非线性反应（Mansour，2016）。混凝土本构模型采用OpenSEES中提供的Concrete02线性拉伸强化材料模型（基于Kent-Scott-Park模型）；钢筋本构采用OpenSEES中提供的Steel02各项同性强化材料模型（基于Menegotto-Pinto模型）；截面对象采用纤维模型；结构梁、柱单元对象分别采用基于有限单元柔度法的非线性梁柱单元（Nonlinear Beam Column Element）和塑性铰梁柱单元（Beam With Hinges）。

对该非正规设计框架结构的外部支撑形式，采用最简单、常见的杆件斜向支撑形式，如图 2所示。考虑到农村地区最常见的材料为木材，本文中所用的杆件为木质杆件，截面为200mm×300mm。木杆件与地面成45°，底部与地面的连接简化成铰支，与结构的连接节点简化为固接。木材的材料参数，参照欧洲胶合木规范等级为GL24h的材料选取，并将木材简化为理想弹塑性模型。为了研究不同支撑位置的支撑效果，将易损框架结构分为3部分，分别为边框1、中框、边框2（图 1（b））。按照边框和中框的区别，设计了2个不同支撑位置的支撑方案，见表 1。

图 2 结构支撑方式示意

Figure 2. Schematic of support form

下载: 全尺寸图片幻灯片

表 1 方案设计

Table 1. Designing schemes

方案	支撑位置
方案	边框1	中框	边框2
非正规设计框架结构	无	无	无
支撑中框	无	有	无
支撑边框	有	无	无

下载: 导出CSV

| 显示表格

2. 数值计算方法

增量动力分析（IDA）是近年来发展起来的一种用于评价结构抗震性能的动力参数分析方法，现已广泛用于评估结构在遭遇不同强度地震时的动力响应和结构的抗倒塌性能。分析结果用IDA曲线体现，表现出地震动强度IM与结构响应DM间的关系（Vamvatsikos等，2002）。本文用峰值加速度（PGA）代表地震动强度IM，用最大层间位移角代表结构响应DM。

结构倒塌是结构的动态不稳定问题，包含竖向不稳定和侧向不稳定。本文将结构侧向不稳定视为结构倒塌。侧向倒塌是由结构侧向塑性变形过大而引起的侧向承载力丧失（吕大刚等，2011），在IDA曲线上的体现是较小地震强度的增加，导致了较大的结构横向变形，即IDA曲线变平。理论上，IDA曲线变平不可实现，本文的结构倒塌判断标准是Vamvatsikos等（2002）提到的基于地震强度（IM-based）的判断方法。该方法能够在IDA曲线上判断出1个点，定义为倒塌能力点（Capacity Point）。倒塌能力点清楚地将IDA曲线划分为2个区域，较低地震强度区为非倒塌区域，较高地震强度区为倒塌区域。本文根据基于结构响应的方法来定义结构的倒塌能力点，即IDA曲线斜率与结构弹性阶段初始斜率的比值为λ（λ＜1，联邦急救管理局推荐的20%）时，曲线上这一点所对应的地震强度为结构的倒塌能力点，当IDA曲线上有多个这样的点时，取地震强度最大的点作为结构的倒塌能力点。

结构倒塌易损性表现为结构在不同的地震强度下结构发生倒塌的概率，通常用对数正态分布函数$ \Phi $表示（吕大刚等，2011；Baker，2015），IM＝x时，结构的倒塌概率可表示为：

$$ P\left({C|IM = x} \right) = \Phi \left({\frac{{{\rm{ln}}\left({x{\rm{/}}\theta } \right)}}{\beta }} \right) $$

(1)

其中，θ、β为易损性函数参数，θ为均值参数，代表结构50%倒塌概率对应的PGA值，β为方差参数，可分别通过式（2）、（3）进行参数拟合：

$$ {\rm{ln}}\hat \theta = \frac{1}{n}\sum\limits_{i = 1}^n {{\rm{ln}}{{\left({IM} \right)}_i}} $$

(2)

$$ \hat \beta = \sqrt {\frac{1}{{n - {\rm{1}}}}\sum\limits_{i = 1}^n {{{({\rm{ln}}({{\left( {IM} \right)}_i}/\hat \theta ))}^2}} } $$

(3)

其中，(IM)_i为第i条地震记录的结构倒塌能力点，n表示用于计算的地震记录总数，$ \hat \theta $为倒塌能力点的对数平均值，$ \hat \beta $为对应的对数标准差。

3. 计算结果分析

IDA计算时，需要选取多条地震记录对结构进行非线性分析。地震记录可以是人工合成地震动或实际记录的地震动，本文从美国太平洋地震工程研究中心（Pacific Earthquake Engineering Research Center，PEER）选取了12个台站的24条地震记录作为地震输入，台站信息见表 2。

表 2 地震记录台站基本信息

Table 2. Basic information of the ground motion records

地震	年份	台站	M_W	V_S30/m·s^－1
帕克菲尔德	1966	Cholame - Shandon Array #12	6.19	408.93
帕克菲尔德	1966	San Luis Obispo	6.19	493.5
科林加	1983	Parkfield - Cholame 2E	6.36	522.74
科林加	1983	Parkfield - Cholame 3E	6.36	397.36
科林加	1983	Parkfield - Fault Zone 10	6.36	372.73
科林加	1983	Parkfield - Fault Zone 9	6.36	372.26
科林加	1983	Parkfield - Gold Hill 3E	6.36	450.61
科林加	1983	Parkfield - Gold Hill 4W	6.36	421.2
科林加	1983	Parkfield - Gold Hill 5W	6.36	441.37
科林加	1983	Parkfield - Stone Corral 2E	6.36	566.33
科林加	1983	Parkfield - Vineyard Cany 4W	6.36	386.19
科林加	1983	Parkfield - Vineyard Cany 6W	6.36	392.24

下载: 导出CSV

| 显示表格

3.1 易损性分析对比

根据IDA计算结果和结构的倒塌判断方法，2种支撑方案和该非正规设计框架结构的IDA曲线及对应的倒塌能力点如图 3所示。从图中可以看出，结构的IDA曲线走势基本类似，不同地震记录计算得到的结构倒塌点的分布离散性较大。

图 3 IDA曲线与对应的倒塌能力点

Figure 3. IDA curves and corresponding capacity points of 3 cases

下载: 全尺寸图片幻灯片

为保证结构在地震作用下的安全性，基本要求是结构在大震作用下的倒塌概率足够低（Haselton等，2011）。FEMA推荐的可接受的倒塌概率是10%或者更低（Federal Emergency Management Agency，2009），即结构在倒塌概率低于10%（或者更低）时视为结构安全。本文选取韩小雷等（2015）提到的5%作为结构可接受的倒塌概率。通过IDA计算结果拟合结构易损性曲线，得到2种支撑方案作用下的结构和非正规设计框架结构的倒塌易损性曲线，如图 4所示。本文建筑结构位于Ⅷ度设防区，罕遇地震等级PGA＝0.4g，从图中可以看出，易损框架在罕遇地震作用下的倒塌概率为17.5%，远超过可接受的倒塌概率，属于危险建筑。通过2种方案的支撑，易损框架结构倒塌概率降至6%和2.4%。其中，支撑中框的方案已将非正规设计框架结构的倒塌概率降低到了可接受的倒塌概率以下，可以视为安全。

图 4 3个工况倒塌易损性曲线

Figure 4. Collapse fragility curves of 3 cases

下载: 全尺寸图片幻灯片

对2种支撑方案的对比说明，支撑结构中间框架，对非正规设计框架结构在地震作用下抗倒塌能力的提升效果明显优于支撑边框。因为支撑中间框架时，保证了结构的对称性。因此，条件允许的情况下，应该尽量支撑中间框架。

3.2 结构地震响应的对比分析

为了进一步研究结构的抗震性能，本文研究了结构在不同地震强度（常遇、设防以及罕遇地震）作用下的结构响应。图 5为结构在选取的地震波作用下，结构平均层间位移角的对比（不予考虑0.4g作用下导致结构倒塌的地震波）。由图可以看出，结构通过支撑后，地震作用下结构响应形式与非正规设计框架结构类似，最大层间位移角都出现在底层。从图中还可以看出，支撑后的结构底层层间位移角均比非正规设计框架结构小，且支撑中间框架的底层结构响应也比支撑边框时小，说明支撑中间框架对结构的地震响应比支撑边框的效果好，这与上述支撑方案对结构地震作用下的抗倒塌性能的影响一致。

图 5 3种不同PGA下结构平均层间位移角

Figure 5. Average of structural seismic response under 3 different ground motion intensities

下载: 全尺寸图片幻灯片

4. 结论

本文针对农村地区易损框架结构的简单支撑，研究了外部斜向支撑位置对易损框架结构抗倒塌性能和抗震性能的影响。基于OpenSEES平台，通过IDA计算和分析结构易损性，量化评估不同支撑位置对易损框架结构的抗倒塌性能，并对比了结构在3种强度地震作用下的结构响应，主要得到以下结论：

（1）简单的支撑形式对易损框架结构地震作用下的抗倒塌能力提升明显，且支撑位置对支撑效果的影响十分明显，支撑中框的效果比支撑边框显著。

（2）合理安排支撑位置对结构抗震性能的提高十分有必要，支撑中间框架对地震作用下的结构响应也十分有利，因此在条件允许情况下，建议优先支撑中间框架。

图 1 研究区

Figure 1. Case study area of this work

下载: 全尺寸图片幻灯片

图 2 建筑物结构类型分布

Figure 2. Distribution of building structure types

下载: 全尺寸图片幻灯片

图 3 建筑物使用类型分布

Figure 3. Distribution of building usage types

下载: 全尺寸图片幻灯片

图 4 各使用类型建筑物在室人口密度变化

Figure 4. In-door population density of different building usage

下载: 全尺寸图片幻灯片

图 5 建筑结构地震易损性曲线

Figure 5. Seismic fragility curves of building structures

下载: 全尺寸图片幻灯片

图 6 研究区工作日人口分布情况

Figure 6. Population distribution of weekdays in the case study area

下载: 全尺寸图片幻灯片

图 7 情景1发生在10点钟时的埋压人员分布

Figure 7. Distribution of buried people in earthquake scenario 1 at 10 o'clock

下载: 全尺寸图片幻灯片

图 8 情景1发生在24点钟时的埋压人员分布

Figure 8. Distribution of buried people in earthquake scenario 1 at 24 o'clock

下载: 全尺寸图片幻灯片

图 9 情景2发生在10点钟时的埋压人员分布

Figure 9. Distribution of buried people in earthquake scenario 2 at 10 o'clock

下载: 全尺寸图片幻灯片

图 10 情景2发生在24点钟时的埋压人员分布

Figure 10. Distribution of buried people in earthquake scenario 2 at 24 o'clock

下载: 全尺寸图片幻灯片

参考文献(23)

安立强, 2022. 基于遥感数据的地震极灾区快速评估. 哈尔滨: 中国地震局工程力学研究所.

An L. Q. , 2022. Rapid damage assessment of extreme earthquake disaster areas based on remote sensing data. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese)

白仙富, 戴雨芡, 叶燎原等, 2018. 基于高精度房屋、人口数据的地震压埋人员和紧急搜救区评估模型研究. 防灾减灾学报, 34（4）: 1—12

Bai X. F. , Dai Y. Q. , Ye L. Y. , 2018. Assessment models of trapped-victim in earthquake & emergency rescue area based on high-precision building and population data. Journal of Disaster Prevention and Reduction, 34(4): 1—12. (in Chinese)

谷国梁, 安立强, 朱宏等, 2021. 城市地震压埋人员分布评估研究——以天津市区为例. 地震工程学报, 43（6）: 1352—1360

Gu G. L. , An L. Q. , Zhu H. , et al. , 2021. Assessment of seismic buried personnel in urban area: a case study of Tianjin urban area, China. China Earthquake Engineering Journal, 43(6): 1352—1360. (in Chinese)

李海华, 1987. 震灾最大伤亡率的初步分析. 灾害学, （2）: 41—47

Li H. H. , 1987. Preliminary analysis on the highest rate of casualty in earthquake disasters. Journal of Catastrophology, 2(2): 41—47. (in Chinese)

马玉宏, 谢礼立, 2000. 地震人员伤亡估算方法研究. 地震工程与工程振动, 20（4）: 140—147.

Ma Y. H., Xie L. L., 2000. Methodologies for assessment of earthquake casualty. Earthquake Engineering and Engineering Dynamics, 20(4): 140—147. (in Chinese)

任静, 张方浩, 李志强等, 2021.2021年5月21日云南漾濞6.4级地震灾害损失快速评估结果分析. 震灾防御技术, 16（3）: 454—466

Ren J. , Zhang F. H. , Li Z. Q. , et al. , 2021. May 21, 2021, Yangbi, Yunnan M_S6.4 analysis of the results of rapid assessment of earthquake disaster losses. Technology for Earthquake Disaster Prevention, 16(3): 454—466. (in Chinese)

田丽莉, 2012. 地震灾害人员伤亡影响因素分析及人员伤亡估算公式研究. 北京: 首都经济贸易大学.

王东明, 李永佳, 陈洪富等, 2012. 汶川地震与玉树地震自救互救调查情况比较研究. 国际地震动态, （5）: 19—25

Wang D. M. , Li Y. J. , Chen H. F. , et al. , 2012. Comparative studies of status of Self Rescue and mutual aid in Wenchuan and Yushu earthquake. Recent Developments in World Seismology, (5): 19—25. (in Chinese)

王紫荆, 陈文凯, 苏浩然等, 2020. 基于区域特征的地震埋压人员评估模型研究. 地震工程学报, 42（5）: 1236—1247

Wang Z. J. , Chen W. K. , Su H. R. , et al. , 2020. Evaluation model of earthquake-buried personnel based on the regional characteristics. China Earthquake Engineering Journal, 42(5): 1236—1247. (in Chinese)

王紫荆, 2021. 面向地震应急的人员压埋模型研究. 兰州: 中国地震局兰州地震研究所.

Wang Z. J. , 2021. Research on personnel burying model for earthquakes emergency. Lanzhou: Lanzhou Institute of Seismology. (in Chinese)

吴兵, 李志强, 齐文华等, 2017. 以乡镇居民地为单元的震后埋压人员分布与救援对策——以乌鲁木齐市为例. 华北地震科学, 35（3）: 13—19

Wu B. , Li Z. Q. , Qi W. H. , et al. , 2017. Trapped population distribution in earthquake and rescue policy based on township residential Area scale: a case study of Urumqi Area. North China Earthquake Sciences, 35(3): 13—19. (in Chinese)

吴晨, 殷海涛, 孙庆文等, 2011. 基于ArcGIS的地震搜救信息系统可行性研究. 高原地震, 23（2）: 54—59

Wu C. , Yin H. T. , Sun Q. W. , et al. , 2011. Study on the feasibility of the search and rescue information system based on ArcGIS. Plateau Earthquake Research, 23(2): 54—59. (in Chinese)

肖东升, 黄丁发, 陈维锋等, 2009. 地震压埋人员压埋率预估模型. 西南交通大学学报, 44（4）: 574—579

Xiao D. S. , Huang D. F. , Chen W. F. , et al. , 2009. Prediction model for buried personnel probabilty in earthquake. Journal of Southwest Jiaotong University, 44(4): 574—579. (in Chinese)

谢江丽, 姚远, 阿布都瓦里斯等, 2021. 城市人口发展与地震救援的思考——以乌鲁木齐市为例. 内陆地震, 35（1）: 100—108.

许建东, 危福泉, 张来泉等, 2008. 地震人员伤亡与压埋人员评估方法的初步研究——以福建省漳州市区为例. 地震研究, 31（4）: 382—387

Xu J. D. , Wei F. Q. , Zhang L. Q. , et al. , 2008. Preliminary study on evaluating the number of casualties and trapped victims by an earthquake——A case study of Zhangzhou City, Fujian Province. Journal of Seismological Research, 31(4): 382—387. (in Chinese)

杨明儒, 2014. 地震压埋人员分布评估方法研究. 哈尔滨: 中国地震局工程力学研究所.

Yang M. R. , 2014. Study on the assessment method of buried personnel distribution in earthquake. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese)

余世舟, 张令心, 杨明儒, 2015. 地震压埋人员分布评估方法初探. 地震工程与工程振动, 35（2）: 138—143

Yu S. Z. , Zhang L. X. , Yang M. R. , 2015. Assessment of buried people distribution after earthquake. Earthquake Engineering and Engineering Dynamics, 35(2): 138—143. (in Chinese)

郑山锁, 孙龙飞, 龙立等, 2019 a. 城市地震灾害损失评估: 理论方法、系统开发与应用. 北京: 科学出版社.

郑山锁, 张睿明, 陈飞等, 2019 b. 地震人员伤亡评估理论及应用研究. 世界地震工程, 35（1）: 87—96

Zheng S. S. , Zhang R. M. , Chen F. , et al. , 2019 b. Research on theory and application of earthquake casualty estimates. World Earthquake Engineering, 35(1): 87—96. (in Chinese)

周中红, 陈文凯, 何少林等, 2019. 基于人口公里格网的地震压埋人员分布评估方法的应用与评价——以甘肃岷县漳县6.6级地震为例. 地震研究, 42（2）: 288—294

Zhou 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)

An L. Q. , Zhang J. F. , 2022. Impact of urbanization on seismic risk: a study based on remote sensing data. Sustainability, 14(10): 6132. doi: 10.3390/su14106132

Shinozuka M. , Feng M. Q. , Lee J. , et al. , 2000. Statistical analysis of fragility curves. Journal of Engineering Mechanics, 126(12): 1224—1231. doi: 10.1061/(ASCE)0733-9399(2000)126:12(1224)

Wei B. Y., Nie G. Z., Su G. W., et al., 2017. Risk assessment of people trapped in earthquake based on km grid: a case study of the 2014 Ludian earthquake, China. Geomatics, Natural Hazards and Risk, 8(2): 1289—1305.

施引文献