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

地震动峰值加速度衰减关系对中国西部地区的适用性分析

周典 李小军

杨钢, 吴建超, 蔡永建, 乔岳强, 胡庆, 汤勇. 地震现场房屋裂缝宽度快速检测[J]. 震灾防御技术, 2019, 14(1): 220-230. doi: 10.11899/zzfy20190121
引用本文: 周典, 李小军. 地震动峰值加速度衰减关系对中国西部地区的适用性分析[J]. 震灾防御技术, 2019, 14(3): 513-523. doi: 10.11899/zzfy20190305
Yang Gang, Wu Jianchao, Cai Yongjian, Qiao Yueqiang, Hu Qing, Tang Yong. A Rapid Width Detection Method of House Crack in Earthquake Field[J]. Technology for Earthquake Disaster Prevention, 2019, 14(1): 220-230. doi: 10.11899/zzfy20190121
Citation: Zhou Dian, Li Xiaojun. Analysis and Application of Attenuation Relationships for Peak Accelerations of Ground Motion in Western China[J]. Technology for Earthquake Disaster Prevention, 2019, 14(3): 513-523. doi: 10.11899/zzfy20190305

地震动峰值加速度衰减关系对中国西部地区的适用性分析

doi: 10.11899/zzfy20190305
基金项目: 

国家重点研发计划 2018YFC1506400

国家自然科学基金 U1839202

国家自然科学基金 51421005

详细信息
    作者简介:

    周典, 女, 生于1994年。硕士研究生。主要从事地震动特性与衰减规律研究工作。E-mail:njzhoudian@sina.com

    通讯作者:

    李小军, 男, 生于1965年。研究员。主要从事地震波动数值模拟研究工作。E-mail:beerli@vip.sina.com

Analysis and Application of Attenuation Relationships for Peak Accelerations of Ground Motion in Western China

  • 摘要: 地震动参数衰减关系描述了地震动参数随震级和距离等因素的变化情况,选取合适的地震动衰减关系是地震危险性分析中确定工程场地地震动参数的关键环节。本文收集、整理和分析了中国西部地区2012年3月9日—2017年9月16日间发生的42次地震事件中获取的强震动记录资料,并选择了3个国际上基于不同地区强震动记录数据建立的地震动加速度衰减关系(ASB11、SCEMY97和LLCS11),开展了加速度衰减关系计算结果与中国西部强震动记录值的对比研究。比较了衰减关系峰值加速度的预测中值与实际记录值之间的差异,并分析了其残差随震级和距离的变化。得到了以下结论:①相比于SCEMY97和LLCS11,ASB11更适合于中国西部地区;②研究中国西部地区的衰减关系时,有必要考虑高频地震动的震级饱和现象。同时,对于如何修正已有衰减关系以使其更适合于中国西部地区给出了相应的建议。
  • 地震动衰减关系是表征地震动参数随震级、距离和场地等因素变化规律的函数。地震动衰减关系作为工程地震学中最常用的估计或预测地震动的方法,被广泛用于工程地震学领域,包括地震危险性分析、地震区划、地震小区划和设计地震动确定等多方面(袁一凡等,2012)。地震动衰减关系通常采用相对固定的函数形式,大致由震源影响函数、距离影响函数、场地影响函数和修正项4个部分构成。衰减关系的系数值一般利用大量强震动记录资料统计确定。每个地区的地震构造环境、地壳介质以及地表和近地表局部场地条件等因素具有很大的差异性,因此,地震动衰减关系具有很强的地区性。欧洲(Akkar等,2014)、美国(Sadigh等,1997)、中国台湾(Lin等,2011)等国家和地区都获得基于强震动记录的适用于相应地区的地震动衰减关系。

    何种地震动衰减关系更为适合1个特定的地区始终是研究的热点。Campbell(2016)利用NGA-West2中的数据,研究了3个针对欧洲和中东地区的地区性衰减关系和1个NGA衰减关系,分别比较了地震动衰减关系的预测中值和标准偏差,并进行了残差分析。认为Campbell-Bozorgnia衰减关系在更具有工程实用价值的大震和近场上优于另外3个衰减关系。Zafarani等(2017)为了找到更适合伊朗高原的小震和中震的地震动预测方程(GMPE),使用3种不同的统计学方法,对8个GMPE进行了比较,其中只有Zafarani于2015年提出的衰减关系在全部的频率范围内均与数据较吻合。

    虽然中国有很长的历史地震记载,但已获取强震动记录的地震事件不多,总体的强震动记录资料较少。因此,长久以来中国研究地震影响的分布往往采用烈度衰减关系。然而,地震烈度的空间分布虽然也能体现地震动强度的空间变化,反映其衰减特征,但它并非能直接反映地震动特性的相关物理量,甚至可能受到一些主观因素的干扰。经过多年来的发展,目前中国已经积累了一定量的强震动加速度记录,这或将使得我们有可能直接建立地震动加速度衰减关系(王继等,2008)。此外,也可以使用中国的强震动加速度记录,与目前国际上已有的衰减关系进行比较,找出其中较为适用的衰减关系,并进而研究对已有衰减关系进行修正处理的方法。司宏俊等(2014)比较了汶川地震和芦山地震中基岩和土层场地观测台站得到的加速度反应谱和基于日本强震动观测数据得出的基岩加速度反应谱衰减关系,发现两者有着相似的衰减特征,并通过拟合强震动记录加速度反应谱及其与衰减关系的残差,分析得到汶川地震和芦山地震的地震动记录台站的场地影响,与日本衰减关系中平均场地参数比较接近,VS30约为500m/s。Dangkua等(2018)考察了NGA-West2项目给出的地震动衰减关系对于中国地区的适用性,认为NGA-West2中的衰减关系对短周期地震动适用性较好,但在长周期范围结果偏高;此外,还利用对数最大似然法,为这些衰减关系分配了权重。

    在地震动衰减关系中,峰值加速度(PGA)是主要被考虑的参数。地震动峰值加速度是指在一次特定地震事件中取得的1条加速度记录的绝对最大值,常被用于描述地震动的强度以及对结构反应的影响程度。虽然峰值加速度并非能完全控制结构反应的地震动参数,特别是对于较长自振周期的结构,但因为其是强震动记录中最容易取得的参数,且结构抗震设计规范中多采用标准形式的地震动反应谱与峰值加速度结合使用,因此,人们在进行衰减关系的研究时仍然十分关心且首先关注峰值加速度(Douglas,2003)。

    在地震动衰减关系研究中,为了关注有显著工程意义的地震动记录,同时也为了减小独立变量的误差,可以选取MW>5地震的地震动记录进行研究(Douglas,2003)。因此,本文采用中国西部地区2012年3月9日—2017年9月16日间发生的42次MW>5地震事件中共916条有效水平方向强震动记录(峰值加速度绝对值超过10gal)。震中位置和发震时间来自地震简报,震源机制和震级参考GCMT(Dziewonski等,1981Ekström等,2012),详细情况如表 1所示。研究中涉及的地震事件和地震台站的位置分别见图 2图 3

    表 1  地震事件基本信息
    Table 1.  Basic information of earthquake events
    事件编号 地震时间 纬度/°N 经度/°E 矩震级MW
    1 2012年3月9日6:50 39.4 81.3 5.9
    2 2012年9月7日11:19 27.5 104 5.6
    3 2013年1月18日20:42 31 99.4 5.6
    4 2013年3月3日13:41 25.9 99.7 5.4
    5 2013年3月11日11:01 40.2 77.5 5.2
    6 2013年3月29日13:01 43.4 86.8 5.4
    7 2013年4月17日9:45 26.9 99.8 5.3
    8 2013年4月20日8:02 30.3 103 6.6
    9 2013年7月22日7:45 34.5 104.2 6.0
    10 2013年8月30日13:27 43.8 87.6 5.3
    11 2013年9月20日5:37 37.7 101.5 5.1
    12 2014年5月24日4:49 25 97.8 5.7
    13 2014年5月30日9:20 25 97.8 5.9
    14 2014年7月9日5:52 39.3 78.3 5.0
    15 2014年8月3日16:30 27.1 103.3 6.2
    16 2014年8月17日6:07 28.1 103.5 5.1
    17 2014年10月2日23:56 36.4 97.8 5.2
    18 2014年10月7日21:49 23.4 100.5 6.1
    19 2014年11月22日16:55 30.26 101.69 6.1
    20 2014年11月25日23:19 30.18 101.73 5.7
    21 2014年12月6日2:43 23.31 100.5 5.6
    22 2014年12月6日18:20 23.31 100.5 5.5
    23 2015年1月10日13:21 40.2 77.3 5.1
    24 2015年3月1日18:24 23.5 98.9 5.3
    25 2015年4月15日15:39 39.8 106.3 5.5
    26 2015年7月3日9:07 37.6 78.2 6.4
    27 2015年10月26日17:09 36.5 70.8 7.5
    28 2015年11月23日5:02 38 100.4 5.2
    29 2015年12月7日15:50 38.2 72.9 7.2
    30 2016年1月14日5:18 42.19 84.12 5.1
    31 2016年1月21日1:13 37.71 101.7 5.9
    32 2016年5月18日0:48 26.1 99.53 5.0
    33 2016年6月26日16:17 39.43 73.4 6.4
    34 2016年9月23日1:23 30.1 99.61 5.2
    35 2016年11月25日22:24 39.27 74.04 6.6
    36 2016年12月8日13:15 43.83 86.35 6.0
    37 2017年3月27日7:55 25.89 99.8 5.1
    38 2017年6月3日18:11 37.99 103.56 5
    39 2017年8月8日21:19 33.2 103.82 6.5
    40 2017年8月9日7:27 44.27 82.89 6.5
    41 2017年9月16日18:11 42.11 83.43 5.5
    42 2017年9月30日14:14 32.27 105 5.1
    下载: 导出CSV 
    | 显示表格
    图 2  地震事件分布
    Figure 2.  The spatial distribution of earthquake events
    图 3  强震动台站分布
    Figure 3.  The location distribution of strong ground motion stations

    VS30即场地地下30m深度内岩土层的平均剪切波速,是用于表征场地条件的重要物理参数。研究中涉及的强震动观测台站,其VS30值需要借助于已有的研究成果确定。因此,本文多数观测台站的VS30值采用了姜治军(2017)俞畑等(2015)的研究成果,而对于已有研究成果中没有涉及的台站,其VS30值借助100km范围内最近台站的VS30或同类场地的VS30均值进行近似处理来确定。此外,本次研究还利用观测场地的VS30值对观测台站的场地类型进行划分。传统上将VS30=750m/s的场地作为基岩场地和土层场地的分界波速值,但因为VS30>750m/s的台站十分缺乏,为此,本文以VS30=360m/s将场地划分两类,一类包括基岩和硬土场地,另一类包括中硬土和软土场地。这一场地分类处理与Sadigh等(1997)的做法一致。

    文中选择的强震动记录大部分来自基岩和硬土场地,仅有少部分数据来自中硬土和软土场地(图 1)。

    图 1  强震动记录的震级-震中距分布
    Figure 1.  The magnitude v.s. epic-center distance of strong ground motion records used in this study

    选择来自欧洲(Akkar等,2014)、美国(Sadigh等,1997)和中国台湾(Lin等,2011)的3个较有代表性的地区性地震动衰减关系。这些衰减关系的建立均未利用中国西部地区的强震动记录数据,在确定衰减关系系数时,是以当地的强震动记录数据为主,并加入其它地区发生的大地震的强震动记录进行补充。选择的衰减关系均采用浅层地壳地震的强震动记录进行推导,与本文采用的中国西部地区强震动记录的地震具有可比性。所选地震动衰减关系的具体情况见表 2

    表 2  所用地震动峰值加速度衰减关系及相关参数
    Table 2.  Peak acceleration attenuation relations and their parameters adopted in this study
    衰减关系编号 地区 震级MW 震中距/km 参数
    ASB14(Akkar等,2014) 泛欧和中东 4.0—7.6 <200 RJBRepiRhypMWVS30FNFR
    SCEMY97(Sadigh等,1997) 美国西部 3.8—7.4 <200 Rrup(Repi),MWFNFS
    LLCS11(Lin等,2011) 中国台湾 3.5—7.6 <240 Rrup(Repi),MW
    下载: 导出CSV 
    | 显示表格

    表 2中,RJB为Joyner-Boore距离,Rhyp为震源距,Repi为震中距,Rrup为台站到破裂平面的距离(震级不太大时可近似为Repi),MW为矩震级,VS30为30m处剪切波速,FNFRFS分别为正断层、逆冲断层和走滑断层的指标,一般为0;当且仅当震源机制为正断层时,FN=1;当且仅当震源机制为逆冲时,FR=1;当且仅当震源机制为走滑时,FS=1。

    针对本文所用的地震动峰值加速度衰减关系的震级与震中距的适用范围,进一步将所采用的地震动记录限制在震级MW 5.0—7.5、震中距小于200km范围内,即只利用了震中距小于200km范围内的地震动记录,而不考虑震中距更大的那些地震动记录。

    衰减关系LLCS11还考虑了断层上、下盘效应对地震动的影响,并分别为上、下盘给出了不同的系数。考虑到当地的地震构造背景,这种划分也许是有益的。但是,因为本文获得的地震动数据难以区分断层上、下盘,故本文将LLCS11中同种场地的上、下盘取平均使用。

    本文采用的3个地震动衰减关系的峰值加速度衰减计算公式如下(峰值加速度单位为g)。

    2.1.1   衰减关系ASB14
    $$ \text{ln}(Y)=\text{ln}[{{Y}_{\text{REF}}}({{M}_{w}}, R, SoF)]+\text{ln}[S({{V}_{S30}}\text{, PG}{{\text{A}}_{\text{REF}}})]+\varepsilon \sigma $$ (1)
    $$ {\rm{ln}}({Y_{{\rm{REF}}}}) = \\ \left\{ \begin{array}{l} {a_1} + {a_2}({M_{\rm{W}}} - {c_1}) + {a_3}{(8.5 - {M_{\rm{W}}})^2} + [{a_4} + {a_5}({M_{\rm{W}}} - {c_1})]\ln (\sqrt {{R^2} + a_6^2}) + {a_8}{F_N} + {a_9}{F_R}, {M_{\rm{W}}} \le {c_{\rm{1}}}\\ {a_1} + {a_7}({M_{\rm{W}}} - {c_1}) + {a_3}{(8.5 - {M_{\rm{W}}})^2} + [{a_4} + {a_5}({M_{\rm{W}}} - {c_1})]\ln (\sqrt {{R^2} + a_6^2}) + {a_8}{F_N} + {a_9}{F_R}, {M_{\rm{W}}}{c_{\rm{1}}} \end{array} \right. $$ (2)
    $$ \ln \left( S \right) = \left\{ \begin{array}{l} {b_1}\ln \left( {{V_{{\rm{S}}30}}/{V_{{\rm{REF}}}}} \right) + {b_2}\ln \left[ {\frac{{{\rm{PG}}{{\rm{A}}_{{\rm{REF}}}} + c{{\left( {{V_{{\rm{S}}30}}/{V_{{\rm{REF}}}}} \right)}^n}}}{{\left( {{\rm{PG}}{{\rm{A}}_{{\rm{REF}}}} + c} \right){{\left( {{V_{{\rm{S}}30}}/{V_{{\rm{REF}}}}} \right)}^n}}}} \right],{V_{{\rm{S30}}}} \le {V_{{\rm{REF}}}}\\ {b_1}{\rm{ln}}\left[ {\frac{{{\rm{min}}\left( {{V_{{\rm{S}}30}},{V_{{\rm{CON}}}}} \right)}}{{{V_{{\rm{REF}}}}}}} \right],{V_{{\rm{S}}30}} > {V_{{\rm{REF}}}} \end{array} \right. $$ (3)

    其中,PGAREFVS30REF为参考场地的峰值地震动和VS30,参考场地的VS30为750m/s。YREF表示由衰减关系计算的参考场地的地震动变量,当计算参考场地的峰值地震动时,YREF=PGAREF。各系数取值为a1=2.52977,a2=0.0029,a3=-0.05496,a4=-1.31001,a5=0.2529,a6=7.5,a7=-0.5096,a8=-0.1091,a9=0.937,b1=-0.41997,b2=-0.28846,c=2.5,c1=6.75,n=3.2。

    2.1.2   衰减关系SCEMY97

    对于基岩场地:

    $$ {\rm{ln}}(y) = {C_1} + {C_2}M + {C_3}{(8.5M)^{2.5}} + {C_4}\ln ({R_{rup}} + \exp ({C_5} + {C_6}M)) + {C_7}\ln ({R_{rup}} + 2) $$ (4)

    其中,MW≤6.5:C1=-0.524,C2=1.0,C3=0.000,C4=-2.100,C5=1.29649,C6=0.250,C7=0.0;MW>6.5:C1=-1.274,C2=1.1,C3=0.000,C4=-2.100,C5=-0.48451,C6=0.524,C7=0.0。

    对于土层场地:

    $$ {\rm{ln}}(y) = {C_1} + {C_2}M - {C_3}\ln ({R_{rup}} + {C_4}{e^{{C_5}M}}) + {C_6} + {C_7}{(8.5 - M)^{2.5}} $$ (5)

    其中,对于走滑型地震,C1=-2.17;对于正断层地震和逆冲型地震,C1=-1.92;当MW≤6.5,C4=2.1861,C5=0.32;当MW>6.5,C4=0.3825,C5=0.5882;C2=1.0,C3=1.70,C6=0.0,C7=0.0。

    2.1.3   衰减关系LLCS11

    对于基岩场地:

    $$ \left\{ \begin{array}{l} {\rm{ln}}({\rm{PGA}}) = - 3.279 + 1.035M - 1.651\ln (R + 0.152{e^{0.623M}}), 上盘\\ {\rm{ln}}({\rm{PGA}}) = - 3.232 + 1.047M - 1.662\ln (R + 0.192{e^{0.630M}}), 下盘 \end{array} \right. $$ (6)

    对于土层场地:

    $$ \left\{ \begin{array}{l} {\rm{ln}}({\rm{PGA}}) = - 3.248 + 0.943M - 1.471\ln (R + 0.100{e^{0.648M}}), 上盘\\ {\rm{ln}}({\rm{PGA}}) = - 3.218 + 0.935M - 1.464\ln (R + 0.125{e^{0.650M}}), 下盘 \end{array} \right. $$ (7)

    分别根据上述3个衰减关系公式计算峰值加速度,并与实际强震动记录中的峰值加速度进行比较,计算残差(R)。同时计算了残差在不同震中距和震级上的标准偏差(σ),得到了残差随距离与震级的分布图像(图 4图 5图 6)。本文残差为实际记录中的峰值加速度减去衰减关系计算值,标准偏差为方差的平方根。

    图 4  ASB14计算峰值加速度与实际强震动记录的残差随震级和距离的分布
    Figure 4.  The distribution of residuals between peak ground acceleration calculated using ASB14 and recorded against magnitude and distance
    图 5  SCEMY97计算峰值加速度与实际强震动记录的残差随距离和震级的分布
    Figure 5.  The distribution of residuals between peak ground acceleration calculated using SCEMY97 and recorded against magnitude and distance
    图 6  LLCS11计算峰值加速度与实际强震动记录的残差随距离和震级的分布
    Figure 6.  The distribution of residuals between peak ground acceleration calculated using LLCS11 and recorded against magnitude and distance

    图 4图 5图 6可以看出,大多数数据点的残差虽然不为0,但均位于[R-σR+σ]范围内,且残差的标准偏差与震级、距离关系不大,说明可能存在与距离和震级无关的误差。

    SCEMY97和LLCS11整体上存在正的残差,即根据这2个衰减关系计算得到的峰值加速度小于中国西部地区的强震动记录值。这种趋势随距离和震级基本不变,可能意味着2个衰减关系需要引入1个常数项进行修正。ASB14虽然也有整体为正的残差,但残差随距离增大而增大,说明该衰减关系需要引入1个随距离增大而增大的项进行修正。总体来说,ASB14的残差小于SCEMY97和LLCS11。

    此外,分别计算了ASB14、SCEMY97和LLCS11的残差与震级、距离的相关系数(表 3),发现SCEMY97和LLCS11的残差与距离成中等相关,两者之间可能有中等程度的线性关系。这意味着若要对这些衰减关系进行修改,可以尝试引入线性项。考虑到相关系数的特点,这种线性关系可能需要更多的观测记录才能确定。

    表 3  残差与震级、距离的相关系数
    Table 3.  Residual-magnitude and residual-distance correlation coefficients
    衰减关系编号 残差与震级的相关系数 残差与距离的相关系数
    ASB14 -0.121235 0.495025
    SCEMY97 0.138204 0.64689
    LLCS11 0.0942092 0.582699
    下载: 导出CSV 
    | 显示表格

    观察上述3个地震动衰减(ASB14、SCEMY97和LLCS11)的残差与震级的关系,可以发现在地震动峰值加速度衰减关系中,考虑震级的平方项有重要的意义。衰减关系中震级的平方项与高频地震动的震级饱和现象有关,而这一现象在震级较大时表现更为明显。因为实际工程中更应关注的是震级较大的地震,所以在建立地震动峰值加速度衰减关系时,应该尽量采用具有震级的平方项的形式。

    此外,上述3个地震动衰减关系的残差关于距离均存在一定程度的偏离。观察残差关于距离的分布,发现若对ASB14进行常数项的修正,对SCEMY97和LLCS11进行距离的线性项和常数项的修正,可能会使3个地震动衰减关系更适合于中国西部地区。

    对于中国西部强震动记录数据,ASB14的表现优于SCEMY97和LLCS11,如果直接利用国外的地震动衰减关系,ASB14更适合于中国西部地区。然而总体来看,无论是随着震级还是距离的变化,采用3个衰减关系预测中国西部地区的地震动将会明显低估地震动峰值加速度。因此,发展强震动观测以积累更多的强震动记录,并基于中国不同区域的强震动记录资料直接统计确定相应的地震动衰减关系,仍然是需要努力的目标。

  • 图  2  地震事件分布

    Figure  2.  The spatial distribution of earthquake events

    图  3  强震动台站分布

    Figure  3.  The location distribution of strong ground motion stations

    图  1  强震动记录的震级-震中距分布

    Figure  1.  The magnitude v.s. epic-center distance of strong ground motion records used in this study

    图  4  ASB14计算峰值加速度与实际强震动记录的残差随震级和距离的分布

    Figure  4.  The distribution of residuals between peak ground acceleration calculated using ASB14 and recorded against magnitude and distance

    图  5  SCEMY97计算峰值加速度与实际强震动记录的残差随距离和震级的分布

    Figure  5.  The distribution of residuals between peak ground acceleration calculated using SCEMY97 and recorded against magnitude and distance

    图  6  LLCS11计算峰值加速度与实际强震动记录的残差随距离和震级的分布

    Figure  6.  The distribution of residuals between peak ground acceleration calculated using LLCS11 and recorded against magnitude and distance

    表  1  地震事件基本信息

    Table  1.   Basic information of earthquake events

    事件编号 地震时间 纬度/°N 经度/°E 矩震级MW
    1 2012年3月9日6:50 39.4 81.3 5.9
    2 2012年9月7日11:19 27.5 104 5.6
    3 2013年1月18日20:42 31 99.4 5.6
    4 2013年3月3日13:41 25.9 99.7 5.4
    5 2013年3月11日11:01 40.2 77.5 5.2
    6 2013年3月29日13:01 43.4 86.8 5.4
    7 2013年4月17日9:45 26.9 99.8 5.3
    8 2013年4月20日8:02 30.3 103 6.6
    9 2013年7月22日7:45 34.5 104.2 6.0
    10 2013年8月30日13:27 43.8 87.6 5.3
    11 2013年9月20日5:37 37.7 101.5 5.1
    12 2014年5月24日4:49 25 97.8 5.7
    13 2014年5月30日9:20 25 97.8 5.9
    14 2014年7月9日5:52 39.3 78.3 5.0
    15 2014年8月3日16:30 27.1 103.3 6.2
    16 2014年8月17日6:07 28.1 103.5 5.1
    17 2014年10月2日23:56 36.4 97.8 5.2
    18 2014年10月7日21:49 23.4 100.5 6.1
    19 2014年11月22日16:55 30.26 101.69 6.1
    20 2014年11月25日23:19 30.18 101.73 5.7
    21 2014年12月6日2:43 23.31 100.5 5.6
    22 2014年12月6日18:20 23.31 100.5 5.5
    23 2015年1月10日13:21 40.2 77.3 5.1
    24 2015年3月1日18:24 23.5 98.9 5.3
    25 2015年4月15日15:39 39.8 106.3 5.5
    26 2015年7月3日9:07 37.6 78.2 6.4
    27 2015年10月26日17:09 36.5 70.8 7.5
    28 2015年11月23日5:02 38 100.4 5.2
    29 2015年12月7日15:50 38.2 72.9 7.2
    30 2016年1月14日5:18 42.19 84.12 5.1
    31 2016年1月21日1:13 37.71 101.7 5.9
    32 2016年5月18日0:48 26.1 99.53 5.0
    33 2016年6月26日16:17 39.43 73.4 6.4
    34 2016年9月23日1:23 30.1 99.61 5.2
    35 2016年11月25日22:24 39.27 74.04 6.6
    36 2016年12月8日13:15 43.83 86.35 6.0
    37 2017年3月27日7:55 25.89 99.8 5.1
    38 2017年6月3日18:11 37.99 103.56 5
    39 2017年8月8日21:19 33.2 103.82 6.5
    40 2017年8月9日7:27 44.27 82.89 6.5
    41 2017年9月16日18:11 42.11 83.43 5.5
    42 2017年9月30日14:14 32.27 105 5.1
    下载: 导出CSV

    表  2  所用地震动峰值加速度衰减关系及相关参数

    Table  2.   Peak acceleration attenuation relations and their parameters adopted in this study

    衰减关系编号 地区 震级MW 震中距/km 参数
    ASB14(Akkar等,2014) 泛欧和中东 4.0—7.6 <200 RJBRepiRhypMWVS30FNFR
    SCEMY97(Sadigh等,1997) 美国西部 3.8—7.4 <200 Rrup(Repi),MWFNFS
    LLCS11(Lin等,2011) 中国台湾 3.5—7.6 <240 Rrup(Repi),MW
    下载: 导出CSV

    表  3  残差与震级、距离的相关系数

    Table  3.   Residual-magnitude and residual-distance correlation coefficients

    衰减关系编号 残差与震级的相关系数 残差与距离的相关系数
    ASB14 -0.121235 0.495025
    SCEMY97 0.138204 0.64689
    LLCS11 0.0942092 0.582699
    下载: 导出CSV
  • 姜治军, 2017.强震动台站的场地VS30估计及场地放大系数模型研究.哈尔滨: 中国地震局工程力学研究所.
    司宏俊, 纐缬一起, 三宅弘惠等, 2014.基于基岩衰减公式的汶川·芦山地震强震动评价.地震工程与工程振动, 34(4):93-100. http://d.old.wanfangdata.com.cn/Conference/8421966
    王继, 俞言祥, 龚飞等, 2008.中国西部地震动衰减关系的适用性分析——以部分新疆地震为例.中国地震, 24(2):126-133. doi: 10.3969/j.issn.1001-4683.2008.02.004
    喻畑, 李小军, 2015.四川、甘肃地区VS30经验估计研究.地震工程学报, 37(2):525-533. doi: 10.3969/j.issn.1000-0844.2015.02.0525
    袁一凡, 田启文, 2012.工程地震学.北京:地震出版社.
    Akkar S., Sandıkkaya M. A., Bommer J. J., 2014. Empirical ground-motion models for point- and extended-source crustal earthquake scenarios in Europe and the Middle East. Bulletin of Earthquake Engineering, 12(1):359-387. doi: 10.1007/s10518-013-9461-4
    Campbell K. W., 2016. Comprehensive comparison among the Campbell-Bozorgnia NGA-West2 GMPE and three GMPEs from Europe and the Middle East. Bulletin of the Seismological Society of America, 106(5):2081-2103. doi: 10.1785/0120160024
    Dangkua D. T., Rong Y. F., Harold M., 2018. Evaluation of NGA-West2 and Chinese Ground-motion prediction equations for developing seismic hazard maps of Mainland China. Bulletin of the Seismological Society of America, 108(5A):2422-2443. doi: 10.1785/0120170186
    Douglas J., 2003. Earthquake ground motion estimation using strong-motion records:a review of equations for the estimation of peak ground acceleration and response spectral ordinates. Earth-Science Reviews, 61(1-2):43-104. doi: 10.1016/S0012-8252(02)00112-5
    Dziewonski A. M., Chou T. A., Woodhouse J. H., 1981. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity. Journal of Geophysical Research:Solid Earth, 86(B4):2825-2852. doi: 10.1029/JB086iB04p02825
    Ekström G., Nettles M., Dziewoński A. M., 2012. The global CMT project 2004-2010:Centroid-moment tensors for 13, 017 earthquakes. Physics of the Earth and Planetary Interiors, 200-201:1-9. doi: 10.1016/j.pepi.2012.04.002
    Lin P. S., Lee C. T., Cheng C. T., et al., 2011. Response spectral attenuation relations for shallow crustal earthquakes in Taiwan. Engineering Geology, 121(3-4):150-164. doi: 10.1016/j.enggeo.2011.04.019
    Sadigh K., Chang C. Y., Egan J. A., et al., 1997. Attenuation relationships for shallow crustal earthquakes based on California strong motion data. Seismological Research Letters, 68(1):180-189. doi: 10.1785/gssrl.68.1.180
    Zafarani H., Farhadi A., 2017. Testing ground-motion prediction equations against small-to-moderate magnitude data in Iran. Bulletin of the Seismological Society of America, 107(2):912-933. doi: 10.1785/0120160046
  • 期刊类型引用(0)

    其他类型引用(6)

  • 加载中
图(6) / 表(3)
计量
  • 文章访问数:  162
  • HTML全文浏览量:  41
  • PDF下载量:  31
  • 被引次数: 6
出版历程
  • 收稿日期:  2019-04-11
  • 刊出日期:  2019-09-01

目录

/

返回文章
返回