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

基于地震动反应谱与傅里叶谱确定场地卓越频率的差异性

李小军 钱玉 荣棉水 孔小山

李小军,钱玉,荣棉水,孔小山,2024. 基于地震动反应谱与傅里叶谱确定场地卓越频率的差异性. 震灾防御技术,19(3):421−435. doi:10.11899/zzfy20240301. doi: 10.11899/zzfy20240301
引用本文: 李小军,钱玉,荣棉水,孔小山,2024. 基于地震动反应谱与傅里叶谱确定场地卓越频率的差异性. 震灾防御技术,19(3):421−435. doi:10.11899/zzfy20240301. doi: 10.11899/zzfy20240301
Li Xiaojun, Qian Yu, Rong Mianshui, Kong Xiaoshan. Difference Analysis of Site Dominant Frequencies Obtained from Response Spectra and Fourier Spectra of Earthquake Motion[J]. Technology for Earthquake Disaster Prevention, 2024, 19(3): 421-435. doi: 10.11899/zzfy20240301
Citation: Li Xiaojun, Qian Yu, Rong Mianshui, Kong Xiaoshan. Difference Analysis of Site Dominant Frequencies Obtained from Response Spectra and Fourier Spectra of Earthquake Motion[J]. Technology for Earthquake Disaster Prevention, 2024, 19(3): 421-435. doi: 10.11899/zzfy20240301

基于地震动反应谱与傅里叶谱确定场地卓越频率的差异性

doi: 10.11899/zzfy20240301
基金项目: 国家自然科学基金资助项目(52192675);国家重点研发计划项目课题(2022YFC3003503)
详细信息
    作者简介:

    李小军,男,博士。教授,博士生导师。主要从事地震工程方面的研究。E-mail:beerli@vip.sina.com

    通讯作者:

    荣棉水,男,博士。教授,博士生导师。主要从事地震工程方面的研究。E-mail:waltrong@126.com

Difference Analysis of Site Dominant Frequencies Obtained from Response Spectra and Fourier Spectra of Earthquake Motion

  • 摘要: 强震动记录的HVSR法常用于评估场地卓越频率,计算强震动记录HVSR时常采用加速度反应谱或加速度傅里叶谱,但两者会给出不同的评估值。为揭示反应谱比和傅里叶谱比评估场地卓越频率的差异,本文选取日本KiK-net台网中场地条件可近似为一维场地模型的16个台站,以其获取的强震动记录开展场地卓越频率研究。首先提出了评估场地卓越频率的数据处理方法,主要包括S波截取、Taper预处理、基于高斯拟合的自动寻峰。探讨并给出了阻尼比、平滑的带宽系数取值对场地卓越频率评估的影响规律;对反应谱阻尼比取10%,对傅里叶谱平滑的带宽系数取20~40之间获取的场地卓越频率较为准确。然后对比分析了利用地震动加速度反应谱比和傅里叶谱比得到的场地卓越频率与场地土层模型计算得到的基于传递函数的自振频率。研究结果表明,对大多数台站而言,采用傅里叶谱比计算场地卓越频率具有明显的优势,对于Ⅱ、Ⅲ、Ⅳ类场地上的台站均有如此结论,只有对少数特定台站,采用反应谱比方法效果更好。
  • 地震引起工程场地地震动,导致工程破坏和其他地震灾害,震级、断层规模、场地效应和距离是影响地震动的最重要因素(Chen等,2022a),其中地表局部场地条件会导致地震动的显著变化。因此,探讨影响地震动的场地动力特征和表征参数一直是关注和重点研究的内容,场地卓越周期是其中反映场地土地震动力特性的重要参数(高广运等,2000)。目前确定场地卓越周期的方法主要有3种(童广才等,2000):①当场地内有强震动记录时,通过频谱分析确定;②由常时微动测试分析确定;③根据场地分层剪切波速等勘测资料计算确定。然而,利用强震动观测记录分析能获取更准确的场地影响特征(陈永新等,2016Chen等,2022b),强震动记录得到的卓越周期真实反映了地震动的卓越周期(陈鹏等,2009)。陈永新等(2016)整理分析了3种基于强震动记录获取场地卓越周期的方法:①地表记录的傅里叶谱分析法;②地表水平/垂直傅里叶谱比法(任叶飞等,2013);③地表/地下傅里叶谱比法,研究表明3种方法基本可以得到较为准确地结果,但对于不同场地不同方法所得场地卓越周期值存在一定的差异。Nakamura(1989)提出了谱比法以评估近地表场地动力特性,该方法早期利用地脉动记录及微震记录,后由Yamazaki等(1997)推广到了利用强震动加速度记录。罗桂纯等(2014)、荣棉水等(2016)开展了HVSR法适应性研究,都认为HVSR方法能较准确的确定场地卓越周期。在基于HVSR法得到的场地卓越周期研究中,如何降低场地卓越周期值评估的不确定性成为被关注问题。利用地震动加速度记录谱比方法计算场地卓越周期还存在另一方面的问题,即利用地震动加速度反应谱比是否比利用傅里叶谱比更有效。姚鑫鑫(2017)以傅里叶谱比曲线确定卓越频率研究了汶川地震强余震场地非线性反应;刘宇实等(2018)、李文倩等(2019)通过傅里叶谱计算HVSR的卓越频率研究了其与覆盖层厚度、等效剪切波速等场地特征参数之间的相关性。传统谱比法更多的是采用傅里叶谱研究场地效应,而Zhao等(2006)认为5%阻尼比的反应谱比也可用于评估场地响应特征;张立宝(2018)、丁毅等(2021)等利用地表/井下地震动记录的反应谱比值得到的卓越频率变化来研究表征场地土非线性的参数。袁一凡等(2012)阐明了傅里叶谱与反应谱在谱值曲线上的相似性,但没有进一步探究傅里叶谱与反应谱计算场地卓越频率的差异;罗桂纯等(2019)基于HVSR法研究场地反应非线性特征时,浅析了傅里叶谱和反应谱在识别场地卓越频率方面存在的问题,但并没有对如何获取准确的场地卓越频率及两者获取途径的差异性进行深入探究。继Thompson等(2012)基于日本KiK-net强震动数据统计识别了16个场地可以简化为一维土层模型的观测台站,一些学者提出通过这些台站计算场地的传递函数,所得的自振频率可以作为场地实际卓越频率的参考值。研究还发现,基于强震动数据计算HVSR的过程中数据处理的很多环节会影响HVSR的计算结果,姚鑫鑫等(2019) 探究了HVSR法中对地震动记录进行S波截取和Taper预处理的必要性。针对谱比曲线的有效峰值识别及峰值点准确定位问题,Zhu等(2021)定义了一个有效峰选取标准,并通过自动寻峰程序对谱比图进行有效峰值的选取,但其选取的显著峰会受到谱比数据局部误差的影响,而出现多个显著峰,从而影响卓越峰值的确定。Hassani等(2016)利用高斯函数的优势进行曲线拟合提高了峰值点的准确性,但并未深入探究高斯拟合阶数以及拟合区段的划分对峰值点准确性的影响。

    从上述研究来看,数据处理环节对场地卓越频率评价结果有较大影响,需要进一步探讨场地卓越频率评价中数据处理环节的相关问题并完善处理方法。本文针对日本KiK-net台网中16个可被近似作为一维土层模型场地的观测台站开展强震动记录分析,首先对这些观测台站的强震动记录进行有效的预处理,探讨并给出了阻尼比(Damp ratio)、平滑的带宽系数 B 的取值对场地卓越频率识别结果的影响,而后基于在MATLAB平台对台站场地地表观测记录的HVSR曲线进行自动寻峰与高斯曲线拟合处理,探究高斯拟合的阶数和区段选取的影响。在进行一套完整的处理后,分别得到HVSR反应谱比、傅里叶谱比的卓越频率,再与一维场地土层模型得到的自振频率做比较分析,并论证两者在不同场地分类台站中的差异性。本文研究旨在从强震动数据处理、评估方法等角度获取较为准确的场地卓越周期值,进而在实际应用上更好地考虑地表地震动受场地条件影响的程度以及研究土体场地的放大效应。

    研究选取Thompson等(2012)确定的日本KiK-net台网中可被近似作为一维土层模型场地的16个台站,依次计算每个台站场地30 m厚度的平均剪切波速,并利用美国场地分类标准对不同钻孔深度的场地类型进行划分,同时结合覆盖土层厚度,依据GB 50011—2010《建筑抗震设计规范》、进行划分。台站详细信息及划分结果如表1所示,台站分布如图1所示。考虑到地表地震动记录的峰值加速度(PGA)大于10 Gal时记录受噪声影响较小(Wen等,1994),而地震动强度过大会引起明显场地土非线性效应等因素,最终选取PGA范围在10~20 Gal的16个台站,共459条强震动记录(表1)。

    表 1  16个选定台站详细信息
    Table 1.  Details of the 16 selected stations
    台站名称 台站编号 纬度 经度 覆盖土层厚度/m 钻孔深度/m 场地分类 中国场地
    分类类别
    选取强震动
    记录数量/个
    VS30/(m·s−1) 类别
    TAKAHAGI IBRH13 66°33'N 140°57'E 24 100 144 E 61
    KASUMIGAURA IBRH17 36°08'N 140°31'E 235 510 335 D 38
    TAMAYAMA IWTH02 39°82'N 141°38'E 19 102 168 E 55
    KUJI-N IWTH08 40°26'N 141°78'E 20 100 301 D 27
    KANEGASAKI IWTH24 39°19'N 141°01'E 56 150 390 C 12
    RIKUZENTAKATA IWTH27 39°03'N 141°53'E 4 100 375 C 50
    TSURUI-E KSRH06 43°22'N 144°42'E 70 237 240 D 21
    HAMANAKA KSRH10 43°20'N 145°11'E 36 255 486 C 9
    KAWANISHI NIGH11 37°17'N 138°74'E 56 205 237 D 4
    UJIIE TCGH12 36°69'N 139°98'E 50 120 670 C 39
    TAIKI TKCH08 42°48'N 143°15'E 36 100 326 D 21
    YABUKI FKSH11 37°20'N 140°33'E 86 115 204 D 41
    TSURUI-S KSRH07 43°13'N 144°32'E 82 222 305 D 22
    IWAKI-E FKSH14 37°02'N 140°97'E 52 147 213 D 24
    BEKKAI-E NMRH04 43°39'N 145°12'E 186 216 353 D 16
    ISHIGE IBRH10 36°11'N 139°98'E 190 900 344 D 19
    下载: 导出CSV 
    | 显示表格
    图 1  16个近似一维土层模型台站分布情况
    Figure 1.  Distribution of 16 approximate 1D stations

    强震动数据预处理有多个重要环节,限于篇幅本文仅介绍S波截取、Taper预处理、自动寻峰与高斯拟合结合处理。

    考虑到利用地表强震动记录的全时程计算HVSR,可能引起不合理的场地效应分析结果,针对强震动记录进行S波提取处理。同时强震动记录的HVSR计算过程中,截取S波时会出现数据首尾两端不为零且数值较大的现象。为此我们在截取S波窗口后的强震动记录首尾两端各加 5% S 波时长的余弦窗进行Taper预处理,能有效消除边瓣效应的影响。具体方法参考姚鑫鑫等(2019)的研究结果。S波截取、Taper预处理如图2所示,强震动数据处理流程如图3所示。

    图 2  S波的截取与 Taper预处理
    Figure 2.  Taper pre-processing and S-wave interception
    图 3  强震动数据处理流程图
    Figure 3.  Flow chart of strong vibration data processing

    针对HVSR曲线,可分别从归一化振幅、显著度、突出度及锐度4种统计特性进行显著峰筛选。Zhu等(2021)通过这些特性设计了一种新的显著峰识别方案(图4),该方案中,当某个局部峰通过3个选择标准,才能被视为1个显著峰:(1)局部某峰值ai应超过全局阈值$ \text{max}({C}_{1},{C}_{2}+{\text{log}}_{10}\overline{HVSR}) $(系数C1C2分别设定为lg(2.2)和lg(1.4));(2)突出度Pi超过设定的阈值C3(阈值被设置为$ {\text{lg}}1.8 $,相当于75%突出值的平均值);(3)锐度(突出度/宽度wi)应超过设定的阈值C4,此处依据跟踪误差法将C4设置为 0.5。依据显著峰识别方案,可以筛选出HVSR曲线的众多显著峰,再将选取的显著峰中最大的显著峰值所对应的频率作为卓越频率。但Zhu等(2021)在HVSR曲线的显著峰筛选过程中由于数据处理误差并不能完全通过方案要求。为避免尖而窄的显著峰对卓越频率确定的影响,本文提出通过自动显著峰选取与峰曲线的高斯拟合结合的处理方法。高斯函数的曲线拟合属于多项式拟合类,其表达式如下:

    图 4  显著峰选取标准示意图
    Figure 4.  Significant peak selection criteria
    $$ Y = a \times \exp \left( { - {{\left( {\frac{{X - b}}{c}} \right)}^2}} \right) $$ (1)

    式中,a初始幅值;b中心横坐标;c标准差。

    高斯函数拟合可以对强震动时程因个别数据点的错误而产生的尖而窄的峰进行误差处理,利用拟合后的曲线以函数形式确定显著峰值对应的频率,可以提高数据可信度。高斯曲线拟合的准确度也是需要关注的问题,应尽可能准确控制高斯拟合的参数。对于不同数据,高斯函数的参数abc不是定值,只能在95%置信区间内取中值,如图5所示。由于地震动数据存在多段显著峰,需要提高相应拟合的阶数,本文考虑的最高阶数为8阶,图5表明不同阶数的拟合程度不同。通常认为在进行高斯拟合时选取高阶数的拟合可使拟合更加准确,但分析发现并不全然这样,针对不同数量峰的情况采取不同阶数的高斯拟合匹配度是不同的,需要通过拟合优度中的SSE(残差平方和)参数大小进行判断并选取相应的阶数,这样得到的高斯拟合效果会更加符合要求,这里SSE值越趋于0拟合程度越好,但面对多峰的数据,由于数据量不够或者峰的可捕捉拟合情况复杂,也会导致SSE值偏大从而影响拟合效果。另一方面,为使拟合效果更优,本研究放弃了原本的全区段拟合,改为采用局部高斯拟合。图6为分段拟合与全区段拟合结果的对比。考虑到相应峰局部区域选取以及数据点数量都会影响拟合效果,采用将高斯拟合SSE值尽可能地控制在1以内,以越趋于0方向越优的处理思路。

    图 5  不同阶数的高斯拟合的效果
    Figure 5.  Effect of Gaussian fitting with different orders
    图 6  全区域高斯拟合和分段局部高斯拟合对比
    Figure 6.  Comparison between full-area Gaussian fitting and local Gaussian fitting

    反应谱与傅里叶谱的处理步骤不同,需要分别进行处理。针对反应谱,首先进行不同阻尼比取值的计算,阻尼比取值越大相应的反应谱曲线平滑程度越高;然后根据HVSR法绘制谱比图,再基于高斯拟合的自动寻峰程序找出显著峰值并标记。针对傅里叶谱,先进行傅里叶变换,并进行Konno-Ohmachi平滑处理;然后根据HVSR法绘制谱比图,在傅里叶谱平滑处理过程中,平滑函数的带宽系数在2~100之间取值,带宽系数过小会导致平滑过度,将影响谱比的幅值,而带宽系数过大则平滑太弱,容易出现极端峰、多相似显著峰,从而影响卓越频率的识别,所以针对16个台站进行不同阻尼比以及不同的平滑带宽系数取值论证,反应谱阻尼比按照1%、5%、10%依次选取,傅里叶谱平滑带宽系数在20~70范围以10为间隔依次选取。不同台站的谱比图如图7图8所示,峰值变化数据如表2表3所示。地震事件用灰细实线表示,均值用粗实线表示,均值加减方差($ \mathrm{A}\mathrm{v}\mathrm{g}\pm \mathrm{S}\mathrm{D} $)用橙色虚线表示,场地土层模型计算得到的传递函数用紫色实线表示,1分段高斯拟合曲线用蓝色实线表示,2分段高斯拟合曲线用红色实线表示。就反应谱的不同阻尼比取值看,不同的阻尼比会使得峰值数据成一定趋势的变化,提高阻尼比取值时对应的反应谱比峰值的对应频率大多呈上升趋势,个别台站呈下降趋势,但几乎都是朝着场地模型计算传递函数相应的自振频率靠近。从数值接近程度看,阻尼比值宜采用10%,从趋势看或许会有疑问是否阻尼比取值越大的峰值对应频率越接近场地模型计算传递函数相应的自振频率。研究中也选取15%的阻尼比进行分析,发现阻尼比为15%时所绘出的反应谱比图峰值平滑过度,显现出与传递函数结果相关性较差。综上分析,在后期研究中均采用10%以内阻尼比值。针对傅里叶谱比曲线平滑的带宽系数取值问题,在了解众多学者平滑的带宽系数取值后,我们把平滑的带宽系数取值范围定在20~70之间。平滑的带宽系数取值的影响分析显示,提高平滑的带宽系数取值时对应的傅里叶谱比峰值对应频率大多呈上升趋势,均朝场地模型计算传递函数相应的自振频率靠近,个别台站呈下降趋势,有远离场地模型计算传递函数相应的自振频率的趋势。从谱比图中谱比峰看,平滑的带宽系数值增加往往会使得谱比图中原本1个峰被分化为2个峰,且2个峰的数值以中心为轴向相反趋势变化。也有平滑的带宽系数值增加但峰的数量并没有变化的情况,此时平滑的带宽系数值不管在20~70之间取多少均相差无几。为此,平滑的带宽系数取峰有分化和未分化的界限值时,平滑曲线与场地模型计算传递函数相关性最好。

    图 7  典型台站不同阻尼比取值的谱比图
    Figure 7.  Spectral ratios of typical stations with different damping ratio values
    图 8  典型台站不同平滑的带宽系数取值谱比图
    Figure 8.  Spectrogram of different values of smooth bandwidth factor for typical stations
    表 2  典型台站不同阻尼比取值峰值频率变化
    Table 2.  Peak frequency variation of typical stations with different damping ratio values
    台站名 阻尼比 峰值频率/Hz 传递函数峰值频率/Hz
    IWTH02 0.01 5.3413 4.91
    0.05 5.1903
    0.1 5.0582
    KSRH07 0.01 2.4020 2.51
    0.05 2.4066
    0.1 2.4128
    IWTH08 0.01 2.4883 3.37
    0.05 2.4916
    0.1 2.5366
    下载: 导出CSV 
    | 显示表格
    表 3  典型台站不同平滑的带宽系数取值峰值频率变化
    Table 3.  Peak frequency variation of typical stations with different values of smooth bandwidth factor
    台站名 平滑的带宽系数 峰值频率/Hz 传递函数峰值频率/Hz
    IWTH02205.93104.91
    306.0793
    406.0942
    506.1148
    606.1354
    706.1580
    KSRH07202.57782.51
    302.5740
    402.5759
    502.5791
    602.5745
    702.5781
    IWTH08202.73083.37
    302.7280
    402.7157
    502.7268
    602.6971
    702.7232
    下载: 导出CSV 
    | 显示表格

    从研究可知,每个台站强震动数据绘出的反应谱比图与傅里叶谱比图峰的数量不同,反应谱取不同的阻尼比、傅里叶谱取不同的平滑的带宽系数都会导致台站谱比曲线峰的数量发生变化,本文对16个一维台站相应峰的数量进行了统计、归类(表4图9)。

    表 4  台站峰值数量统计
    Table 4.  Peak number statistics of 16 one-dimensional stations
    台站名(场地类别)显著峰数量/个峰的数量变化
    阻尼比平滑的带宽系数
    0.010.050.1203040506070
    IWTH02(Ⅱ)1111111122
    IBRH13(Ⅱ)1211112333
    IWTH27(Ⅱ)1111111111
    KSRH06(Ⅱ)1111122233
    KSRH10(Ⅱ)2432223344
    TCGH12(Ⅱ)2211233444
    IWTH08(Ⅱ)3333334445
    TKCH08(Ⅱ)3553344667
    KSRH07(Ⅲ)2222222222
    FKSH11 (Ⅲ)3432335555
    FKSH14(Ⅲ)5544456678
    IBRH10(Ⅳ)4432445555
    FKSH14(Ⅲ)5544456678
    NMRH04(Ⅳ)5742556778
    IBRH17(Ⅱ)6922689999
    NIGH11(Ⅱ)64435667811
    下载: 导出CSV 
    | 显示表格
    图 9  不同数量峰的示意图
    Figure 9.  Classification of different peak numbers

    表4可知,16个台站记录均显现出不同的显著峰数量,Ⅱ类场地的主要呈现为单峰形式,同时也包含了双峰和多峰形式;Ⅲ类与Ⅳ类场地仅显现出多峰形式。已有研究表明,场地条件对场地卓越频率及附近范围频率成分的影响较为显著,基于场地卓越频率(第一自振频率)间接评估场地地震反应的方法具有其合理性,但场地第二阶及高阶自振频率对地震动的影响也不能忽视,对于某些特定场地,高阶自振频率的作用较第一自振频率更大;而且考虑到较软弱场地的第一自振频率较小,可能位于基岩地震动主要频率成分范围外,为此在很多情况下,对场地表面地震动特性起主要作用的是第二、三阶自振频率。因此以场地第一自振频率作为场地影响的主要频率并不具有普遍性,如图9所示。从强震动记录的反应谱比或傅里叶谱比图不难看出,不同场地的显著峰数量不同,相对应的场地影响的主要频率也不唯一,可以粗略估计Ⅱ类场地大多数成单峰形式主要是第一自振频率作为场地影响的主要频率;个别双峰及多峰形式依然以第一自振频率作为场地的主要影响频率,唯独台站TCGH08呈多峰形式但是以第三自振频率作为场地影响的主要频率。Ⅲ类场地通常以多峰形式出现,从FKSH14可看到其以第一自振频率作为场地影响的主要频率。Ⅳ场地也通常以多峰形式出现,但其以第二自振频率作为场地影响的主要频率。本文的强震动统计分析结果进一步证实了作者早年的一些研究结果(李小军,1992)。

    本文最关键的工作是论证基于加速度反应谱(SA)与加速度傅里叶谱(FAS)的场地卓越周期评估方法的差异性和优劣。从目前研究现状看,暂无学者论证2种强震动数据处理方法获取场地卓越周期的具体差异,同时也没有确切比较二者优劣的准则。本研究选用的16个台站均为LG观测台站。LG由Thompson等(2012)在其台站特征分类研究中提出,用以区分不同台站场地可近似为一维场地的程度,L(Low)代表同一台站地震事件之间的低差异性,G(Good)代表ETF(经验传递函数)与SH波入射下的一维场地土层模型(SH1 D)计算TTF(理论传递函数)高匹配度。Thompson等(2012)指出当观测台站场地符合一维场地模型SH波传播假设,理论与观测的数据之间会显现出较好的相关性。Lermo等(1993)利用强震动记录对墨西哥3个场地分别进行了HVSR与经验传递函数的比较研究,结果表明二者吻合较好,因此推断强震动记录中S波部分的HVSR可以被用于估计场地的经验传递函数。Theodulidis等(1996)在现场波速测试与场地19 m深度内钻孔取样波速测定的基础上,给出了Garner Valley Downhole Array(GVDA)场地的一维模型,并将该模型的理论传递函数与28次地震事件的经验传递函数均值曲线进行了比较,二者吻合良好。传递函数是研究场地放大效应公认的合理手段,其峰值频率即为场地卓越频率,学者们也证实HVSR方法能准确地揭示场地卓越频率。鉴于理论传递函数、经验传递函数、HVSR三者之间在一维场地模型下的良好对应性,同时传递函数与HVSR方法都能较准确地获取场地卓越频率且结果大体一致,本文提出将由场地钻孔数据获取的理论传递函数作为衡量不同途径HVSR方法得到场地卓越频率优劣的准则。这一准则下的方法比较也只是一种探讨性研究,希望这一探讨性研究能为相关工作提供参考价值。作为比较SA与FAS方法之间的差异性准则,就准则的可信度是否满足要求,本研究将通过台站一维场地模型进行分析论证。选择皮尔逊样本相关系数r,以衡量主要测算峰值相关程度。谱比曲线的峰值频率能与传递函数对应较好的台站如图10所示(r代表全区域显著峰的相关系数,r1、r23分别代表显著峰1、显著峰峰23的相关系数),且相应台站的场地类别也包含了3种场地类别,其中Ⅱ类场地数量最多,符合实际的场地类别分布情况,说明以传递函数作为标准具有一定的普遍性。部分台站(IBRH13、KSRH07、FKSH14台站)多个峰值频率与传递函数高度对应,说明整体与传递函数有着高度相关性。有些台站只在个别峰值频率上有高度相关性,如KSRH10台站在首峰频率对应良好,IBRH10台站第二、三峰频率对应良好,这些研究结果为后续SA与FAS方法的比较提供了较为有力的支撑。

    图 10  不同台站谱比曲线与传递函数对应状况
    Figure 10.  Correspondence between spectral ratio curves and transfer functions at different stations

    利用前述比较准则,针对16个台站进行SA与FAS的全段或局部比较,相应结果如表5所示。对表5的结果分析可得到以下结论:

    表 5  16个一维观测台站峰值数据统计
    Table 5.  Peak number statistics of 16 one-dimensional stations
    场地类别 显著峰数量/个 反应谱最优参数对应的峰值频率 傅里叶谱最优参数对应的峰值频率 理论传递函数对应的
    峰值频率/Hz
    阻尼比 峰值频率/Hz 带宽系数 峰值频率/Hz
    IBRH13(Ⅱ) 1 0.1 2.5910 20 3.0665 3.11
    IWTH02(Ⅱ) 1 0.1 5.0582 20 5.9310 4.91
    IWTH08(Ⅱ) 1 0.01 2.5883 20 2.7308 3.37
    2 0.1 6.4100 40 8.0736 9.23
    3 0.1 10.9237 40 15.3087 15.01
    IWTH27(Ⅱ) 1 0.01 6.1773 20 7.7455 9.4
    KSRH10(Ⅱ) 1 0.1 1.8353 20 1.9810 1.99
    IBRH17(Ⅱ) 1 0.01 0.3534 20 0.3504 0.39
    2 0.01 0.9096 20 0.9047 1.06
    3 0.01 1.3833 20 1.4935 1.5
    IWTH24(Ⅱ) 1 0.1 2.5140 20 2.7376 3.37
    2 0.1 6.3388 40 8.1129 9.23
    3 0.1 10.8505 20 15.3724 15.01
    KSEH06(Ⅱ) 1 0.1 5.7812 20 6.3617 6.65
    NIGH11(Ⅱ) 2 0.1 4.5118 20 4.6040 3.9
    3 0.1 6.8949 20 8.4756 6.09
    4 0.1 10.01 20 14.5878 9.36
    TCGH12(Ⅱ) 1 0.1 5.2806 20 6.6765 7.22
    2 0.05 7.2384 20 9.3454 9.25
    TKCH08(Ⅱ) 3 0.1 6.7590 20 8.1318 8.35
    FKSH11 (Ⅲ) 1 0.1 1.5472 60 1.6730 2.12
    2 0.1 4.5552 70 5.2341 5.86
    3 0.01 6.8218 50 8.0801 9.61
    KSRH07(Ⅲ) 1 0.01 2.4128 60 2.5635 2.51
    2 0.01 6.9196 20 8.8538 7.02
    FKSH14(Ⅲ) 1 0.1 1.1760 40 1.1875 1.35
    2 0.1 3.6658 20 3.9909 4.01
    3 0.1 5.2186 40 5.9271 6.4
    4 0.1 6.3233 40 7.7832 8.83
    NMRH04(Ⅳ) 3 0.1 0.4053 20 0.4089 0.61
    4 0.1 1.7428 30 1.9726 1.61
    5 0.1 2.3917 20 4.3434 2.04
    IBRH10(Ⅳ) 1 0.01 0.7869 40 0.7931 0.78
    2 0.01 1.1389 20 1.1897 1.27
    下载: 导出CSV 
    | 显示表格

    (1)单峰的场地卓越频率比较中,IBRH13(Ⅱ)、IWTH27(Ⅱ)、KSRH10(Ⅱ)台站的f1(基阶频率)比较显示FAS较有优势;IWTH02(Ⅱ)台站的f1比较显示SA较有优势。

    (2)KSRH07(Ⅲ)双峰的场地卓越频率比较中,KSRH07(Ⅲ)台站的f1比较显示FAS有优势,f2(第二阶频率)比较显示SA有优势;IBRH10(Ⅳ)台站的f1比较显示SA较有优势,f2比较显示FAS较有优势。

    (3)IWTH08(Ⅱ)、FKSH14(Ⅲ)、FKSH11(Ⅲ)台站三峰的场地卓越频率比较显示,FAS与传递函数相关性更好,FAS的f1f2f3(第三阶频率)也更加合理准确。

    (4)从场地类别看,除KSRH07(Ⅲ)台站的f2、IWTH02(Ⅱ)和IBRH10(Ⅳ)的f1场地外,SA的结果占据了优势;其余Ⅱ、Ⅲ、Ⅳ类场地台站均是FAS占据绝对优势。

    综上,在实际工程中只要规避一些特定场地,FAS更适用于评估场地卓越频率或其范围,在Ⅱ类场地上效果最佳。

    基于SA与FAS的场地卓越频率评估方法差异性分析表明,FAS方法较SA方法更具优势,已有研究也说明HVSR法是目前公认的评估场地卓越频率方法中更切实际的方法。一些学者认为卓越频率并不单指HVSR图中幅值最高的显著峰对应的频率,卓越频率也可能存在于其相邻的显著峰中。针对上述情况,本文按照峰的数量分析,结合传递函数的对应标准,对只有单峰形式且与传递函数高度相关的台站获取其场地卓越频率;对双峰或多峰形式,按照其对应的多个显著峰给出一定可信度的场地卓越频率范围(表6)。由表6可知,单峰形式均是Ⅱ类场地,在工程运用中遇到单峰的II类场地则可通过FAS法确定一个较为准确的场地卓越频率;针对双峰及多峰形式的场地情况,可给出一定可信度的场地卓越频率范围。

    表 6  场地卓越频率或其范围统计
    Table 6.  Frequency of excellence or its range statistics
    场地类别 显著峰数量/个 卓越频率值/Hz 卓越频率值范围/Hz 选取方法
    IWTH02(Ⅱ) 1 5.05 SA
    KSRH10(Ⅱ) 1 1.98 FAS
    IWTH27(Ⅱ) 1 7.74 FAS
    IBRH13(Ⅱ) 1 3.06 FAS
    KSRH06(Ⅱ) 1 6.36 FAS
    KSRH07(Ⅲ) 2 2.560~8.850 FAS
    FKSH14(Ⅲ) 2 1.180~7.780 FAS
    IBRH10(Ⅳ) 2 0.790~1.180 FAS
    TCGH12(Ⅱ) 2 0.107~0.381 FAS
    IWTH08(Ⅱ) 3 2.730~8.070 FAS
    TKCH08(Ⅱ) 3 0.123~0.551 FAS
    IWTH24(Ⅱ) 3 0.084~0.593 FAS
    FKSH11 (Ⅲ) 3 1.670~8.00 FAS
    NMRH04(Ⅳ) 5 0.263~2.467 SA
    NIGH11(Ⅱ) 6 0.097~0.355 SA
    IBRH17(Ⅱ) 6 0.163~0.297 FAS
    下载: 导出CSV 
    | 显示表格

    本文针对不同途径的HVSR法获取场地卓越频率的方法开展了研究,以日本KiK-net台网中16个可近似为一维场地模型的台站数据为研究对象,重点关注了数据处理过程中的影响因素控制,规划了一套数据处理流程,在此基础上提出了评估场地卓越频率或其范围的方法,并得到以下结论:

    (1)HVSR方法可用以评估场地卓越频率,但获取的场地卓越频率准确性有待提高。采用S波截取、Taper预处理、高斯拟合的自动寻峰程序完善数据处理流程十分必要。

    (2)多数情况下HVSR方法中反应谱计算取10%阻尼比较为合适;Ⅱ、Ⅲ、Ⅳ类场地下,傅里叶谱计算平滑的带宽系数在20~40之间取值较为合适,由此得到的谱比曲线能更容易识别和准确获取场地卓越频率。

    (3)从KiK-net强震动数据的显著峰数量看,影响地震动的场地自振频率不唯一。对于Ⅱ类场地台站,谱比曲线多数以单峰形式出现,场地影响的主要频率基本上是第一阶自振频率;对于Ⅲ类场地台站,谱比曲线以多峰形式出现,场地影响的主要频率基本上还是第一阶自振频率;对于Ⅳ类场地台站,谱比曲线常以多峰形式出现,但场地影响的主要频率是第二阶自振频率。

    (4)基于SA与FAS的场地卓越频率评估方法比较发现,FAS方法在大多数场地条件下优于SA方法,SA方法只在少数台站场地上具有优势。建议在评估场地卓越频率或了解其范围时,优先选择FAS方法。

    本文研究采用的强震动记录峰值加速度均介于10~20 Gal范围,以避免场地土体非线性特性的影响。如果要进一步采用峰值加速度大于20 Gal的强震动记录进行研究,则需考虑场地土体的动力非线性影响问题,这一研究将更为复杂,这也是我们后续要关注工作。

    致谢 感谢日本防灾科学技术研究所(NIED)所属的KiK-net台网提供的强震动记录数据与场地钻孔数据。

  • 图  1  16个近似一维土层模型台站分布情况

    Figure  1.  Distribution of 16 approximate 1D stations

    图  2  S波的截取与 Taper预处理

    Figure  2.  Taper pre-processing and S-wave interception

    图  3  强震动数据处理流程图

    Figure  3.  Flow chart of strong vibration data processing

    图  4  显著峰选取标准示意图

    Figure  4.  Significant peak selection criteria

    图  5  不同阶数的高斯拟合的效果

    Figure  5.  Effect of Gaussian fitting with different orders

    图  6  全区域高斯拟合和分段局部高斯拟合对比

    Figure  6.  Comparison between full-area Gaussian fitting and local Gaussian fitting

    图  7  典型台站不同阻尼比取值的谱比图

    Figure  7.  Spectral ratios of typical stations with different damping ratio values

    图  8  典型台站不同平滑的带宽系数取值谱比图

    Figure  8.  Spectrogram of different values of smooth bandwidth factor for typical stations

    图  9  不同数量峰的示意图

    Figure  9.  Classification of different peak numbers

    图  10  不同台站谱比曲线与传递函数对应状况

    Figure  10.  Correspondence between spectral ratio curves and transfer functions at different stations

    表  1  16个选定台站详细信息

    Table  1.   Details of the 16 selected stations

    台站名称 台站编号 纬度 经度 覆盖土层厚度/m 钻孔深度/m 场地分类 中国场地
    分类类别
    选取强震动
    记录数量/个
    VS30/(m·s−1) 类别
    TAKAHAGI IBRH13 66°33'N 140°57'E 24 100 144 E 61
    KASUMIGAURA IBRH17 36°08'N 140°31'E 235 510 335 D 38
    TAMAYAMA IWTH02 39°82'N 141°38'E 19 102 168 E 55
    KUJI-N IWTH08 40°26'N 141°78'E 20 100 301 D 27
    KANEGASAKI IWTH24 39°19'N 141°01'E 56 150 390 C 12
    RIKUZENTAKATA IWTH27 39°03'N 141°53'E 4 100 375 C 50
    TSURUI-E KSRH06 43°22'N 144°42'E 70 237 240 D 21
    HAMANAKA KSRH10 43°20'N 145°11'E 36 255 486 C 9
    KAWANISHI NIGH11 37°17'N 138°74'E 56 205 237 D 4
    UJIIE TCGH12 36°69'N 139°98'E 50 120 670 C 39
    TAIKI TKCH08 42°48'N 143°15'E 36 100 326 D 21
    YABUKI FKSH11 37°20'N 140°33'E 86 115 204 D 41
    TSURUI-S KSRH07 43°13'N 144°32'E 82 222 305 D 22
    IWAKI-E FKSH14 37°02'N 140°97'E 52 147 213 D 24
    BEKKAI-E NMRH04 43°39'N 145°12'E 186 216 353 D 16
    ISHIGE IBRH10 36°11'N 139°98'E 190 900 344 D 19
    下载: 导出CSV

    表  2  典型台站不同阻尼比取值峰值频率变化

    Table  2.   Peak frequency variation of typical stations with different damping ratio values

    台站名 阻尼比 峰值频率/Hz 传递函数峰值频率/Hz
    IWTH02 0.01 5.3413 4.91
    0.05 5.1903
    0.1 5.0582
    KSRH07 0.01 2.4020 2.51
    0.05 2.4066
    0.1 2.4128
    IWTH08 0.01 2.4883 3.37
    0.05 2.4916
    0.1 2.5366
    下载: 导出CSV

    表  3  典型台站不同平滑的带宽系数取值峰值频率变化

    Table  3.   Peak frequency variation of typical stations with different values of smooth bandwidth factor

    台站名 平滑的带宽系数 峰值频率/Hz 传递函数峰值频率/Hz
    IWTH02205.93104.91
    306.0793
    406.0942
    506.1148
    606.1354
    706.1580
    KSRH07202.57782.51
    302.5740
    402.5759
    502.5791
    602.5745
    702.5781
    IWTH08202.73083.37
    302.7280
    402.7157
    502.7268
    602.6971
    702.7232
    下载: 导出CSV

    表  4  台站峰值数量统计

    Table  4.   Peak number statistics of 16 one-dimensional stations

    台站名(场地类别)显著峰数量/个峰的数量变化
    阻尼比平滑的带宽系数
    0.010.050.1203040506070
    IWTH02(Ⅱ)1111111122
    IBRH13(Ⅱ)1211112333
    IWTH27(Ⅱ)1111111111
    KSRH06(Ⅱ)1111122233
    KSRH10(Ⅱ)2432223344
    TCGH12(Ⅱ)2211233444
    IWTH08(Ⅱ)3333334445
    TKCH08(Ⅱ)3553344667
    KSRH07(Ⅲ)2222222222
    FKSH11 (Ⅲ)3432335555
    FKSH14(Ⅲ)5544456678
    IBRH10(Ⅳ)4432445555
    FKSH14(Ⅲ)5544456678
    NMRH04(Ⅳ)5742556778
    IBRH17(Ⅱ)6922689999
    NIGH11(Ⅱ)64435667811
    下载: 导出CSV

    表  5  16个一维观测台站峰值数据统计

    Table  5.   Peak number statistics of 16 one-dimensional stations

    场地类别 显著峰数量/个 反应谱最优参数对应的峰值频率 傅里叶谱最优参数对应的峰值频率 理论传递函数对应的
    峰值频率/Hz
    阻尼比 峰值频率/Hz 带宽系数 峰值频率/Hz
    IBRH13(Ⅱ) 1 0.1 2.5910 20 3.0665 3.11
    IWTH02(Ⅱ) 1 0.1 5.0582 20 5.9310 4.91
    IWTH08(Ⅱ) 1 0.01 2.5883 20 2.7308 3.37
    2 0.1 6.4100 40 8.0736 9.23
    3 0.1 10.9237 40 15.3087 15.01
    IWTH27(Ⅱ) 1 0.01 6.1773 20 7.7455 9.4
    KSRH10(Ⅱ) 1 0.1 1.8353 20 1.9810 1.99
    IBRH17(Ⅱ) 1 0.01 0.3534 20 0.3504 0.39
    2 0.01 0.9096 20 0.9047 1.06
    3 0.01 1.3833 20 1.4935 1.5
    IWTH24(Ⅱ) 1 0.1 2.5140 20 2.7376 3.37
    2 0.1 6.3388 40 8.1129 9.23
    3 0.1 10.8505 20 15.3724 15.01
    KSEH06(Ⅱ) 1 0.1 5.7812 20 6.3617 6.65
    NIGH11(Ⅱ) 2 0.1 4.5118 20 4.6040 3.9
    3 0.1 6.8949 20 8.4756 6.09
    4 0.1 10.01 20 14.5878 9.36
    TCGH12(Ⅱ) 1 0.1 5.2806 20 6.6765 7.22
    2 0.05 7.2384 20 9.3454 9.25
    TKCH08(Ⅱ) 3 0.1 6.7590 20 8.1318 8.35
    FKSH11 (Ⅲ) 1 0.1 1.5472 60 1.6730 2.12
    2 0.1 4.5552 70 5.2341 5.86
    3 0.01 6.8218 50 8.0801 9.61
    KSRH07(Ⅲ) 1 0.01 2.4128 60 2.5635 2.51
    2 0.01 6.9196 20 8.8538 7.02
    FKSH14(Ⅲ) 1 0.1 1.1760 40 1.1875 1.35
    2 0.1 3.6658 20 3.9909 4.01
    3 0.1 5.2186 40 5.9271 6.4
    4 0.1 6.3233 40 7.7832 8.83
    NMRH04(Ⅳ) 3 0.1 0.4053 20 0.4089 0.61
    4 0.1 1.7428 30 1.9726 1.61
    5 0.1 2.3917 20 4.3434 2.04
    IBRH10(Ⅳ) 1 0.01 0.7869 40 0.7931 0.78
    2 0.01 1.1389 20 1.1897 1.27
    下载: 导出CSV

    表  6  场地卓越频率或其范围统计

    Table  6.   Frequency of excellence or its range statistics

    场地类别 显著峰数量/个 卓越频率值/Hz 卓越频率值范围/Hz 选取方法
    IWTH02(Ⅱ) 1 5.05 SA
    KSRH10(Ⅱ) 1 1.98 FAS
    IWTH27(Ⅱ) 1 7.74 FAS
    IBRH13(Ⅱ) 1 3.06 FAS
    KSRH06(Ⅱ) 1 6.36 FAS
    KSRH07(Ⅲ) 2 2.560~8.850 FAS
    FKSH14(Ⅲ) 2 1.180~7.780 FAS
    IBRH10(Ⅳ) 2 0.790~1.180 FAS
    TCGH12(Ⅱ) 2 0.107~0.381 FAS
    IWTH08(Ⅱ) 3 2.730~8.070 FAS
    TKCH08(Ⅱ) 3 0.123~0.551 FAS
    IWTH24(Ⅱ) 3 0.084~0.593 FAS
    FKSH11 (Ⅲ) 3 1.670~8.00 FAS
    NMRH04(Ⅳ) 5 0.263~2.467 SA
    NIGH11(Ⅱ) 6 0.097~0.355 SA
    IBRH17(Ⅱ) 6 0.163~0.297 FAS
    下载: 导出CSV
  • 陈鹏,刘文锋,付兴潘,2009. 关于场地卓越周期和特征周期的若干讨论. 青岛理工大学学报,30(6):30−35. doi: 10.3969/j.issn.1673-4602.2009.06.006

    Chen P., Liu W. F., Fu X. P., 2009. Discussions on site predominant period and characteristic period. Journal of Qingdao Technological University, 30(6): 30−35. (in Chinese) doi: 10.3969/j.issn.1673-4602.2009.06.006
    陈永新,迟明杰,李小军,2016. 基于强震动记录确定的场地卓越周期. 地震学报,38(1):138−145. doi: 10.11939/jass.2016.01.014

    Chen Y. X., Chi M. J., Li X. J., 2016. Determination of site dominant period based on strong motion records. Acta Seismologica Sinica, 38(1): 138−145. (in Chinese) doi: 10.11939/jass.2016.01.014
    丁毅,王玉石,王宁等,2021. 地表/井下反应谱比值非线性统计特征与影响因素研究. 震灾防御技术,16(2):362−370. doi: 10.11899/zzfy20210215

    Ding Y., Wang Y. S., Wang N., et al., 2021. Study on nonlinear statistical characteristics of surface/downhole response spectrum ratio and influencing factors. Technology for Earthquake Disaster Prevention, 16(2): 362−370. (in Chinese) doi: 10.11899/zzfy20210215
    高广运,吴世明,周健等,2000. 场地卓越周期的讨论与测定. 工程勘察,(5):29−31.

    Gao G. Y., Wu S. M., Zhou J., et al., 2000. Discussion and measurement of site predominant period. Geotechnical Investigation & Surveying, (5): 29−31. (in Chinese)
    李文倩,何金刚,朱皓清,2019. 基于 H/V 谱比法的场地卓越频率研究. 内陆地震,33(4):314−320.

    Li W. Q., He J. G., Zhu H. Q., 2019. Study on site predominant frequency based on H/V spectral ratio method. Inland Earthquake, 33(4): 314−320. (in Chinese)
    李小军,1992. 场地土层对地震地面运动影响的分析方法. 世界地震工程,(2):49−60.
    刘宇实,师黎静,2018. 基于地脉动谱比法的场地特征参数快速测定. 振动与冲击,37(13):235−242.

    Liu Y. S., Shi L. J., 2018. Site characteristic parameters’ quick measurement based on micro-tremor’s H/V spectra. Journal of Vibration and Shock, 37(13): 235−242. (in Chinese)
    罗桂纯,李小军,王玉石等,2014. 关于中村(Nakamura)方法分析结构响应有效性的讨论. 地震学报,36(3):491−499. doi: 10.3969/j.issn.0253-3782.2014.03.014

    Luo G. C., Li X. J., Wang Y. S., et al., 2014. Discussion on validity of structural response by Nakamura’s technique. Acta Seismologica Sinica, 36(3): 491−499. (in Chinese) doi: 10.3969/j.issn.0253-3782.2014.03.014
    罗桂纯,李小军,傅磊等,2019. 基于 HVSR 谱比法研究场地反应非线性特征. 地震研究,42(4):546−554. doi: 10.3969/j.issn.1000-0666.2019.04.013

    Luo G. C., Li X. J., Fu L., et al., 2019. Study on nonlinearity of site effect with the HVSR spectral ratio. Journal of Seismological Research, 42(4): 546−554. (in Chinese) doi: 10.3969/j.issn.1000-0666.2019.04.013
    任叶飞,温瑞智,山中浩明等,2013. 运用广义反演法研究汶川地震场地效应. 土木工程学,46(S2):146−151.

    Ren Y. F., Wen R. Z., Hiroaki Y., et al., 2013. Research on site effect of Wenchuan earthquake by using generalized inversion technique. China Civil Engineering Journal, 46(S2): 146−151. (in Chinese)
    荣棉水,李小军,王振明等,2016. HVSR方法用于地震作用下场地效应分析的适用性研究. 地球物理学报,59(8):2878−2891. doi: 10.6038/cjg20160814

    Rong M. S., Li X. J., Wang Z. M., et al., 2016. Applicability of HVSR in analysis of site-effects caused by earthquakes. Chinese Journal of Geophysics, 59(8): 2878−2891. (in Chinese) doi: 10.6038/cjg20160814
    童广才,刘康和,2000. 场地卓越周期的确定. 电力勘测,(2):43−46.

    Tong G. C., Liu K. H., 2000. Determine predominant period of site. Electric Power Survey & Design, (2): 43−46. (in Chinese)
    姚鑫鑫,2017. 利用强震动记录的谱比分析场地非线性反应. 北京:中国地震局工程力学研究所.

    Yao X. X., 2017. Identification of the nonlinear site response using strong motion records based on spectral ratio method. Beijing:Institute of Engineering Mechanics,China Earthquake Administration. (in Chinese)
    姚鑫鑫,任叶飞,温瑞智等,2019. 强震动记录 H/V 谱比法计算处理的若干关键环节. 震灾防御技术,14(4):719−730. doi: 10.11899/zzfy20190403

    Yao X. X., Ren Y. F., Wen R. Z., et al., 2019. Some technical notes on the data processing of the spectral ratio based on the strong-motion records. Technology for Earthquake Disaster Prevention, 14(4): 719−730. (in Chinese) doi: 10.11899/zzfy20190403
    袁一凡,田启文,2012. 关于反应谱比与傅里叶谱比曲线的比较. 见:工程地震学. 北京:地震出版社.
    张立宝,2018. 基于钻井台阵数据的场地非线性效应研究. 北京:中国地震局地球物理研究所.

    Zhang L. B., 2018. Study on nonlinearity of site effects on ground motion based on vertical borehole array data. Beijing:Institute of Geophysics,China Earthquake Administration. (in Chinese)
    Chen J., Tang H., Chen W. K., et al., 2022a. A prediction method of ground motion for regions without available observation data (LGB-FS) and its application to both Yangbi and Maduo earthquakes in 2021. Journal of Earth Science, 33(4): 869−884. doi: 10.1007/s12583-021-1560-6
    Chen W. K. , Wang D. , Zhang C. , et al. , 2022b. Estimating seismic intensity maps of the 2021 MW 7.3 Madoi, Qinghai and MW 6.1 Yangbi, Yunnan, China Earthquakes. Journal of Earth Science, 33 (4): 839−846.
    Hassani B., Atkinson G. M., 2016. Applicability of the site fundamental frequency as a V S30 proxy for central and eastern north America. Bulletin of the Seismological Society of America, 106(2): 653−664. doi: 10.1785/0120150259
    Lermo J., Chávez-García F. J., 1993. Site effect evaluation using spectral ratios with only one station. Bulletin of the Seismological Society of America, 83(5): 1574−1594. doi: 10.1785/BSSA0830051574
    Nakamura Y. , 1989. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Railway Technical Research Institute, Quarterly Reports, 30 (1).
    Theodulidis N., Bard P. Y., Archuleta R., et al., 1996. Horizontal-to-vertical spectral ratio and geological conditions: the case of Garner Valley downhole array in southern California. Bulletin of the Seismological Society of America, 86(2): 306−319. doi: 10.1785/BSSA0860020306
    Thompson E. M., Baise L. G., Tanaka Y., et al., 2012. A taxonomy of site response complexity. Soil Dynamics and Earthquake Engineering, 41: 32−43. doi: 10.1016/j.soildyn.2012.04.005
    Wen K. L., Beresnev I. A., Yeh Y. T., 1994. Nonlinear soil amplification inferred from downhole strong seismic motion data. Geophysical Research Letters, 21(24): 2625−2628. doi: 10.1029/94GL02407
    Yamazaki F., Ansary M. A., 1997. Horizontal-to-vertical spectrum ratio of earthquake ground motion for site characterization. Earthquake Engineering & Structural Dynamics, 26(7): 671−689.
    Zhao J. X., Irikura K., Zhang J., et al., 2006. An empirical site-classification method for strong-motion stations in Japan using H/V response spectral ratio. Bulletin of the Seismological Society of America, 96(3): 914−925. doi: 10.1785/0120050124
    Zhu C. B., Weatherill G., Cotton F., et al., 2021. An open-source site database of strong-motion stations in Japan: K-NET and KiK-net (v1.0. 0). Earthquake Spectra, 37(3): 2126−2149. doi: 10.1177/8755293020988028
  • 加载中
图(10) / 表(6)
计量
  • 文章访问数:  142
  • HTML全文浏览量:  41
  • PDF下载量:  48
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-05-17
  • 网络出版日期:  2024-10-15
  • 刊出日期:  2024-09-01

目录

/

返回文章
返回