CN109408965A - 基于地震动参数的房屋震害矩阵曲线化分析系统及方法 - Google Patents
基于地震动参数的房屋震害矩阵曲线化分析系统及方法 Download PDFInfo
- Publication number
- CN109408965A CN109408965A CN201811259732.7A CN201811259732A CN109408965A CN 109408965 A CN109408965 A CN 109408965A CN 201811259732 A CN201811259732 A CN 201811259732A CN 109408965 A CN109408965 A CN 109408965A
- Authority
- CN
- China
- Prior art keywords
- vulnerability
- motion parameter
- ground motion
- earthquake
- function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- E—FIXED CONSTRUCTIONS
- E04—BUILDING
- E04H—BUILDINGS OR LIKE STRUCTURES FOR PARTICULAR PURPOSES; SWIMMING OR SPLASH BATHS OR POOLS; MASTS; FENCING; TENTS OR CANOPIES, IN GENERAL
- E04H9/00—Buildings, groups of buildings or shelters adapted to withstand or provide protection against abnormal external influences, e.g. war-like action, earthquake or extreme climate
- E04H9/02—Buildings, groups of buildings or shelters adapted to withstand or provide protection against abnormal external influences, e.g. war-like action, earthquake or extreme climate withstanding earthquake or sinking of ground
Landscapes
- Engineering & Computer Science (AREA)
- Architecture (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Structural Engineering (AREA)
- Environmental & Geological Engineering (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Theoretical Computer Science (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于地震动参数的房屋震害矩阵曲线化分析方法及相应的系统，基于烈度与地震动参数之间的对应关系，采用最大似然估计方法，将正确的房屋易损性矩阵或实际震害统计的破坏比结果转化为双参数易损性曲线，弥补了基于烈度易损性曲线的不足，特别是针对土木、砖石等难以进行模拟分析的结构类型提供了更为简便的获得易损性曲线的途径；基于烈度与地震动参数的关系，给出各类房屋结构的地震易损性曲线特征参数，可为基于地震动参数的房屋建筑震害评估提供基础数据。该方法和系统可靠性高、适用广泛，可用于对已有的易损性矩阵或曲线进行校验和完善，以便更好地服务于地震灾害快速评估、风险评价和保险估算等领域。
Description
技术领域
本发明属于建筑信息处理技术领域，特别涉及一种基于地震动参数的房屋震害矩阵曲线化分析系统及方法。
背景技术
我国是世界上陆地地震灾害最严重的的国家之一，特别是随着经济建设和城市化进程的不断加快、地震灾害潜在风险水平不断提高。易损性分析作为灾害风险分析的3个主要部分(地震危险性分析、承灾体的易损性分析、灾情损失评估)之一，其核心是建立建筑结构出现或超越各种破坏状态的概率与地震强度之间的关系。开展建筑结构的地震易损性分析对于地震破坏和损失预测以及减轻地震灾害具有十分重要的意义。
在长期实际地震现场和震害预测工作中，研究人员积累了大量的基于烈度的建筑结构震害经验关系，即用烈度-震害程度的比例或概率构成的易损性矩阵，然而烈度的获得是基于人的感觉、器物反应、房屋为主的工程结构破坏以及地表破坏等4类宏观现象来评定的。这也意味着随着我国结构物抗震设防水平的提高，不同时期烈度本身的度量尺度也会发生变化，所以基于烈度的易损性矩阵具有较大的不确定性。然而以往大多数震害预测中关于房屋建筑易损性的评估，都是基于烈度进行，即表示结构在不同烈度下各破坏状态的概率。为了准确预测建筑结构的抗震性能，更好地进行结构抗震设计，完成对现有建筑物的加固和维修，加快地震保险在全国范围内的推行，基于烈度的易损性矩阵很难满足实际需要，特别是对于土木结构、砖石结构等简易房屋，难以通过结构模拟分析获取其易损性曲线。因此，完善现有的易损性矩阵，并将基于烈度的易损性矩阵转化为更为合理的基于地震动参数的易损性矩阵，成为本领域技术人员急需解决的技术问题。
发明内容
为了解决上述技术问题，本发明提供了一种基于地震动参数的房屋震害矩阵曲线化分析系统及方法。
本发明具体技术方案如下：
本发明提供了一种基于地震动参数的房屋震害矩阵曲线化分析系统，包括如下步骤：
资料收集统计模块，用于收集待研究区域内的历史震害数据和房屋建筑调查信息，并据此为各类建筑结构分别生成基于烈度的易损性矩阵；
易损性超越概率第一分析模块，用于根据所述基于烈度的易损性矩阵，为各类建筑结构分别生成基于烈度的易损性超越概率矩阵；
地震动参数获取模块，用于获取与每一级地震烈度对应的地震动参数，所述地震动参数包括峰值加速度和峰值速度中的至少一种；
易损性曲线双参数计算模块，用于采用对数累计正态分布函数构建基于地震动参数的易损性函数，通过最大似然估计算法，计算所述易损性函数的中位数θ和地震动强度的对数标准差β；
易损性超越概率第二分析模块，用于根据所述易损性函数中的中位数θ和对数标准差β，为各类建筑结构分别生成基于地震动参数的易损性超越概率曲线；
结构易损性曲线生成模块，用于根据所述基于地震动参数的易损性超越概率曲线，计算不同建筑结构在不同破坏状态下的损害概率，并据此为各类建筑结构分别生成基于地震动参数的结构易损性曲线。
进一步地，所述易损性曲线双参数计算模块包括如下部分：
易损性函数定义单元，用于采用对数累计正态分布函数构建基于地震动参数的易损性函数
式中：P(C|IM＝x)是地震动强度达到IM＝x时，建筑结构到达某一极限状态的概率；是标准累积正态分布函数；θ是易损性函数的中位数；β为lnIM的标准差；
密度函数定义单元，用于构建似然函数的密度函数；
对于超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
对于未超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
似然函数定义单元，用于根据不同的破坏状态有不同的中位数θ和对数标准差β，定义似然函数为：
式中：N为统计的建筑结构总数；xi表示第i栋建筑结构对应的加速度PGA；ti为伯努利事件，达到破坏状态时为1，否则为0；
双参数计算单元，用于通过公式计算所述易损性函数中的中位数θ和对数标准差β：
进一步地，所述易损性超越概率第二分析模块包括如下部分：
超越概率曲线绘制单元，用于根据所述易损性函数的中位数θ和地震动强度的对数标准差β，为各类建筑结构分别绘制基于地震动参数的易损性超越概率曲线；
校验单元，用于通过最大似然估计法拟合出易损性超越概率校核曲线，并对所述基于地震动参数的易损性超越概率曲线进行校验，如果曲线形状一致，则校验通过，认为所述基于地震动参数的易损性超越概率曲线合理、准确。
进一步地，所述结构易损性曲线生成模块包括如下部分：
损害概率计算单元，用于创建不同建筑结构在不同破坏状态下的损害概率的计算公式
式中：表示房屋建筑在DSk状态下的概率；表示极限破坏状态的概率；k＝1、2或3；
根据所述结构易损性超越概率曲线结果与所述损害概率的计算公式，得到不同破坏状态下的基于地震动参数的结构易损性矩阵；
结构易损性曲线绘制单元，用于根据所述基于地震动参数的结构易损性矩阵，绘制基于地震动参数的结构易损性曲线。
进一步地，所述系统还包括如下部分：
图表生成模块，用于根据各个模块生成的数据结果制作相应的图表。
本发明另一方面提供了一种基于地震动参数的房屋震害矩阵曲线化分析方法，包括如下步骤：
S1：收集待研究区域内的历史震害数据和房屋建筑调查信息，并据此为各类建筑结构分别生成基于烈度的易损性矩阵；
S2：根据所述基于烈度的易损性矩阵，为各类建筑结构分别生成基于基于烈度的易损性超越概率矩阵；
S3：获取与每一级地震烈度对应的地震动参数，所述地震动参数包括峰值加速度和峰值速度中的至少一种；
S4：采用对数累计正态分布函数构建基于地震动参数的易损性函数，通过最大似然估计算法，计算所述易损性函数的中位数θ和地震动强度的对数标准差β；
S5：根据计算得到的θ和β，为各类建筑结构分别生成基于地震动参数的易损性超越概率曲线；
S6：根据所述基于地震动参数的易损性超越概率曲线，计算不同建筑结构在不同破坏状态下的损害概率，并据此为各类建筑结构分别生成基于地震动参数的结构易损性曲线。
进一步地，步骤S4的具体方法如下：
S4.1：采用对数累计正态分布函数构建基于地震动参数的易损性函数
式中：P(C|IM＝x)是地震动强度达到IM＝x时，建筑结构到达某一极限状态的概率；是标准累积正态分布函数；θ是易损性函数的中位数；β为lnIM的标准差；
S4.2：对于超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
S4.3：对于未超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
S4.4：根据不同的破坏状态有不同的中位数θ和对数标准差β，定义似然函数为：
式中：N为统计的建筑结构总数；xi表示第i栋建筑结构对应的加速度PGA；ti为伯努利事件，达到破坏状态时为1，否则为0；
S4.5：通过公式计算所述易损性函数中的中位数θ和对数标准差β：
进一步地，步骤S5的具体方法如下：
根据所述易损性函数的中位数θ和地震动强度的对数标准差β，为各类建筑结构分别绘制基于地震动参数的易损性超越概率曲线；通过最大似然估计法拟合出易损性超越概率校核曲线，并对所述基于地震动参数的易损性超越概率曲线进行校验，如果曲线形状一致，则校验通过，认为所述基于地震动参数的易损性超越概率曲线合理、准确。
进一步地，步骤S6的具体方法如下：
S6.1：创建不同建筑结构在不同破坏状态下的损害概率的计算公式
式中：表示房屋建筑在DSk状态下的概率；表示极限破坏状态的概率；k＝1、2或3；
S6.2：根据所述结构易损性超越概率曲线结果，结合S6.1中的公式，得到不同破坏状态下的基于地震动参数的结构易损性矩阵，并据此绘制基于地震动参数的结构易损性曲线。
本发明的有益效果如下：本发明提供了一种基于地震动参数的房屋震害矩阵曲线化分析方法及相应的系统，基于烈度与地震动参数之间的对应关系，采用最大似然估计方法，将正确的房屋易损性矩阵或实际震害统计的破坏比结果转化为双参数易损性曲线，具有如下优点：
(1)弥补了基于烈度易损性曲线的不足，特别是针对难以进行模拟分析的结构类型，如土木结构、砖石结构等简易房屋，提供了更为简便的获得易损性曲线的途径；
(2)基于中国地震烈度表中烈度与地震动参数的关系，给出研究区域内各类房屋结构的地震易损性曲线特征参数，分析结果可为研究区域基于地震动参数的房屋建筑震害评估提供基础数据；
(3)该方法和系统具有较高的可靠性和适用性，在进一步研究工作中可应用该方法对已有的易损性矩阵或曲线进行校验和完善，以便更好地服务于地震灾害快速评估、风险评价和保险估算等领域。
附图说明
图1为实施例1所述的一种基于地震动参数的房屋震害矩阵曲线化分析方法的流程图；
图2为实施例1中Ⅲ区土木结构震害矩阵散点图；
图3为实施例1中Ⅲ区土木结构震害矩阵直方图；
图4为实施例1中基于烈度的易损性超越概率曲线不同破坏状态示意图；
图5为实施例1中基于地震动参数的Ⅲ区土木结构易损性超越概率曲线；
图6为实施例1中基于地震动参数的Ⅲ区土木结构震害矩阵拟合曲线；
图7为基于实施例1所述的基于地震动参数的房屋震害矩阵曲线化分析平台绘制的四川省分区域各类型房屋结构震害矩阵拟合曲线；
图8为实施例2所述的一种基于地震动参数的房屋震害矩阵曲线化分析系统的结构示意图。
具体实施方式
下面结合附图和以下实施例对本发明作进一步详细说明。
四川省地质构造复杂、地震活动频繁、房屋结构类型多样，历史上发生过多次7级以上特大地震。因此，以下实施例中以四川省为例，将其分为三个区：Ⅰ区为西部地区，Ⅱ区为中北部地区，Ⅲ区为中南部地区。采用本申请提供的系统和方法，对四川省各个区域内各类房屋建筑的分布特征和抗震能力进行研究，生成各区域各类建筑结构的基于地震动参数的易损性矩阵。
实施例1
如图1所示，本发明实施例1提供了一种基于地震动参数的房屋震害矩阵曲线化分析方法，包括如下步骤：
S1：收集待研究区域内的历史震害数据和房屋建筑调查信息，并据此为各类建筑结构分别生成基于烈度的易损性矩阵；
在四川省内选取204个调查点，共调查了23946栋房屋建筑，根据GB/T24335-2009的规定，将建筑物破坏等级分为5类：DSk(k＝0，1，2，3，4)分别表示“基本完好”、“轻微破坏”、“中等破坏”、“严重破坏”、“毁坏”五种破坏状态；LSs(s＝1，2，3，4)为5种破坏状态之间的阈值，分别代表“轻微破坏”、“中等破坏”、“严重破坏”、“毁坏”四种破坏状态的极限值；以土木结构为例，根据上述调查信息，计算土木结构建筑出现不同程度破坏的概率，生成基于烈度的易损性矩阵，并生成图表。
S2：根据基于烈度的易损性矩阵，为各类建筑结构分别生成基于基于烈度的易损性超越概率矩阵；
以Ⅲ区为例，通过土木结构建筑基于烈度的易损性矩阵，从破坏状态DSi的破坏概率计算出超越每个极限破坏状态LSi的累计概率，得到建筑结构的易损性超越概率矩阵(如表1所示)，其中i＝1、2、3、4或5，同时绘制生成基于烈度的易损性超越概率曲线。
表1 Ⅲ区土木结构超越概率矩阵
S3：获取与每一级地震烈度对应的地震动参数，地震动参数包括峰值加速度(PGA)和峰值速度(PGA)中的至少一种；
本实施例直接采用《中国地震烈度表》(GB/T17742-2008)给出的烈度与峰值加速度(PGA)和峰值速度(PGA)之间的关系，具体如表2所示；
表2 不同地震烈度对应的峰值加速度PGA和峰值速度PGV
S4：采用对数累计正态分布函数构建基于地震动参数的易损性函数，通过最大似然估计算法，计算易损性函数的中位数θ和地震动强度的对数标准差β，以峰值加速度PGA为例，具体过程如下：
S4.1：采用对数累计正态分布函数构建基于PGA的易损性函数
式中：P(C|IM＝x)是地震动强度达到IM＝x时，结构物到达某一极限状态的概率；是标准累积正态分布函数；θ是易损性函数的中位数；β为lnIM的标准差；
S4.2：对于超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
S4.3：对于未超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
S4.4：根据不同的破坏状态有不同的中位数θ和对数标准差β，定义似然函数为：
式中：N为统计的建筑结构总数；xi表示第i栋建筑结构对应的加速度PGA；ti为伯努利事件，达到破坏状态时为1，否则为0；
S4.5：通过公式计算易损性函数中的中位数θ和对数标准差β，结果如表3所示：
表3 Ⅲ区土木结构不同极限状态下超越概率易损性曲线中位数θ和对数标准差β
通过上述方法，可进一步计算出四川省内各个区域所有房屋建筑类型的基于PGA的θ和β，具体如表4所示。基于PGV的θ和β也可以通过同样方法计算得出，具体如表5所示。
表4 基于PGA的不同结构不同破坏状态下的易损性曲线中位数θ和对数标准差β
S5：根据易损性函数的中位数θ和地震动强度的对数标准差β，生成土木结构建筑的基于PGA的易损性超越概率矩阵，如表6所示，并绘制基于PGA的易损性超越概率曲线；通过最大似然估计法拟合出易损性超越概率校核曲线，并对基于地震动参数的易损性超越概率曲线进行校验，如果曲线形状一致，则校验通过，认为基于地震动参数的易损性超越概率曲线合理、准确。
S6：根据基于地震动参数的易损性超越概率曲线，计算不同建筑结构在不同破坏状态下的损害概率，生成土木结构易损性矩阵，如表7所示，并据此为各类建筑结构分别生成基于地震动参数的结构易损性曲线，如图6所示。
如图7所示，通过上述方法，可进一步生成四川省所有区域内所有类型房屋建筑的基于PGA的震害拟合曲线。
表5 基于PGV的不同结构不同破坏状态下的易损性曲线中位数θ和对数标准差β
表6 基于PGA的Ⅲ区土木结构超越概率易损性矩阵(％)
表7 基于PGA的Ⅲ区土木结构易损性矩阵(％)
实施例2
如图2所示，本发明实施例2提供了一种实现上述计算方法的基于地震动参数的房屋震害矩阵曲线化分析系统，包括如下部分：
资料收集统计模块1，用于收集待研究区域内的历史震害数据和房屋建筑调查信息，并据此为各类建筑结构分别生成基于烈度的易损性矩阵；
易损性超越概率第一分析模块2，用于根据基于烈度的易损性矩阵，为各类建筑结构分别生成基于烈度的易损性超越概率矩阵；
地震动参数获取模块3，用于获取与每一级地震烈度对应的地震动参数，地震动参数包括峰值加速度和峰值速度中的至少一种；
易损性曲线双参数计算模块4，用于采用对数累计正态分布函数构建基于地震动参数的易损性函数，通过最大似然估计算法，计算易损性函数的中位数θ和地震动强度的对数标准差β；
易损性超越概率第二分析模块5，用于根据易损性函数中的中位数θ和对数标准差β，为各类建筑结构分别生成基于地震动参数的易损性超越概率曲线；
结构易损性曲线生成模块6，用于根据基于地震动参数的易损性超越概率曲线，计算不同建筑结构在不同破坏状态下的损害概率，并据此为各类建筑结构分别生成基于地震动参数的结构易损性曲线；
图表生成模块7，用于根据各个模块生成的数据结果制作相应的图表。
易损性曲线双参数计算模块4包括如下部分：
易损性函数定义单元41，用于采用对数累计正态分布函数构建基于地震动参数的易损性函数
式中：P(C|IM＝x)是地震动强度达到IM＝x时，结构物到达某一极限状态的概率；是标准累积正态分布函数；θ是易损性函数的中位数；β为lnIM的标准差；
密度函数定义单元42，用于构建似然函数的密度函数；
对于超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
对于未超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
似然函数定义单元43，用于根据不同的破坏状态有不同的中位数θ和对数标准差β，定义似然函数为：
式中：N为统计的建筑结构总数；xi表示第i栋建筑结构对应的加速度PGA；ti为伯努利事件，达到破坏状态时为1，否则为0；
双参数计算单元44，用于通过公式计算易损性函数中的中位数θ和对数标准差β：
易损性超越概率第二分析模块5包括如下部分：
超越概率曲线绘制单元51，用于根据易损性函数的中位数θ和地震动强度的对数标准差β，为各类建筑结构分别绘制基于地震动参数的易损性超越概率曲线；
校验单元52，用于通过最大似然估计法拟合出易损性超越概率校核曲线，并对基于地震动参数的易损性超越概率曲线进行校验，如果曲线形状一致，则校验通过，认为基于地震动参数的易损性超越概率曲线合理、准确。
结构易损性曲线生成模块6包括如下部分：
损害概率计算单元61，用于创建不同建筑结构在不同破坏状态下的损害概率的计算公式
式中：表示房屋建筑在DSk状态下的概率；表示极限破坏状态的概率；k＝1、2或3；
根据结构易损性超越概率曲线结果与损害概率的计算公式，得到不同破坏状态下的基于地震动参数的结构易损性矩阵；
结构易损性曲线绘制单元62，用于根据基于地震动参数的结构易损性矩阵，绘制基于地震动参数的结构易损性曲线。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (9)
1.一种基于地震动参数的房屋震害矩阵曲线化分析系统，其特征在于，包括如下部分：
资料收集统计模块(1)，用于收集待研究区域内的历史震害数据和房屋建筑调查信息，并据此为各类建筑结构分别生成基于烈度的易损性矩阵；
易损性超越概率第一分析模块(2)，用于根据所述基于烈度的易损性矩阵，为各类建筑结构分别生成基于烈度的易损性超越概率矩阵；
地震动参数获取模块(3)，用于获取与每一级地震烈度对应的地震动参数，所述地震动参数包括峰值加速度和峰值速度中的至少一种；
易损性曲线双参数计算模块(4)，用于采用对数累计正态分布函数构建基于地震动参数的易损性函数，通过最大似然估计算法，计算所述易损性函数的中位数θ和地震动强度的对数标准差β；
易损性超越概率第二分析模块(5)，用于根据所述易损性函数中的中位数θ和对数标准差β，为各类建筑结构分别生成基于地震动参数的易损性超越概率曲线；
结构易损性曲线生成模块(6)，用于根据所述基于地震动参数的易损性超越概率曲线，计算不同建筑结构在不同破坏状态下的损害概率，并据此为各类建筑结构分别生成基于地震动参数的结构易损性曲线。
2.如权利要求1所述的基于地震动参数的房屋震害矩阵曲线化分析系统，其特征在于，所述易损性曲线双参数计算模块(4)包括如下部分：
易损性函数定义单元(41)，用于采用对数累计正态分布函数构建基于地震动参数的易损性函数
式中：P(C|IM＝x)是地震动强度达到IM＝x时，结构物到达某一极限状态的概率；是标准累积正态分布函数；θ是易损性函数的中位数；β为lnIM的标准差；
密度函数定义单元(42)，用于构建似然函数的密度函数；
对于超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
对于未超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
似然函数定义单元(43)，用于根据不同的破坏状态有不同的中位数θ和对数标准差β，定义似然函数为：
式中：N为统计的建筑结构总数；xi表示第i栋建筑结构对应的加速度PGA；ti为伯努利事件，达到破坏状态时为1，否则为0；
双参数计算单元(44)，用于通过公式计算所述易损性函数中的中位数θ和对数标准差β：
3.如权利要求2所述的基于地震动参数的房屋震害矩阵曲线化分析系统，其特征在于，所述易损性超越概率第二分析模块(5)包括如下部分：
超越概率曲线绘制单元(51)，用于根据所述易损性函数的中位数θ和地震动强度的对数标准差β，为各类建筑结构分别绘制基于地震动参数的易损性超越概率曲线；
校验单元(52)，用于通过最大似然估计法拟合出易损性超越概率校核曲线，并对所述基于地震动参数的易损性超越概率曲线进行校验，如果曲线形状一致，则校验通过，认为所述基于地震动参数的易损性超越概率曲线合理、准确。
4.如权利要求3所述的基于地震动参数的房屋震害矩阵曲线化分析系统，其特征在于，所述结构易损性曲线生成模块(6)包括如下部分：损害概率计算单元(61)，用于创建不同建筑结构在不同破坏状态下的损害概率的计算公式
式中：表示房屋建筑在DSk状态下的概率；表示极限破坏状态的概率；k＝1、2或3；
根据所述结构易损性超越概率曲线结果与所述损害概率的计算公式，得到不同破坏状态下的基于地震动参数的结构易损性矩阵；
结构易损性曲线绘制单元(62)，用于根据所述基于地震动参数的结构易损性矩阵，绘制基于地震动参数的结构易损性曲线。
5.如权利要求1～4中任一项所述的基于地震动参数的房屋震害矩阵曲线化分析系统，其特征在于，所述系统还包括如下部分：
图表生成模块(7)，用于根据各个模块生成的数据结果制作相应的图表。
6.一种基于地震动参数的房屋震害矩阵曲线化分析方法，其特征在于，包括如下步骤：
S1：收集待研究区域内的历史震害数据和房屋建筑调查信息，并据此为各类建筑结构分别生成基于烈度的易损性矩阵；
S2：根据所述基于烈度的易损性矩阵，为各类建筑结构分别生成基于基于烈度的易损性超越概率矩阵；
S3：获取与每一级地震烈度对应的地震动参数，所述地震动参数包括峰值加速度和峰值速度中的至少一种；
S4：采用对数累计正态分布函数构建基于地震动参数的易损性函数，通过最大似然估计算法，计算所述易损性函数的中位数θ和地震动强度的对数标准差β；
S5：根据计算得到的θ和β，为各类建筑结构分别生成基于地震动参数的易损性超越概率曲线；
S6：根据所述基于地震动参数的易损性超越概率曲线，计算不同建筑结构在不同破坏状态下的损害概率，并据此为各类建筑结构分别生成基于地震动参数的结构易损性曲线。
7.如权利要求6所述的基于地震动参数的房屋震害矩阵曲线化分析平台，其特征在于，步骤S4的具体方法如下：
S4.1：采用对数累计正态分布函数构建基于地震动参数的易损性函数
式中：P(C|IM＝x)是地震动强度达到IM＝x时，结构物到达某一极限状态的概率；是标准累积正态分布函数；θ是易损性函数的中位数；β为lnIM的标准差；
S4.2：对于超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
S4.3：对于未超越破坏状态的结构，构建地震动强度达到IM＝x处的似然函数的密度函数：
S4.4：根据不同的破坏状态有不同的中位数θ和对数标准差β，定义似然函数为：
式中：N为统计的建筑结构总数；xi表示第i栋建筑结构对应的加速度PGA；ti为伯努利事件，达到破坏状态时为1，否则为0；
S4.5：通过公式计算所述易损性函数中的中位数θ和对数标准差β：
8.如权利要求7所述的基于地震动参数的房屋震害矩阵曲线化分析平台，其特征在于，步骤S5的具体方法如下：
根据所述易损性函数的中位数θ和地震动强度的对数标准差β，为各类建筑结构分别绘制基于地震动参数的易损性超越概率曲线；通过最大似然估计法拟合出易损性超越概率校核曲线，并对所述基于地震动参数的易损性超越概率曲线进行校验，如果曲线形状一致，则校验通过，认为所述基于地震动参数的易损性超越概率曲线合理、准确。
9.如权利要求8所述的基于地震动参数的房屋震害矩阵曲线化分析平台，其特征在于，步骤S6的具体方法如下：
S6.1：创建不同建筑结构在不同破坏状态下的损害概率的计算公式
式中：表示房屋建筑在DSk状态下的概率；表示极限破坏状态的概率；k＝1、2或3；
S6.2：根据所述结构易损性超越概率曲线结果，结合S6.1中的公式，得到不同破坏状态下的基于地震动参数的结构易损性矩阵，并据此绘制基于地震动参数的结构易损性曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811259732.7A CN109408965B (zh) | 2018-10-26 | 2018-10-26 | 基于地震动参数的房屋震害矩阵曲线化分析系统及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811259732.7A CN109408965B (zh) | 2018-10-26 | 2018-10-26 | 基于地震动参数的房屋震害矩阵曲线化分析系统及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109408965A true CN109408965A (zh) | 2019-03-01 |
| CN109408965B CN109408965B (zh) | 2023-04-28 |
Family
ID=65469411
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811259732.7A Active CN109408965B (zh) | 2018-10-26 | 2018-10-26 | 基于地震动参数的房屋震害矩阵曲线化分析系统及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109408965B (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110210650A (zh) * | 2019-05-06 | 2019-09-06 | 同济大学 | 一种应急避难场所避难人口需求的预测方法及系统 |
| CN110309781A (zh) * | 2019-07-01 | 2019-10-08 | 中国科学院遥感与数字地球研究所 | 基于多尺度光谱纹理自适应融合的房屋损毁遥感识别方法 |
| CN110472350A (zh) * | 2019-08-20 | 2019-11-19 | 广西科技大学 | 一种基于虚拟仿真的地震震害评估方法 |
| CN112270125A (zh) * | 2020-10-23 | 2021-01-26 | 同济大学 | 一种基于有向图逻辑模型和蒙特卡洛模拟的系统地震易损性分析方法 |
| CN112329376A (zh) * | 2020-11-02 | 2021-02-05 | 同济大学 | 一种基于蒙特卡洛模拟的变电站系统抗震韧性量化评估算法 |
| CN112541666A (zh) * | 2020-12-08 | 2021-03-23 | 同济大学 | 考虑地震易损性模型不确定性的盾构隧道风险评估方法 |
| CN112613189A (zh) * | 2020-12-31 | 2021-04-06 | 中国再保险(集团)股份有限公司 | 一种考虑地震动幅值不确定性的建筑物易损性曲线建立方法及装置 |
| CN113434538A (zh) * | 2021-06-24 | 2021-09-24 | 北京市地震局 | 一种基于人口位置数据的地震强度评估方法及系统 |
| WO2023052570A1 (en) | 2021-09-29 | 2023-04-06 | Swiss Reinsurance Company Ltd. | Aerial and/or satellite imagery-based, optical sensory system and method for quantitative measurements and recognition of property damage after an occurred natural catastrophe event |
| CN117951921A (zh) * | 2024-03-26 | 2024-04-30 | 云南农业大学 | 一种边坡地震易损性评估方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003132296A (ja) * | 2001-10-29 | 2003-05-09 | Mitsui Constr Co Ltd | 建物の地震被害予測方法，それに用いる最適条件決定処理方法およびそれらのプログラム |
| CN106649954A (zh) * | 2016-10-08 | 2017-05-10 | 中冶华天工程技术有限公司 | 一种基于扩展pbee2理论框架的地震易损性分析方法 |
| CN108256141A (zh) * | 2017-12-11 | 2018-07-06 | 哈尔滨工业大学 | 一种基于Copula理论的主余震联合易损性分析方法 |
-
2018
- 2018-10-26 CN CN201811259732.7A patent/CN109408965B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003132296A (ja) * | 2001-10-29 | 2003-05-09 | Mitsui Constr Co Ltd | 建物の地震被害予測方法，それに用いる最適条件決定処理方法およびそれらのプログラム |
| CN106649954A (zh) * | 2016-10-08 | 2017-05-10 | 中冶华天工程技术有限公司 | 一种基于扩展pbee2理论框架的地震易损性分析方法 |
| CN108256141A (zh) * | 2017-12-11 | 2018-07-06 | 哈尔滨工业大学 | 一种基于Copula理论的主余震联合易损性分析方法 |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110210650A (zh) * | 2019-05-06 | 2019-09-06 | 同济大学 | 一种应急避难场所避难人口需求的预测方法及系统 |
| CN110210650B (zh) * | 2019-05-06 | 2021-07-02 | 同济大学 | 一种应急避难场所避难人口需求的预测方法及系统 |
| CN110309781A (zh) * | 2019-07-01 | 2019-10-08 | 中国科学院遥感与数字地球研究所 | 基于多尺度光谱纹理自适应融合的房屋损毁遥感识别方法 |
| CN110472350A (zh) * | 2019-08-20 | 2019-11-19 | 广西科技大学 | 一种基于虚拟仿真的地震震害评估方法 |
| CN112270125A (zh) * | 2020-10-23 | 2021-01-26 | 同济大学 | 一种基于有向图逻辑模型和蒙特卡洛模拟的系统地震易损性分析方法 |
| CN112329376B (zh) * | 2020-11-02 | 2022-06-07 | 同济大学 | 一种基于蒙特卡洛模拟的变电站系统抗震韧性量化评估方法 |
| CN112329376A (zh) * | 2020-11-02 | 2021-02-05 | 同济大学 | 一种基于蒙特卡洛模拟的变电站系统抗震韧性量化评估算法 |
| CN112541666A (zh) * | 2020-12-08 | 2021-03-23 | 同济大学 | 考虑地震易损性模型不确定性的盾构隧道风险评估方法 |
| CN112541666B (zh) * | 2020-12-08 | 2022-09-13 | 同济大学 | 考虑地震易损性模型不确定性的盾构隧道风险评估方法 |
| CN112613189B (zh) * | 2020-12-31 | 2023-12-26 | 中国再保险(集团)股份有限公司 | 一种考虑地震动幅值不确定性的建筑物易损性曲线建立方法及装置 |
| CN112613189A (zh) * | 2020-12-31 | 2021-04-06 | 中国再保险(集团)股份有限公司 | 一种考虑地震动幅值不确定性的建筑物易损性曲线建立方法及装置 |
| CN113434538A (zh) * | 2021-06-24 | 2021-09-24 | 北京市地震局 | 一种基于人口位置数据的地震强度评估方法及系统 |
| CN113434538B (zh) * | 2021-06-24 | 2024-03-08 | 北京市地震局 | 一种基于人口位置数据的地震强度评估方法及系统 |
| WO2023052570A1 (en) | 2021-09-29 | 2023-04-06 | Swiss Reinsurance Company Ltd. | Aerial and/or satellite imagery-based, optical sensory system and method for quantitative measurements and recognition of property damage after an occurred natural catastrophe event |
| CN117951921A (zh) * | 2024-03-26 | 2024-04-30 | 云南农业大学 | 一种边坡地震易损性评估方法 |
| CN117951921B (zh) * | 2024-03-26 | 2024-06-04 | 云南农业大学 | 一种边坡地震易损性评估方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109408965B (zh) | 2023-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109408965A (zh) | 基于地震动参数的房屋震害矩阵曲线化分析系统及方法 | |
| AU2007207834B2 (en) | A statistical-deterministic approach to natural disaster prediction | |
| CN106897510A (zh) | 一种桥梁结构三维地震易损性分析方法 | |
| CN108256141B (zh) | 一种基于Copula理论的主余震联合易损性分析方法 | |
| Lin et al. | Windborne debris risk assessment | |
| Awadallah | Selecting optimum locations of rainfall stations using kriging and entropy | |
| Mahmoud et al. | Framework for lifecycle cost assessment of steel buildings under seismic and wind hazards | |
| CN108898247A (zh) | 一种电网暴雨灾害风险预测方法、系统及存储介质 | |
| Arnalds et al. | Avalanche hazard zoning in Iceland based on individual risk | |
| Ciftci et al. | Analysis of the probability of failure for open-grown trees during wind storms | |
| CN106250590A (zh) | 一种基于cfd数值模拟的高空露天平台行人风环境评估方法 | |
| Liu et al. | Numerical study on factors influencing typhoon-induced storm surge distribution in Zhanjiang Harbor | |
| CN114218835A (zh) | 一种考虑风致疲劳效应的输电塔结构全寿命抗多灾性能评估方法 | |
| Takemi et al. | Large-eddy-simulation analysis of airflows and strong wind hazards in urban areas | |
| Yum et al. | Development of vulnerability curves of buildings to windstorms using insurance data: An empirical study in South Korea | |
| Wang et al. | Dynamic modeling of sea-level rise impact on coastal flood hazard and vulnerability in New York City's built environment | |
| CN105069299A (zh) | 一种事故时放射性核素大气扩散轨迹集合预测计算方法 | |
| CN107403004A (zh) | 一种基于地形数据的遥测雨量站点可疑数值检验方法 | |
| CN112232689A (zh) | 一种茶叶连片种植区的生态地质脆弱性评价方法 | |
| CN109784602A (zh) | 一种基于ptva模型的多灾种耦合物理脆弱性评估方法 | |
| Gilleland et al. | Evaluating NARCCAP model performance for frequencies of severe-storm environments | |
| Pan | Damage prediction of low-rise buildings under hurricane winds | |
| Mansouri et al. | Seismic human loss estimation in Tehran using GIS | |
| Song et al. | WARP^ 2: Wind assessment of roofs to pullout & pullover for priority cultural heritage structures in the Philippines | |
| Kellens et al. | A GIS tool for flood risk analysis in Flanders (Belgium) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |