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之间获取的场地卓越频率较为准确。然后对比分析了利用地震动加速度反应谱比和傅里叶谱比得到的场地卓越频率与场地土层模型计算得到的基于传递函数的自振频率。研究结果表明，对大多数台站而言，采用傅里叶谱比计算场地卓越频率具有明显的优势，对于Ⅱ、Ⅲ、Ⅳ类场地上的台站均有如此结论，只有对少数特定台站，采用反应谱比方法效果更好。Abstract: The Horizontal-to-Vertical Spectral Ratio (HVSR) method for strong motion records is commonly used to evaluate a site’s predominant frequency. Both the acceleration response spectrum and the acceleration Fourier spectrum are frequently employed to calculate the HVSR of strong motion records, but these two approaches often yield different results. To investigate the differences in predominant site frequency between the response spectrum ratio and the Fourier spectrum ratio, this study analyzes 16 stations from Japan’s KiK-net network, where site conditions can be approximated using a one-dimensional site model. First, data processing methods for evaluating site predominant frequency were developed. These include S-wave interception, taper preprocessing, and automatic peak identification based on Gaussian fitting. The effects of damping ratio and smoothed bandwidth ratio on site frequency evaluation were examined. It was found that using a 10% damping ratio for the response spectrum and a smoothing bandwidth ratio of 20~40 for the Fourier spectrum yielded more accurate results. Next, the predominant site frequencies obtained from the ground motion acceleration response spectrum and Fourier spectrum were compared with the natural frequency calculated from the site soil model using the transfer function. The results indicate that the Fourier spectrum ratio provides a more accurate estimate of the predominant site frequency compared to the response spectrum ratio for most stations, particularly those located on Class II, III, and IV sites. However, a few specific stations showed better results when using the response spectrum method.
-
引言
国内外地震灾害分析表明，隔震技术在提高结构抗震性能、经济效益和社会效益等方面具有显著优势，近几十年隔震技术在建筑结构与桥梁工程领域得到一定应用。但《建筑抗震设计规范》(GB 50011—2010)(中华人民共和国住房和城乡建设部等，2016)中第12.2.9条第2款规定隔震层应满足嵌固刚度比要求，此要求使隔震支墩长细比较小，须单独增设1个隔震层，增加了结构造价，在一定程度上限制了隔震技术的应用。
目前对隔震支墩长细比的研究较少，主要集中于对底部二层框架上部多塔楼底隔震结构的数值模拟与试验研究(徐忠根等，2005)、对首层薄弱层框架结构柱顶隔震性能的分析(吴应雄等，2011)、对层间隔震减震结构的理论分析和振动台试验研究(黄襄云，2008)、对近断层脉冲型地震动作用下高层建筑组合隔震结构的减震性能研究(潘钦锋等，2019)等，上述研究中模型均设置了单独的隔震层，未涉及隔震支墩长细比的研究。另外，对屋盖结构柱顶隔震技术的研究仅针对大跨度结构，得出柱顶隔震技术可应用于屋盖结构的结论(唐柏鉴等，2005)，对高举架立式圆筒型储液容器柱顶隔震地震动响应的研究仅针对构筑物(孙建刚等，2018)。本文使用大型通用有限元分析软件ANSYS分别对无隔震支座、常规隔震、柱顶隔震、柱底隔震4种方案的同一建筑进行地震时程反应分析，研究隔震支墩长细比对隔震效果的影响，进而证实柱顶隔震、柱底隔震方案的可行性，从而降低隔震建筑造价，为医院、学校等多层建筑隔震设计提供经济可行性，同时也为《建筑抗震设计规范》(GB 50011—2010)的修订提供一定理论依据。
1. 工程概况
本建筑为混凝土框架结构民用建筑，共8层，层高3m，纵向5跨，跨度均为6m，横向3跨，跨度分别为6.3m、2.7m、6.3m。本建筑处于8度设防区，设计地震加速度为0.20g，地震分组为第二组，场地类别为Ⅱ类，根据《建筑抗震设计规范》(GB 50011—2010)，得该场地特征周期为0.4s，属于重点设防类、乙类建筑，结构平面图如图 1所示。
2. 数值模拟分析
2.1 有限元模型
本文使用大型通用有限元分析软件ANSYS进行地震时程反应分析。本工程结构为钢筋混凝土框架结构，通过实体单元Solid65模拟。为更好地模拟地震动，用3个相互垂直的弹簧单元Combin14分别模拟普通橡胶隔震支座x、y、z方向的分量，用2个相互垂直的弹簧单元Combin40分别模拟铅芯橡胶隔震支座x、y方向的分量，用Combin14单元模拟铅芯橡胶隔震支座z方向的分量(王新敏等，2011)。建筑物有限元模型如图 2，普通隔震支座有限元模拟示意图如图 3。
图 3 普通隔震支座有限元模拟示意图Figure 3. The finite element simulation diagram of ordinary isolation bearing无隔震支座方案为常规建筑底部无隔震支座；常规隔震方案为在底部单独增设1个隔震层，在隔震层底部安装隔震支座；柱顶隔震方案为不在首层、底部单独增设隔震层，在首层柱柱顶安装隔震支座；柱底隔震方案为不在首层、底部单独增设隔震层，在首层柱柱底安装隔震支座。
2.2 PKPM柱底压力计算结果
本文利用PKPM软件计算无隔震支座方案建筑在重力荷载代表值作用下的首层柱柱底压力，结果见表 1。
表 1 首层柱柱底压力计算结果(kN)Table 1. The bottom pressure result of the first layer's coloum on PKPM (kN)横轴 纵轴 Ⓐ Ⓑ Ⓒ Ⓓ Ⓔ Ⓕ ① 1524.6 2065.7 2083.7 2083.7 2065.7 1524.6 ② 1768.7 2301.2 2322.3 2322.3 2301.2 1768.7 ③ 1768.7 2301.2 2322.3 2322.3 2301.2 1768.7 ④ 1524.6 2065.7 2083.7 2083.7 2065.7 1524.6 2.3 有限元模型参数
柱顶隔震方案中，首层和2层柱截面尺寸为800mm×800mm，柱底隔震方案中只有底层柱截面尺寸为800mm×800mm，其他柱截面尺寸均为600mm×600mm，梁截面尺寸为300mm×500mm，板厚100mm。柱顶隔震方案中，首层和2层柱Solid65单元弹性模量取3.80×104N/mm2，柱底隔震方案中只有底层柱Solid65单元弹性模量取3.80×104N/mm2，其他柱Solid65单元弹性模量取均3.45×104N/mm2，泊松比取0.2，密度取5300kg/m3。根据重力荷载代表值作用下的柱底压力计算结果(表 1)，按隔震橡胶支座最大竖向受力及外部支座扭转位移较大的原则，使用3类隔震橡胶支座，支座1为LNR(普通橡胶支座，下同)600，支座2为LRB(铅芯橡胶支座，下同)500，支座3为LRB600，支座力学参数和结构尺寸见表 2和表 3，其中，Ⓐ、Ⓕ轴选用LRB500支座，①、④轴(除与Ⓐ、Ⓕ轴相交的支座外)选用LRB600支座，其余支座选用LNR600支座(图 4)，无隔震支座、常规隔震、柱顶隔震、柱底隔震4种方案单元数分别为49328、54144、49400、49400，模型底部采用固定约束。
表 2 隔震支座力学参数Table 2. The mechanics parameter of isolation bearing序号 类型 设计荷载/kN 竖向刚度/kN·mm-1 等效水平刚度/kN·mm-1 屈服前刚度/kN·mm-1 屈服后刚度/kN·mm-1 屈服力/kN 等效阻尼比/% 1 LNR600 4230 2803 0.917 2 LRB500 2945 2451 1.186 10.139 12.121 40 20.9 3 LRB600 4241 2917 1.459 0.780 0.932 63 21.9 表 3 隔震支座结构尺寸Table 3. The structure size of isolation bearing序号 类型 支座外径/mm 橡胶保护层厚度/mm 内部橡胶厚度/mm 有效直径/mm 铅芯直径/mm 不含连板高度/mm 第一形状系数 第二形状系数 连接板外形尺寸/mm 连接板厚度/mm 螺栓直径/mm 螺栓个数 支座高度/mm 1 LNR600 620 10 118.9 600 241.3 34.8 5.0 700 25.0 M36 4 291.3 2 LRB500 520 10 98.6 500 80 213.0 36.8 5.1 600 22.0 M30 4 257.0 3 LRB600 620 10 118.9 600 100 241.3 36.6 5.0 700 25.0 M36 4 291.3 3. 有限元分析
3.1 地震动时程反应分析
用时程分析法分析结构地震响应时，除给出必要的结构参数外，还应确定相应的输入地震动时程。地震动时程影响因素较多，在相同烈度下，同一场地类别的观测点观测到的地震加速度记录在峰值、波形、频谱和持续时间上不同，即使是同一震源先后两次相同震级的地震，同一观测点观测到的地震加速度记录也不同(杨溥等，2000；王丽娟，2013)。但只要合理选择地震动主要参数(幅值、频谱、持时)，时程分析结果较可靠。因此，选择合适的地震动参数，并调整计算结果尤为重要(兰雁，2012；李建亮等，2011)。
《建筑抗震设计规范》(GB 50011—2010)中第5.1.2条规定，采用时程分析法时，应按建筑场地类别和设计地震分组选用实际强震记录和人工模拟的加速度时程曲线，其中实际强震记录数量不应少于总数的2/3。当取3组加速度时程曲线输入时，计算结果宜取时程法包络值和振型分解反应谱法的较大值。《建筑抗震设计规范》(GB 50011—2010)第5.1.2条规定，当结构采用三维空间模型需要双向(2个水平方向)或三向(2个水平方向和1个竖向)地震动输入时，其加速度最大值通常按1(水平方向1):0.85(水平方向2):0.65(竖向)的比例调整。
3.2 地震动时程反应分析结果
本工程分别采用罕遇地震水平和设防地震水平3条不同地震动时程进行时程反应分析，研究建筑遭遇罕遇地震和设防地震时的反应情况，输入的地震动时程分别为由设计反应谱得到的人工地震动时程、相同场地类别并经过调幅的El-Centro地震动时程及唐山地震动时程。
调幅后罕遇地震水平SPECTRUM地震动、唐山地震动、El-Centro地震动时程曲线及3条地震动时程反应谱曲线、平均反应谱曲线和规范反应谱曲线的对比如图 5。《建筑抗震设计规范》(GB 50011—2010)中第5.1.2条规定，输入的地震加速度时程曲线有效持续时间一般从首次达到该时程曲线最大峰值的10%那一刻算起，到最后一刻达到最大峰值的10%为止；无论是实际的强震记录还是人工模拟波形，有效持续时间一般为结构基本周期的5—10倍。由图 5(a)、5(b)、5(c)可知，3条地震动时程有效持续时间均达到建筑自身结构基本周期(2.32s)的5—10倍，满足规定。由图 5(d)可知，各地震动时程平均反应谱曲线与规范反应谱曲线较接近，满足要求。由此可知，这3条地震动时程可用于本工程地震动时程输入。
图 5 水平罕遇地震作用下3条地震动时程及反应谱对比Figure 5. Three seismic wave time history of horizontal rare occurrence earthquake level and comparison of response spectrum输入调幅后的罕遇地震水平唐山地震动时程后无隔震支座、常规隔震、柱顶隔震、柱底隔震4种方案建筑物纵向第3跨最高处中部纵向位移时程曲线如图 6。由图 6可知，遭遇地震时，无隔震支座方案建筑物摆动频率较大，使用隔震支座后，隔震支座使整个建筑物刚度变小，摆动频率降低，周期变大，减小了地震对结构、非结构构件、内部附属物品的损坏；隔震建筑物顶部位移小于无隔震支座建筑物，有利于减小地震对结构、非结构构件、内部附属物品的损坏；常规隔震方案与柱顶隔震方案、柱底隔震方案建筑物顶部位移反应基本相同，说明对于建筑物顶部位移反应而言，隔震支墩长细比基本能达到普通层柱长细比，可根据建筑物自身需求决定采用柱顶隔震方案或柱底隔震方案。
图 6 输入调幅唐山地震动时程后建筑物纵向第3跨最高处中部纵向位移反应Figure 6. The longitudinal displacement response on the middle part of the construction's third longitudinal span after inputting amplitude adjusted tangshan seismic wave time history输入调幅后的罕遇地震水平El-Centro地震动时程后无隔震支座、常规隔震、柱顶隔震、柱底隔震4种方案建筑物纵向第3跨最高处中部最大纵向位移时建筑的纵向位移云图如图 7，由图 7可知，遭遇地震时，无隔震支座方案建筑物运动方式为左右摆动，使用隔震支座后，支座上部建筑物运动方式变为平动，大部分位移由隔震支座承担，建筑自身层间相对位移很小，隔震支座消耗了大部分地震能量，可减小地震对结构、非结构构件、内部附属物品的损坏；常规隔震、柱顶隔震、柱底隔震3种方案建筑物位移反应相差较小，其中柱顶隔震方案建筑物位移相对较小，仅柱顶隔震方案建筑物首层柱位移较大，不利于装修和非结构构件的安装，加大了地震发生时首层顶部非结构构件和附属构件脱落致人伤亡的危险，综合比较得知，隔震支墩长细比基本能达到普通层柱长细比，出于安全考虑，应优选柱底隔震方案。
图 7 输入调幅El-Centro地震动时程后建筑物最大纵向位移云图(mm)Figure 7. The max longitudinal displacement cloud diagram of the construct after inputting amplitude adjusted El-Centro seismic wave time history (mm)输入罕遇地震和设防地震水平地震动时程后，柱顶隔震方案和柱底隔震方案与常规隔震方案建筑物隔震支座纵向和横向最大位移比见表 4，由表 4可知，对于铅芯橡胶支座，遭遇罕遇地震和设防地震水平地震动时，柱顶隔震方案和柱底隔震方案建筑物隔震支座纵向和横向位移均小于常规隔震方案，比值为0.73—0.99；对于普通橡胶支座，遭遇罕遇地震和设防地震水平地震动时，柱顶隔震方案和柱底隔震方案建筑物隔震支座纵向和横向位移均大于常规隔震方案，比值为1.15—1.43。可考虑将普通橡胶支座全部换成铅芯橡胶支座或更换大一号的普通橡胶支座避免遭遇地震时可能出现的破坏。综上所述，隔震支墩长细比基本能达到普通层柱长细比。
表 4 柱顶隔震方案和柱底隔震方案与常规隔震方案建筑物隔震支座位移比Table 4. The isolation bearing's displacement ratio of the column's top and bottom isolation project with the conventional isolation project横轴 工况 支座位置 方向 纵轴 Ⓐ Ⓑ Ⓒ Ⓓ Ⓔ Ⓕ ① 罕遇地震 柱顶 纵向 0.91 0.87 0.87 0.87 0.87 0.92 横向 0.93 0.88 0.88 0.88 0.88 0.92 柱底 纵向 0.98 0.99 0.99 0.99 0.99 0.91 横向 0.99 0.92 0.93 0.92 0.91 0.89 设防地震 柱顶 纵向 0.91 0.87 0.87 0.87 0.87 0.92 横向 0.93 0.88 0.88 0.88 0.88 0.92 柱底 纵向 0.97 0.99 0.98 0.99 0.99 0.98 横向 0.98 0.91 0.92 0.91 0.90 0.88 ② 罕遇地震 柱顶 纵向 0.91 1.15 1.15 1.15 1.15 0.92 横向 0.91 1.17 1.17 1.17 1.17 0.91 柱底 纵向 0.90 1.33 1.33 1.33 1.33 0.92 横向 0.99 1.44 1.43 1.42 1.42 0.96 设防地震 柱顶 纵向 0.91 1.15 1.15 1.15 1.15 0.92 横向 0.91 1.17 1.17 1.17 1.17 0.91 柱底 纵向 0.89 1.33 1.33 1.33 1.33 0.91 横向 0.99 1.43 1.43 1.42 1.41 0.95 ③ 罕遇地震 柱顶 纵向 0.91 1.15 1.15 1.15 1.15 0.91 横向 0.91 1.16 1.16 1.16 1.16 0.91 柱底 纵向 0.90 1.33 1.33 1.33 1.33 0.91 横向 0.99 1.43 1.43 1.43 1.41 0.95 设防地震 柱顶 纵向 0.91 1.15 1.15 1.15 1.15 0.91 横向 0.91 1.16 1.16 1.16 1.16 0.91 柱底 纵向 0.89 1.32 1.32 1.33 1.33 0.91 横向 0.99 1.42 1.42 1.42 1.41 0.95 ④ 罕遇地震 柱顶 纵向 0.91 0.87 0.86 0.86 0.86 0.91 横向 0.93 0.89 0.89 0.89 0.89 0.93 柱底 纵向 0.88 0.88 0.88 0.88 0.88 0.90 横向 0.97 0.92 0.91 0.91 0.90 0.89 设防地震 柱顶 纵向 0.84 0.87 0.86 0.86 0.86 0.91 横向 0.93 0.89 0.89 0.89 0.89 0.93 柱底 纵向 0.80 0.87 0.87 0.87 0.87 0.89 横向 0.96 0.91 0.90 0.90 0.89 0.88 输入罕遇地震和设防地震水平地震动时程后，柱顶隔震方案和柱底隔震方案与常规隔震方案建筑物每层纵向和横向最大层间位移角比见表 5，由表 5可知，遭遇罕遇地震和设防地震水平地震动时，柱顶隔震方案和柱底隔震方案与常规隔震方案建筑物相比，横向和纵向层间位移角增幅较小，层间位移角比为0.82—1.06，这是因为在采用柱顶隔震方案和柱底隔震方案的建筑物中，对隔震支座相邻上下层柱混凝土尺寸、强度等级和钢筋等进行了一定优化，使其抗变形能力增强，在一定程度上抵消了由于隔震支墩柔性、位移变形较大造成的支座相邻上下层层间位移角增大，使其与常规隔震方案建筑物层间位移角相比未发生太大变化；另外，纵向层间位移角比值基本小于横向层间位移角比值，这是由于纵向层间刚度大于横向层间刚度，使其受隔震方案变化的影响偏小。综上所述，隔震支墩长细比基本能达到普通层柱长细比。
表 5 柱顶隔震方案和柱底隔震方案与常规隔震方案建筑物层间位移角比Table 5. The displacement angle's ratio between layers of the column's top and bottom isolation project with the conventional isolation project层数 罕遇地震 设防地震 柱底 柱顶 柱底 柱顶 纵向 横向 纵向 横向 纵向 横向 纵向 横向 1 1.00 1.05 1.06 1.03 1.04 1.00 1.01 1.03 2 1.04 1.05 1.03 1.05 1.02 1.01 1.02 1.02 3 0.87 1.00 0.92 1.01 0.96 1.01 0.91 1.01 4 0.83 1.02 0.87 0.98 0.97 0.98 0.85 0.92 5 0.82 1.01 0.84 0.97 0.99 0.97 0.89 0.95 6 0.82 1.04 0.85 1.04 1.00 0.99 0.87 0.96 7 0.83 1.01 0.90 1.02 1.00 0.89 0.88 1.01 8 0.82 1.01 0.95 1.01 1.02 1.02 0.82 0.97 输入罕遇地震和设防地震水平地震动时程后，柱顶隔震方案和柱底隔震方案与常规隔震方案建筑物每层纵向和横向最大层间剪力比见表 6，由表 6可知，遭遇罕遇地震和设防地震水平地震动时，柱顶隔震方案建筑物首层位移基本由隔震支座承担，隔震支墩位移较小，基本不受剪力作用，因此剪力比为0，2层层间剪力比为1.01—1.03；柱底隔震方案建筑物首层层间剪力增幅较小，层间剪力比为1.01—1.05，这是因为柱顶隔震方案和柱底隔震方案建筑物中，对隔震支座相邻上下层柱混凝土尺寸、强度等级和钢筋等进行了一定优化，使其抗剪能力增强，在一定程度上抵消了由于隔震支墩柔性、位移变形较大造成的隔震支座相邻上下层层间剪力增大，使其与常规隔震方案建筑物层间剪力相比未发生太大变化；其他层层间剪力由于远离隔震支座所在层，在未进行柱相关参数优化的基础上，层间剪力比为0.85—1.05。综上所述，隔震支墩长细比基本能达到普通层柱长细比。
表 6 柱顶隔震方案和柱底隔震方案与常规隔震方案建筑物层间剪力比Table 6. The shear's ratio between layers of the column's top and bottom isolation project with the conventional isolation project层数 罕遇地震 设防地震 柱底 柱顶 柱底 柱顶 纵向 横向 纵向 横向 纵向 横向 纵向 横向 1 1.02 1.05 0.00 0.00 1.01 1.04 0.00 0.00 2 0.89 0.96 1.01 1.03 1.05 0.90 1.02 1.03 3 0.93 0.87 0.98 1.03 0.96 0.85 1.01 0.92 4 1.05 1.04 1.02 0.91 1.04 0.87 0.89 0.85 5 0.95 0.90 0.94 0.94 0.95 0.98 0.96 0.87 6 0.86 1.04 0.98 0.88 0.98 0.96 0.91 0.86 7 1.00 0.91 0.94 1.01 1.03 0.92 0.85 0.85 8 0.99 1.01 1.01 0.87 1.01 0.87 0.98 0.99 4. 结论
利用ANSYS软件分别对采用无隔震支座、常规隔震、柱顶隔震、柱底隔震4种方案的相同建筑进行罕遇地震和设防地震水平地震动时程反应分析，得出以下结论：
(1) 由于地震随机性很强，地面运动影响因素较多，同一建筑在输入相同场地条件和幅值的不同地震动时程后，其位移、应力等差别较大。
(2) 遭遇地震时，无隔震支座建筑物摆动频率较大，增设隔震支座后使整个建筑物摆动频率降低，周期变大；隔震建筑物大部分位移由隔震支座承担，建筑自身层间相对位移较小，减小了地震对结构、非结构构件、内部附属物品的损坏；3种隔震方案建筑物位移反应相差较小，柱顶隔震方案建筑物位移相对较小。
(3) 遭遇罕遇地震和设防地震水平地震动时，对于铅芯橡胶支座，柱顶隔震方案和柱底隔震方案建筑物隔震支座位移均小于常规隔震方案；对于普通橡胶支座，柱顶隔震方案和柱底隔震方案建筑物隔震支座位移均大于常规隔震方案，可考虑将普通橡胶支座全部换成铅芯橡胶支座或更换大一号的普通橡胶支座避免遭遇地震时可能出现的破坏。
(4) 遭遇罕遇地震和设防地震水平地震动时，柱顶隔震方案和柱底隔震方案建筑物层间位移角和层间剪力与常规隔震方案相比无实质性增大，这是因为对隔震支座相邻上下层柱混凝土尺寸、强度等级和钢筋等进行了一定优化，使其抗变形和抗剪能力增强，使其与常规隔震方案建筑物层间位移角和层间剪力相比未发生太大变化，使隔震支墩长细比基本能达到普通层柱长细比。
(5) 当隔震支墩长细比不能达到普通层柱长细比时，可通过更换大一号隔震橡胶支座或加装阻尼器增加地震发生时消耗的地震能量，以此减少隔震支座位移、层间位移角和层间剪力，使隔震支墩长细比能达到普通层柱长细比。由此可见，柱顶隔震、柱底隔震方案可行，同时推动隔震技术在多层建筑特别是学校、医院等建筑中的应用，提高隔震建筑经济性，也为《建筑抗震设计规范》(GB 50011—2010)的修订提供一定理论依据。
-
图 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
图 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 IWTH02 20 5.9310 4.91 30 6.0793 40 6.0942 50 6.1148 60 6.1354 70 6.1580 KSRH07 20 2.5778 2.51 30 2.5740 40 2.5759 50 2.5791 60 2.5745 70 2.5781 IWTH08 20 2.7308 3.37 30 2.7280 40 2.7157 50 2.7268 60 2.6971 70 2.7232
下载: 导出CSV
表 4 台站峰值数量统计
Table 4. Peak number statistics of 16 one-dimensional stations
台站名（场地类别） 显著峰数量/个 峰的数量变化 阻尼比 平滑的带宽系数 0.01 0.05 0.1 20 30 40 50 60 70 IWTH02（Ⅱ） 1 1 1 1 1 1 1 1 2 2 IBRH13（Ⅱ） 1 2 1 1 1 1 2 3 3 3 IWTH27（Ⅱ） 1 1 1 1 1 1 1 1 1 1 KSRH06（Ⅱ） 1 1 1 1 1 2 2 2 3 3 KSRH10（Ⅱ） 2 4 3 2 2 2 3 3 4 4 TCGH12（Ⅱ） 2 2 1 1 2 3 3 4 4 4 IWTH08（Ⅱ） 3 3 3 3 3 3 4 4 4 5 TKCH08（Ⅱ） 3 5 5 3 3 4 4 6 6 7 KSRH07（Ⅲ） 2 2 2 2 2 2 2 2 2 2 FKSH11 （Ⅲ） 3 4 3 2 3 3 5 5 5 5 FKSH14（Ⅲ） 5 5 4 4 4 5 6 6 7 8 IBRH10（Ⅳ） 4 4 3 2 4 4 5 5 5 5 FKSH14（Ⅲ） 5 5 4 4 4 5 6 6 7 8 NMRH04（Ⅳ） 5 7 4 2 5 5 6 7 7 8 IBRH17（Ⅱ） 6 9 2 2 6 8 9 9 9 9 NIGH11（Ⅱ） 6 4 4 3 5 6 6 7 8 11
下载: 导出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.006Chen 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.014Chen 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/zzfy20210215Ding 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.014Luo 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.013Luo 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/cjg20160814Rong 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/zzfy20190403Yao 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 -
下载:
点击查看大图
计量
- 文章访问数: 145
- HTML全文浏览量: 41
- PDF下载量: 49
- 被引次数: 0
