Energy magnitude determination of moderate-strong earthquakes
-
摘要: 本文根据地震波衰减特性,开展了利用宽频带地震波形数据测定地震波能量ES和能量震级Me的方法研究。利用震中距处于20°—98°范围内的宽频带远震P波波形数据,测定了4次国外和4次国内中强震的能量震级Me,并对其面波震级MS、矩震级MW及能量震级Me进行了分析对比。结果表明:面波震级MS表示的是地震在某一固定频率所辐射的地震波能量大小;矩震级MW与地震所产生的断层长度、断层宽度、震源破裂的平均位错量等静态构造效应密切相关;而能量震级Me反映的是震源动态特征,与地震震源的动力学特性密切相关。由于地震是以地震波形式辐射,能量主要集中在震源谱的拐角频率附近,因此能量震级Me更适合描述地震的破坏性。由此可见,联合测定面波震级MS,矩震级MW和能量震级Me对于地震定量研究以及地震灾害与风险评估具有重要作用。
-
关键词:
- 地震波能量 /
- 能量震级 /
- 宽频带远震P波波形数据 /
- 震源动态特征 /
- 震源谱
Abstract: This paper, based on the seismic wave attenuation characteristics, explores a method for utilizing broadband seismic waveform data to measure seismic energy ES and energy magnitude Me. The energy magnitudes Me of four overseas and four domestic moderate-strong earthquakes are measured by using broadband teleseismic P waveform data with epicentral distance in the range of 20°−98°. The comparison of surface wave magnitude MS, moment magnitude MW and energy magnitude Me of several earthquakes shows that: surface wave magnitude MS represents the magnitude of the energy of seismic wave radiating at a fixed frequency; the moment magnitude MW is closely related to the static structural effects such as the length of fault, the width of fault and the average dislocation amount of the source rupture produced by the earthquake; the energy magnitude Me reflects source dynamic characteristics and is closely related to the dynamic characteristics of the seismic source. The energy radiated by the earthquake in the form of seismic waves is mainly concentrated near the corner frequency of the source spectrum, so the energy magnitude Me is more suitable to describe the destructiveness of the earthquake. Therefore, joint determination of surface wave magnitude MS, moment magnitude MW and energy magnitude Me are of great significance for quantifying and evaluating the disasters caused by earthquakes or tsunamis. -
引言
场地条件对地震动的影响起决定作用,场地条件主要考虑以下几个方面(薄景山等,2021):近地表的工程地质条件(迟明杰等,2021)、水文地质条件(李伟华,赵成刚,2015)、局部地形效应(王伟等,2016;李平等,2018)及断层场地效应(周正华等,2003)。近地表的工程地质条件可从岩土类型、覆盖层厚度及土层结构等三个方面来考察。一般来说,在相同的地震动作用下,基岩场地上的地震动较小,硬土场地次之,软土场地地震动最大。已有分析表明:土层结构对地震动也具有显著的影响(王海云,2011;李平等,2012),主要表现为:硬夹层(张海等,2016;周正华等,2019,2020; Su et al,2020)埋深增加时,地表峰值加速度(peak ground acceleration,缩写为PGA)、反应谱谱值增大;厚度增加时,地表峰值加速度先逐渐减小、继而逐渐增大,反应谱谱值增大。软夹层(王亚红等,2019)埋深或厚度增加时,地表峰值加速度、反应谱峰值减小。
近年来,含软弱土层的场地对场地地震反应的影响已引起地震工程领域人员的重视,一些学者从不同的角度对其展开了深入的研究。许建聪等(2005)对福州软土场地的地震破坏机理进行了研究,其结果表明,软土在一定程度上增大了场地卓越周期,且卓越周期与地基土的刚度密切相关。Yao等(2022)认为局部软夹层的存在可以显著放大或衰减地震动,从而可能影响空间变化与地震动之间的滞后空间一致性。曹志翔(2006)分析了软夹层埋深对地表地震动的影响,结果表明,随软弱夹层位置的加深,场地的放大作用减小,场地卓越频率减小。王亚红等(2019)通过等效线性方法的场地地震反应分析,研究软夹层埋深对地表地震动参数的影响,结果显示:随着软夹层埋深由浅到深的变化,加速度峰值、速度峰值和反应谱峰值也会随着减小,特征周期和反应谱最大值所对应的周期也会随着软夹层的埋深增加而呈现递增的变化趋势。田守岐(2013)提出了相同Ⅲ类场地条件下,含弱软土层场地的地表地震动参数值与不含软弱土层的场地差异很大,软弱土层在土层结构中的埋深不同,地震动参数值的差异也很大。闫孔明等(2017)重点通过加速度放大效应和地震波波动机制探讨软弱夹层和斜坡效应对斜坡场地动力响应的影响,并结合传统傅里叶谱和希尔伯特(Hilbert)边际谱方法,在频域内展示了场地的频谱变化特性。李美娟和夏雄(2017)针对软土夹层在埋深不变的情况下,计算软夹层厚度分别为3 m,5 m和9 m的三个剖面来分析软土夹层的埋深和厚度对地表地震动参数的影响。王竞等(2022)采用理想场地作为基础研究模型,改变软土层所在位置得到相关计算模型,通过土层反应分析来探究不同位置软夹层对场地地表峰值加速度、放大系数、等效剪切波速计算值的影响,得到了软弱夹层对场地的峰值加速度影响为先放大后缩小、四类场地的传播能力有一定范围、等效剪切波速不能反映土层结构等相关结论。
综上所述,软弱土层对场地地震反应具有显著影响,该影响尤其表现为特征周期显著增大。经规准后的反应谱特征周期与我国GB 50011−2010 《 建筑抗震设计规范 》 (2016版)(中国建筑科学研究院,2016)相比,含软弱土层的场地地震动加速度反应谱特征周期远大于规范规定的取值。目前虽有针对软弱土层对场地地震反应影响的分析,但鲜有关于含软弱土层场地反应谱特征周期的调整方法研究。鉴于此,本文将在前人现有研究基础上,建立含淤泥土层的软弱场地模型,针对软弱土层对场地地震反应的影响进行分析,并提出场地地震动加速度反应谱特征周期的调整方法,以期为含软弱土层场地地震动反应谱特征周期的确定提供理论依据。
1. 地震动输入
含淤泥层的软弱场地地震反应分析中,输入地震动通过人工合成(Thráinsson,Kiremidjian,2002;徐国栋等,2010;刘帅等,2018)得到,峰值加速度分别为50 cm/s2,100 cm/s2和200 cm/s2,时间间隔为0.02 s,离散点数为2 048。对加速度时程幅值缩小50% (25 cm/s2,50 cm/s2,100 cm/s2)作为一维土层地震反应分析模型的计算基底输入地震动。输入地震动加速度时程及加速度反应谱如图1所示,其反应谱特征周期分别为0.30 s,0.35 s,0.40 s。
2. 场地地震反应分析
2.1 场地地震反应计算模型及动力参数的确定
为了分析含软弱土层场地地震反应,依据特定场地工程地震条件勘察资料,分别建立了含一层淤泥不同埋深的六个分析模型和含两层淤泥不同埋深的六个分析模型,以分析软弱土层厚度及埋深对场地地震反应的影响特征。
场地土层地震动反应分析,需要土层剖面的土层分层厚度及土层土体性状描述资料,同时也需要土层中土体的力学特性资料,包括土体的剪切波速值、土体的密度值及土体动力非线性特性参数值。依据实际现场勘查资料及实验数据,并考虑不同淤泥层厚度及不同埋深条件,建立了含一层与两层淤泥且埋深不同的六个分析模型。模型1的剖面与力学特性参数列于表1,模型2—模型6是在模型1的基础上通过改变淤泥层上覆粉质黏土的层数得到的,逐渐将模型1的淤泥层下伏粉质黏土逐层上移至淤泥层以上,如模型2上覆有土类号为3 (周正华等,2019)的一层粉质黏土,模型3上覆有土类号为3,4的两层粉质黏土,模型4上覆有土类号为3,4,5的三层粉质黏土,模型5上覆有土类号为3,4,5,6的四层粉质黏土,模型6上覆有土类号为3,4,5,6和7的五层粉质黏土。模型7的剖面与力学特性参数列于表1,模型8—12是在模型7的基础上亦通过改变淤泥层上覆粉质黏土的层数得到的,即模型8上覆有土类号为3的1层粉质黏土,模型9上覆有土类号为3,4的两层粉质黏土,模型10上覆有土类号为3,4和5的三层粉质黏土,模型11上覆有土类号为3,4,5,6的四层粉质黏土,模型12上覆有土类号为3,4,5,6,7的五层粉质黏土。
表 1 分析模型1和模型7的剖面和力学特性参数Table 1. Profile and mechanical characteristic parameters of Analysis model 1模型 土层
序号岩土名称 土类号 层底深度
/m层厚
/m剪切波速
/(m·s−1)密度
/(kg·m−3)模型1 1 淤泥 1 5.0 5.0 112 1580 2 粉质黏土 3 9.5 4.5 160 1860 3 粉质黏土 4 13.0 3.5 165 1870 4 粉质黏土 5 17.0 4.0 199 1880 5 粉质黏土 6 21.0 4.0 212 1960 6 粉质黏土 7 24.0 3.0 242 1980 7 圆砾 8 27.0 3.0 258 2200 8 全风化安山岩 8 30.0 3.0 393 2250 9 计算基底 9 516 2650 模型7 1 淤泥 1 5.0 5.0 112 1580 2 淤泥 2 10.0 5.0 112 1660 3 粉质黏土 3 14.5 4.5 160 1860 4 粉质黏土 4 18.0 3.5 165 1870 5 粉质黏土 5 22.0 4.0 199 1880 6 粉质黏土 6 26.0 4.0 212 1960 7 粉质黏土 7 29.0 3.0 242 1980 8 圆砾 8 32.0 3.0 258 2200 9 全风化安山岩 8 35.0 3.0 393 2250 10 计算基底 9 516 2650 不同计算场地模型各土层土体动力非线性参数及土体密度值源于实际工程场地地震安全性评价项目的试验结果,其各土层土体动力非线性参数列于表2。
表 2 各土层不同剪应变水平下的动力剪切非线性参数Table 2. Nonlinear parameters of dynamic shear of all soils under different shear strain levels土类号 土层名称 参数 剪应变/(10−4) 0.05 0.1 0.5 1 5 10 50 100 1 淤泥 模量比G/Gmax 0.990 2 0.9808 60.910 5 0.835 8 0.504 5 0.337 4 0.092 3 0.048 3 阻尼比ζ 0.017 3 0.024 4 0.052 5 0.071 1 0.123 6 0.142 9 0.167 2 0.171 2 2 淤泥 模量比G/Gmax 0.991 3 0.982 7 0.918 9 0.850 0 0.531 3 0.361 7 0.101 8 0.053 6 阻尼比ζ 0.008 8 0.013 5 0.035 6 0.052 5 0.107 3 0.130 3 0.161 5 0.166 9 3 粉质黏土 模量比G/Gmax 0.991 8 0.983 8 0.924 1 0.858 8 0.548 9 0.378 3 0.108 5 0.057 3 阻尼比ζ 0.013 8 0.019 9 0.045 9 0.064 1 0.120 1 0.142 8 0.173 5 0.178 8 4 粉质黏土 模量比G/Gmax 0.992 5 0.985 1 0.929 6 0.868 4 0.568 9 0.397 5 0.116 6 0.061 9 阻尼比ζ 0.012 3 0.017 6 0.040 2 0.056 1 0.105 3 0.125 8 0.154 2 0.159 2 5 粉质黏土 模量比G/Gmax 0.993 9 0.987 8 0.941 9 0.890 3 0.618 7 0.447 9 0.139 6 0.075 0 阻尼比ζ 0.015 7 0.021 8 0.046 1 0.062 6 0.113 6 0.135 6 0.167 7 0.173 6 6 粉质黏土 模量比G/Gmax 0.994 3 0.988 7 0.946 0 0.897 5 0.636 5 0.466 8 0.149 0 0.080 5 阻尼比ζ 0.018 1 0.024 9 0.051 2 0.068 8 0.123 4 0.147 3 0.182 7 0.189 4 7 粉质黏土 模量比G/Gmax 0.995 0 0.990 1 0.952 4 0.909 2 0.666 9 0.500 3 0.166 8 0.091 0 阻尼比ζ 0.010 6 0.015 2 0.034 2 0.047 8 0.093 6 0.115 4 0.150 4 0.157 0 8 圆砾及卵石 模量比G/Gmax 0.990 0.970 0.900 0.850 0.700 0.550 0.320 0.200 阻尼比ζ 0.004 0.006 0.019 0.030 0.075 0.090 0.110 0.120 9 基岩 模量比G/Gmax 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 阻尼比ζ 0.004 0.008 0.010 0.015 0.021 0.030 0.036 0.046 2.2 计算结果分析
对各分析模型采用等效线性化方法(高武平等,2012;荣棉水等,2013)进行场地地震反应计算时,分别以三个峰值加速度水平(25 cm/s2,50 cm/s2,100 cm/s2)对应的输入地震动加速度时程作为一维土层反应分析模型的基底入射地震动,通过计算各分析模型进行水平向地震反应,得到各模型地表水平向地震反应的地震动加速度时程和反应谱值。
各分析模型不同地震动作用下地表水平向地震反应峰值加速度及地表水平向地震反应动力放大系数列于表3。由表3可以看出,输入相同峰值加速度,软弱土层厚度越厚,地表峰值加速度越小;软弱土层埋深越深,地表峰值加速度亦越小。由表3还可以看出,输入相同峰值加速度,软弱土层厚度越厚,地表峰值加速度动力放大系数越小;软弱土层埋深越深,地表峰值加速度动力放大系数亦越小,且随埋深增加,动力放大系数衰减变慢;在同一软弱土层厚度和相同埋深条件下,随输入峰值加速度的增加,地表峰值加速度动力放大系数逐渐减小,说明场地土具有显著的非线性。
表 3 各分析模型地表峰值加速度和地震反应动力放大系数Table 3. The peak acceleration and dynamic amplification coefficient for surface seismic response of each analysis model模型 输入不同峰值加速度的地表水平向峰值加速度/(cm·s−2) 模型 输入不同峰值加速度的地表水平向地震反应动力放大系数 PGA=25 cm/s2 PGA=50 cm/s2 PGA=100 cm/s2 PGA=25 cm/s2 PGA=50 cm/s2 PGA=100 cm/s2 模型1 53.1 97.1 193.2 模型1 2.124 1.942 1.932 模型2 43.8 86.5 152.3 模型2 1.752 1.73 1.523 模型3 41.2 72.4 118.7 模型3 1.648 1.448 1.187 模型4 38.7 58.2 109.9 模型4 1.548 1.164 1.099 模型5 36.4 52.6 100.6 模型5 1.456 1.052 1.006 模型6 30.6 45 87.4 模型6 1.224 0.900 0.874 模型7 48.9 96 174.2 模型7 1.956 1.92 1.742 模型8 38.7 62.6 110.1 模型8 1.548 1.252 1.101 模型9 33.7 55.2 102.8 模型9 1.348 1.104 1.028 模型10 31.1 53.9 83.9 模型10 1.244 1.078 0.839 模型11 30.4 48.2 69.4 模型11 1.216 0.964 0.694 模型12 29.6 46.4 64.6 模型12 1.184 0.928 0.646 依据表3给出了地表峰值加速度随软弱土层不同埋深下的变化,如图2所示。同样,依据表3给出了地表峰值加速度动力放大系数随软弱土层不同埋深下的变化,如图3所示。
由图2可以看出:输入地震动水平相同情况下,随软弱土层埋深加大,地表峰值加速度越小,且近地表衰减较快,埋深增加,衰减变慢;输入地震动水平不同时,输入峰值加速度越大,地表峰值加速度衰减越快;不同软弱土层厚度对地表峰值加速度有影响,土层越厚,地表峰值加速度越小,浅地表衰减越快,且随输入地震动峰值加速度增加,差异越明显。
由图3可以看出:输入地震动相同时,随软弱土层埋深加大,地表峰值加速度动力放大系数越小,且近地表衰减较快,随埋深增加,衰减变慢;输入地震动水平不同时,输入峰值加速度越大,动力放大系数衰减越快;不同软弱土层厚度对地表峰值加速度有影响,厚度越厚,动力放大系数越小,浅地表衰减越快,且随输入地震动峰值加速度的增加,差异越明显。
各模型地震反应分析中亦得到了阻尼比5%的场地相关反应谱,如图4所示。
由图4a可以看出:对于含一层淤泥的分析模型,相同输入地震动峰值加速度水平下,在场地相关反应谱起伏变化频带内,总体表现为随软弱土层埋深越深,加速度反应谱值越低,且场地相关反应谱衰减段起始频率越来越低,反应谱优势频带越来越宽;输入地震动峰值加速度不同时,各模型反应谱变化趋势大致一致,但输入峰值加速度越大,反应谱值越大;周期在1 s以上的加速度反应谱,各模型结果基本接近,表明软弱土层埋深对长周期地震动影响较小。
由图4b可以看出:对于含两层淤泥层的计算模型,输入相同地震动峰值加速度,在场地相关反应谱起伏变化频带内,总体表现为随软弱土层埋深加大,加速度反应谱值越低,且场地相关反应谱衰减段起始频率越来越低,反应谱优势频带越来越宽;不同输入地震动峰值加速度水平下,各模型反应谱变化趋势大致一致,但输入峰值加速度越大,反应谱值越大。周期在1 s以上的加速度反应谱,各模型结果基本接近,说明软弱土层埋深对长周期地震动影响较小。
依据工程场地自由地表水平向加速度反应谱,计算得到了阻尼比5%的场地相关归一化反应谱,如图5所示。
由图5a可以看出:对于含一层淤泥层的计算模型,相同输入地震动水平下,场地相关归一化反应谱均随软弱层埋深变化而变化,归一化反应谱起伏变化频带内,总体表现为随软弱土层埋深加大,归一化反应谱值越低,且场地相关归一化反应谱衰减段起始频率越来越低,归一化反应谱优势频带越来越宽;不同输入地震动下,各模型归一化反应谱变化趋势大致一致,但输入峰值加速度越大,归一化反应谱越小;周期在0.1 s及以下时,各模型的归一化反应谱基本接近,说明软弱土层埋深对短周期归一化反应谱影响较小,但在0.1 s以上周期段内,软弱土层埋深对归一化反应谱影响较大,总体表现为,软弱土层埋深越深,归一化反应谱值越大,彼此差异越明显。
对于含两层淤泥层的计算模型(图5b):输入地震动相同时,场地相关归一化反应谱均随软弱层埋深变化而变化,归一化反应谱起伏变化频带内,总体表现为随软弱土层埋深越深,归一化反应谱值越低,且场地相关归一化反应谱衰减段起始频率越来越低,归一化反应谱优势频带越来越宽;不同输入地震动水平下,各模型归一化反应谱变化趋势大致一致,但输入峰值加速度越大,归一化反应谱越小;在0.1 s及以下周期段内,各模型归一化反应谱基本接近,说明软弱土层埋深对短周期的归一化反应谱影响较小,但在0.1 s以上周期段内,软弱土层埋深对归一化反应谱影响较大,总体表现为,软弱土层埋深越深,归一化反应谱值越大,彼此差异越明显。
3. 反应谱特征周期调整方法
根据 《 建筑抗震设计规范 》 (中国建筑科学研究院,2016)的场地类别确定方法划分各分析模型场地类别,结果如表4所示。
表 4 分析模型的场地参数及场地类别Table 4. Site parameters and categories of analysis models分析
模型覆盖层
厚度/m等效剪切
波速/m·s−1场地
类别分析
模型覆盖层
厚度/m等效剪切
波速/m·s−1场地
类别模型1 30 156.0 Ⅱ 模型7 35 134.5 Ⅲ 模型2 30 156.0 Ⅱ 模型8 35 134.5 Ⅲ 模型3 30 156.0 Ⅱ 模型9 35 134.5 Ⅲ 模型4 30 156.0 Ⅱ 模型10 35 142.0 Ⅲ 模型5 30 161.3 Ⅱ 模型11 35 161.3 Ⅱ 模型6 30 182.4 Ⅱ 模型12 35 182.4 Ⅱ 由表4可以看出本文分析模型场地类别为Ⅱ 和 Ⅲ类,结合 《 建筑抗震设计规范 》 表5.1.4-2,若按第三设计地震分组考虑,其场地特征周期为0.45 s (对应于Ⅱ类场地)或0.65 s (对应于Ⅲ类场地)。
根据 《 建筑抗震设计规范 》 规定的规整反应谱格式(中国建筑科学研究院,2016)得到场地规准反应谱,并得到各模型场地地震动加速度反应谱特征周期(图6),如表5所示。
表 5 各分析模型场地规准反应谱特征周期TgTable 5. The characteristic periodic values of site standard response spectrum of each analysis model模型 输入不同峰值加速度各模型的反应谱特征周期 模型 输入不同峰值加速度各模型的反应谱特征周期 PGA=25 cm/s2 PGA=50 cm/s2 PGA=100 cm/s2 PGA=25 cm/s PGA=50 cm/s2 PGA=100 cm/s2 模型1 0.5 0.55 0.6 模型7 0.7 0.75 0.95 模型2 0.6 0.65 0.8 模型8 0.9 0.95 1.2 模型3 0.7 0.75 0.85 模型9 1 1.05 1.35 模型4 0.75 0.8 1.05 模型10 1.1 1.2 1.45 模型5 0.8 0.9 1.1 模型11 1.15 1.3 1.5 模型6 0.85 0.95 1.2 模型12 1.25 1.35 1.65 图 9 不同输入地震动下补充模型1、2与分析模型2场地相关加速度反应谱(a) 输入PGA=25 cm/s2;(b) 输入PGA=50 cm/s2;(c) 输入PGA=100 cm/s2Figure 9. The site-related acceleration response spectra of supplementary model 1,2 and analysis model 2 under different input ground motions(a) Input PGA=25 cm/s2;(b) Input PGA=50 cm/s2;(c) Input PGA=100 cm/s2由图6和表5可知:相同输入峰值加速度水平下,软弱土层厚度越厚,场地反应谱特征周期越大;软弱土层埋深越深,反应谱特征周期亦越大,且随埋深增加,反应谱特征周期增幅逐渐变慢,浅地表增幅较快;同一软弱土层厚度和相同埋深条件下,随输入峰值加速度的增加,反应谱特征周期逐渐增大。
为了进一步分析输入峰值加速度水平对含软弱土层场地地震反应的影响,在进行模型1—6的一维场地地震反应分析时,还计算了对应峰值加速度为25 cm/s2地震动时程调幅至50 cm/s2,100 cm/s2,150 cm/s2,200 cm/s2,300 cm/s2输入下的地震反应,其场地相关加速度反应谱见图7。在此基础上,按照前述方法对计算得到的场地相关加速度反应谱规准化,得到了规准反应谱和相关参数,规准谱见图7,规准谱特征周期见表6。
表 6 输入不同地震动水平下模型1—6的反应谱特征周期TgTable 6. Characteristic periods Tg of model 1−6 under different input ground motion分析模型 不同地震动输入水平下的反应谱特征周期 PGA=25 cm/s2 PGA=50 cm/s2 PGA=100 cm/s2 PGA=150 cm/s2 PGA=200 cm/s2 PGA=300 cm/s2 模型1 0.50 0.55 0.60 0.65 0.70 0.75 模型2 0.60 0.65 0.70 0.75 0.80 0.90 模型3 0.70 0.75 0.80 0.85 0.90 0.95 模型4 0.75 0.80 0.85 0.90 0.95 1.00 模型5 0.80 0.85 0.90 0.95 1.00 1.05 模型6 0.85 0.90 0.95 1.00 1.05 1.20 注:输入峰值加速度25 cm/s2的各模型特征周期来自于表5 由图7和表6可知,在不同输入地震动水平下,随输入地震动水平不断增加,各分析模型反应谱特征周期不断增大,例如:模型1反应谱特征周期从0.50 s逐渐增大到0.75 s,模型2从0.6 s增大到0.9 s,模型3从0.7 s增大到0.95 s,模型4从0.75 s增大到1.0 s,模型5从0.8 s增大到1.05 s,模型6从0.85 s增大到1.2 s。
综上所述,软弱土层对场地加速度反应谱特征周期具有显著影响,经规准后的反应谱特征周期与我国抗震设计规范(中国建筑科学研究院,2016)相比,含软弱土层的场地反应谱特征周期远大于规范取值。本文将以建立的12个场地模型为分析模型,依据软弱土层厚度和埋深对场地加速度反应谱特征周期的影响特征,提出含软弱土层场地加速度反应谱特征周期调整方法,以期为含软弱土层场地地震反应谱特征周期的确定提供理论依据。
从含软弱土层场地取高为h的单位面积土柱(如图8所示),鉴于软弱土层刚度相对于其上覆及下伏土层较低,土柱表面的变形可视为由土柱中的软弱层决定,而忽略上、下部分土体变形。基于这一假定,土柱可简化为如图8所示的弹簧−质量单自由度体系,则这一体系的刚度k和质量m为
图 8 弹簧−质量单自由度体系分析示意图图中h,v和ρ分别表示土层厚度、剪切波速、密度;下标su,s,sd分别表示上覆土层、软弱夹层、下伏土层Figure 8. Schematic diagram of spring mass single degree of freedom system analysisIn the figure,h,v and ρ respectively represent soil layer thickness,shear wave velocity, and density;subscriptsu,s and sd represent the overlying soil layer,weak interlayer,and underlying soil layer respectively$$ k=\frac{{\rho }_{\mathrm{s}}\times {v }_{\mathrm{s}}^{2}}{{h}_{\mathrm{s}}} \text{,} $$ (1) $$ m={h}_{\mathrm{s}}{\rho }_{\mathrm{s}} + {h}_{\mathrm{u}}{\rho }_{\mathrm{u}} \text{,} $$ (2) 式中:ρs,νs和hs分别为软弱土层的密度、剪切波速、厚度;ρu和hu分别为上覆土层密度与厚度。
不难得到这一弹簧−质量单自由度体系的固有周期为:
$$ {T}_{0}=2 \text{π} \sqrt{\frac{{m}}{{k}}} .$$ (3) 结合式(1)—(3)可看出,随软弱土层厚度或上覆土层厚度的增加,弹簧−质量单自由度体系的刚度减小、质量增加,体系的固有周期变大。如此势必影响场地地震反应的频谱特征,继而使场地地震反应的特征周期增大。因此随地震动输入强度的增加,含软弱土层场地特征周期也增加。输入地震动强度越大,浅地表软弱土层应变水平越高,软弱土体非线性越强,其剪切模量比衰减越快,值越小,继而造成剪切模量降低,以致软弱土层剪切刚度越小,如图8所示的弹簧−质量单自由度体系固有周期变大,从而引起含软弱土层场地的特征周期变大。这一定性解释亦适用于上覆土层、软弱夹层和下伏土层的剪切波速(vsu,vs,vsd)差异不是很大的情形下,通常含软弱土层场地地震反应分析结果表明,应变最大值发生在软弱土层中,且软弱土体的剪切模量比衰减较快,由此导致了软弱土层的剪切模量比小于上覆与下伏土层,从而使三者的剪切模量差异变大,且随输入地震动强度的增加,其差异会更显著,使上述假设基本成立。
为了进一步说明三个土层间波速差异不是很大时的合理性,在原模型2 (见表7)的基础上,将淤泥层分别用上覆和下伏粉质黏土替换,构建了两个分析模型(见表8,表9),并计算得到了相应的加速度反应谱,如图9所示。
表 7 原分析模型剖面和力学特性参数Table 7. The original analysis model序号 土层名称 土类号 层底深度
/m层厚
/m剪切波速
/(m·s−1)密度
/(kg·m−3)1 粉质黏土 3 4.5 4.5 160 1860 2 淤泥 1 9.5 5.0 112 1580 3 粉质黏土 4 13.0 3.5 165 1870 4 粉质黏土 5 17.0 4.0 199 1880 5 粉质黏土 6 21.0 4.0 212 1960 6 粉质黏土 7 24.0 3.0 242 1980 7 圆砾 8 27.0 3.0 258 2200 8 全风化安山岩 8 30.0 3.0 393 2250 9 计算基底 9 516 2650 表 8 补充分析模型1Table 8. The supplementary analysis model 1土层序号 土层名称 土类号 层底深度
/m层厚
/m剪切波速
/(m·s−1)密度
/(kg·m−3)1 粉质黏土 3 4.5 4.5 160 1860 2 粉质黏土 3 9.5 5.0 160 1860 3 粉质黏土 4 13.0 3.5 165 1870 4 粉质黏土 5 17.0 4.0 199 1880 5 粉质黏土 6 21.0 4.0 212 1960 6 粉质黏土 7 24.0 3.0 242 1980 7 圆砾 8 27.0 3.0 258 2200 8 全风化安山岩 8 30.0 3.0 393 2250 9 计算基底 9 516 2650 表 9 补充分析模型2Table 9. The supplementary analysis model 2土层序号 土层名称 土类号 层底深度
/m层厚
/m剪切波速
/(m·s−1)密度
/(kg·m−3)1 粉质黏土 3 4.5 4.5 160 1860 2 粉质黏土 4 9.5 5.0 165 1870 3 粉质黏土 4 13.0 3.5 165 1870 4 粉质黏土 5 17.0 4.0 199 1880 5 粉质黏土 6 21.0 4.0 212 1960 6 粉质黏土 7 24.0 3.0 242 1980 7 圆砾 8 27.0 3.0 258 2200 8 全风化安山岩 8 30.0 3.0 393 2250 9 计算基底 9 516 2650 由图9可以看出,在淤泥层波速与上下土层波速差异不明显的情况下,反应谱特征周期具有明显的差异,且随输入地震动强度增加,差异越明显。
上述给出了含软弱土层使场地地震反应特征周期变大的定性解释,为了进一步定量刻画软弱土层厚度及埋深对场地特征周期的影响特征,基于上述各模型分析数据,通过统计回归方法,得到含软弱土层厚度及埋深对场地特征周期的影响规律,同时也取得了输入地震动强度对含软弱土层场地特征周期的影响规律。
由上述定性分析可知,对于特定的含软弱土层场地,其特征周期与场地卓越周期是相关的,为此假定含软弱土层场地特征周期Tg由软土层的固有周期T0、输入地震动特征周期$ {T'_{\mathrm{g}} }$、下伏土体卓越周期T确定,即
$$ {T_{\mathrm{g}}} = a {T_0} + {T'_{\mathrm{g}}} +b T \text{,} $$ (4) 式中,假设$ T ={{4{h_{\mathrm{d}}}}}/{{{v _{{\mathrm{sd}}}}}} $,而${v _{{\mathrm{Sd}}}} $为软弱土层下土体等效剪切波速,hd为软弱土层下伏土体厚度;a和b为回归系数。
由式(4)可知,场地特征周期考虑了输入地震动频谱特性、软弱土层上下伏土体和软弱土层的影响。以表5和表6的数据为基础,通过最小二乘法(高秋英等,2021)回归得到了式(4)中的回归系数a和b (表10),回归结果如图10所示。
表 10 不同输入地震动水平下场地反应谱特征周期拟合结果Table 10. Fitting results of characteristic period of site response spectrum under different input ground motions输入地震动/(cm·s−2) 模型1—6 输入地震动/(cm·s−2) 模型7—12 a b R2 a b R2 25 0.838 −0.066 0.999 80 25 0.972 −0.297 0.999 12 50 0.941 −0.166 0.999 05 50 1.050 −0.410 0.999 64 100 1.231 −0.270 0.999 72 100 1.246 −0.223 0.999 38 图 10 不同输入地震动下含一层 (a) 和含两层 (b) 淤泥的各分析模型场地反应谱特征周期随软弱土层厚度及埋深的变化关系Figure 10. Variation of the characteristic period of the site response spectrum with the thickness and burial depth of the weak soil layer for the analysis models with one-layer (a) and two-layer (b) silt under different input ground motions在统计学中对变量进行线行回归分析,采用最小二乘法进行参数估计时,R2为回归平方和与总离差平方和的比值,表示总离差平方和中可以由回归平方和解释的比例,这一比例越大,模型越精确,回归效果越显著。R2介于0−1之间,越接近1,回归拟合效果越好,一般认为超过0.8的模型拟合优度比较高。在不同输入地震动水平下,软弱土层厚度相同情况下,含软弱土层场地特征周期,随软弱土层埋深的统计关系R2约为1,表明在不同输入地震动水平下,场地地震动加速度反应谱特征周期与软弱土层埋深的相关性较好。模型1—6在不同输入地震动水平下场地地震动加速度反应谱特征周期随软弱土层埋深的变化曲线近似于直线,模型7—12在不同输入地震动水平下场地反应谱特征周期随软弱土层埋深的变化曲线较模型1—6的曲线曲率大,且软弱土层厚度越厚,其固有周期对场地地震动加速度反应谱特征周期的贡献值越大。
此外,依据表6的数据得到了含特定厚度和埋深软弱土层场地的反应谱特征周期Tg与输入地震动峰值加速度A的变化趋势,如图11中的散点所示。由该图可看出,不同模型场地加速度反应谱特征周期Tg随输入峰值加速度增加的变化特征基本相近,且与输入峰值加速度A近似成线性关系。因此,拟通过式Tg=α + βA来拟合场地地震动加速度反应谱特征周期Tg随输入峰值加速度A的变化规律。基于最小二乘法线性回归得到了式中回归系数α和β,如表11所示,回归结果如图11中直线所示。
表 11 不同输入水平下模型1−模型6的反应谱特征周期拟合结果Table 11. Fitting results of characteristic periods of model 1−6 at different input ground motions分析模型 α β R2 模型1 0.501 5 0.000 9 0.962 4 模型2 0.588 4 0.001 0 0.993 7 模型3 0.701 5 0.000 9 0.962 4 模型4 0.751 5 0.000 9 0.962 4 模型5 0.801 5 0.000 9 0.962 4 模型6 0.825 4 0.001 2 0.995 3 回归分析结果R2均接近1,表明场地地震动加速度反应谱特征周期与不同输入地震动峰值加速度具有较好的相关性,对于含特定厚度和埋深软土层场地的反应谱特征周期随输入峰值加速度大致呈线性增长,表明输入地震动峰值加速度越大,场地地震动加速度反应谱特征周期越长。
通过考虑土层结构对场地地震动加速度反应谱特征周期的影响,提出对含软弱土层的场地地震动加速度反应谱特征周期的调整公式,该公式形式简单,物理意义明确,不仅考虑软弱土层厚度与埋深的影响,还考虑了输入不同地震动的影响。利用典型土层结构场地模型进行场地地震反应分析,结果表明,该调整方法不仅合理反映了软弱土层在剖面中的埋深对场地地震动加速度反应谱特征周期的影响规律,且估计值近似等于实际工程场地反应谱特征周期。
4. 讨论与结论
基于建立的12个含软弱土层的场地模型,在不同输入地震动水平下,进行了一维等效线性化场地地震反应分析,讨论了软弱土层厚度和埋深,以及输入地震动强度对场地地震反应的影响及对场地地震动加速度反应谱特征周期的影响特征,提出了含软弱土层场地地震动加速度反应谱特征周期调整方法,获得的主要结论如下:
1) 相同输入地震动水平下,随软弱土层埋深越深,地表峰值加速度越小,且近地表衰减较快,埋深增加,衰减变慢;不同输入地震动水平下,输入峰值加速度越大,地表峰值加速度衰减越快;不同软弱土层厚度对地表峰值加速度亦有影响,厚度越厚,地表峰值加速度越小,浅地表衰减越快,且随输入地震动峰值加速度增加,差异越明显。
2) 相同输入地震动水平下,随软弱土层埋深越深,地表峰值加速度动力放大系数越小,且近地表衰减较快,随埋深增加,衰减变慢;不同输入地震动水平下,输入峰值加速度越大,动力放大系数衰减越快;软弱土层厚度对地表峰值加速度有影响,土层越厚,动力放大系数越小,浅地表衰减越快,且随输入地震动峰值加速度的增加,差异越明显。
3) 随软弱土层的埋深或厚度增加,场地地震动加速度反应谱特征周期逐渐增大;输入地震动峰值增加,场地地震动加速度反应谱特征周期亦增大。继而,根据分析软弱土层厚度、埋深,以及输入地震动强度对场地地震动加速度反应谱特征周期的影响特征,提出了含软弱土层场地地震动加速度反应谱特征周期调整方法。
-
表 1 计算九寨沟地震能量震级Me的单台数据
Table 1 Single station data for calculating energy magnitude Me of Jiuzhaigou earthquake
台站名 震中距/° Me 台站名 震中距/° Me 台站名 震中距/° Me IU.MAKZ 21.46 7.02 IU.MBWA 56.42 6.22 IU.PAB 81.55 6.18 BK.CMB 97.75 6.66 IU.WAKE 57.19 6.47 IU.FUNA 82.63 6.21 II.AAK 24.86 6.70 IU.KBS 57.87 6.24 IU.KIP 84.85 6.31 II.KURK 25.44 6.45 IU.ADK 58.32 6.46 II.TAU 85.86 6.34 IU.MAJO 28.23 6.36 IU.PMG 59.25 6.16 IU.LSZ 86.81 6.12 II.SIMI 28.57 6.50 II.MSEY 59.38 6.17 II.MSVF 87.12 6.66 II.BRVK 31.06 6.19 II.WRAB 60.58 6.06 IU.POHA 87.71 6.53 II.ERM 32.05 6.08 II.KWAJ 63.29 6.58 II.FFC 89.31 6.47 IU.YSS 32.37 6.12 IU.KONO 63.83 6.15 IU.COR 91.20 6.34 IU.DAV 32.98 6.76 II.ALE 64.10 6.08 IU.AFI 92.97 6.75 IU.YAK 33.20 6.30 IU.FURI 64.42 6.24 GT.LBTB 94.43 6.58 II.ARU 38.59 6.09 IU.MIDW 66.27 6.45 II.CMLA 94.85 6.48 IU.MA2 40.40 6.09 IU.CTAO 66.90 6.05 IU.XMAS 96.14 6.28 II.KAPI 41.07 6.31 IU.NWAO 67.31 6.56 IU.KOWA 96.41 6.95 IU.GUMO 42.06 6.53 IU.COLA 67.41 6.12 IU.RAO 96.80 6.22 II.UOSS 42.10 6.25 IU.HNR 68.24 6.57 IU.TSUM 97.39 7.07 IU.PET 43.32 5.89 II.KDAK 68.99 6.36 碾子山 20.30 6.63 IU.GNI 47.05 6.22 II.BFO 69.18 5.94 五常 21.50 6.55 II.KIV 47.89 6.29 IU.KMBO 71.25 5.92 讷河 21.77 6.01 IU.BILL 49.33 6.47 IU.TARA 71.83 6.85 宾县 21.97 6.35 II.DGAR 50.38 6.52 II.ESK 71.91 6.02 宝清 21.98 6.47 II.OBN 50.66 6.18 II.BORG 73.18 6.28 靖宇 21.06 6.14 II.RAYN 51.60 6.23 II.ABPO 75.12 6.18 柳河 20.95 6.21 II.LVZ 51.76 6.03 II.MBAR 76.27 6.63 丰满 21.55 6.42 IU.KEV 54.58 5.97 IU.SFJD 77.95 6.41 抚松 20.65 6.31 IU.ANTO 55.99 6.19 IU.JOHN 78.17 6.01 牡丹江 22.91 6.55 平均能量震级 6.3 注:表中10个中国台站的地震引自国家测震台网数据备份中心(2007)和郑秀芬等(2009)。 表 2 不同机构测定的4次国外和4次国内地震的各类震级
Table 2 Various magnitudes measured by different institutes for four overseas and four domestic earthquakes
发震时间 发震地点 震源深度
/km震源破裂
时间/s震源机制
类型MS
(CENC)MW
(GCMT)Me
(IRIS)Me
(本文)1991-04-22 哥斯达黎加 10.0 81 T 8.1 7.6 7.4 7.4 1992-09-02 尼加拉瓜 15.0 172 T 7.7 7.6 6.7 6.6 2006-11-15 昆士兰 39.0 133 T 8.0 8.3 7.7 7.7 2007-01-13 昆士兰 10.0 101 N 7.9 8.1 8.6 8.1 2008-05-12 四川汶川 12.8 126 TS 8.0 7.9 8.1 8.0 2013-04-20 四川芦山 12.2 39 T 7.0 6.6 6.7 6.6 2014-08-03 云南鲁甸 10.0 46 S 6.5 6.2 6.4 6.4 2017-08-08 四川九寨沟 9.0 43 S 7.0 6.5 6.4 6.3 注:T为逆断层,N为正断层,S为走滑断层,TS为包含走滑成分的逆断层。 -
陈培善. 1990. 地震定量的国际现状[J]. 地震地磁观测与研究, 11(3): 33-39. Chen P S. 1990. The international status of earthquake quantification[J]. Seismological and Geomagnetic Observation and Research, 11(3): 33-39 (in Chinese).
国家测震台网数据备份中心. 2007. 国家测震台网地震波形数据[DB/OL]. 中国地震局地球物理研究所. http://www. seisdmc.ac.cn. doi: 10.11998/SeisDmc/SN. Data Management Centre of China National Seismic Network. 2007. Waveform data of China national seismic network[DB/OL]. Institute of Geophysics, China Earthquake Administration. http://www.seisdmc.ac.cn. doi: 10.11998/SeisDmc/SN (in Chinese).
李志强, 侯建盛, 李洋, 苏桂武, 姜立新, 孙柏涛, 袁一凡, 林均岐. 2013. 2013年4月20日四川芦山MS7.0地震灾害特点分析[J]. 地震地质, 35(2): 398-410. doi: 10.3969/j.issn.0253-4967.2013.02.018 Li Z Q, Hou J S, Li Y, Su G W, Jiang L X, Sun B T, Yuan Y F, Lin J Q. 2013. Analysis on the Characteristics of the MS7.0 Lushan, Sichuan province, earthquake hazard on April 20, 2013[J]. Seismology and Geology, 35(2): 398-410 (in Chinese).
赵仲和. 2014. 2014年8月3日云南鲁甸M6.5地震的地震波能量[J]. 国际地震动态, (9): 24-28, 34. doi: 10.3969/j.issn.0235-4975.2014.09.007 Zhao Z H. 2014. Seismic energy of Ludian M6.5 earthquake, August 3, 2014[J]. Recent Developments in World Seismology, (9): 24-28, 34 (in Chinese).
郑秀芬, 欧阳飚, 张东宁, 姚志祥, 梁建宏, 郑洁.2009. " 国家测震台网数据备份中心”技术系统建设及其对汶川大地震研究的数据支撑[J].地球物理学报, 52(5): 1412—1417. doi: 10.3969/j.issn.0001-5733.2009.05.031 Zheng X F, Ouyang B, Zhang D N, Yao Z X, Liang J H, Zheng J. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake. Chinese Journal of Geophysics, 52(5):1412—1417(in Chinese).
中国地震局. 2013. 中国地震局发布四川省芦山" 4·20”7.0级强烈地震烈度图[EB/OL]. [2018−11−15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5195229/index.html. China Earthquake Administration. 2013. China Earthquake Administration release the seismic intensity map of Lushan Sichuan M7.0 strong earthquake on April 20, 2017[EB/OL]. [2018-11-15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5195229/index.html (in Chinese).
中国地震局. 2017. 中国地震局发布四川九寨沟7.0级地震烈度图[EB/OL]. [2018−11−15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5206273/index.html. China Earthquake Administration. 2017. China Earthquake Administration release the seismic intensity map of Jiuzhaigou Sichuan M7.0 earthquake[EB/OL]. [2018-11-15]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5206273/index.html (in Chinese).
Aki K, Richards P G. 1980. Quantitative Seismology: Theory and Methods[M]. San Francisco: W. H. Freeman: 932.
Ammon C J, Kanamori H, Lay T. 2008. A great earthquake doublet and seismic stress transfer cycle in the central Kuril islands[J]. Nature, 451(7178): 561-565. doi: 10.1038/nature06521
Boatwright J, Choy G L. 1986. Teleseismic estimates of the energy radiated by shallow earthquakes[J]. J Geophys Res, 91(B2): 2095-2112. doi: 10.1029/JB091iB02p02095
Bormann P, Liu R F, Ren X, Gutdeutsch R, Kaiser D, Castellaro S. 2007. Chinese National Network magnitudes, their relation to NEIC Magnitudes, and recommendations for new IASPEI magnitude standards[J]. Bull Seismol Soc Am, 97(1B): 114-127. doi: 10.1785/0120060078
Bormann P, Saul J. 2008. The new IASPEI standard broadband magnitude mB[J]. Seismol Res Lett, 79(5): 698-705. doi: 10.1785/gssrl.79.5.698
Bormann P, Saul J. 2009a. Earthquake magnitude[G]//Encyclopedia of Complexity and Systems Science. Heidelberg: Springer: 2473−2496.
Bormann P, Saul J. 2009b. A fast, non-saturating magnitude estimator for great earthquakes[J]. Seism Res Lett, 80(5): 808-816. doi: 10.1785/gssrl.80.5.808
Choy G L, Boatwright J L. 1995. Global patterns of radiated seismic energy and apparent stress[J]. J Geophys Res, 100(B9): 18205-18228. doi: 10.1029/95JB01969
Choy G L, Boatwright J. 2007. The energy radiated by the 26 December 2004 Sumatra-Andaman earthquake estimated from 10-minute P-wave windows[J]. Bull Seismol Soc Am, 97(1A): S18-S24. doi: 10.1785/0120050623
Convers J A, Newman A V. 2011. Global evaluation of large earthquake energy from 1997 through mid-2010[J]. J Geophys Res, 116(B8): B08304.
Di Bona M, Rovelli A. 1988. Effects of bandwidth limitation on stress drop estimated from integrals of the ground motion[J]. Bull Seismol Soc Am, 78: 1818-1825.
Di Giacomo D, Grosser H, Parolai S, Bormann P, Wang R J. 2008. Rapid determination of Me for strong to great shallow earthquakes[J]. Geophys Res Lett, 35(10): L10308.
Di Giacomo D. 2010. Determination of the Energy Magnitude ME: Application to Rapid Response Purposes and Insights to Regional/Local Variabilities[D]. Potsdam: Universität Potsdam: 49−57.
Duda S J, Yanovskaya T B. 1993. Spectral amplitude-distance curves for P-waves: effects of velocity and Q-distribution[J]. Tectonophysics, 217(3/4): 255-265.
Goes S D B, Velasco A A, Schwartz S Y, Lay T. 1993. The April 22, 1991, Valle de la Estrella, Costa Rica (MW=7.7) earthquake and its tectonic implications: A broadband seismic study[J]. J Geophys Res, 98(B5): 8127-8142. doi: 10.1029/93JB00019
Gutenberg B, Richter C F. 1954. Seismicity of the Earth and Associated Phenomena[M]. 2nd ed. Princeton: Princeton University Press: 1−15.
Haskell N A. 1964. Total energy spectra density of elastic waves for propagating faults[J]. Bull Seismol Soc Am, 54(6): 1811-1841.
Ide S, Beroza G C. 2001. Does apparent stress vary with earthquake size?[J]. Geophys Res Lett, 28(17): 3349-3352. doi: 10.1029/2001GL013106
Kanamori H, Mori J, Hauksson E, Heaton T H, Hutton L K, Jones L M. 1993. Determination of earthquake energy-release and ML using TERRAscope[J]. Bull Seismol Soc Am, 83(2): 330-346.
Kanamori H, Kikuchi M. 1993. The 1992 Nicaragua earthquake: a slow tsunami earthquake associated with subducted sediments[J]. Nature, 361(6414): 714-716. doi: 10.1038/361714a0
Kennett B L N, Engdahl E R, Buland R. 1995. Constraints on seismic velocities in the Earth from traveltimes[J]. Geophys J Int, 122(1): 108-124. doi: 10.1111/gji.1995.122.issue-1
Lomax A, Michelini A, Piatanesi A. 2007. An energy-duration procedure for rapid determination of earthquake magnitude and tsunamigenic potential[J]. Geophys J Int, 170(3): 1195-1209. doi: 10.1111/gji.2007.170.issue-3
Montagner J P, Kennett B L N. 1996. How to reconcile body-wave and normal-mode reference Earth models[J]. Geophys J Int, 125(1): 229-248. doi: 10.1111/gji.1996.125.issue-1
NEIC. 1991. Today in earthquake history[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/learn/today/index.php? month=4&day=22.
NEIC. 1992. M7.7−near the coast of Nicaragua[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/earthquakes/eventpage/ usp0005ddn/executive.
NEIC. 2006. M8.3−Kuril Islands[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/earthquakes/eventpage/usp000exfn/executive.
NEIC. 2007. M8.1−east of the Kuril Islands[EB/OL]. [2018−11−15]. https://earthquake.usgs.gov/earthquakes/eventpage/ usp000f2ab/executive.
Newman A V, Okal E A. 1998. Teleseismic estimates of radiated seismic energy: the E/M0 discriminant for tsunami earthquakes[J]. J Geophys Res, 103(B11): 26885-26898. doi: 10.1029/98JB02236
Parolai S. 2009. Denoising of seismograms using the S transform[J]. Bull Seismol Soc Am, 99(1): 226-234. doi: 10.1785/0120080001
Richter C F. 1935. An instrumental earthquake magnitude scale[J]. Bull Seismol Soc Am, 25(1): 1-32.
Rudnicki J W, Freund L B. 1981. On energy radiation from seismic sources[J]. Bull Seismol Soc Am, 71(3): 583-595.
Singh S K, Ordaz M. 1994. Seismic energy release in Mexican subduction zone earthquakes[J]. Bull Seismol Soc Am, 84(5): 1533-1550.
Stockwell R G, Mansinha L, Lowe R P. 1996. Localization of the complex spectrum: the S transform[J]. IEEE Trans Signal Process, 44(4): 998-1001. doi: 10.1109/78.492555
Venkataraman A, Kanamori H. 2004a. Effect of directivity on estimates of radiated seismic energy[J]. J Geophys Res, 109(B4): B04301.
Wang R J. 1999. A simple orthonormalization method for stable and efficient computation of Green’s functions[J]. Bull Seismol Soc Am, 89(3): 733-741.
-
期刊类型引用(1)
1. 李祥秀,范世凯,李小军,刘爱文. 场地地震动特征周期对高层建筑结构工程材料用量和破坏状态影响的研究. 地震科学进展. 2024(08): 497-505 . 百度学术
其他类型引用(0)