四川雅安地区场地地震动影响特征分析1

闫静茹 张郁山 邓 菲 李 伟

（中国地震灾害防御中心, 北京 100029）

2022 年6 月1 日四川雅安地区发生6.1 级地震，现场震害调查发现地震造成道路中断，隧道及公路破坏严重。这种灾害情况一方面与地形结构较复杂、山区落石较多有关，另一方面还应关注场地条件导致的震害表现，尤其是峰值加速度、峰值速度、峰值位移等参数的放大效应影响。大量研究表明，场地的地震效应与场地结构特征密切相关，场地条件不同对基岩运动的放大缩小效应不同，引起的地震灾害形式不同。近年来，越来越多的学者关注近地表土壤条件对地面运动的影响，冯希杰等（2001）通过实际土层地震反应结果的统计分析和强震加速度观测结果的对比，讨论了不同场地条件对基岩峰值加速度的放大效应及特点；
薄景山等（2003a，2003b）通过构造场地计算剖面，研究了不同土层结构对地表加速度峰值及反应谱平台值的影响；
陈党民等（2013）基于西安地区大量钻孔资料，进行了典型场地剖面的地震反应计算，得到场地峰值放大系数随等效剪切波速、覆盖层厚度及基岩输入地震动强度的变化关系。综上所述，学者们关于场地效应的研究多集中于峰值加速度及加速度反应谱，对峰值速度、峰值位移等其他地震动参数的研究较少。

通过工程实践及相关研究表明（O’Rourke 等，2001；
陈波，2013），峰值速度、峰值位移可表征地震动的中低频信息，对中长周期结构影响较大，尤其是峰值速度与地震震级有较好的相关性，可提供有关地震动频率含量和强震持时的有关信息，相较于其他强度指标，峰值速度与变形需求的相关性更好（Akkar 等，2005），同时，O’Rourk 等（2001）在生命线工程岩土性能研究中发现峰值速度与埋深管道的地震破坏有密切关系，推测本次隧道与公路的严重破坏及峰值速度的变化关系有关。本文在全面分析雅安地区场地地震动特征的基础上，以峰值速度作为地震动参数指标，研究四川雅安地区场地条件引起的峰值速度放大效应，分析等效剪切波速与覆盖土层厚度及峰值速度的变化关系，研究雅安地区峰值速度引发震动反应的一般特征。本文研究成果可用于雅安地区场地峰值速度放大倍数指标的确定，为地震备灾工作提供技术支持。同时，本文明确指出本地区软土场地更易引起场地峰值速度的放大，可为生命线工程的抗震选址等工作提供参考。

雅安地区钻孔基础资料的获取来源于已通过评审的雅安地区重大工程场地地震安全性评价项目，通过收集整理，获取雅安地区钻孔有效数据27 个，基本可描述雅安地区的场地情况。根据对钻孔数据的土层分布、剪切波速、覆盖层厚度、成孔信息进行梳理后，建立相应钻孔的土层反应模型，同时，根据《建筑抗震设计规范》（GB 50011−2010）（中华人民共和国住房和城乡建设部等，2010）中场地类别划分标准，建立的27 个模型均属于Ⅱ类场地，钻孔分布情况如表1 所示。

表1 场地钻孔资料分布情况Table 1 Site borehole data distribution map

为更清楚地认识雅安地区的地形条件及场地土层情况，图1 给出了本次地震震中及芦山县地形地貌，图2 给出了所选雅安地区场地钻孔土层的剪切波速随埋深的变化曲线。由图1 和图2 可知，该地区地形有一定起伏变化，周边山区较多，选取的场地土层较坚硬，大部分土层等效结构剪切波速大于350 m/s，且场地上覆盖土层较浅。

图1 雅安及周边芦山县地形地貌Fig. 1 Topographic and geomorphological map of Ya"an and surrounding Lushan

图2 场地剪切波速随埋深变化Fig. 2 Site shear wave velocity varies with burial depth

由于发震时间距今较短，未搜集到本次地震较有效的强震观测记录信息，但地震备灾及灾后重建需建立在基于场地条件的地震灾害影响分析基础上。本文为全面分析雅安地区场地条件对强震输入峰值速度的影响效应，为本地区今后的地震灾害防御工作提供技术支撑，本文采用可搜集到的雅安芦山县7.0 级地震的强震记录作为基岩输入时程（表2、图3）。同时，为全面考虑当地震输入峰值加速度一致时，峰值速度及频谱分布特征不同造成的结果变化，从NGA 数据库中选取震级和震源深度较相近的地震事件作为补充，NGA 强震数据如表3、图3 所示。

图3 基岩反应谱Fig. 3 Bedrock response spectrum

表2 输入基岩时程信息（芦山地震强震数据）Table 2 Input bedrock time-histories information（Data from Lushan earthquake）

为分析雅安地区场地地震动反应特征，本文设计如下计算分析方案：

（1）分别建立场地计算模型；

（2）利用表2、表3 中的时程记录信息作为输入进行土层反应分析；

表3 输入基岩时程信息（NGA 强震数据）Table 3 Input bedrock time-histories information（Data from NGA）

（3）利用等效线性化方法进行场地地震反应计算。

3.1 计算结果

3.1.1 场地地表加速度反应谱

图4 为不同基岩输入下场地地表地震反应谱。由图4 可知，地表反应谱谱形受基岩输入谱谱形影响，不同钻孔反应谱谱形较一致，说明选取的雅安地区钻孔场地基岩上覆土层结构变化较小，土类较一致。对比芦山地震及NGA 强震记录输入下的场地反应计算结果，在相同PGA 基岩输入的基础上，由于输入谱谱形的不一致，导致场地加速度峰值及反应谱变化不同，尤其是在记录0789 与HYT-0201 基岩输入下，计算反应谱和峰值加速度均出现明显差异，这说明场地的震动破坏不仅与震级和震源深度有关，而且与地震发生时产生的基岩地震动频谱分布特征密切相关。此外，除记录0789 基岩输入外，随着震动强度的增加，谱加速度在短周期范围内出现显著的发散，这是由于记录0789 基岩输入长周期成分较丰富，能量较强，对地震峰值反应起主导作用，而大量的长周期成分不易被土层过滤，受土层放大效应的影响较小，短周期反应谱谱值发散不明显。同时，不同地震作用下的场地反应谱峰值多处于0.2 s 周期范围，说明雅安地区场地卓越周期处于此区间，可导致自振周期为0.2 s 的建筑物发生较严重的破坏。

图4 场地地表加速度反应谱Fig. 4 Site ground acceleration response spectrum

3.1.2 场地地表峰值加速度分布情况

图5 和图6 分别给出了芦山地震及NGA 强震输入下的场地地表峰值加速度分布情况。由图5 和图6 可知，随着基岩峰值加速度的增大，地表峰值加速度呈增大趋势，但这种增大趋势在高震级地震动输入时趋于放缓。芦山地震与NGA 强震输入下，除记录 HYT-0201 和0789 外，场地地表峰值加速度变化区间较一致。同时，弱地震记录PXZ-0202 和3390 输入下得到的场地地表峰值加速度为12.7～25.5 Gal，强地震记录TQL-0202 和0994 输入下得到的场地地表峰值加速度为291.7～591.2 Gal。

图5 场地地表峰值加速度分布（芦山地震输入）Fig. 5 Site ground peak acceleration distribution（Lushan earthquake）

图6 场地地表峰值加速度分布（NGA 地震输入）Fig. 6 Site ground peak acceleration distribution（NGA）

3.1.3 场地地表峰值速度分布情况

图7 和图8 分别给出了芦山地震及NGA 强震输入下场地地表峰值速度分布情况。由图7 和图8 可知，随着基岩峰值速度的增大，地表峰值速度呈显著增大趋势，且随着地震强度的提高，地表峰值速度呈更大的增长趋势。与峰值加速度不同的是，场地峰值速度在芦山地震及NGA 强震输入下呈不同的变化情况，普遍来说，NGA 强震会引起较大的地表峰值速度，地震强度越大，这种影响越明显。同时，相比峰值加速度，峰值速度未呈现同样的增大趋势，如记录PXZ-0202 及0789 的峰值速度大于记录MNC-0202 及1021，这与峰值加速度的变化不同。峰值速度反映了场地的变形指标，与地下管道和隧道公路的破坏关系较大，因此在以往研究中单一进行场地峰值加速度分析易忽略峰值速度对场地变形的影响，从而影响分析结果的全面性。

图7 场地地表峰值速度分布（芦山地震输入）Fig. 7 Site ground peak velocity distribution（Lushan earthquake）

图8 场地地表峰值速度分布（NGA 地震输入）Fig. 8 Site ground peak velocity distribution（NGA）

3.2 峰值速度放大效应分析

为进一步分析峰值速度变化特征，引入场地速度放大系数kvi：

式中，kvi为编号为i的基岩地震输入速度放大系数；
vsi为编号为i的地震输入地表峰值速度；
vpi为编号为i的基岩地震输入峰值速度。

3.2.1 峰值速度放大系数随基岩输入的变化

图9 和图10 分别给出了芦山地震及NGA 强震输入下峰值速度放大系数随基岩输入的变化曲线。由图9和图10 可知，不同地震输入下峰值速度放大系数为1～1.8，除记录PXZ-0202、0789 和1091 外，其他地震输入下的地表峰值速度放大系数均较分散。观察基岩反应谱，发现记录PXZ-0202 及0789 均存在中低频成分丰富，而高频成分较少的现象，说明场地峰值速度放大效应变化与加速度反应谱频谱特性分布有一定的相关性。芦山地震下，高强度地震输入（如记录TQL13-0202）引起的场地放大效应程度更大，但此规律在NGA输入时并不明显。

图9 场地地表峰值速度放大系数分布（芦山地震输入）Fig. 9 Site ground peak velocity magnification factor distribution

图10 场地地表峰值速度放大系数分布（NGA 地震输入）Fig. 10 Site ground peak velocity magnification factor distribution

3.2.2 峰值速度放大系数随场地特征的变化

为比较芦山地震和NGA 强震输入下，雅安地区场地结构对峰值速度的影响效应，给出地表峰值速度放大系数随场地等效剪切波速和覆盖层厚度的变化曲线，如图11 所示，其中散点为各场地钻孔的峰值速度放大系数，直线为根据散点图绘制的拟合曲线。由图11 可知，除极个别情况外，峰值速度放大系数kv在不同强度基岩输入下，大致呈随等效剪切波速的增大而减小的趋势，而kv随覆盖层厚度的变化规律不明显，说明场地剪切波速对峰值速度的影响较大，软土场地更易引起场地效应的放大。由峰值速度放大系数kv随等效剪切波速的变化特征可知，芦山地震与NGA 强震引起的场地峰值速度随等效剪切波速的放大趋势不同，以记录PJD-0201 和1091 为例，可看出场地剪切波速对记录1091 的放大效应影响不明显。总体来说，雅安地区软土场地（等效剪切波速较小）在芦山地震作用下呈现出更大的放大效应，这种情况尤其在强震输入时更显著。

图11 场地峰值速度放大系数随场地特征变化情况Fig. 11 Site ground peak velocity magnification factor varies with site characteristics

本文采用等效线性化方法进行雅安地区场地反应计算，得到该地区各类参数的场地反应特征，重点开展峰值速度对场地的放大效应分析，得出以下结论：

（1）雅安地区钻孔反应谱谱形较一致，土层结构变化较小、土类分布均匀。该地区场地卓越周期约为0.2 s，本次地震引起的自振周期为0.2 s 的建筑物破坏相对严重。

（2）对比芦山地震及NGA 强震记录输入下的场地反应计算结果，在相同PGA 基岩输入的基础上，由于输入谱谱形的不同，导致场地加速度峰值及反应谱变化不同。

（3）不同地震动输入下峰值速度放大系数为1～1.8。随着基岩峰值加速度的增大，地表峰值加速度和峰值速度均呈增大趋势，同时峰值加速度在高震级地震动输入时增大趋势减缓。普遍来说，NGA 强震会引起较大的地表峰值速度，地震强度越大，影响越明显。

（4）峰值速度放大系数kv在 不同强度基岩输入下，大致随等效剪切波速的增大而减小，而kv随覆盖层厚度的变化规律不明显，说明场地剪切波速对峰值速度的影响较大，软土场地更易引起场地峰值速度的放大。同时，这种放大效应在芦山地震强震动作用下更为显著。

致谢 感谢四川省地震局提供的大量场地钻孔资料，感谢中国地震局工程力学研究所工程技术研究中心提供的强震动观测数据，感谢中国地震灾害防御中心张郁山研究员提供的建议和帮助。