申请/专利权人：中铁二院贵阳勘察设计研究院有限责任公司

申请日：2017-08-23

公开（公告）日：2020-05-26

公开（公告）号：CN107292545B

主分类号：G06Q10/06(20120101)

分类号：G06Q10/06(20120101);G06Q50/26(20120101)

优先权：

专利状态码：有效-授权

法律状态：2020.05.26#授权;2017.11.24#实质审查的生效;2017.10.24#公开

摘要：本发明主要采用地震动强度指标中谱加速度Sa，引入PupIMi为桥梁结构在强度为IM时的地震的损伤概率，ΔλH,i为各地震设防水准下地震危险性贡献参数，利用PupIMi和ΔλH,i计算出的桥梁结构的地震风险值λUP，对桥梁结构的抗震性能进行评估。通过本方法求解出桥梁各危险构件的地震风险值，根据结构风险评估相关规范，研究桥梁各构件的抗震性能，还可以对比研究不同桥梁减隔震措施下各构件的地震风险值，评估不同桥梁减隔震措施的可行性和适用性。

主权项：1.一种基于地震风险评估的桥梁抗震分析方法，其特征在于，包括如下步骤：1根据桥址处的历史地震的记录数据和对应的地震参数，建立地震动强度预测方程组，计算出地震动谱加速度，结合地震设防水准，每种地震水准对应一个地震动强度和这种强度地震的年平均发生概率，将所有水准下的值连成线得到桥址处的地震危险性曲线；2根据所述地震危险性曲线，计算不同地震设防水准下的地震危险性贡献参数；3建立地震动-桥梁分析样本库并开展时程分析，提取结构最大动力响应数据，确定构件损伤指标及其量化方法，计算特定强度地震作用下的桥梁结构的地震损伤破坏概率，建立桥梁结构损伤破坏概率随地震动强度变化的关系曲线；4根据桥梁结构的地震损伤破坏概率和地震危险性贡献参数计算桥梁结构地震风险值。

全文数据：一种基于地震风险评估的桥梁抗震分析方法技术领域[0001]本发明涉及桥梁结构抗震分析领域，特别涉及一种基于地震风险评估的桥梁抗震分析方法。背景技术[0002]我国是世界上地震活动最强烈的国家之一，地震灾害历来是我国最大的自然灾害之一，最近几十年来，全球发生的多次破坏性地震，都造成了非常严重的生命财产损失，尤其是作为公路、铁路交通网络系统关键点的桥梁，是交通生命线系统中的重要枢纽结构，也是地震作用下相对容易损伤的基础设施，桥梁结构在地震中遭到了严重破坏，切断了震区交通生命线，造成救灾工作的巨大困难，使次生灾害加重，从而导致了非常大的损失。[0003]地震风险分析就是利用地震危险性的分析结果，结合对受灾物体和有关社会累计的知识，对特定区域、特定结构物可能遭受的地震灾害程度进行量化分析，地震风险一般是通过结构地震损伤曲线和工程场地危险曲线积分求解得到，但这两种曲线的具体表达式往往非常复杂，直接积分求解往往会存在很多的问题，为了将地震风险评估应用于桥梁抗震性能分析中，因此急需提出新的分析方法。发明内容[0004]本发明所要解决的技术问题在于针对现有技术的不足，提供一种基于地震风险评估的桥梁结构抗震分析方法，实现桥梁单个构件或整体结构在地震作用下地震风险评估、桥梁减隔震装置的适应性评估。[0005]为了实现上述发明目的，本发明提供了以下技术方案：[0006]—种基于地震风险评估的桥梁抗震分析方法，包括如下步骤：[0007]1根据桥址处的历史地震的记录数据和对应的地震参数，建立地震动强度预测方程组，计算出地震动谱加速度，结合地震设防水准，每种地震水准对应一个地震动强度和这种强度地震的年平均发生概率，将所有水准下的值连成线得到桥址处的地震危险性曲线；[0009]其中，SaTj是工程场地的地震动在特定周期Tj的谱加速度预测值;M为地震矩等级;rrup为是距离断层面的直线距离，单位为km;F表示断层类型的系数，其中正断层为1，逆断层为〇.5，其他为0;Hff是地震动上盘效应的系数，其中1为表示为上盘效应，0为其它效应；S是地质条件系数，其中0为岩石土体或浅层土体，1为深层土体；是岩体峰值地面加速度PGA的预期值;〇为预测方程的标准差，函数M，rrup是地震动强度预测函数的基本函数，断层类型函数f2M描述断层类型与滑移效益之间地震动的差异性，函数f3M，r_用于考虑上盘效应，场地响应函数用于考虑土层的非线性响应。[0010]进一步的，桥梁构件在特定地震作用下的损伤破坏概率的计算方法为：[0011][0012]其中，C表示桥梁地震响应的实际需求值;D表示桥梁构件抵抗地震作用的能力值；頂表示地震动强度值；頂i表示第i条地震动对应的强度值;im表示地震动强度的具体值;hM表示地震动的带宽参数;H表示带宽矩阵，IhI表示对η矩阵求行列式“·）τ表示矩阵的转置；ΐ中，X表示[0013]进一步的，步骤⑷包含以下步骤：[00Μ]4-1计算各地震设防水准下的桥梁结构的地震损伤破坏概率和地震危险性贡献参数的乘积；[0015]4-2将各地震设防水准下的乘积值按公式求和，得到桥梁结构综合地震风险值，其中，λυΡ为桥梁结构的地震风险值;M1为第i组设防水准下结构特征周期对应的地震动强度指标，PlipM1表示强度为IM1的地震作用下桥梁结构发生某种破坏状态所对应的发生概率，AλΗ1*各地震设防水准下地震危险性贡献参数。[0016]与现有技术相比，本发明的有益效果:本发明基于国内外桥梁抗震设计理论、地震风险评估理论、地震工程学、概率学与数理统计理论，建立了工程场地地震强度预测函数，提出了工程场地地震危险性贡献参数计算方法，将核密度估计引入到桥梁结构地震损伤分析中，建立了桥梁结构地震损伤破坏概率计算理论，在此基础上，建立了一种新的地震风险评估方法，进来评估了地震作用下桥梁结构的抗震性能，实现桥梁单个构件或整体结构在地震作用下地震风险评估，还可以对比研究不同桥梁减隔震措施下各构件的地震风险值，评估不同桥梁减隔震措施的可行性和适用性，具有较好的应用价值。附图说明[0017]图1为一种基于地震风险评估的桥梁抗震分析方法的流程图。[0018]图2为实施例1中的地震危险性曲线示意图。[0019]图3为计算特定强度地震作用下的桥梁结构地震损伤破坏率流程图。[0020]图4为地震动强度与桥梁结构损伤破坏率的曲线图。具体实施方式[0021]下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范围。[0022]实施例1[0023]—种基于地震风险评估的桥梁抗震分析方法，如图1，包括如下步骤：[0024]1根据桥址处的历史地震的记录数据和对应的地震参数，建立地震动强度预测方程组，计算出地震动谱加速度，结合地震设防水准，每种地震水准对应一个地震动强度和这种强度地震的年平均发生概率，将所有水准下的值连成线得到桥址处的地震危险性曲线；[0025]2根据所述地震危险性曲线，计算不同地震设防水准下的地震危险性贡献参数；[0026]3建立地震动-桥梁分析样本库并开展时程分析，提取结构最大动力响应数据，确定构件损伤指标及其量化方法，计算特定强度地震作用下的桥梁结构的地震损伤破坏概率，建立桥梁结构损伤破坏概率随地震动强度变化的关系曲线；[0027]4根据桥梁结构的地震损伤破坏概率和地震危险性贡献参数计算桥梁结构地震风险值。[0028]首先建立地震动强度预测方程组，即利用非参数估计的方法，将实测场地的地震动强度与地震震源的特性，如震级、断层距、断层类型、地震动上盘效应和地质条件系数关联起来，通过公式：[0029][0030]计算桥址处在特征周期为Tj条件下的地震动谱加速度SaTj，其中，M为地震矩等级;rrup为是距离断层面的直线距离，单位为km;F表示断层类型的系数正断层为1，逆断层为0.5,其他为0;HW是地震动上盘效应的系数（1为表示为上盘效应，0为其它效应）；S是地质条件系数〇为岩石土体或浅层土体，1为深层土体），σ为预测方程的标准差;hM，rrup是地震动强度预测函数的基本函数：[0031][0032]其中时，a=a2;M6.4时，a=a4;[0033][0034]f2M用于描述断层类型与滑移效益之间地震动的差异性，它与震级和周期有关；[0035]f3M，rrup=ftwMXfWrrup，函数f3M，rrup是用来考虑上盘效应，其中，[0036][0037]场地响应函数f4表示土层的非线性响应，是岩体峰值地面加速度PGA的预期值；[0038][0039]其中：上述各函数中的相关系数ai，...，a6，a9,...，313，〇4，〇5，11，5和6由桥址处所在地区周边的实测地震动数据及对应的地震参数如震级、断层距、断层类型等带入上述公式中进行拟合来确定，然后根据桥址实际的地震参数，预测桥址处发生某种地震等级时，地震动谱加速度Sa的值。[0040]地震设防水准共有6种，分别为50年超越概率为63%的小地震、100年超越概率为63%的中小地震、50年超越概率为10%的中震、100年超越概率为10%的中大地震、50年超越概率为2%的大震和100年超越概率为2%的特大地震，每种地震水准对应一个地震动强度和这种强度地震的年发生概率λΗ1，将所有水准下的值连成线便得到了地震危险曲线，如图2,根据公式可计算出各种发生概率下的地震水准对应的地震动强度，其中，Φ·代表标准正太分布涵教,P表示某种地震设防水准下的发生概率为预测出的InSa⑴）的平均值由前述公式可得，再根据前述相关公式，反算桥址处各种地震水准下的地震动强度僧，基于各地震水准下的数据点，j=l，2···，6,利用回归分析法便得到了桥址处的地震危险性曲线，地震危险性曲线的横坐标为地震动谱加速度，纵坐标为地震危险年平均超越概率，在地震危险性曲线上找出所有相邻两地震水准等级Sa1的平均谱加速度值，然后计算出相邻两个平均谱加速度值对应的地震年超越概率的差值，即为地震危险性贡献参数△λΗ4。[0041]如图3,运用拉丁超立方体抽样法对桥梁结构参数抽样，建立桥梁结构参数样本库，研究地震动输入的随机性建立地震动样本库，桥梁结构参数样本和地震动样本随机配对形成地震动-桥梁分析样本库，并开展时程分析，提取结构最大动力响应数据并确定构件损伤指标及其量化方法，利用时程分析结果，计算最优带宽参数h，求解带宽矩阵Η，然后根据公式，[0042][0043]计算特定强度地震作用下的桥梁结构的地震损伤破坏概率，其中，C表示桥梁地震响应的实际需求值；D表示桥梁构件抵抗地震作用的能力值；IM表示地震动强度值；为D与IM的联合概率密度分布函数，fIM·为地震强度的边缘分布函数，引入核密度估计法，概率密度分布函数的核密度估计为，其中，件表示高斯核函数，，〇为乂的标准差，h表示带宽参数，最优带宽参数取，〇为X的标准差，可求解地震强度边缘分布函数，即通过η组UM1J=1，...，η}地震动强度估计出，其中，已知η元变量X={χι,·χη}，各xii=l，2,·η变量之间相互独立，且各Xi变量中有m个样本点，则,„的核密度估计为将多元参数变量的标准正态核函数带入的核密度估计公式，则核密度估计按求解，H表示带宽矩阵，属于对称正定矩阵，IHI表示对H矩阵求行列式，（·）τ表示矩阵的转置，而后对桥梁开展η组地震作用下的时程分析，提取结构最大动力响应及地震动强度值便可计算，将和带入再通过积分求解，最终得到桥梁构件在特定地震作用下的损伤破坏概率为：[0044]:建立桥梁结构损伤破坏概率随地震动强度变化的关系曲线，如图4。[0045]最后，根据桥梁结构损伤破坏概率PupIMi和各地震设防水准下的地震危险性贡献参数△λΗ4,计算各水准地震（即地震谱加速度为Sai作用下，桥梁地震损伤概率PupMi和地震危险性贡献AλΗ,4^乘积值，然后将各水准地震的乘积值进行求和，得到桥梁结构地震风险概率值λυΡ，计算公式为，其中，IM1为第i组设防水准下结构特征周期对应的地震动强度指标，PlipM1表示强度为IM1的地震作用下桥梁结构发生某种破坏状态所对应的发生概率，通过桥梁结构地震风险概率值λυΡ来评估地震作用下桥梁结构的抗震性能，实现桥梁单个构件或整体结构在地震作用下地震风险评估，还可以对比研究不同桥梁减隔震措施下各构件的地震风险值，评估不同桥梁减隔震措施的可行性和适用性。[0046]通过上述方法，求解出桥梁各危险构件的地震风险值，根据结构风险评估相关规范，研究桥梁各构件的抗震性能，还可以对比研究不同桥梁减隔震措施下各构件的地震风险值，评估不同桥梁减隔震措施的可行性和适用性。

权利要求：1.一种基于地震风险评估的桥梁抗震分析方法，其特征在于，包括如下步骤：1根据桥址处的历史地震的记录数据和对应的地震参数，建立地震动强度预测方程组，计算出地震动谱加速度，结合地震设防水准，每种地震水准对应一个地震动强度和这种强度地震的年平均发生概率，将所有水准下的值连成线得到桥址处的地震危险性曲线；⑵根据所述地震危险性曲线，计算不同地震设防水准下的地震危险性贡献参数；3建立地震动-桥梁分析样本库并开展时程分析，提取结构最大动力响应数据，确定构件损伤指标及其量化方法，计算特定强度地震作用下的桥梁结构的地震损伤破坏概率，建立桥梁结构损伤破坏概率随地震动强度变化的关系曲线；⑷根据桥梁结构的地震损伤破坏概率和地震危险性贡献参数计算桥梁结构地震风险值。2.如权利要求1所述的一种基于地震风险评估的桥梁抗震分析方法，其特征在于，所述地震设防水准共有6种，分别为50年超越概率为63%的小地震、100年超越概率为63%的中小地震、50年超越概率为10%的中震、100年超越概率为10%的中大地震、50年超越概率为2%的大震和100年超越概率为2%的特大地震。3.如权利要求1所述的一种基于地震风险评估的桥梁抗震分析方法，其特征在于，所述地震危险性贡献参数为所述地震危险性曲线上相邻两个地震年超越概率的差值。4.如权利要求1所述的一种基于地震风险评估的桥梁抗震分析方法，其特征在于，在步骤⑶中，分别建立桥梁结构参数样本库和地震动样本库，桥梁结构参数样本和地震动样本随机配对形成地震动-桥梁分析样本库。5.如权利要求1所述的一种基于地震风险评估的桥梁抗震分析方法，其特征在于，在步骤3中，引入核密度估计法计算桥梁结构的地震损伤破坏概率。6.如权利要求1-5任一所述的一种基于地震风险评估的桥梁抗震分析方法，其特征在于，步骤⑴包含以下步骤：1-1首先搜集实际桥址处周边的实测地震动数据及对应的地震参数，包括震级、断层距、断层类型参数；1-2然后利用收集的实测数据进行拟合，确定各公式的系数；1-3根据待评估桥址附近实际存在的地震断层距离、断层类型参数，预测桥址处发生某种地震等级时，地震动谱加速度的值;其中工程场地特征周期为Tj条件下地震动谱加速度的计算公式为：其中，SaTj是工程场地的地震动在特定周期Tj的谱加速度预测值;M为地震矩等级；rrup为是距离断层面的直线距离，单位为km;F表示断层类型的系数，其中正断层为1，逆断层为〇.5，其他为0;HW是地震动上盘效应的系数，其中1为表示为上盘效应，0为其它效应;S是地质条件系数，其中〇为岩石土体或浅层土体，1为深层土体；是岩体峰值地面加速度PGA的预期值;〇为预测方程的标准差。7.如权利要求1-5任一所述的一种基于地震风险评估的桥梁抗震分析方法，其特征在于，桥梁构件在特定地震作用下的损伤破坏概率的计算方法为：其中，C表示桥梁地震响应的实际需求值;D表示桥梁构件抵抗地震作用的能力值;IM表示地震动强度值;Mi表示第i条地震动对应的强度值;im表示地震动强度的具体值;hM表示地震动的带宽参数;H表示带宽矩阵，|H|表示对H矩阵求行列式；（·）T表示矩阵的转置；8.如权利要求1-5任一所述的一种基于地震风险评估的桥梁抗震分析方法，其特征在于，步骤⑷包含以下步骤：4-1计算各地震设防水准下的桥梁结构的地震损伤破坏概率和地震危险性贡献参数的乘积；4-2将各地震设防水准下的乘积值按公式：，求和，得到桥梁结构综合地震风险值，其中，λυΡ为桥梁结构的地震风险值;为第i组设防水准下结构特征周期对应的地震动强度指标，PlipM1表示强度为IM1的地震作用下桥梁结构发生某种破坏状态所对应的发生概率，△λΗ1*各地震设防水准下地震危险性贡献参数。

百度查询： 中铁二院贵阳勘察设计研究院有限责任公司 一种基于地震风险评估的桥梁抗震分析方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD