基于连续小波变换的斜坡动力响应特征分析——以四川长宁Ms6.0级地震为例

金 刚，王运生，何先龙，史丙新，周宇航

（1.成都理工大学地质灾害防治与地质环境保护国家重点实验室，四川 成都 610059；2.中国地震局工程力学研究所，河北 廊坊 065000；3.四川省地震局，四川 成都 610041）

0 引言

强震是人类所面临的最危险的自然灾害之一[1]。在强切割山区，山体在地震作用下往往发生震裂松动，其完整性、强度大幅度降低，诱发一系列的崩、滑灾害[2-3]。崩滑灾害堆积物又会为泥石流的发生提供物源[4-5]。上述“崩滑流”灾害造成的人员伤亡和财产损失往往不亚于地震本身[6]。随着对地震滑坡研究的深入，一些学者认为，山体在地震作用的失稳机理与山体的地震放大效应存在耦合关系。罗永红等[7]对汶川地震次生地质灾害展开调查后，结合灾害区域地貌单元，认为特殊地形尺寸与地震波长的耦合作用导致地震放大效应。范刚[8]基于自贡西山公园台站记录的地震数据，分析了地形对峰值加速度、峰值位移等地震动参数的影响规律，认为地形因素仅对水平向的地震响施加影响，而在垂直向不明显。VALAGUSSA 等[9]通过对1993年Papua New Guinea 地震、1999年ChiChi 地震、1994年Northridge地震、2004年Niigata-Chuetsu 地震、2008年Iwate-Miyagi Nairiku 地震和2008年汶川地震后造成的地貌改观展开分析，探讨了地震震级和震中距与山体失稳规模的耦合关系。CROISSANT 等[10]通过2D 数值模拟，发现地震作用下断层作为一种不稳定结构面会对山体的稳定性产生不利影响。GUAN 等[11]通过物理模拟，发现隧道深处即山体内部的水平向峰值加速度要小于同一高程山体表面处。HARZTEL 等[12]对1989年Loma Prieta 地震展开研究后表明，体波在山脊内的多向反射和散射以及瑞利波与勒夫波的相互作用是导致Robinwood山脊受震灾破坏严重的重要原因。CELEBI[13]通过对1976年Guatemala 地震和1989年Loma Prieta 地震数据展开分析，发现山脊处的地震放大效应与地震波入射角度存在耦合关系。SU 等[14]通过对地震波进行小波变换，讨论了地震波不同频率段的地面运动属性。2014年康定地震后，王运生等[15]对冷竹关地震台阵所记录到的数据进行分析，发现山体内部的地震动放大效应要弱于同高程山体表面，且山体表面傅立页谱图较山体内部，成分更为复杂。刘峡等[16]利用小波变换的方法，对地震信号进行了时域特征上的分析。

地震波是一种复杂的非稳定波型，其携带的能量在时域上和频率域上分布极不均匀[16]，一旦某时刻能量突然集中会对受震物体的结构稳定产生不利影响。目前针对地震动响应的研究多侧重于傅立页谱分析，而傅立页谱仅能表明频率域上能量分布，而忽略掉了相应频率出现的时间。小波分析作为数学放大镜，突破了傅立页变换和短时傅立页变换局限，可以清楚地表明地震波能量、时域、频率域之间的关系，以及地震波成分的演变。因此基于连续小波变换深入探讨地震在特殊山地单元下的地震动响应规律，对揭示地震致灾机理和开展防灾减灾工作意义重大。

1 监测数据来源与地震动响应分析

1.1 强震监测点概况

据中国地震台网，2019年6月17日22 时55 分，四川省宜宾市长宁县发生Ms6.0 级地震，震中位于长宁县双河镇附近，震源深度16 km。云南，四川和重庆等多地均有明显震感。成都理工大学地质灾害防治与地质环境保护国家重点实验室依托中国地震局“川西深切河谷斜坡地震动评价技术研究”项目，在雅安市石棉县南桠河两岸设立的3 个地震监测点记录到了此次Ms6.0 级地震。其中南桠河右岸鸡公山上放置有两台强震监测仪，左岸陈家包山上放置有一台强震监测仪器，其监测点分布位置如图1所示。鸡公山上两台站置于前震旦纪晋宁期具碎裂结构的弱风化花岗岩平硐内，陈家包台站置于坡残积厚覆盖层之上，下伏前震旦纪晋宁期花岗岩，在南桠河左岸坡脚有部分出露。红色剖面线所代表的地质剖面图见图2。三台站仪器放置场地稳定，无外界因素干扰。其场地属性如表1所示。用于地震监测的仪器则由中国地震局工程力学研究所研制的G01NET-3 结构与斜坡震动监测仪，采样时间间隔为0.005 12 s。同时，四川省地震局也提供了石棉县先锋台阵记录到Ms6.0 级地震监测数据作为参照基准。

图1 南桠河两岸监测仪器分布图（方位角118°）Fig.1 Distribution of monitoring sites on both banks of Nanya River

图2 南桠河两岸红线地质剖面图Fig.2 Geological profile of both banks at Nanya River

表1 各监测点所在位置场地属性Table 1 Properties of monitoring sites

1.2 石棉县南桠河两岸山体地震动响应综述

山体的地震动响应主要表现在峰值加速度（Peak Ground Acceleration, PGA）和阿里亚斯强度（Arias Intensity,AI）的放大[17]。对3 个监测点以及1 个参照点所记录到Ms6.0 级地震的数据进行低通30 Hz 滤波处理后，通过软件SeismoSignal 处理后可以读出其地震动响应的基本参数属性。其4 个点的峰值加速度、阿里亚斯强度以及峰值加速度放大倍数如表2所示。

表2 各监测点地震动响应参数Table 2 Ground motion response parameters at each monitoring site

罗永红[17]、贺建先等[18]、祁生文等[19]结合大量实测地震数据分析发现，山体的山脊处，山体转折部位以及第四纪覆盖层会对山体的地震动响应起控制作用，具体表现为上述部位的峰值加速度放大，在地震来临时候更加强烈的震动响应会导致山体发生失稳。由表2可见，同位于鸡公山上的1#及2#监测点，位于山脊处的2#监测点其峰值加速度要大于1#监测点。而3#监测点场地类型为厚覆盖层，其水平向峰值加速度远远大于1#和2#监测点，而垂直向峰值加速度放大效应则弱于水平向，但也有一定程度的放大。可以明显看出，峰值加速度的放大效应在水平和垂直方向上具有显著的差异。

2 地震信号的时频分析

傅立页变换作为一种良好的信号分析手段，在世界范围内被广泛应用[16]。但对于非稳定信号，傅立页变换只能体现出信号的频域特性，而隐略了随机信号中的时域特性。小波变换自1986年被发明以来，因其具有良好的时频分析能力，在信号处理领域内发展迅速[20]。小波变换分析不仅可以描述非稳定信号的局部变化特征，而且可以通过尺度变换，从不同频域上观察信号的演变特征。因此，相较于使用传统手段傅立页变换对地震信号频率域上展开解析，利用小波变换分析地震信号，可以反映地震信号的时频特性，明确地震时频演变规律。

2.1 连续小波变换

将四台强震监测仪记录到的地震数据做小波变换处理。所选用的变换方式为连续小波变换（Continuous Wavelet Transform），小波基函数为Mexh，其中心频率（Centre Frequency）为0.25。选择矢量尺度为1 至140，由140 至1 作公差为-0.5 的递减。由尺度-频率换算公式计算可得，其覆盖频率范围0.35～50 Hz。

地震数据进行连续小波变换后，可以三维的方式呈现，其中x轴表示时间，单位为s；y轴表示频率，单位为Hz；z轴表示地震加速度幅值，单位为20 gal。3 个监测点的地震时频图如表3所示。地震信号经过连续小波变换之后，透过传统时程图增加了频率信息，做到了时域和频率域上的结合。

表3 监测点3 分量连续小波分解图Table 3 Three-component continuous wavelet decomposition

2.2 地震信号时频成分统计

地震信号经过连续小波变换后，其时域频域的结合有助与对地震信号成分的分析。4 处监测点的3 分量地震信号时频图中，随着时间增加，有两个波峰，其中峰值较低先到达的为P 波，而后到达的峰值较大的则为S 波。在S 波后面则为一段由体波激发的面波。S 波在到达受震地表后，分异为高频和低频两种成分：形状尖锐、变化剧烈的低频成分与光滑平顺、平稳降低的高频成分。地震发生时，受震物体的破坏主要依靠S 波和体波交会激发的面波[21]。P 波主要为垂直方向上的振动，且幅值小于S 波，对受震物体结构破坏较小[6]，因此本文分析重点为S 波成分变化。现对各监测点三分量上P 波和S 波波峰空间坐标进行测量，其结果如表4所示，其中第一行表示波峰出现的时间位置，第二行表示波峰出现的频率位置，第三行表示波峰幅值绝对值。

3 放大效应分析

3.1 山脊部位放大效应分析

历史地震监测表明，地震作用下山脊部位往往放大效应更为显著[17-18]。根据表2的时程数据分析结果来看，2#监测点所处的山脊部位，其水平向上PGA和AI数值要高于同在山腰部位的1#监测点。其中PGA放大系数（2#/1#）在东西南北和垂直向上分别为1.74、1.23和1.05；AI放大系数（2#/1#）在三个方向上则分别为1.52、1.47 和0.91。这说明2#监测点所在山体部位地震动力响应更为强烈。由表4可以看到地震成分的演变，1#与2#监测点的S 波低频部分，其波峰在频域上的发生位置完全相同，这可能与两处监测点位于一座山体有关。相较于1#监测点，2#监测点的东西向、南北向和垂直向S 波高频成分，波峰幅值放大倍数分别为1.73、1.53 和1.16，而S 波低频部分的放大倍数则分别为1.14、1.04 和0.76。2#监测点三分量S 波高频波峰幅值全部高于低频部分，而1#监测点则相反。即在基岩场地下，山脊相较于山腰，地震放大效应主要体现为高频成分能量的增加，低频部分有小幅增大，但远不如高频部分。此外，放大效应还存在显著的方向性，无论是在PGA、AI放大系数上还是地震波峰成分演变上，水平方向的放大效应要更为显著。受地震作用，在2#点所在的山脊处，落石与斜坡变形较其他地方更为发育（图3）。

表4 监测点波峰成分统计表Table 4 Statistical table of signal peak at each site

图3 2#监测点附近山脊部位的落石以及斜坡变形Fig.3 Rock falls and slope deformation around monitoring site 2# in Shimian County

3.2 覆盖层场地放大效应分析

3#监测点的场地条件为厚覆盖层，下伏同晋宁期花岗岩，仅在南桠河左岸边坡支护剖面上有所出露。通过表3可以看到，3#监测点水平向时频三维图上高频成分衰减剧烈，其低频波峰幅值在东西向和南北向更是分别高达0.630 和0.656，相较于同在山脊部位的2#监测点S 波最大幅值，其放大倍数达8.08 和3.53。但在垂直向，依然可以清楚地看到S 波的高频和低频成分，低频成分在时间序列上出现的次数更多。此外，参照点场地条件为薄覆盖层，下伏同套岩组。同样可以看到，薄覆盖层场地的低频部分相较于基岩场地，其低频部分已经有所发育。而PGA放大系数（3#/2#）在东西、南北和垂直向上的放大系数为3.48、2.79 和1.25；AI放大系数（3#/2#）则分别为27.4、7.36 和3.39。地震时，震感更加强烈。

现有研究表明，当成层厚覆盖层接收来自基岩传来的剪切波时，在地震波向地表传播的过程中，会发生多次反射与折射，波能量发生叠加而进一步增强，尤其为长周期波更为显著，这与表征能量的AI显著放大相吻合[21-22]。此外，覆盖层土体对入射剪切波具有吸收和放大作用。当剪切波由基岩入射到覆盖层后，其中短周期成分被吸收，而长周期成分由于与土体自身振动周期相近，极易发生共振作用[22]。为了探究3#厚覆盖层场地的S 波低频部分剧烈变化是否为共振效应引起，现在对3#监测数据进行加速度反应谱分析。加速度反应谱即在某一阻尼比的作用下，反应一系列单质点在振动时，其最大响应绝对值与场地结构周期的关系的频谱图。现采用0.05，0.1，0.2 的阻尼比[20]，对3#监测点三分量加速度反应谱展开计算，用于评价场地的特征周期。其三分量反应谱图如表5所示。加速度反应谱中有多个波峰，表明覆盖层成分不均匀。其最高点所代表的时间即为场地的特征周期[20]，分析后其特征周期如表6。同时，由同样方法求得基岩场地特征周期如表7。

表5 3#监测点场地3 分量反应图谱Table 5 Three-component response spectrum of monitoring site No.3

由表6和表7可知，基岩场地特征周期远小于厚覆盖层场地。3#监测点东西向和南北向特征频率分别为1.35 Hz 和1.13 Hz，即场地的特征频率十分接近于S 波低频成分波峰峰值处频率，频率差值分别为0.05 Hz 和0.03 Hz。而垂直向特征频率则与水平向有较大差异。在水平方向上，由于频率接近，地震发生时，S 波中低频成分振动与场地自振交汇，产生强烈的共振效应，导致频率域上其特征周期附近区间内振动峰值显著增加。加之厚覆盖层土体对短周期波成分的吸收效应和长周期波成分的折射叠加效应，3#监测点的S 波高频部分衰减显著[21-22]。厚覆盖层对地震波显著放大，地震时其震感更加强烈，加重震害损失，这一现象在震中长宁灾害调查中得以体现。3#点所在村落多为刚度大，结合性差的低矮砖混结构房屋和土坯房，此次地震由于震中距较大，震级较小，未对3#点所在村落建筑物造成严重损害。但如果震级加大，受到山脊放大效应和覆盖层放大效应的双重影响，震感强烈，相较于基岩场地这里会遭受更加严重地表震害。

表6 3#监测点三分量特征周期/频率表Table 6 Dominant period / frequency of site 3#

表7 1#与2#监测点三分量特征周期Table 7 Dominant period of site 1# and 2#

4 结论

基于连续小波变换的信号处理手段，对石棉县城南桠河两岸强震监测点所记录到的长宁Ms6.0 级地震数据，作三维时频图谱，并结合时域和频域对地震信号成分进行分析。现有结论如下：

（1）通过连续小波变换所呈的三维时频图，可以清楚看到S 波到达地表时分为高频和低频两种成分。其中高频成分主频在3.5 Hz 附近，而低频成分主频在1.1 Hz附近，两种成分在时域上近乎同时出现。

（2）地震时山脊部分的地震动响应更加强烈。山脊部位场地类型为基岩，通过时频分析结果可知，放大效应主要体现为S 波高频成分能量的增加，而低频部分则无明显变化。

（3）覆盖层场地对S 波具有显著的低频放大和高频衰减效应。S 波在厚覆盖层的折射和反射作用，以及3#覆盖层场地水平方向的自振与S 波低频成分交汇产生共振效应，是导致覆盖层场地低频放大效应显著的原因。地震发生时，3#场地震感强烈，地表震害将更加严重。

（4）地震动响应规律具有极强的方向差异性，从时程分析结果和时频分析结果可以看出，其差异性主要体现在水平向和垂直向上，其水平向的放大效应相较于垂直向要更为显著。