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

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

李小军 钱玉 荣棉水 孔小山

白建方, 马立龙. Rayleigh波场的数值模拟及其应用[J]. 震灾防御技术, 2019, 14(2): 328-340. doi: 10.11899/zzfy20190207
引用本文: 李小军,钱玉,荣棉水,孔小山,2024. 基于地震动反应谱与傅里叶谱确定场地卓越频率的差异性. 震灾防御技术,19(3):421−435. doi:10.11899/zzfy20240301. doi: 10.11899/zzfy20240301
Bai Jianfang, Ma Lilong. Numerical Modeling Techniques of Rayleigh Wave Field and Its Application[J]. Technology for Earthquake Disaster Prevention, 2019, 14(2): 328-340. doi: 10.11899/zzfy20190207
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之间获取的场地卓越频率较为准确。然后对比分析了利用地震动加速度反应谱比和傅里叶谱比得到的场地卓越频率与场地土层模型计算得到的基于传递函数的自振频率。研究结果表明,对大多数台站而言,采用傅里叶谱比计算场地卓越频率具有明显的优势,对于Ⅱ、Ⅲ、Ⅳ类场地上的台站均有如此结论,只有对少数特定台站,采用反应谱比方法效果更好。
  • 关于场地地震反应的分析已有大量研究成果,研究表明土壤在地震作用下会表现出材料非线性效应ADDIN EN.CITE.DATA(Joyner等,1975Huang等,2001Arslan等,2006Hosseini等,2012)。等效线性化方法ADDIN EN.CITE.DATA(Schnabel等,1972Idriss等,1992Bardet等,2000王笃国等,2016)是一种频域方法,通过在不同土体应变条件下选择等效阻尼比和剪切模量,将非线性问题转化为线性问题。当采用材料非线性本构模型描述土体非线性时,需采用时间积分算法求解非线性动力有限元方程。时间积分算法可分为隐式方法和显式方法。隐式算法每时刻需求解线性代数方程组,计算效率相对较低,如Wilson-θ法和Newmark法等。显式算法无需求解线性代数方程组,适合于强非线性和自由度数目较大的问题。研究者已提出多种显式时间积分算法ADDIN EN.CITE.DATA(Chung等,1994王进廷等,2002Belytschko等,2014)。作者近期提出一种二阶精度的单步显式算法,该算法适合变时步问题,在线弹性范围内稳定性较好。本文将该算法推广至求解非线性动力有限元方程中,并将其应用于地震波垂直入射时非线性地震反应分析。

    设已知非线性体系第${t_i}$时步的受力状态,求解第${t_{i + 1}}$时步的非线性结构动力学方程:

    $${\boldsymbol{M}}{{\boldsymbol{\ddot u}}_{i + 1}}{\boldsymbol{ + C}}{{\boldsymbol{\dot u}}_{i + 1}} + {\boldsymbol{f}}_{i + 1}^S{\boldsymbol{ = }}{{\boldsymbol{f}}_{i + 1}}$$ (1)

    式中MC、${{\boldsymbol{f}}^S}$和${\boldsymbol{f}}$分别表示非线性体系的质量矩阵、阻尼矩阵、内力向量和外荷载向量;u表示位移,点号对时间t求导,i+1表示第${t_{i + 1}}$时刻。第i+1时刻时间步长为:

    $${\boldsymbol{\Delta }}{t_i} = {t_{i + 1}} - {t_i}$$ (2)

    文献显式方法求解非线性方程(1)的过程如下,第i+1时刻位移${{\boldsymbol{u}}_{i + 1}}$为:

    $${{\boldsymbol{u}}_{i + 1}} = {{\boldsymbol{u}}_i} + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}{{\boldsymbol{\dot u}}_i} + \frac{{\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}^2}}{2}{{\boldsymbol{\ddot u}}_i}$$ (3)

    i+1时刻位移增量$\mathit{\Delta }{{\boldsymbol{u}}_i}$、内力增量$\mathit{\Delta }{\boldsymbol{f}}_i^S$和内力全量${\boldsymbol{f}}_{i + 1}^S$分别为:

    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{{\boldsymbol{u}}_i} = {{\boldsymbol{u}}_{i + 1}} - {{\boldsymbol{u}}_i}$$ (4)
    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{f}}_i^S = {\boldsymbol{f}}(\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{{\boldsymbol{u}}_i})$$ (5)
    $${\boldsymbol{f}}_{i + 1}^S = {\boldsymbol{f}}_i^S + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{f}}_i^S$$ (6)

    i+1时刻预估速度${{\boldsymbol{\dot {\tilde u}}}_{i + 1}}$、预估加速度${{\boldsymbol{\ddot {\tilde u}}}_{i + 1}}$、速度${{\boldsymbol{\dot u}}_{i + 1}}$和加速度${{\boldsymbol{\ddot u}}_{i + 1}}$分别为

    $${{\boldsymbol{\dot {\tilde u}}}_{i + 1}} = {{\boldsymbol{\dot u}}_i} + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}{{\boldsymbol{\ddot u}}_i}$$ (7)
    $${{\boldsymbol{\ddot {\tilde u}}}_{i + 1}} = {{\boldsymbol{M}}^{ - 1}}({{\boldsymbol{f}}_{i + 1}} - {\boldsymbol{C\dot {\tilde u}}}_{i + 1}^{} - {\boldsymbol{f}}_{i + 1}^S)$$ (8)
    $${{\boldsymbol{\dot u}}_{i + 1}} = {{\boldsymbol{\dot u}}_i} + \frac{{\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}}}{2}({{\boldsymbol{\ddot u}}_i} + {{\boldsymbol{\ddot {\tilde u}}}_{i + 1}})$$ (9)
    $${{\boldsymbol{\ddot u}}_{i + 1}} = {{\boldsymbol{M}}^{ - 1}}({{\boldsymbol{f}}_{i + 1}} - {\boldsymbol{C\dot u}}_{i + 1}^{} - {\boldsymbol{f}}_{i + 1}^S)$$ (10)

    式(3)—式(10)为求解式(1)的显式算法。算法中需由位移增量计算内力增量,目前常用的应力计算方法包括向前欧拉法、向后欧拉法和完全隐式计算法等ADDIN EN.CITE.DATA(Sloan等,19922001Ahadi等,2003)。下面给出式(5)由位移增量计算内力增量的过程,即一种带误差控制的修正欧拉算法。

    对于每个有限单元,由位移增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{u}}_i^e$计算应变增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e$的表达式为:

    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e = {{\boldsymbol{B}}^e}\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{u}}_i^e$$ (11)

    式中Be为应变矩阵。将ti时刻单元应变增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e$赋值给子步应变增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$,ti时刻单元应力${\boldsymbol{ \pmb{\mathit{ σ}} }}_i^e$赋值给${\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e$,初始化子步应变增量和应力状态分别为:

    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e \leftarrow \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e$$ (12)
    $${\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e \leftarrow {\boldsymbol{ \pmb{\mathit{ σ}} }}_i^e$$ (13)

    每个子步中应力增量计算思路见图 1,具体计算公式如下:

    $${\boldsymbol{D}}_1^e = {\boldsymbol{D}}({\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e)$$ (14)
    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e = {\boldsymbol{D}}_1^e\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (15)
    $${\boldsymbol{D}}_2^e = {\boldsymbol{D}}({\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e)$$ (16)
    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_2^e = {\boldsymbol{D}}_2^e\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (17)
    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_s^e = \frac{{\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_2^e}}{2}$$ (18)
    图 1  修正欧拉算法计算应力增量
    Figure 1.  Modified Euler algorithm to calculate stress increment

    式中${{\boldsymbol{D}}^e}$为单元应力-应变关系矩阵。判断每个子步中应力增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{{\boldsymbol{ \pmb{\mathit{ σ}} }}_s}$是否符合精度要求的误差判断式为:

    $${e_r} = \frac{{\left\| {\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e - \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_2^e} \right\|}}{{\left\| {{\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_s^e} \right\|}}$$ (19)

    判断误差er是否小于预先给定的判断值st,条件不满足时,缩小子步应变增量为:

    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e \leftarrow A\sqrt {{{{s_t}} \mathord{\left/ {\vphantom {{{s_t}} {{e_r}}}} \right. } {{e_r}}}} \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (20)

    式中A为误差峰值系数。采用缩小的子步应变增量重新进行式(14)—式(19)的计算与判断,循环直至满足精度要求,更新剩余应变增量和应力状态分别为:

    $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e \leftarrow \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e - \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (21)
    $${\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e \leftarrow {\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_s^e$$ (22)

    利用更新剩余应变增量和应力状态循环执行式(14)—式(20),直至剩余应变增量小于等于零结束。

    利用求得的第i+1时刻单元应力可得到单元应力增量和内力增量分别为:

    $$ \Delta \boldsymbol{\sigma }_i^e = \boldsymbol{\sigma }_{i + 1}^e - \boldsymbol{\sigma }_i^e $$ (23)
    $$ \Delta {\boldsymbol{f}}_i^S{\rm{ = }}\sum\limits_e {\int {{{\boldsymbol{B}}^{e{\rm{T}}}}\boldsymbol{\Delta }{\boldsymbol{\sigma }}_i^e{\bf{d}}A} } $$ (24)

    本节将上述非线性有限元方程的显式时间积分算法应用于地震波垂直入射时场地非线性地震反应分析中。假定基岩为线弹性半空间,考虑基岩上覆土层的材料非线性,不考虑土体阻尼。在土层下部设置黏性边界条件模拟半空间基岩的辐射阻尼,并在该处以等效结点力的方式实现地震动输入。

    计算模型见图 2,选取A点作为观测点。土体非线性材料本构模型选取邓肯-张模型,土体线弹性参数见表 1,未给出配套的非线性参数,故算例中的非线性参数参考实际情况选取,后续研究中将使用更真实表现土体非线性行为的本构模型及真实工程场地参数。算例中的大气压参数取100kPa,内摩擦角增量取0°。入射地震动分别选取狄拉克脉冲和实测地震动(Gilroy Array #3,Coyote Lake, 1979)。入射狄拉克脉冲见图 3,观测点结果见图 4,实测地震动见图 5,观测点结果见图 6图 4图 6中给出采用中心差分法的计算结果作为参考解,由图 4图 6可知,本文算法与中心差分法计算结果吻合较好,说明本文算法的有效性。

    图 2  大开车站沿线土层纵断面构造
    Figure 2.  Site condition of the Daikai subway station in vertical direction
    表 1  土层参数
    Table 1.  Parameters of soils
    土质 深度/
    m
    $\rho $/
    (g/cm3
    cs/
    (m/s)
    v
    -
    EN
    -
    Rf
    -
    c/
    (MPa)
    θ/(°) D
    -
    F
    -
    人工填土 0—1.0 1.9 140 0.33 0.33 0.758 0.084 26.9 1.06 0.021
    全新世砂土 1.0—5.1 1.9 140 0.32 0.33 0.758 0.084 26.9 1.06 0.021
    全新世砂土 5.1—8.3 1.9 170 0.32 0.36 0.768 0.120 31.0 1.11 0.015
    更新世粘土 8.3—11.4 1.9 190 0.40 0.44 0.822 0.188 28.4 1.01 0.012
    更新世粘土 11.4—17.2 1.9 240 0.30 0.44 0.822 0.188 28.4 1.01 0.012
    更新世砂土 17.2—22.2 2.0 330 0.26 0.51 0.840 0.300 30.0 1.02 0.011
    基岩 >22.2 2.0 330 0.26 - - - - - -
    下载: 导出CSV 
    | 显示表格
    图 3  狄拉克脉冲速度和加速度时程图
    Figure 3.  Velocity and acceleration time history of the Dirac pulse
    图 4  狄拉克脉冲入射时场地反应分析结果
    Figure 4.  Results of site analysis under the incident of Dirac pulse
    图 5  实测地震动速度和加速度时程图
    Figure 5.  Velocity and acceleration time history of the seismic motion
    图 6  实测地震动入射时场地反应分析结果
    Figure 6.  Results of site reaction analysis under the incident of the seismic motion

    表 1ρcsvENRfcθ为模型参数,分别表示密度、剪切波速、泊松比、无量纲幂次、破坏比、土的内聚力、土的摩擦角。DF为试验常数。

    本文发展一种求解材料非线性结构动力学方程的显式时间积分算法,并应用于地震波竖直入射时非线性地震反应分析中,通过算例验证了该方法的有效性。该显式算法具有无需对角阻尼矩阵、单步、稳定性良好等优点。本文考虑了邓肯-张非线性弹性本构模型,下步研究可考虑将该显式算法扩展到弹塑性本构模型及更能反映土层真实变形的本构模型中。

  • 图  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
  • 期刊类型引用(2)

    1. 沈银澜,周敬轩,王利辉,刘辉,周海宾,吴震东,宫逸飞. 北方村落木结构民居榫卯节点抗震性能试验研究. 震灾防御技术. 2021(01): 165-175 . 本站查看
    2. 郭光玲,郭瑞,徐乾. 略阳县农居房屋抗震现状调查与评价研究. 震灾防御技术. 2021(03): 492-500 . 本站查看

    其他类型引用(4)

  • 加载中
图(10) / 表(6)
计量
  • 文章访问数:  145
  • HTML全文浏览量:  41
  • PDF下载量:  49
  • 被引次数: 6
出版历程
  • 收稿日期:  2023-05-17
  • 网络出版日期:  2024-10-15
  • 刊出日期:  2024-09-01

目录

/

返回文章
返回