• ISSN 1673-5722
  • CN 11-5429/P

基于多租户的高性能地震保险损失评估SaaS云平台设计与实现

王小东 粟春晨 宋圣楠 冯键 翟亮亮 熊政辉 马腾飞

杨贝贝, 王志涛, 张秀彦. 基于IDA方法的框架结构震害风险评估[J]. 震灾防御技术, 2020, 15(1): 21-32. doi: 10.11899/zzfy20200103
引用本文: 王小东,粟春晨,宋圣楠,冯键,翟亮亮,熊政辉,马腾飞,2021. 基于多租户的高性能地震保险损失评估SaaS云平台设计与实现. 震灾防御技术,16(2):253−262. doi:10.11899/zzfy20210204. doi: 10.11899/zzfy20210204
Yang Beibei, Wang Zhitao, Zhang Xiuyan. Seismic Risk Assessment of Frame Structures Based on IDA Method[J]. Technology for Earthquake Disaster Prevention, 2020, 15(1): 21-32. doi: 10.11899/zzfy20200103
Citation: Wang Xiaodong, Su Chunchen, Song Shengnan, Feng Jian, Zhai Liangliang, Xiong Zhenghui, Ma Tengfei. Design and Implementation of Multi-tenant SaaS Cloud Platform for High-Performance Seismic Insurance Loss Assessment[J]. Technology for Earthquake Disaster Prevention, 2021, 16(2): 253-262. doi: 10.11899/zzfy20210204

基于多租户的高性能地震保险损失评估SaaS云平台设计与实现

doi: 10.11899/zzfy20210204
基金项目: 国家重点研发计划(2018YFC1504604)
详细信息
    作者简介:

    王小东,男,生于1988年。硕士。主要从事巨灾损失评估模型研发方面的工作。E-mail:wxd_mail@163.com

    通讯作者:

    冯键,男,生于1972年。博士。主要从事企业数字化转型方面的工作。E-mail:fengj@chinare.com.cn

Design and Implementation of Multi-tenant SaaS Cloud Platform for High-Performance Seismic Insurance Loss Assessment

  • 摘要: 地震巨灾风险的特点是低频高损,历史震害数据缺乏、风险暴露快速变迁等因素导致基于大数定理的费率厘定方法无法针对各区域不同建筑类型的风险暴露进行精细化定价。本文基于“五代图”潜在震源区模型的随机事件集解决观测数据不足的问题;并使用“五代图”所采用的地震动参数衰减关系模型与工程易损性方法计算地震事件对风险暴露造成的损失,从而计算费率厘定、地震风险管理需要的必备参数。本文采用云计算平台的弹性伸缩计算技术,实现动态按需分配计算资源,满足多用户并发使用的业务需求;同时采用以业务数据为单元的数据隔离方案,构建支持多租户的高性能地震保险损失评估SaaS云平台。
  • 地震是一种极具破坏性的灾难,全世界地震多发的国家都因地震造成不可弥补的损失,人民生命与财产安全受到严重威胁。对建筑结构采取抗震设防措施是减轻地震灾害损失的根本手段,而地震区划是工程抗震设防的基础和前提条件。现阶段超过设防水平的大地震的发生不能忽视,如唐山地震(1978年)、汶川地震(2008年)等多次地震中,地震实际烈度往往超出抗震设计规范的设防烈度(陈肇元,2008),造成大面积楼房倒塌,数百万人流离失所。因此,灾前对既有建筑结构开展震害评估,并采取相应防灾措施尤为重要。传统的地震预测分析方法通常指对结构在特定地震作用下发生不同破坏状态的可能性进行分析,由于未充分考虑地震活动的不确定性,不能较好地反映建筑结构潜在的地震风险水平。Preciado等(2015)进行建筑结构震害风险评估时联合结构易损性和场地危险性,并以二者卷积为表征,表示结构在未来某一段时间内超过或发生某种地震破坏状态的概率,从而明确度量建筑结构未来可能发生的风险水平,为做好灾前预防、地震应急准备和灾后救援工作提供科学的决策依据,也是贯彻“从减少灾害损失向减轻灾害风险转变”(罗华春,2017)思想的体现。

    在长期的实践探索中,国内外学者在结构震害风险研究方面取得了一系列有价值的成果。Melani等(2016)采用IDA方法对某低层钢筋混凝土框架结构开展易损性分析,并基于框架结构概率风险四次积分公式,对其进行地震风险分析和年损失估计。Berto等(2012)基于强震环境对钢筋混凝土构件的影响研究,提出概率风险管理方法,计算与建筑结构长期性能相关的地震风险。谷音等(2012)针对某典型矮塔斜拉桥,采用LHS-MC方法对其地震风险概率进行评估。尹超等(2016)以实际路堤结构为例,采用蒙特卡罗抽样方法评价有无挡土墙路堤的震害风险,提出基于风险可接受度的路堤抗震设计与管理方法。然而,城市中存在的大量钢筋混凝土框架结构震害风险尚未可知,对不同抗震设防烈度区框架风险等级的划分仍处于研究阶段。

    基于上述认识,综合考虑我国地域辽阔且国土范围内地震危险性差异大的实际情况(马玉宏等,2002)和目前国内最主要的钢筋混凝土建筑结构形式,本文以钢筋混凝土框架结构为例,按照不同设防烈度要求设计3个结构模型,采用IDA方法对模型进行地震易损性分析和不同地震作用的可能性分析,评估结构震害风险,为抗震防灾对策的制定提供参考。

    根据地震影响区域、地震损伤数据采集方法和计算方法的不同,地震易损性评价的方法主要有经验方法和理论方法(Wu等,2012)。理论方法内容主要包括确定损伤判别条件、确定结构模型、选择和输入地震动、IDA分析、绘制易损性曲线等,主要通过计算数值反映结构各损伤状态与地震动强度之间的概率关系(孙柏涛等,2018)。IDA方法是目前性能化地震工程(Performance-Based Earthquake Engineering,PBEE)中评价结构地震反应、预测结构震害最常用的方法之一。该方法可对结构进行从弹性状态到弹塑性状态直至倒塌的全过程分析,能较好地反映结构在未来不同地震作用下的抗震能力变化规律(侯炜,2013)。因此,本文采用IDA方法进行结构地震易损性分析。

    IDA分析结果由地震动强度指标IM(Intensity Measures)与结构性能参数指标DM(Damage Measures)之间的关系表示,目前国内外常用的地震动强度指标主要包括第一模态谱加速度Sa(T1,5%)和地震加速度峰值PGA(g),为了与地震危险性关联,本文选取PGA作为地震动强度指标。最大层间位移角${\theta _{\max }}$通过体现层间变形的综合结果在整体层面反映结构抗震性能(蒋欢军等,2009),因此,本文选取${\theta _{\max }}$作为结构性能参数指标。

    建筑结构地震破坏等级随着破坏程度的加深分为基本完好、轻微破坏、中等破坏、严重破坏和倒塌。框架结构最大层间位移角${\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
    下载: 导出CSV 
    | 显示表格

    本文以住宅中常用的内廊式结构为例进行研究,考虑所建模型应能从一定程度上反映按照不同设防标准设计和建造的一般框架结构抗震性能,采用较规则的平面布置,如图 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

    图 1  结构平面布置(单位:mm)
    Figure 1.  The plan layout of structure(unit: mm)
    表 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
    抗震等级 三级 二级 一级
    下载: 导出CSV 
    | 显示表格

    不同地震动输入得到的地震作用效应值差异较大,有时可达数十倍(王亚勇等,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
    下载: 导出CSV 
    | 显示表格

    统计和实际工程中认为结构能力参数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)

    式中,abAB为回归系数,可根据结构IDA分析结果拟合得到。

    易损性曲线表示不同地震作用下结构反应D超出某一破坏等级所定义的结构能力C的条件概率(贾晗曦等,2019),即结构失效概率${P_{\rm{f}}}$:

    $${P_{\rm{f}}} = P(C/D \leqslant 1)$$ (5)

    根据式(1)和式(2),已知DC均服从对数正态分布,则结构失效概率可表示为:

    $${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曲线簇,如图 24所示。

    图 2  6度设防烈度框架结构IDA分析结果
    Figure 2.  IDA analysis results of the 6 intensity frame structure
    图 3  7度设防烈度框架结构IDA分析结果
    Figure 3.  IDA analysis results of the 7 intensity frame structure
    图 4  8度设防烈度框架结构IDA分析结果
    Figure 4.  IDA analysis results of the 8 intensity frame structure

    对结构模型进行IDA分析后,将横纵坐标取对数,拟合线性关系,并根据上述计算方法分别得到结构易损性曲线,如图 57所示。

    图 5  6度设防烈度框架结构易损性曲线
    Figure 5.  Fragility curves of the 6 intensity frame structure
    图 6  7度设防烈度框架结构易损性曲线
    Figure 6.  Fragility curves of the 7 intensity frame structure
    图 7  8度设防烈度框架结构易损性曲线
    Figure 7.  Fragility curves of the 8 intensity frame structure

    图 57可知,结构发生各种破坏状态的累积超越概率曲线均随着PGA的增大呈上升趋势,PGA越大,结构震害越严重;对于同一破坏状态的累积超越概率曲线,8度设防烈度结构上升趋势最缓,同一地震动强度作用下8度设防烈度结构发生严重破坏和倒塌的概率更小,说明8度设防烈度结构较6度和7度设防烈度结构具有更好的抗震能力。

    根据高小旺等(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
    下载: 导出CSV 
    | 显示表格
    图 8  地震烈度概率密度曲线
    Figure 8.  Curves of probabilistic densities of seismic intensities

    结构震害风险可定义为:在考虑场地地震危险性的基础上,结构发生各种破坏状态的可能性,在数值上等于场地危险性和地震易损性的卷积。结构在未来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)$一致。

    蒙特卡罗方法是以概率和统计理论方法为基础的数值计算方法,也称统计模拟方法或随机抽样技术,通过使用随机数(或伪随机数)解决计算问题。本文根据蒙特卡罗思想产生大量符合地震概率密度的随机数,模拟未来场地可能遭遇的地震烈度。考虑积分计算的精确性,将地震烈度作为连续变量,以式(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)

    依据上述分析模型,结合本文结构易损性分析结果,计算得到不同设防烈度框架结构50年震害风险结果,如表 57所示。

    表 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
    下载: 导出CSV 
    | 显示表格
    表 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
    下载: 导出CSV 
    | 显示表格
    表 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
    下载: 导出CSV 
    | 显示表格

    表 57可知,在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%。综合结构地震易损性结果可以看出,造成结构震害风险的因素除结构抗震能力外,还与其所处抗震设防烈度地区有关。

    为综合考虑结构不同震害等级的发生风险,并实现对不同设防烈度地区结构震害风险的排序,本文采用综合震害风险指数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
    下载: 导出CSV 
    | 显示表格

    结构综合震害风险指数计算公式如下:

    $$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
    震害风险等级 Ⅰ级 Ⅱ级 Ⅲ级 Ⅳ级
    下载: 导出CSV 
    | 显示表格

    结合表 57中不同设防烈度框架结构发生各震害等级的风险,按照式(15)进行计算,得到结构综合震害风险指数,进而根据表 9中综合震害风险指数范围判断结构震害风险等级,计算和判断结果如表 10所示。

    表 10  结构综合震害风险指数和震害风险等级
    Table 10.  Composite seismic risk indexes and risk grades of buildings
    结构 6度设防烈度框架结构 7度设防烈度框架结构 8度设防烈度框架结构
    综合震害风险指数 0.18 0.22 0.28
    震害风险等级 Ⅰ级 Ⅰ级 Ⅱ级
    下载: 导出CSV 
    | 显示表格

    经过计算,6度、7度、8度设防烈度框架结构综合震害风险指数分别为0.18、0.22、0.28,可知8度设防烈度地区框架结构震害风险最高,7度设防烈度地区框架结构震害风险次之,6度设防烈度地区框架结构震害风险最低。根据表 9可知,6度、7度设防烈度框架结构综合震害风险指数为0—0.25,震害风险等级为Ⅰ级;8度设防烈度框架结构震害风险指数为0.25—0.5,震害风险等级为Ⅱ级。

    1)进行建筑结构震害风险分析时,综合考虑结构地震易损性和地震危险性,通过对风险水平进行量化与评估,更全面和直观地反映建筑结构未来一段时间存在潜在破坏可能,对制定科学合理的抗震防灾对策更有意义。

    2)在抗震能力方面,按照6度、7度、8度设防烈度设计的框架结构随着设防烈度的升高,抗震性能有所提高。但由于不同设防烈度地区地震危险性的差异,导致综合考虑地震危险性后结构震害风险趋势发生变化,如8度设防烈度地区结构震害风险等级为Ⅱ级,而6度和7度设防烈度地区结构震害风险等级为Ⅰ级。因此地震危险性在震害风险中起着不容忽视的作用,在城市抗震防灾规划中应注意高烈度地区建筑结构震害风险的防范。

    3)由于建筑结构使用功能不同,使其平立面布置等具有较大差异,必然导致其在结构设计方面存在差异,难以通过单一模型表示。本文仅以3个3跨10层内廊式规则框架结构为例,评估其在不同设防烈度地区的震害风险,存在一定局限性。建筑结构模型的不确定性、地震活动的不确定性、认知不确定性及其他因素对风险的影响有待进一步研究。

  • 图  1  地震保险损失评估SaaS云平台技术架构

    Figure  1.  Technical architecture

    图  2  任务队列的计算引擎调度策略

    Figure  2.  Engine scheduling strategy based on task queue

    图  3  地震保险损失评估SaaS云平台功能架构

    Figure  3.  Functional architecture

    图  4  多种损失评估模式

    Figure  4.  Multiple options for loss assessment

    图  5  损失汇总及展示

    Figure  5.  Loss aggregate and display

    图  6  费率厘定模块

    Figure  6.  Module of rate making

    图  7  地震保险损失评估SaaS云平台业务架构示意图

    Figure  7.  Business architecture

    表  1  系统内置的行业风险暴露数据示例数据

    Table  1.   Sample data for built-in industial exposure

    编号经度/°E纬度/°N行政区划代码区/县住宅建筑价值/元商业建筑价值/元工业建筑价值/元
    1101.3521.825532823云南省西双版纳傣族自治州勐腊县300000
    2107.521.825451422广西壮族自治区崇左市宁明县160000
    3109.27521.825450521广西壮族自治区北海市合浦县3500000
    4110.821.825440981广东省茂名市高州市50000016000005000000
    5100.02521.85532822云南省西双版纳傣族自治州勐海县400000
    6100.22521.85532822云南省西双版纳傣族自治州勐海县1300000
    7101.0521.85532801云南省西双版纳傣族自治州景洪市500000
    8101.27521.85532823云南省西双版纳傣族自治州勐腊县900000
    9107.121.85451422广西壮族自治区崇左市宁明县100000
    下载: 导出CSV

    表  2  随机地震事件集示例数据

    Table  2.   Sample data for random earthquake event set

    事件编号经度/°E纬度/°N深度/km走向/°面波震级衰减分区
    100 000 000 401136121.97523.45856705.960
    100 000 000 405140100.15822.167171206.221
    100 000 000 496113183.61843.670181506.592
    100 000 000 531116689.59746.13451256.072
    100 000 000 5641199122.67423.53318706.050
    100 000 000 581121679.01330.368121546.11
    100 000 000 6491284121.81423.90531606.230
    100 000 000 707134279.60530.388301356.411
    100 000 000 756226120.93424.66913206.140
    100 000 000 76823879.20137.062281705.981
    100 000 000 854212478.15435.72681016.51
    100 000 000 8602130105.80738.45081606.423
    下载: 导出CSV

    表  3  事件损失表示例数据

    Table  3.   Sample data for ELT

    随机事件集编号总损失/元毛损失/元
    100 000 299 22335350 829 985.24280 663 988.19
    100 000 303 19735358 336 146 996.12286 668 917 596.90
    100 000 553 049642 642 638 483.252 114 110 786.60
    100 000 582 342671 774 458.681 419 566.94
    100 000 724 75383425 193.08340 154.46
    100 000 895 924103141 852 462 611.31113 481 970 089.05
    100 001 000 898115462 850 303.77370 280 243.02
    100 001 284 622147265 095 628.27212 076 502.62
    100 001 469 99216829 291 919.4423 433 535.55
    ……………………
    下载: 导出CSV

    表  4  年累积损失表示例数据

    Table  4.   Sample data for YLT

    总损失/元毛损失/元
    35358 686 976 981.36286 949 581 585.09
    642 642 638 483.252 114 110 786.60
    671 774 458.681 419 566.94
    83425 193.08340 154.46
    103141 852 462 611.31113 481 970 089.05
    115462 850 303.77370 280 243.02
    147265 095 628.27212 076 502.62
    16829 291 919.4423 433 535.55
    ………………
    下载: 导出CSV
  • [1] 陈洪富, 孙柏涛, 陈相兆等, 2013. HAZ-China地震灾害损失评估系统研究. 土木工程学报, 46(S2): 294—300

    Chen 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/zzfy20200301

    Lin 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—102

    Meng 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—366

    Niu 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.1024

    Xu 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
  • 加载中
图(7) / 表(4)
计量
  • 文章访问数:  193
  • HTML全文浏览量:  47
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-04-12
  • 刊出日期:  2021-06-30

目录

/

返回文章
返回