Design and Implementation of Multi-tenant SaaS Cloud Platform for High-Performance Seismic Insurance Loss Assessment
-
摘要: 地震巨灾风险的特点是低频高损,历史震害数据缺乏、风险暴露快速变迁等因素导致基于大数定理的费率厘定方法无法针对各区域不同建筑类型的风险暴露进行精细化定价。本文基于“五代图”潜在震源区模型的随机事件集解决观测数据不足的问题;并使用“五代图”所采用的地震动参数衰减关系模型与工程易损性方法计算地震事件对风险暴露造成的损失,从而计算费率厘定、地震风险管理需要的必备参数。本文采用云计算平台的弹性伸缩计算技术,实现动态按需分配计算资源,满足多用户并发使用的业务需求;同时采用以业务数据为单元的数据隔离方案,构建支持多租户的高性能地震保险损失评估SaaS云平台。Abstract: Earthquake catastrophe risk is characterized by low frequency and extreme severity. Due to the lack of historical earthquake loss data and the rapid changes of exposure, it is difficult to price the various exposure in different regions precisely based on the law of large numbers. And in this paper, the problem of insufficient data is solved by the application of the random event set based on the potential seismic source zones of the Fifth Generation National Seismic Zoning Map. The losses caused by those earthquakes are calculated by utilizing the Ground Motion Prediction Equation (GMPE) adopted by the Fifth Generation National Seismic Zoning Map and classic engineering vulnerability methods. Based on the output, we calculate the necessary parameters for insurance pricing and risk managing of seismic catastrophe. To meet the concurrent use needs of a large number of users, we implement the dynamic on-demand management of computing resources by adopting the elastic scaling computing technology based on the cloud computing platform. And to ensure the separation and security of user data, we isolate and split the data by the exposures uploaded by users. By the methods above, we implement a multi-tenant high-performance seismic insurance loss assessment SaaS cloud platform.
-
Key words:
- Earthquake Insurance Loss /
- Risk Assessment /
- Catastrophe Model /
- Pricing
-
引言
地震是一种极具破坏性的灾难,全世界地震多发的国家都因地震造成不可弥补的损失,人民生命与财产安全受到严重威胁。对建筑结构采取抗震设防措施是减轻地震灾害损失的根本手段,而地震区划是工程抗震设防的基础和前提条件。现阶段超过设防水平的大地震的发生不能忽视,如唐山地震(1978年)、汶川地震(2008年)等多次地震中,地震实际烈度往往超出抗震设计规范的设防烈度(陈肇元,2008),造成大面积楼房倒塌,数百万人流离失所。因此,灾前对既有建筑结构开展震害评估,并采取相应防灾措施尤为重要。传统的地震预测分析方法通常指对结构在特定地震作用下发生不同破坏状态的可能性进行分析,由于未充分考虑地震活动的不确定性,不能较好地反映建筑结构潜在的地震风险水平。Preciado等(2015)进行建筑结构震害风险评估时联合结构易损性和场地危险性,并以二者卷积为表征,表示结构在未来某一段时间内超过或发生某种地震破坏状态的概率,从而明确度量建筑结构未来可能发生的风险水平,为做好灾前预防、地震应急准备和灾后救援工作提供科学的决策依据,也是贯彻“从减少灾害损失向减轻灾害风险转变”(罗华春,2017)思想的体现。
在长期的实践探索中,国内外学者在结构震害风险研究方面取得了一系列有价值的成果。Melani等(2016)采用IDA方法对某低层钢筋混凝土框架结构开展易损性分析,并基于框架结构概率风险四次积分公式,对其进行地震风险分析和年损失估计。Berto等(2012)基于强震环境对钢筋混凝土构件的影响研究,提出概率风险管理方法,计算与建筑结构长期性能相关的地震风险。谷音等(2012)针对某典型矮塔斜拉桥,采用LHS-MC方法对其地震风险概率进行评估。尹超等(2016)以实际路堤结构为例,采用蒙特卡罗抽样方法评价有无挡土墙路堤的震害风险,提出基于风险可接受度的路堤抗震设计与管理方法。然而,城市中存在的大量钢筋混凝土框架结构震害风险尚未可知,对不同抗震设防烈度区框架风险等级的划分仍处于研究阶段。
基于上述认识,综合考虑我国地域辽阔且国土范围内地震危险性差异大的实际情况(马玉宏等,2002)和目前国内最主要的钢筋混凝土建筑结构形式,本文以钢筋混凝土框架结构为例,按照不同设防烈度要求设计3个结构模型,采用IDA方法对模型进行地震易损性分析和不同地震作用的可能性分析,评估结构震害风险,为抗震防灾对策的制定提供参考。
1. 结构易损性分析
根据地震影响区域、地震损伤数据采集方法和计算方法的不同,地震易损性评价的方法主要有经验方法和理论方法(Wu等,2012)。理论方法内容主要包括确定损伤判别条件、确定结构模型、选择和输入地震动、IDA分析、绘制易损性曲线等,主要通过计算数值反映结构各损伤状态与地震动强度之间的概率关系(孙柏涛等,2018)。IDA方法是目前性能化地震工程(Performance-Based Earthquake Engineering,PBEE)中评价结构地震反应、预测结构震害最常用的方法之一。该方法可对结构进行从弹性状态到弹塑性状态直至倒塌的全过程分析,能较好地反映结构在未来不同地震作用下的抗震能力变化规律(侯炜,2013)。因此,本文采用IDA方法进行结构地震易损性分析。
1.1 IDA指标选取
IDA分析结果由地震动强度指标IM(Intensity Measures)与结构性能参数指标DM(Damage Measures)之间的关系表示,目前国内外常用的地震动强度指标主要包括第一模态谱加速度Sa(T1,5%)和地震加速度峰值PGA(g),为了与地震危险性关联,本文选取PGA作为地震动强度指标。最大层间位移角${\theta _{\max }}$通过体现层间变形的综合结果在整体层面反映结构抗震性能(蒋欢军等,2009),因此,本文选取${\theta _{\max }}$作为结构性能参数指标。
1.2 结构破坏等级判别
建筑结构地震破坏等级随着破坏程度的加深分为基本完好、轻微破坏、中等破坏、严重破坏和倒塌。框架结构最大层间位移角${\theta _{\max }}$与地震破坏等级之间的对应关系如表 1所示(马东辉等,2007)。
表 1 最大层间位移角与地震破坏等级的关系Table 1. The relationship between the maximum inter-story drift ratios and damage degrees破坏等级 基本完好 轻微破坏 中等破坏 严重破坏 倒塌 ${\theta _{\max }}$ ${\theta _{\max }}$<1/400 1/400≤${\theta _{\max }}$<1/250 1/250≤${\theta _{\max }}$<1/125 1/125≤${\theta _{\max }}$<1/50 ${\theta _{\max }}$≥1/50 1.3 结构模型信息
本文以住宅中常用的内廊式结构为例进行研究,考虑所建模型应能从一定程度上反映按照不同设防标准设计和建造的一般框架结构抗震性能,采用较规则的平面布置,如图 1所示。按照相关规范(中华人民共和国住房和城乡建设部等,2010;中华人民共和国住房和城乡建设部,2012)要求,采用PKPM软件建立3个抗震设防烈度分别为6度(0.05g)、7度(0.1g)、8度(0.2g)的3跨10层框架结构模型。场地类别为Ⅱ类,设计地震分组为第二组,场地特征周期为0.40s。结构首层层高4.5m,其他层层高3.6m,建筑高度和抗震等级如表 2所示。恒荷载标准值取4.6kN/m2,活荷载标准值不上人屋面取0.5kN/m2、楼面取2.5kN/m2。
表 2 结构抗震等级Table 2. Seismic grade of structures项目 模型1 模型2 模型3 设防烈度 6度(0.05g) 7度(0.1g) 8度(0.2g) 建筑高度/m 36.9 36.9 36.9 抗震等级 三级 二级 一级 1.4 地震动选取和输入
不同地震动输入得到的地震作用效应值差异较大,有时可达数十倍(王亚勇等,1991)。如果建筑物所在场地具有实际强震记录,则为最优选择,但大部分情况下并不具备满足条件的强震记录。与美国地质勘测中心(United States Geological Survey,USGS)划分的S2场地基本相似,为使选用的地震动频谱特性尽可能与建筑物所在场地保持一致,从S2场地中选取30条具有代表性的强震记录(FEMA建议20条以上),将PGA调幅后作为地震动输入,如表 3所示。
表 3 地震动基本信息Table 3. Basic information of ground motion编号 名称 时间/年 台站 PGA/g 1 Livermore-02 1980 San Ramon-Eastman Kodak 0.191 2 Westmorland 1981 Parachute Test Site 0.219 3 Imperial Valley-06 1979 Delta 0.284 4 Northridge-06 1994 LA-Century City CC North 0.123 5 Whittier Narrows-01 1987 Downey-Co Maint Bldg 0.177 6 Whittier Narrows-01 1987 Brea Dam(Downstream) 0.231 7 Whittier Narrows-02 1987 LA-116th St School 0.157 8 Imperial Valley-02 1940 El Centro Array #9 0.342 9 Imperial Valley-06 1979 El Centro Array #7 0.42 10 Northridge-01 1994 LA-Century City CC North 0.223 11 Northridge-01 1994 LA-Obregon Park 0.467 12 Sierra Madre 1991 LA-Obregon Park 0.224 13 Superstition Hills-02 1987 Parachute Test Site 0.451 14 Loma Prieta 1989 Capitola 0.480 15 Victoria,Mexico 1980 LA-Obregon Par 0.118 16 Morgan Hill 1984 Capitola 0.117 17 San Fernando 1971 Carbon Canyon Dam 0.070 18 San Fernando 1971 Gormon-Oso Pump Plant 0.087 19 Imperial Valley-07 1979 Bonds Corner 0.092 20 Imperial Valley-07 1979 El Centro Array #7 0.161 21 Livermore-01 1980 San Ramon-Eastman Kodak 0.106 22 Victoria,Mexico 1980 Chihuahua 0.118 23 Whittier Narrows-02 1987 Inglewood-Union Oil 0.133 24 Landers 1992 Downey-Co Maint Bldg 0.046 25 Landers 1992 Inglewood-Union Oil 0.040 26 Northridge-06 1994 LA-Baldwin Hills 0.052 27 Hector Mine 1999 Downey-Co Maint Bldg 0.031 28 Hector Mine 1999 LA-116th St School 0.022 29 CA/Baja Border Area 2002 Calexico Fire Station 0.094 30 CA/Baja Border Area 2002 El Centro Array #7 0.078 1.5 易损性分析
统计和实际工程中认为结构能力参数C服从对数正态分布,假定结构反应参数D也服从对数正态分布(李玉,2015):
$$C = \ln (\mathop C\limits^ \wedge, {\beta _{\rm{c}}})$$ (1) $$D = \ln (\mathop D\limits^ \wedge, {\beta _{\rm{d}}})$$ (2) 式中,$\mathop C\limits^ \wedge $、$\mathop D\limits^ \wedge $分别表示结构能力和结构反应均值;${\beta _{\rm{c}}}$、${\beta _{\rm{d}}}$分别表示结构能力和结构反应对数标准差。
结构进行IDA分析后,可得最大层间位移角${\theta _{\max }}$(结构反应D)和$PGA$的对应关系(何益斌等,2016)为:
$$D = aPG{A^b}$$ (3) 上式两边分别取对数可得:
$$\ln D = \ln a + b\ln PGA = A + B\ln PGA$$ (4) 式中,a、b、A、B为回归系数,可根据结构IDA分析结果拟合得到。
易损性曲线表示不同地震作用下结构反应D超出某一破坏等级所定义的结构能力C的条件概率(贾晗曦等,2019),即结构失效概率${P_{\rm{f}}}$:
$${P_{\rm{f}}} = P(C/D \leqslant 1)$$ (5) 根据式(1)和式(2),已知D和C均服从对数正态分布,则结构失效概率可表示为:
$${P_{\rm{f}}} = \mathit{\Phi} \left[ {\frac{{\ln (\mathop D\limits^ \wedge /\mathop C\limits^ \wedge)}}{{\sqrt {\beta _{\rm{c}}^{\rm{2}} + \beta _{\rm{d}}^2} }}} \right]$$ (6) 式中,$\mathop C\limits^ \wedge $根据表 1中以最大层间位移角表示的结构性能量化指标取值;$\sqrt {\beta _d^{\rm{2}} + \beta _c^{\rm{2}}} $根据规定:当结构易损性曲线以地震峰值加速度$PGA$为自变量时取0.5,本文取0.5;$\mathit{\Phi}( \cdot )$为正态分布函数,其取值可通过查询标准分布表确定。
对30条地震动PGA进行调幅,分别调至0.2g、0.4g、0.6g、0.8g、1.0g、1.2g,得到结构性能参数(θmax)与地震动强度(PGA)之间的关系曲线,即IDA曲线簇,如图 2—4所示。
对结构模型进行IDA分析后,将横纵坐标取对数,拟合线性关系,并根据上述计算方法分别得到结构易损性曲线,如图 5—7所示。
由图 5—7可知,结构发生各种破坏状态的累积超越概率曲线均随着PGA的增大呈上升趋势,PGA越大,结构震害越严重;对于同一破坏状态的累积超越概率曲线,8度设防烈度结构上升趋势最缓,同一地震动强度作用下8度设防烈度结构发生严重破坏和倒塌的概率更小,说明8度设防烈度结构较6度和7度设防烈度结构具有更好的抗震能力。
2. 震害风险评价
2.1 地震烈度概率分布
根据高小旺等(1985)对华北、西北、西南地区的45个城镇地震危险性分析结果,一个地区在未来50年内发生的地震烈度服从极值Ⅲ型分布,分布函数为:
$$F(x) = \exp \left[ { - {{\left({\frac{{\omega - x}}{{\omega - \varepsilon }}} \right)}^K}} \right]$$ (7) 对分布函数求导可得概率密度函数$f(x)$,即某一场地未来50年内发生各地震烈度的概率密度:
$$f(x) = \frac{{K{{(\omega - x)}^{K - 1}}}}{{{{(\omega - \varepsilon)}^K}}}\exp \left[ { - {{\left({\frac{{\omega - x}}{{\omega - \varepsilon }}} \right)}^K}} \right]$$ (8) 式中,$x$表示地震烈度,为1—12的离散变量;$\omega $为烈度上限值,根据目前通用的地震烈度划分方法,$\omega $取12;$\varepsilon $为众值烈度,即未来50年内超越概率为63.2%的地震烈度;K为形状参数,取决于一个地区地震背景的复杂性,根据文献(吕大刚,1999),形状参数K的取值如表 4所示。由式(8)绘制不同设防烈度地区地震概率密度曲线,如图 8所示。
表 4 形状参数K的取值Table 4. Values of shape parameter K基本烈度 6度 7度 8度 K 9.7932 8.3339 6.8713 2.2 震害风险计算方法
结构震害风险可定义为:在考虑场地地震危险性的基础上,结构发生各种破坏状态的可能性,在数值上等于场地危险性和地震易损性的卷积。结构在未来50年内发生震害的概率P计算公式如下:
$$P = \sum\limits_{{\rm{all}}\;{x_i}} {P(C \leqslant D|PGA = {x_i})} P(PGA = {x_i})$$ (9) 式中,${x_i}$为未来可能遭遇的$PGA$。由于$PGA$的取值是连续的,因此式(9)可写为:
$$P{\rm{ = }}\int {P(C \leqslant D|PGA = {x_i})} f({x_i}){\rm{d}}x$$ (10) 式中,$f({x_i})$为$PGA$取值概率密度函数,与地震烈度概率密度函数$f(x)$一致。
2.3 基于蒙特卡罗方法的风险评价
蒙特卡罗方法是以概率和统计理论方法为基础的数值计算方法,也称统计模拟方法或随机抽样技术,通过使用随机数(或伪随机数)解决计算问题。本文根据蒙特卡罗思想产生大量符合地震概率密度的随机数,模拟未来场地可能遭遇的地震烈度。考虑积分计算的精确性,将地震烈度作为连续变量,以式(7)为母体,采用蒙特卡罗方法通过MATLAB软件对地震烈度进行n次抽样,并根据公式将地震烈度转化为PGA,则每个PGA在整体中出现的概率为$1/n$,产生的随机样本n越大,结果越接近真实值。本文取n=50000,50000个随机数即表示该地区可能发生的地震烈度值及其分布特征。当对不同设防烈度地区进行计算时,由于其分布特征不同,因而对结构造成的危险性也不同。
将抽样后得到的地震烈度I根据刘恢先等(1981)提出的公式转换为PGA:
$$ PGA = {10^{\left( {I \cdot \lg 2 - 0.01} \right)}} $$ (11) 根据大数定律,式(10)可写为:
$$P = \frac{1}{n}\sum\limits_{i = 1}^{50000} {P(C \leqslant D|PGA = {x_i})} $$ (12) 结构超越各级破坏的风险概率${P_j}$为:
$${P_j} = \sum\limits_{i = 1}^{50000} {\frac{{P(C \leqslant D|PGA = {x_i})}}{{50000}}} $$ (13) 式中,j=1、2、3、4、5,代表每级震害状态,分别对应基本完好、轻微破坏、中等破坏、严重破坏、倒塌。
同理,可知结构发生各级震害的风险${C_j}$为:
$${C_j} = \left\{ \begin{gathered} {P_j} - {P_{j + 1}}\;\;\;\;\;\;\;j \leqslant 4 \\ {P_j}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;j = 5 \\ \end{gathered} \right.$$ (14) 2.4 震害风险评价
依据上述分析模型,结合本文结构易损性分析结果,计算得到不同设防烈度框架结构50年震害风险结果,如表 5—7所示。
表 5 6度设防烈度框架结构震害风险评价结果(%)Table 5. Seismic risk assessment results of 6 intensity frame structure(%)震害等级 基本完好 轻微破坏 中等破坏 严重破坏 倒塌 超越风险 100 29.04 8.91 1.92 0.14 发生风险 70.96 20.13 6.99 1.78 0.14 表 6 7度设防烈度框架结构震害风险评价结果(%)Table 6. Seismic risk assessment results of 7 intensity frame structure(%)震害等级 基本完好 轻微破坏 中等破坏 严重破坏 倒塌 超越风险 100 42.34 16.09 4.25 0.43 发生风险 57.66 26.25 11.84 3.82 0.43 表 7 8度设防烈度框架震害风险评价结果(%)Table 7. Seismic risk assessment results of 8 intensity frame structure(%)震害等级 基本完好 轻微破坏 中等破坏 严重破坏 倒塌 超越风险 100 55.20 25.64 8.28 1.07 发生风险 44.80 29.56 17.36 7.21 1.07 由表 5—7可知,在3类设防烈度地区结构处于基本完好和轻微破坏状态风险的概率较大,分别为91.09%、83.91%和74.36%,发生严重破坏和倒塌状态风险的概率较低。由6度、7度、8度抗震设防烈度框架结构震害风险对比可知,8度设防烈度框架结构发生严重破坏和倒塌状态风险的概率最大,分别为7.21%、1.07%;7度设防烈度框架结构次之,发生严重破坏和倒塌状态风险的概率分别为3.28%、0.43%;而6度设防烈度框架结构发生严重破坏和倒塌状态风险的概率最低,分别为1.78%、0.14%。综合结构地震易损性结果可以看出,造成结构震害风险的因素除结构抗震能力外,还与其所处抗震设防烈度地区有关。
2.5 震害风险等级划分
为综合考虑结构不同震害等级的发生风险,并实现对不同设防烈度地区结构震害风险的排序,本文采用综合震害风险指数R表征结构震害风险平均程度,其数值范围0—1,结构综合震害风险指数越小,表示结构震害风险水平越低。参考震害指数(尹之潜等,1990)划分方法,建议各震害等级对应的综合震害风险指数范围如表 8所示。
表 8 综合震害风险指数划分Table 8. Composite seismic risk indexes classification震害等级 基本完好 轻微破坏 中等破坏 严重破坏 倒塌 综合震害风险指数范围 0<R1≤0.2 0.2<R2≤0.4 0.4<R3≤0.6 0.6<R4≤0.8 0.8<R5≤1 综合震害风险指数中值 0.1 0.3 0.5 0.7 0.9 结构综合震害风险指数计算公式如下:
$$R = \sum\limits_{i{\rm{ = 1}}}^5 {\sum\limits_{j = 1}^5 {{R_i}{C_j}} } $$ (15) 式中,Cj为结构发生各级震害的风险;Ri为各级震害风险指数中值,i=1、2、3、4、5。
结构综合震害风险指数由式(15)计算,既考虑结构各级震害风险的发生概率,又量化结构震害风险。
为了更直观地对建筑结构震害风险水平进行评定和分类,本文根据结构综合震害风险指数将结构震害风险等级划分为4级,由低到高分别为Ⅰ级、Ⅱ级、Ⅲ级、Ⅳ级,从Ⅰ级到Ⅳ级表示结构震害风险越来越严重,如表 9所示。
表 9 震害风险等级划分Table 9. Seismic risk grades classification综合震害风险指数 0<R≤0.25 0.25<R≤0.5 0.5<R≤0.75 0.75<R≤1 震害风险等级 Ⅰ级 Ⅱ级 Ⅲ级 Ⅳ级 结合表 5—7中不同设防烈度框架结构发生各震害等级的风险,按照式(15)进行计算,得到结构综合震害风险指数,进而根据表 9中综合震害风险指数范围判断结构震害风险等级,计算和判断结果如表 10所示。
表 10 结构综合震害风险指数和震害风险等级Table 10. Composite seismic risk indexes and risk grades of buildings结构 6度设防烈度框架结构 7度设防烈度框架结构 8度设防烈度框架结构 综合震害风险指数 0.18 0.22 0.28 震害风险等级 Ⅰ级 Ⅰ级 Ⅱ级 经过计算,6度、7度、8度设防烈度框架结构综合震害风险指数分别为0.18、0.22、0.28,可知8度设防烈度地区框架结构震害风险最高,7度设防烈度地区框架结构震害风险次之,6度设防烈度地区框架结构震害风险最低。根据表 9可知,6度、7度设防烈度框架结构综合震害风险指数为0—0.25,震害风险等级为Ⅰ级;8度设防烈度框架结构震害风险指数为0.25—0.5,震害风险等级为Ⅱ级。
3. 结论
1)进行建筑结构震害风险分析时,综合考虑结构地震易损性和地震危险性,通过对风险水平进行量化与评估,更全面和直观地反映建筑结构未来一段时间存在潜在破坏可能,对制定科学合理的抗震防灾对策更有意义。
2)在抗震能力方面,按照6度、7度、8度设防烈度设计的框架结构随着设防烈度的升高,抗震性能有所提高。但由于不同设防烈度地区地震危险性的差异,导致综合考虑地震危险性后结构震害风险趋势发生变化,如8度设防烈度地区结构震害风险等级为Ⅱ级,而6度和7度设防烈度地区结构震害风险等级为Ⅰ级。因此地震危险性在震害风险中起着不容忽视的作用,在城市抗震防灾规划中应注意高烈度地区建筑结构震害风险的防范。
3)由于建筑结构使用功能不同,使其平立面布置等具有较大差异,必然导致其在结构设计方面存在差异,难以通过单一模型表示。本文仅以3个3跨10层内廊式规则框架结构为例,评估其在不同设防烈度地区的震害风险,存在一定局限性。建筑结构模型的不确定性、地震活动的不确定性、认知不确定性及其他因素对风险的影响有待进一步研究。
-
表 1 系统内置的行业风险暴露数据示例数据
Table 1. Sample data for built-in industial exposure
编号 经度/°E 纬度/°N 行政区划代码 省 市 区/县 住宅建筑价值/元 商业建筑价值/元 工业建筑价值/元 1 101.35 21.825 532823 云南省 西双版纳傣族自治州 勐腊县 3000 0 0 2 107.5 21.825 451422 广西壮族自治区 崇左市 宁明县 1600 0 0 3 109.275 21.825 450521 广西壮族自治区 北海市 合浦县 35000 0 0 4 110.8 21.825 440981 广东省 茂名市 高州市 500000 1600000 5000000 5 100.025 21.85 532822 云南省 西双版纳傣族自治州 勐海县 4000 0 0 6 100.225 21.85 532822 云南省 西双版纳傣族自治州 勐海县 13000 0 0 7 101.05 21.85 532801 云南省 西双版纳傣族自治州 景洪市 5000 0 0 8 101.275 21.85 532823 云南省 西双版纳傣族自治州 勐腊县 9000 0 0 9 107.1 21.85 451422 广西壮族自治区 崇左市 宁明县 1000 0 0 … … … … … … … … … … 表 2 随机地震事件集示例数据
Table 2. Sample data for random earthquake event set
事件编号 年 日 经度/°E 纬度/°N 深度/km 走向/° 面波震级 衰减分区 100 000 000 401 1 36 121.975 23.458 56 70 5.96 0 100 000 000 405 1 40 100.158 22.167 17 120 6.22 1 100 000 000 496 1 131 83.618 43.670 18 150 6.59 2 100 000 000 531 1 166 89.597 46.134 5 125 6.07 2 100 000 000 564 1 199 122.674 23.533 18 70 6.05 0 100 000 000 581 1 216 79.013 30.368 12 154 6.1 1 100 000 000 649 1 284 121.814 23.905 3 160 6.23 0 100 000 000 707 1 342 79.605 30.388 30 135 6.41 1 100 000 000 756 2 26 120.934 24.669 13 20 6.14 0 100 000 000 768 2 38 79.201 37.062 28 170 5.98 1 100 000 000 854 2 124 78.154 35.726 8 101 6.5 1 100 000 000 860 2 130 105.807 38.450 8 160 6.42 3 … … … … … … … … … 表 3 事件损失表示例数据
Table 3. Sample data for ELT
随机事件集编号 年 总损失/元 毛损失/元 100 000 299 223 35 350 829 985.24 280 663 988.19 100 000 303 197 35 358 336 146 996.12 286 668 917 596.90 100 000 553 049 64 2 642 638 483.25 2 114 110 786.60 100 000 582 342 67 1 774 458.68 1 419 566.94 100 000 724 753 83 425 193.08 340 154.46 100 000 895 924 103 141 852 462 611.31 113 481 970 089.05 100 001 000 898 115 462 850 303.77 370 280 243.02 100 001 284 622 147 265 095 628.27 212 076 502.62 100 001 469 992 168 29 291 919.44 23 433 535.55 …… …… …… …… 表 4 年累积损失表示例数据
Table 4. Sample data for YLT
年 总损失/元 毛损失/元 35 358 686 976 981.36 286 949 581 585.09 64 2 642 638 483.25 2 114 110 786.60 67 1 774 458.68 1 419 566.94 83 425 193.08 340 154.46 103 141 852 462 611.31 113 481 970 089.05 115 462 850 303.77 370 280 243.02 147 265 095 628.27 212 076 502.62 168 29 291 919.44 23 433 535.55 …… …… …… -
[1] 陈洪富, 孙柏涛, 陈相兆等, 2013. HAZ-China地震灾害损失评估系统研究. 土木工程学报, 46(S2): 294—300Chen H. F., Sun B. T., Chen X. Z., et al., 2013. HAZ-China earthquake disaster loss estimation system. China Civil Engineering Journal, 46(S2): 294—300. (in Chinese) [2] 冯键, 2019. 巨灾数字化平台服务国家现代风险治理体系建设. 金融电子化, (5): 32—34 [3] 国家地震局, 国家统计局. 1996. 中国大陆地震灾害损失评估汇编: 1990~1995. 北京: 地震出版社. [4] 林向洋, 文鑫涛, 李华玥等, 2020. 2019年中国大陆地震灾害损失述评. 震灾防御技术, 15(3): 473—483 doi: 10.11899/zzfy20200301Lin X. Y., Wen X. T., Li H. Y., et al., 2020. Review of earthquake damage losses in mainland of China in 2019. Technology for Earthquake Disaster Prevention, 15(3): 473—483. (in Chinese) doi: 10.11899/zzfy20200301 [5] 楼宝棠, 1996. 中国古今地震灾情总汇. 北京: 地震出版社.Lou B. T., 1996. A comprchensive compilatton of historie and recent earthquakes diaster status in China. Beijing: Seismological Press. (in Chinese) [6] 孟生旺, 李政宵, 2018. 地震死亡人数预测与巨灾保险基金测算. 统计研究, 35(10): 89—102Meng S. W., Li Z. X., 2018. Modeling earthquake deaths and calculating catastrophe insurance fund. Statistical Research, 35(10): 89—102. (in Chinese) [7] 牛方曲, 高晓路, 季珏, 2012. 区域中长期房屋震灾损失评估系统. 资源科学, 34(2): 359—366Niu F. Q., Gao X. L., Ji Y., 2012. A house damage evaluation system of regional long-term seism in China. Resources Science, 34(2): 359—366. (in Chinese) [8] 熊政辉, 2019. 面向地震保险的巨灾模型研究与应用. 北京: 中国地震局地球物理研究所.Xiong Z. H., 2019. The study and application of catastrophe model for earthquake insurance. Beijing: Institute of Geophysics, China Earthquake Administration. (in Chinese) [9] 徐伟进, 高孟潭, 左惠强, 2020. 地震巨灾模型中的随机地震事件集模拟. 地震工程学报, 42(4): 1024—1034 doi: 10.3969/j.issn.1000-0844.2020.04.1024Xu W. J., Gao M. T., Zuo H. Q., 2020. Simulation of stochastic seismic event set using the earthquake catastrophe model. China Earthquake Engineering Journal, 42(4): 1024—1034. (in Chinese) doi: 10.3969/j.issn.1000-0844.2020.04.1024 [10] 应急管理部国家减灾中心, 2021. 2020年全国自然灾害基本情况. 北京: 应急管理部. [11] 俞言祥, 2016. 新一代地震区划图地震动参数衰减关系的建立与特点分析. 城市与减灾, (3): 34—38 doi: 10.3969/j.issn.1671-0495.2016.03.009 [12] 中国地震局监测预报司, 2001. 中国大陆地震灾害损失评估汇编: 1996~2000. 北京: 地震出版社. [13] 中国地震局震灾应急救援司, 2010. 2001~2005年中国大陆地震灾害损失评估汇编. 北京: 地震出版社. [14] 中国地震局震灾应急救援司, 2015. 2006~2010年中国大陆地震灾害损失评估汇编. 北京: 地震出版社. [15] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会, 2016. GB 18306—2015 中国地震动参数区划图. 北京: 中国标准出版社.General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, Standardization Administration of the People’s Republic of China, 2016. GB 18306—2015 Seismic ground motion parameters zonation map of China. Beijing: Standards Press of China. (in Chinese) [16] 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 2010. GB 50011—2010 建筑抗震设计规范. 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People’s Republic of China (MOHURD), 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) [17] Li X. J., Jing B. B., Liu C., et al., 2019. Site classification method based on geomorphological and geological characteristics and its application in China. Bulletin of the Seismological Society of America, 109(5): 1843—1854. doi: 10.1785/0120190058 -