俯冲带板缘与板内地震长周期地震动参数预测模型研究1

肖 亮 俞言祥

（中国地震局地球物理研究所，北京 100081）

引言

地震动参数预测模型（常称为衰减关系）的建立是地震动参数区划图编制的关键技术环节之一。我国陆域地震区划图编制已经历数代，逐渐完善发展。然而，对于广袤的海域，地震动参数区划仍需开展进一步研究工作。

目前，国内外多位研究者利用强震记录、数值模拟结果建立了全球性或区域性俯冲带地震动预测模型，如Atkinson 等（2003）利用全世界范围内俯冲带强震记录建立了地震动预测模型（板缘、板内），并通过实际记录对比分析了该模型在Cascadia 区域（华盛顿、北加拿大区域）与其他区域的适用性；Arango 等（2012）对比了世界范围内多组俯冲带地震动预测模型在智利、中美洲俯冲带区域的适用性，并探讨了区域性差异及其在地震危险性分析中的影响；胡进军等(2013，2014)和郝彦春（2015）利用强震资料及数值模拟结果研究了俯冲带地震动衰减特性、局部场地效应等，并探讨了潜在地震对我国海域开发和建设的影响；Zhao 等（2006，2016）利用日本、台湾等地俯冲带强震记录，采用H/V进行场地分类，建立了板内、板缘地震地震动预测模型。但以上研究针对的反应谱频段多限于周期为5～6 s，对长周期地震动的估计存在一定不足。

陆域地震动参数区划主要针对面大量广的民用建筑与一般性重大工程，而海域几乎不涉及民居，其中岛礁建筑、跨海大桥等长周期结构较多。我国海域地震环境存在特殊性，大部分海域远离俯冲带，除短周期地震动参数主要受陆域或近海海域高震级潜在震源区控制外，长周期地震动参数受远场俯冲带深大地震的影响不容忽视。考虑以上原因，本文进行俯冲带地震动参数预测模型建立时，重点关注了对反应谱长周期段的估计。

1 总体思路

近年来，我国强震动观测数据虽得到快速积累，但仍缺乏海域强震资料。因此，根据实际情况，广泛收集与我国海域相似区域（日本东部俯冲带、琉球-冲绳俯冲带）地震记录，即从日本强震台网（K-net 和Kiknet）收集了该区域基岩强震加速度资料，从国际地震学与地球内部物理学联合会（IASPEI）全球数字地震台网收集了该区域宽频带速度测震记录。建立地震动参数预测模型时，综合强震记录近场高频优势与测震记录远场低频优势，可较好地估计反应谱长周期段。

考虑与陆域区划和我国实际工程设计的衔接，建立俯冲带地震动参数预测模型时，使用阻尼比为5%的基岩水平向加速度反应谱（周期至10.0 s），地震震级采用面波震级，考虑断层尺度和震源深度的影响，距离参数采用断层距Rrup（场点至断层破裂面的最短距离）。

2 地震资料与分类

2.1 俯冲带强震资料

除搜集琉球-冲绳俯冲带地震记录外，为弥补高震级地震的缺失，挑选具有较多高质量记录且地震活动性较高的日本东部俯冲带地震记录。为保证数据的可靠性，选用近年强震记录，数据获取时间为2000 年1 月1 日至2020 年7 月1 日。获取强震记录后，对其进行基线校正和截止频率0.025 Hz 的高通滤波，并计算阻尼比为5%的加速度反应谱。采用人工检查的方式剔除信噪比较低或有明显错误的记录。加速度记录主要用于反应谱周期至1.0 s 的短周期段预测模型的确定。

从全球CMT 地震目录获取强震资料地震参数，包括震级、断层性质、震源深度等信息。采用面波震级作为统一震级标度，对于没有面波震级却有矩震级的地震，采用刘瑞丰（2015）震级转换公式计算面波震级。采用断层距作为距离参数，对于可在K-net、Kik-net 网站查询到的有限断层模型地震，直接计算断层距；对于无法查询到有限断层模型的地震（大部分地震震级较低），采用震源距（场点至震源的距离）作为断层距，当距离较远时（大多数情况），二者差别较小。

采用Vs30（地表以下30 m 内折算平均剪切波速）作为场地分类参考依据，采用Vs30 超过500 m/s 的记录作为基岩场地记录。

2.2 俯冲带测震资料

将宽频带数字测震记录作为速度记录，主要来自于IASPEI 全球数字地震台网。该台网台站布设历史较长，记录器中速度平台型宽频带高增益BH 信道具有较多记录，能够提供高质量的俯冲带远场记录，适用于长周期地震动研究。

挑选冲绳俯冲带与日本东部俯冲带宽频带速度记录，剔除井下或具有一定高程的记录，主要关注自由地表。获取记录后，与强震资料处理程序相同，完成基线校正和截止频率0.025 Hz 的高通滤波，对其进行微分得到加速度时程。测震记录具有高增益特性，常存在近场限幅的情况。为保证资料的可靠性，采用人工检查的方式对信噪比较低、限幅或显著有误的记录进行剔除，然后计算阻尼比为5%的加速度反应谱。该部分速度记录主要用于确定反应谱周期超过1.0 s 的长周期段预测模型。

震级与距离参数与强震资料的处理相同，对于场地分类，由于台站放置在基岩或基座上，可不考虑土层的影响。

2.3 板缘与板内地震分类

通过震源深度与震源机制确定板缘、板内地震分类，具体规则如下：

（1）对于震源深度为20～30 km 浅壳地震，震源机制为逆断层，考虑俯冲带上缘与地壳之间无明显界限，将这部分地震归为板缘地震。

（2）对于震源深度为30～50 km 的地震，归为板缘地震。考虑震源深度为50～70 km 处，中深源和深源地震不存在明显界限，一般认为板缘地震震源机制为逆断层，因此，将震源深度为50～70 km 且震源机制为逆或走滑断层的地震归为板缘地震。

（3）对于震源深度为50～70 km 且震源机制为正断层的地震，归为板内地震。将震源深度超过70 km的地震也归为板内地震。

为完成预测模型的回归，选择面波震级Ms≥5.0、Rrup≤500 km 的水平向地震记录进入板缘地震回归数据集；考虑板内地震震源深度较深，适当扩大距离范围，因此选取面波震级Ms≥5.0、Rrup≤600 km 的水平向地震记录进入板内地震回归数据集。为增加统计样本量，同一台站2 条水平向地震记录视为独立记录，强震资料、测震资料地震目录如表1～4 所示，地震震中分布如图1 所示。

图1 地震事件震中分布图Fig. 1 Map of earthquakes epicenter location

表1 强震资料地震目录（板缘地震）Table 1 Catalog of strong earthquakes records (Interface earthquakes)

表2 宽频带资料地震目录（板缘地震）Table 2 Catalog of broad-band records (Interface earthquakes)

表3 强震资料地震目录（板内地震）Table 3 Catalog of strong earthquakes records (Intraslab earthquakes)

表4 测震资料地震目录（板内地震）Table 4 Catalog of broad-band records (Intraslab earthquakes)

板缘地震使用资料震级-断层距分布关系如图2 所示，由图2 可知，震级为6 级以下时，由于断层距较远，导致加速度幅值过小或信噪比较低，多条中小震级地震记录被剔除，对回归产生一定影响。板内地震使用资料震级-断层距分布关系如图3 所示，由图3 可知，板内地震震级上限为7 级，对高震级结果的约束性相对较弱。总体来看，宽频带资料较好地弥补了强震记录在远场的缺失，有利于长周期预测模型的建立。

图2 板缘地震记录震级-断层距分布关系图Fig. 2 Distribution of magnitude-distance of interface earthquakes

图3 板内地震记录震级-断层距分布关系图Fig. 3 Distribution of magnitude-distance of intraslab earthquakes

3 模型与回归方法

借鉴《中国地震动参数区划图》（GB 18306-2015）分段线性模型（高孟潭，2015）建立板缘地震动预测模型，对于板内地震，增加震源深度相关项，回归模型可统一写为：

式中，Y为反应谱值；M为面波震级；R为断层距；T为反应谱周期，当震源深度＜70 km 时，h=0，当震源深度为70～120 km 时，h=震源深度－70 km，当震源深度＞120 km 时，h=50 km；ε为具有标准差σ的随机量，均值为0；其余参数由回归拟合得到。由于板缘地震震源深度＜70 km，震源深度项对于板缘地震不起作用；对于板内地震，震源深度项可在一定程度上反映震级量规函数对震源深度的补偿，由于对深源地震的研究较少，人为设定线性关系在120 km 震源深度处不再增长。

沿用与五代图相配套的分段回归方法（肖亮等，2010）进行预测模型回归，对每个周期点反应谱值进行回归，并对整体谱型进行适当光滑，得到完整的反应谱预测模型。由于俯冲带地震一般远离大陆，近场记录天然性较少，对直接回归近场饱和因子DeEM造成较大困难。考虑DeEM主要反映断层破裂尺度信息，仅对近场区域有影响，因此，直接借用陆域回归结果，即D=0.956，E=0.462，其余参数由回归确定。

4 回归结果

对于阻尼比为5%的基岩水平向加速度反应谱，选取10.0 s 范围内基本等对数间隔的30 个点。强震记录用于短周期段（1.0 s 前）的回归，宽频带速度记录用于长周期段（3.0 s 后）的回归，对于1.5～3.0 s 周期段，采用强震资料与宽频带速度资料混用的方式进行衔接，最终得到周期至10.0 s 的长周期反应谱预测模型。板缘地震和板内地震预测模型相关参数分别如表5、6 所示，相应的反应谱曲线（由下至上Ms=5、6、7、8 级）如图4 所示，周期6.0 s 反应谱值随断层距的衰减曲线如图5 所示（由下至上Ms=5、6、7、8 级）。

图4 基岩水平向加速度反应谱曲线Fig. 4 Spectral acceleration curves for horizontal rock ground motion

图5 周期6.0 s 反应谱值衰减曲线Fig. 5 Attenuation curves for spectrum of period 6.0 s

表5 板缘地震基岩水平向加速度反应谱预测模型相关参数Table 5 Coefficients of horizontal rock spectral acceleration prediction model for interface earthquakes

表6 板内地震基岩水平向加速度反应谱预测模型的参数Table 6 Coefficients of horizontal rock spectral acceleration prediction model for intraslab earthquakes

5 讨论

为观察回归结果的可靠性，给出板缘、板内地震关键周期点（周期6.0 s）反应谱回归残差随震级、断层距的变化关系，如图6、7 所示。由图6，7 可知，随着断层距和震级的变化，回归残差在零值附近水平分布，无明显趋势，说明回归结果真实可靠，且本模型中离散标准差可视为不随断层距和震级变化的常数。为进行合理分析，选择目前国际上常用的俯冲带衰减关系与本文结果进行比较，包括Abrahamson 等（2016）采用全球俯冲带强震记录统计结果和Zhao 等（2006）采用日本俯冲带强震记录统计结果。震级6.5 级、断层距50 km（板缘地震，震源深度30 km），震级7.5 级、断层距200 km（板缘地震，震源深度50 km），震级6.5 级、断层距100 km（板内地震，震源深度70 km），震级7.5 级、断层距200 km（板内地震，震源深度70 km）处反应谱的对比如图8 所示。由图8 可知，在震源深度较浅的板缘地震情况下，对于近场高频部分，本文预测模型得到的结果较高，具有一定保守性，对弈中远场部分，本文预测模型得到的结果居中；在外推较深的板内地震情况下，3 组预测模型得到的结果较接近，对于远场部分，本文预测模型得到的长周期略低。总体来看，3 组预测模型反应谱曲线基本相当，在主要震级-断层距段结果较接近，可在一定程度上验证本文预测模型的可靠性。靠性。Abrahamson 等（2016）的研究使用全球资料，Zhao 等（2006）的研究可用周期截至5.0 s，考虑区域性衰减特性与长周期可用性等，本文预测模型更适用于我国海域及长周期地震动的估计，适用震级为5.0～8.0 级，断层距为0～600 km。

图6 回归残差随震级的变化关系图Fig. 6 Comparison of residuals in relation with magnitude

图7 回归残差随断层距的变化关系图Fig. 7 Comparison of residuals in relation with distance

图8 基岩水平向加速度反应谱曲线对比结果Fig. 8 Comparison of spectral acceleration curves for horizontal rock ground motion

致谢 本文使用日本K-net、Kik-net 地震台网、IASPEI 全球地震台网地震资料，研究和成文过程中得到李小军教授、周正华教授及审稿专家的热心指导与建议，在此一并表示感谢！