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

SH波斜入射下水平成层场地地震动特性研究

陈三红 张郁山 郝明辉

崔恩文,郭恩栋,王倩,孙逊,2021. 改进的标准化累积绝对速度及其在高速铁路地震报警应用中阈值的确定. 震灾防御技术,16(3):527−532. doi:10.11899/zzfy20210312. doi: 10.11899/zzfy20210312
引用本文: 陈三红,张郁山,郝明辉,2024. SH波斜入射下水平成层场地地震动特性研究. 震灾防御技术,19(1):88−95. doi:10.11899/zzfy20240109. doi: 10.11899/zzfy20240109
Cui Enwen, Guo Endong, Wang Qian, Sun Xun. Improved Calculation of Standardized Cumulative Absolute Velocity and Its Seismic Alarm Threshold for High-speed Railway[J]. Technology for Earthquake Disaster Prevention, 2021, 16(3): 527-532. doi: 10.11899/zzfy20210312
Citation: Chen Sanhong, Zhang Yushan, Hao Minghui. Ground Motion Characteristics of Horizontal Layered Site Based on SH Wave Oblique Incidence[J]. Technology for Earthquake Disaster Prevention, 2024, 19(1): 88-95. doi: 10.11899/zzfy20240109

SH波斜入射下水平成层场地地震动特性研究

doi: 10.11899/zzfy20240109
基金项目: 地震科技星火计划( XH21049Y);国家重点研发计划重点专项(2018YFC1504600)
详细信息
    作者简介:

    陈三红,女,生于1994年。硕士研究生,工程师。主要从事复杂场地效应等方面的研究。E-mail:1561980175@qq.com

Ground Motion Characteristics of Horizontal Layered Site Based on SH Wave Oblique Incidence

  • 摘要: 基于位移场理论和传递矩阵法,建立SH波斜入射下水平成层场地非线性地震反应等效线性化求解方法,研究SH波不同角度入射下场地地表峰值加速度放大系数及反应谱谱比变化规律。研究结果表明,地震波入射角度对场地地表峰值加速度影响明显,在大多数情况下,地震波斜入射下地表峰值加速度大于垂直入射;地震波斜入射下反应谱特性变化明显,入射角度对高频段反应谱的放大效应显著,且入射角度对不同类型场地的影响周期范围不同,入射角度对Ⅲ类场地周期影响范围明显宽于Ⅱ类场地。
  • 众所周知,我国是地震频发的国家,地震可能造成高速铁路列车脱轨事故,为此,我国建立了相应的地震报警及紧急处置系统,力争在潜在破坏性地震到达之前发出报警,降低列车发生脱轨或倾覆的概率。目前广泛采用的地震报警参数主要为峰值加速度PGA,然而PGA主要体现地震波幅值,未考虑频谱及持时影响,对于无破坏性的高频小震可能产生误报现象。为避免这种干扰,郭恩栋等(2015)引入谱强度SI作为高速铁路地震报警参数,并给出相应的计算方法,SI从能量方面出发,考虑了幅值及频谱影响,忽略持时的影响,但相比PGA更合理。美国电力研究所EPRI(O'Hara 等,1991)在核电厂地震监测报警中引入累积绝对速度CAV作为报警参数,研究发现低频震动对CAV的贡献明显高于高频震动。我国三代核电站AP1000参考国外标准也引入了CAV作为地震报警参数(胡进军等,2013)。林宜娴(2002)提出CAV、PGA、PGV等参数可作为地震报警参数,其稳定性好,且确定性较高。Kostov(2005)根据欧洲几次强地震动数据,分析得出标准化累积绝对速度CAVSTD在预测地震动破坏能力方面明显优于PGA。伊斯坦布尔在地震预报警系统中利用累积绝对速度快速判别强地震动(Erdik等,2003)。Fahjan等(2011)在地震报警的触警算法中引入CAV,并提出新的用于报警系统的参数。黄俊等(2014)引入CAV作为报警参数,分析其在高速铁路报警中的适用性,研究发现CAV可有效排除破坏性小的高频小震对高速铁路列车地震报警的干扰。目前虽对CAV在地震报警中的应用有了深入研究,但在高速铁路列车地震报警中的研究较少,仅通过相关的地震波数据处理得出,未考虑车轨动力响应关系,且给出的计算参数为CAVSTD,为此,需探讨考虑车轨关系的基于相关起算阈值的CAV计算方法。

    目前,地震工程界主要研究CAV、CAVSTD、CAV5。其中CAV于1988年由美国电力研究所(O'Hara等,1991)在核电厂报警中第2级判别标准中提出的,其概念为记录地震加速度数据的绝对值对记录时间进行积分,记为累积绝对速度CAV:

    $${\rm{ CAV}} = \mathop \int \nolimits_0^{{{{t}}_{{\rm{max}}}}} \left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t$$ (1)

    为提高CAV在核电厂报警中的准确度,EPRI在报告中引入标准化累积绝对速度CAVSTD,其可有效避免长持时记录的地震波尾部对应的幅值小且破坏性小的地震动对其数值的影响,其定义如下:将记录的加速度时程分成以1 s为时间间隔的N个小段,当1 s时间段内存在加速度绝对值≥0.025 g时,对该秒内加速度数值的绝对值进行积分,否则该秒内的积分值记为零,然后将整个记录时程进行相加求和,表达式为:

    $$ {\rm{CA{V_{STD}}}} = \mathop \sum \limits_{{i=1}}^n \mathop \int \nolimits_{{i}}^{{{i}} + 1} {W_i}\left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\;\;\;\;\;{W_i} = \left\{ {\begin{array}{*{20}{c}} {0,\left| {{{a}}\left( {{t}} \right)} \right| < 0.025\;g}\\ {1,\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.025\;{{g}}} \end{array}} \right.$$ (2)

    Kramer等(2002)于2002年提出CAV5,用于场地液化判别,其概念为:将记录的加速度时程中对加速度数值的绝对值≥0.005 g进行积分求和,表达式为:

    $${{\rm{ CAV}}_5} = \displaystyle\int \nolimits_0^\infty \left\langle {\rm{\chi }} \right\rangle \left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\left\langle {\rm{\chi }} \right\rangle = \left\{ {\begin{array}{*{20}{c}} {0\;\left| {{{a}}\left( {{t}} \right)} \right| < 0.005\;g}\\ {1\;\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.005\;{{g}}} \end{array}} \right. $$ (3)

    式中,a(t)为加速度时程,Wi为权重。

    本文引入31自由度的SY97477车辆模型(崔恩文,2014),结合车轨典型非线性接触关系,依据列车安全行驶的3项评判指标(脱轨系数、轮重减载率及横向力),分析车体在轨道不平顺(我国高速铁路无砟轨道标准谱)与简谐波(仅考虑横向输入,与列车行驶方向垂直)共同激励下的动态响应,考虑95%的安全保证,引入2.55倍的动力放大系数(孙汉武等,2007),计算安全极限状态(超出3项指标任意1个限值)下不同车速对应的简谐波幅值随周期的变化关系,如图1所示。具体模型参数、车轨接触关系及评判标准可参考崔恩文(2014)的研究。

    图 1  列车安全运行极限状态下不同车速对应的简谐波幅值随周期的变化关系
    Figure 1.  The relationship between the amplitude of simple harmonic wave corresponding to different speeds and periods under the safe running limit state of the train

    图1可知,车速一定时,即使幅值很小(<0.025 g)的低频简谐波也可能使列车脱轨,因此不能简单地套用标准化累积绝对速度CAVSTD计算公式,有必要研究改进的标准化累积绝对速度应用在高速铁路地震报警阈值中的计算方法。车速一定时,简谐波幅值随周期增加逐渐递减,最终趋于常数,且周期越大,不同车速对应的简谐波幅值之间的差值越小。车速为200 km/h时,当简谐波周期达2.2 s时,其幅值开始趋于常数,该常数约为0.008 g,本文规定此时的周期为简谐波幅值趋于常数对应的起始周期Tst,该常数记为Amin。不同车速下简谐波幅值趋于常数对应的起始周期及该常数的数值如表1所示。

    表 1  不同车速下简谐波幅值趋于常数对应的起始周期及该常数的数值
    Table 1.  The initial period of the simple harmonic wave amplitude tending to be constant at different speeds and the value of the constant
    项目车速/(km·h−1
    200250300350400
    起始周期Tst/s2.22.32.12.01.7
    Amin/g0.0080.0040.0040.0040.004
    下载: 导出CSV 
    | 显示表格

    表1可知,车速为200 km/h与车速≥250 km/h时对应的Amin不同,车速≥250 km/h时Amin相同,但起始周期Tst随车速的增加呈递减趋势,说明车速越大,简谐波趋于常数所需的周期越小。

    车速为200 km/h时,幅值低于0.008 g的简谐波对列车运行安全的影响可忽略不计。车速≥250 km/h时,幅值低于0.004 g的简谐波对列车运行安全的影响可忽略不计。由此可推算出改进的标准化累积绝对速度起算阈值,如表2所示。

    表 2  不同车速下改进的标准化累积绝对速度起算阈值
    Table 2.  Starting threshold of improved calculation of standardized cumulative absolute velocity at different speeds
    项目车速/(km·h−1
    200≥250
    改进的标准化累积绝对速度起算阈值/g0.0080.004
    下载: 导出CSV 
    | 显示表格

    考虑改进的标准化累积绝对速度为累积过程的算法,其值对于加速度超过0.004 g或0.008 g的有限个数的高频小震会很小,因此,可较好地排除高速铁路沿线的高频小震,是较好的地震监测报警参数。

    魏智祥(2015)分析了京武高速铁路线布置的强震仪监测到的列车运行时震动记录数据,发现列车运行产生的震动记录明显高于背景噪声,且发现东西、南北、垂直方向列车振动噪声幅值均≤0.002 g,由此可知本文设定的不同车速下改进的标准化累积绝对速度可有效排除列车振动的影响。

    综上所述,本文给出改进的标准化累积绝对速度用于高速铁路地震报警是可行的。

    参照CAVSTD计算方法,本文给出车速为200 km/h时改进的标准化累积绝对速度CAV8计算公式如下:

    $${{\rm{CAV}}_8} = \mathop \sum \limits_{i=1}^n \mathop \int \nolimits_{{i}}^{{{i}} + 1} {W_i}\left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\;{W_i} = \left\{ {\begin{array}{*{20}{c}} {0,\;\left| {{{a}}\left( {{t}} \right)} \right| < 0.008\;g}\\ {1,\;\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.008\;{{g}}} \end{array}} \right.$$ (4)

    $ a\left(t\right) $$< 0.008\;g$时,${{W}}_{{i}}$=0;当$\left|{a}\left({t}\right)\right|\geqslant 0.008\;g$时,${{W}}_{{i}}$=1。

    本文给出的CAV8定义为:将加速度时程按1 s为时间间隔进行划分,如果在这1 s内加速度绝对值≥0.008 g,则对该秒内的加速度绝对值进行积分,否则忽略该秒内的积分值,然后将整个时程中每秒积分得到的数值进行累积求和。

    车速为250 km/h时改进的标准化累积绝对速度CAV4计算公式如下:

    $${ {\rm{CAV}}_4} = \mathop \sum \limits_{i=1}^{{n}} \mathop \int \nolimits_{{i}}^{{{i}} + 1} {W_i}\left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\;\;{W_i} = \left\{ {\begin{array}{*{20}{c}} {0,\;\left| {{{a}}\left( {{t}} \right)} \right| < 0.004\;g}\\ {1,\;\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.004\;{{g}}} \end{array}} \right.$$ (5)

    ${a}\left({t}\right)$$< 0.004\;g$时,${{W}}_{{i}}=0$;当$\left|{a}\left({t}\right)\right|\geqslant 0.004\;{g}$时,${{W}}_{{i}}$=1。

    本文给出的CAV4定义为:将加速度时程按1 s为时间间隔进行划分,如果在这1 s内加速度绝对值≥0.004 g,则对该秒内的加速度绝对值进行积分,否则忽略该秒内的积分值,然后将整个时程中每秒积分得到的数值进行累积求和。

    为确定基于改进的标准化累积绝对速度高速铁路地震报警阈值,本文从Peer数据库选取6条不同场地地震波(表3),对应的加速度反应谱如图2所示,对其进行归一化处理,分析加速度幅值与评判列车脱轨3项指标之间的关系,确定不同车速下影响行车安全的6条地震波加速度上限值,按照本文给出不同车速下改进的标准化累积绝对速度计算公式,给出6条地震波对应的改进的标准化累积绝对速度上限值,如表4所示。

    表 3  地震波相关信息
    Table 3.  Seismic wave information
    项目地震波台站名称
    Cape MendocinoTCU045CapitolaHectorEL Centro Imp. Co. CentPoe Road
    震级 7.01 7.62 6.93 7.13 6.54 6.54
    峰值加速度/g 0.376 5 0.008 1 0.528 5 0.265 6 0.357 9 0.446 3
    v30/(m·s−1 567.78 704.64 288.62 726.00 192.05 316.64
    加速度记录持时/s 36.0 45.0 39.9 45.3 40.0 22.3
    下载: 导出CSV 
    | 显示表格
    图 2  6条地震波加速度反应谱
    Figure 2.  Acceleration response spectra of six seismic waves
    表 4  不同车速下6条地震波加速度上限值与改进的标准化累积绝对速度上限值
    Table 4.  The upper limit of acceleration amplitudes and upper limit improved calculation of standardized cumulative absolute velocity CAV of six seismic waves at different speeds
    地震波台站名称车速200 km/h车速250 km/h车速300 km/h车速350 km/h车速400 km/h
    加速度
    上限值/
    g
    改进的标准化
    累积绝对速度
    上限值/g-sec
    加速度
    上限值/
    g
    改进的标准化
    累积绝对速度
    上限值/g-sec
    加速度
    上限值/
    g
    改进的标准化
    累积绝对速度
    上限值/g-sec
    加速度
    上限值/
    g
    改进的标准化
    累积绝对速度
    上限值/g-sec
    加速度
    上限值/
    g
    改进的标准化
    累积绝对速度
    上限值/g-sec
    Cape Mendocino0.0740.170.0700.170.0550.130.0400.090.0390.09
    TCU0450.1010.550.0900.490.0780.430.0620.340.0390.21
    Capitola0.1090.340.0780.250.0620.200.0540.170.0390.11
    Hector0.1010.270.0700.190.0620.170.0540.140.0470.12
    EL Centro Imp. Co. Cent0.0780.160.0700.180.0540.130.0390.080.0270.05
    Poe Road0.0740.190.0540.140.0430.110.0390.100.0350.08
    下载: 导出CSV 
    | 显示表格

    表4可知,车速相同时,即使地震波加速度相同,其对应的改进的标准化累积绝对速度不一定相同,如车速为200 km/h时,台站名称为Cape Mendocino与Poe Road的地震波加速度相同,但对应的改进的标准化累积绝对速度不同,且Cape Mendocino台站加速度记录持时较大,说明改进的标准化累积绝对速度与频谱紧密相关。车速≥250 km/h时,对于同一台站的地震波,峰值加速度越大其对应的改进的标准化累积绝对速度越大,如台站名称为Cape Mendocino的地震波,说明对于同一地震波而言,加速度影响累积绝对速度。根据改进的标准化累积绝对速度定义,加速度记录持时影响其数值。综上所述,改进的标准化累积绝对速度数值不仅与加速度幅值有关,还与频谱及加速度记录持时密切相关,可较好地反映地震动三要素。

    选取6条地震波改进的标准化累积绝对速度最小值作为报警阈值,给出不同车速下基于改进的标准化累积绝对速度高速铁路地震报警阈值,如表5所示。由表5可知,改进的标准化累积绝对速度报警阈值与车速整体呈线性递减关系,中间车速对应的改进的标准化累积绝对速度可按线性插入计算。

    表 5  不同车速下改进的标准化累积绝对速度高速铁路地震报警阈值
    Table 5.  High-speed earthquake alarm threshold based on improved calculation of standardized cumulative absolute velocity thresholds at different vehicle speeds
    项目车速/(km·h−1
    200250300350400
    报警阈值/g-sec0.160.140.110.080.05
    下载: 导出CSV 
    | 显示表格

    (1)通过建立车轨模型,分析列车安全运行极限状态下简谐波幅值与周期的关系,得出不同车速下改进的标准化累积绝对速度起算阈值,并论证改进的标准化累积绝对速度在高速铁路地震报警中的可行性。

    (2)参照标准化累积绝对速度CAVSTD计算公式,本文给出不同车速下改进的标准化累积绝对速度计算公式。

    (3)根据本文给出的改进的标准化累积绝对速度计算公式,选取6条不同场地的地震波,分析列车在轨道不平顺及地震波共同激励下的动力响应,计算不同车速下改进的标准化累积绝对速度报警阈值,以供参考。

  • 图  1  SH波不同角度入射下自由表面加速度峰值变化比较

    Figure  1.  Compare of the peak acceleration excited by SH wave with different angles

    图  2  天然地震波加速度时程曲线

    Figure  2.  Acceleration time-history curves of input ground motions

    图  3  天然地震波加速度反应谱曲线

    Figure  3.  Acceleration Spectrum curves of input ground motions

    图  4  场地土剪切模量比和阻尼比随剪应变变化曲线

    Figure  4.  The curve shows change of shear modulus ratio and damping ratio change with shear strain

    图  5  地震波入射角度对峰值加速度放大系数的影响

    Figure  5.  Influence of incident angle on PGA amplification coefficient

    图  6  所有地震波以60°入射下的反应谱谱比曲线

    Figure  6.  Spectral ratio curves under the inputs of natural ground motions at 60°incidence

    图  7  不同强度地震波斜入射下反应谱放大系数

    Figure  7.  Magnification of response spectrum under oblique incidence with different intensities of seismic waves

    表  1  天然地震波基本信息

    Table  1.   Basic information of natural ground motions

    编号地震名称震级/级台站名称震中距/km记录分量
    11979年美国Imperial Valley-06地震6.5Parachute Test Site12.69225
    21999年台湾集集地震7.6TCU04616.74EW
    31978年伊朗Tabas地震7.3Tabas2.05L1
    41999年台湾集集地震7.6TCU0750.89EW
    52004年日本新泻地震6.6TYMH05175.52NS
    61999年美国Hector Mine地震7.1Amboy43.05090
    71999年台湾集集地震7.6TCU1012.11EW
    82010年墨西El Mayor-Cucapah地震7.2El Centro Array #1211.26000
    下载: 导出CSV

    表  2  千叶场地计算模型参数

    Table  2.   Model parameters of the Chiba-ken site

    土层编号土类名称剪切波速/
    (m·s−1
    密度/(g·cm−3层厚/m
    1亚黏土1401.755
    2砂质黏土3201.805
    3细砂3202.0014
    4细砂4202.1016
    5基岩输入5402.20
    下载: 导出CSV

    表  3  地铁17号线某站场计算模型参数

    Table  3.   Model parameters of the Line 17 Subway site

    土层编号土类名称剪切波速/
    (m·s−1
    密度/(g·cm−3层厚/m
    1粉质黏土1921.9010
    2细砂2361.9525
    3粉质黏土3312.0010
    4细砂3552.0026
    5粉质黏土3622.0512
    6粉质黏土4162.1017
    7粉质黏土4672.1528
    8基岩输入5052.20
    下载: 导出CSV
  • 陈维, 2007. 多层土粘弹性边界的二维斜入射SH波时域输入. 山西建筑, 33(25): 84—85

    Chen W. , 2007. A method of inclined-incidence SH waves input on viscous-spring boundary in multilayered soil. Shanxi Architecture, 33(25): 84—85. (in Chinese)
    范留明, 赵钦, 2014. 成层半空间出平面自由波场的界面子波算法. 地震工程学报, 36(4): 765—772

    Fan L. M. , Zhao Q. , 2014. Interfacial wavelet superposition method for out-of-plane wave motion of a free field in a layered half space. China Earthquake Engineering Journal, 36(4): 765—772. (in Chinese)
    高智能, 卓卫东, 谷音, 2017. SH波斜入射时有阻尼成层介质自由场的一维化时域算法. 振动与冲击, 36(16): 37—43, 84

    Gao Z. N. , Zhuo W. D. , Gu Y. , 2017. A 1 D time-domain method for free field motion in layered media with damping under obliquely incident SH wave. Journal of Vibration and Shock, 36(16): 37—43, 84. (in Chinese)
    廖河山, 陈清军, 徐植信, 1994. SH波倾斜入射时土层的非线性响应分析. 同济大学学报, 22(4): 517—522

    Liao H. S. , Chen Q. J. , Xu Z. X. , 1994. Nonlinear responses of layered soils to obliquely incident SH waves. Journal of Tongji University, 22(4): 517—522. (in Chinese)
    刘晶波, 王艳, 2006. 成层半空间出平面自由波场的一维化时域算法. 力学学报, 38(2): 219—225

    Liu J. B. , Wang Y. , 2006. A 1-D time-domain method for 2-D wave motion in elastic layered half-space by antiplane wave oblique incidence. Chinese Journal of Theoretical and Applied Mechanics, 38(2): 219—225. (in Chinese)
    王笃国, 赵成刚, 2016. 地震波斜入射时二维成层介质自由场求解的等效线性化方法. 岩土工程学报, 38(3): 554—561

    Wang D. G. , Zhao C. G. , 2016. Two-dimensional equivalent linear seismic analysis of free field in layered half-space due to oblique incidence. Chinese Journal of Geotechnical Engineering, 38(3): 554—561. (in Chinese)
    尤红兵, 赵凤新, 荣棉水, 2009. 地震波斜入射时水平层状场地的非线性地震反应. 岩土工程学报, 31(2): 234—240

    You H. B. , Zhao F. X. , Rong M. S. , 2009. Nonlinear seismic response of horizontal layered site due to inclined wave. Chinese Journal of Geotechnical Engineering, 31(2): 234—240. (in Chinese)
    Heymsfield E. , 2000. Two-dimensional scattering of SH waves in a soil layer underlain with a sloping bedrock. Soil Dynamics and Earthquake Engineering, 19(7): 489—500. doi: 10.1016/S0267-7261(00)00030-0
    Katayama T. , Yamazaki F. , Nagata S. , et al, 1990. A strong motion database for the Chiba seismometer array and its engineering analysis. Earthquake Engineering & Structural Dynamics, 19(8): 1089—1106.
    Lee V. W. , Luo H. , Liang J. W. , 2006. Antiplane (SH) waves diffraction by a semicircular cylindrical hill revisited: an improved analytic wave series solution. Journal of Engineering Mechanics, 132(10): 1106—1114. doi: 10.1061/(ASCE)0733-9399(2006)132:10(1106)
    Liu G. , Chen H. T. , Liu D. K. et al. , 2010. Surface motion of a half-space with triangular and semicircular hills under incident SH waves. Bulletin of the Seismological Society of America, 100(3): 1306—1319. doi: 10.1785/0120090273
    Tsaur D. H. , Chang K. H. , 2009. Scattering and focusing of SH waves by a convex circular-arc topography. Geophysical Journal International, 177(1): 222—234. doi: 10.1111/j.1365-246X.2008.04080.x
    Tsaur D. H. , 2011. Scattering and focusing of SH waves by a lower semielliptic convex topography. Bulletin of the Seismological Society of America, 101(5): 2212—2219. doi: 10.1785/0120100324
  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  110
  • HTML全文浏览量:  65
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-12-02
  • 刊出日期:  2024-03-31

目录

/

返回文章
返回