鄂尔多斯西缘地区地震动高频衰减特性初步研究1

包金哲 王树波 张 帆

（内蒙古自治区地震局，呼和浩特 010010）

引言

地震动高频衰减在工程地震领域有着重要的应用，常被用于震源、地震动模拟和衰减关系研究等工作中。kappa 值是描述地震动高频衰减参数之一，多位学者对其进行了研究，如Sun 等（2013）选用汶川地震52 个主震加速度数据，得到kappa 值变化范围，认为kappa 值与震源和传播路径有关。Anderson 等（1984，1986，1991）发现在美国加州地区强震数据中，kappa 值与震中距表现出正相关关系。李文倩（2014）基于华北地区近2 000 条数据，利用Anderson 等（1984）的方法得到kappa 值，认为kappa 值与震级和震源深度无关，但离散性较大；朱百慧（2016）采用Anderson 等（1984）经典计算方法、固定应力降计算方法及宽带反演计算方法，以日本东北地区为研究区域，利用K-NET 和KiK-net 强震记录，对kappa 值与传播途径、震级及局部场地特征的关系进行分析，发现kappa 值与震中距具有一定线性关系，随着震中距的增大kappa 值呈增长趋势，kappa 值与震级的线性关系不明显，随着VS30的增大kappa 值呈对数下降趋势。

鄂尔多斯块体是古克拉通块体，其内部结构稳定，但其外部地震频发，尤其是鄂尔多斯块体西缘地区，该区域属于亚洲东部“大三角”地震构造区东边缘（高祥林等，2010），“大三角”地震构造沿华南、鄂尔多斯西边缘向北延伸至贝加尔湖附近，曾被称为中蒙大陆中轴构造带，研究该区域具有深远意义。鄂尔多斯西缘地区由内蒙古强震动台网、宁夏强震动台网属地监控，该区域共有14 个强震动台站，台站积累了真实可靠的地震数据，为研究提供了数据支撑。

本文以鄂尔多斯西缘地区（34°N−42°N，103°E−109°E）14 个固定强震台站（包括8 个内蒙古强震台站和6 个宁夏强震台站）2007−2020 年116 条强震记录（1≤MS≤8）为研究对象，利用Anderson 等（1984）提出的经典算法，通过MATLAB 软件计算kappa0 值，讨论kappa 值与震中距、kappa0 值与VS20、VS30的关系。

1 研究区域台站分布与地质构造特征

鄂尔多斯西缘地区地质条件复杂，北部临河断陷盆地西南方向与潮水−吉兰泰断陷盆地相连。南部有银川地堑、干盐池盆地，东部鄂尔多斯盆地是古老的克拉通地块，地块相对稳定。西北部地区为阿拉善块体，近年中小震频发。区域内有河套平原和宁夏平原，均为第四纪泥沙沉积覆盖的冲击平原。

临河断陷盆地自早第三纪晚期开始发育，新生界沉积厚度最大达10 000 m 以上，第四纪以来断陷活动依然保持了第三纪晚期的活动特征，周边山地继续抬升，断陷继承性沉降，第四纪厚度2 400 m。吉兰泰盆地最初发育于早白垩纪，但晚白垩纪-早第三纪古新世时期抬生为陆地，始新世又重新开始发育，直到现在。该断陷盆地内新生界厚度约2 000 m，其中第四纪最大厚度约400 m，沉降中心靠近巴彦乌拉山东麓断裂。银川地堑是由喜马拉雅造山运动时期构造活跃的贺兰山褶皱带与鄂尔多斯盆地相对上升形成的，其为本区域最大的拉张型活动盆地，北起石嘴山，南至青铜峡，长160 km，最宽55 km，地堑形成于渐新世，其后持续下陷，沉积了巨厚的新生代地层，总厚度达7 000 m，其中古近纪最大厚度3 600 m，新近系最大厚度2 500 m，第四纪最大厚度>1 609 m。干盐池盆地为青藏高原东北边界海原断裂带内最大的拉分盆地，沉积物粒度呈向上逐渐变细的趋势，湖泊沉积为主，干盐池盆地主要受南−西华山与唐家坡−干盐池断层的控制，最深处可能>550 m。第四纪以来，构造运动主要以大面积的周期节奏性震荡上升为主，在西吉县形成了剥蚀构造山地、构造剥蚀红岩丘陵，剥蚀堆积黄土丘陵，侵蚀堆积河谷平原。西吉县大面积被第四纪覆盖。河套平原是冲击平原，河套地区北部、南部断层证明该区域存在断层下陷，黄河流经低洼地区带来泥沙淤积，形成河套平原。宁夏平原也是冲击平原，由于贺兰山与鄂尔多斯高原的抬升，宁夏平原断裂下陷形成地堑带，后由于黄河携带大量泥沙堆积形成宁夏平原，可知宁夏平原是由地壳运动、断裂下陷、河流流水搬运堆积而成。

2 计算分析

2.1 数据选取

研究数据来自于内蒙古自治区地震局和宁夏自治区地震局提供的地震记录数据、钻孔资料、剪切波速测试资料、地质条件资料等，研究区域地震台站信息如表1 所示，台站分布及地质构造如图1 所示。

图1 台站分布及地质构造Fig. 1 The distribution of seismograph station and geologic structure

表1 研究区域地震台站信息Table 1 The research area seismic station information

2.2 计算方法

加速度反应谱法由Anderson 等（1984）提出，是目前工程地震学领域中使用最广泛的方法，本文采用该方法计算kappa 值，开展区域高频衰减特性研究。频率高于特定频率fE的单对数坐标系加速度傅立叶谱幅值衰减规律可表示为：

式中，A0为依赖于震源特征和传播距离的因子，k为谱衰减参数，f为地震动频率，fE为谱幅值对数值近似直线下降的起始点。

一般来说，当频率达fE后，谱幅值随着频率的增加呈指数下降。在线性-对数坐标系（即单对数坐标系）下，高频衰减参数k可用最简单的线性关系式近似得到：

kappa 值计算截止点为傅立叶谱与噪声谱重合的点或幅值平行起始点，利用最小二乘法进行线性拟合，得到的斜率即为kappa 值（图2）。

图2 Anderson 等（1984）经典计算方法计算的kappa 值Fig. 2 Anderson&Hough's the classical method calculates k

为计算和确定kappa 值与震中距R的关系，利用最小二乘法得到线性回归方程：

式中，R(s)为震中距，K0和KR为常数量，K0为由台站附近近地表地下地质构造引起的高频衰减参数值，KR为与区域相关的衰减参数值（Anderson 等，1984）kappa（N-S）和kappa（E-W）分别为kappa 南北向和东西向分量， K¯0表示强震台南北向、东西向kappa0均值。

2.3 数据处理

利用MATLAB 语言编辑计算处理全部程序，包括kappa 值、kappa 值与震中距关系、kappa0 值、kappa0 值与VS20、VS30关系计算程序。

（1）通过编程实现对强震观测记录数据的处理，包括基线调整、采用0.1～20 Hz 的带通滤波（白立新等，2019）、数据初筛及人工复筛等工作（崔建文等，2007；万永革，2012）。

（2）分析单个台站多次地震事件，利用计算机程序获得每次地震每个分量的kappa 值，进而利用每个分量的kappa 值结合震中距，采用最小二乘法进行拟合，得到kappa 值与震中距的线性回归方程（陶正如等，2015）。

2.3.1 kappa 值与震中距的关系

kappa 值与震中距的关系计算结果如表2、图3 所示，由表2、图3 可知，kappa 值与震中距存在线性回归关系；kappa 值在鄂尔多斯西缘大部分地区拟合优度值相对较高，为0.6～0.9，而西吉台站地区南北向、东西向拟合优度值分别为0.47、0.49，可能是由于台站记录地震数据较少、仪器噪声过高等原因导致kappa值相对离散，拟合优度较低。

图3 鄂尔多斯西缘地区14 个台站kappa 值与震中距拟合关系Fig. 3 The fitting relationship between kappa and epicentral distance of 14 stations in the western margin of Erdos region

表2 kappa 值与震中距拟合结果Table 2 The relationship between kappa and epicentral distance

2.3.2 计算结果不确定性分析

通过大量计算得到研究区域116 条地震记录地震动峰值加速度PGA 及震中距，震级与震中距关系分布如图4 所示，震级与峰值加速度PGA 关系分布如图5 所示。由图4、5 可知，震级为Ms2.0～Ms8.0，震中距为0～1 200 km，峰值加速度PGA 为0～180 cm/s2，多数集中在0～60 cm/s2。

图4 震级与震中距关系分布图Fig. 4 Distribution of earthquake magnitude and epicenter distance

图5 震级与峰值加速度PGA 关系分布图Fig. 5 Distribution of earthquake magnitude and PGA(peak acceleration value）

由于研究区域地震较少，覆盖层较厚，台站有效记录地震数据不足，且kappa 值与震中距拟合关系受远震（汶川地震等）影响，故拟合结果存在一定不确定性。针对此问题，剔除了震中距500 km 以上的地震。部分台站拟合结果偏差较大，以巴彦木仁、磴口台站为例进行分析，结果如表3、图6 所示。由表3、图6可知，剔除汶川地震等远震后，拟合优度值明显下降，巴彦木仁台南北向拟合优度值由0.869 11 降至0.618 6，东西向拟合优度值由0.928 0 降至0.708 1；磴口台南北向拟合优度值由0.897 1 降至0.244 7，东西向拟合优度值由0.942 5 降至0.606 4。综上所述，kappa 值与震中距线性拟合受远震的影响较大，因此计算kappa 值时应尽量考虑震中距为500 km 以内的地震数据。

图6 巴彦木仁、磴口台站kappa 值与震中距拟合关系Fig. 6 Fitting relationship between kappa and epicenter distance at Banyanmuren and Dengkou station

表3 巴彦木仁、磴口台站kappa 值与震中距拟合结果Table 3 Fitting results of kappa value and epicenter distance at Bayanmuren and Dengkou station

2.3.3 kappa0 值与剪切波速关系

研究区域有14 个强震动台站，内蒙古强震动台网8 个台站分别为：巴彦木仁台、巴彦浩特台、磴口台、公地台、吉兰泰台、临河台、沙海台、乌海台，宁夏强震动台网6 个台站分别为：宝丰台、干盐池台、灵武台、前进农场台、青铜峡台、西吉台。计算得到研究区域内14 个台站的kappa0 值，结合合站提供的钻孔资料，得到剪切波速VS20、VS30值，进而得出台站kappa0 值与剪切波速VS20拟合关系式为：

kappa0 值与剪切波速VS30拟合关系式为：

kappa0 值与剪切波速VS20、VS30拟合关系曲线如图7 所示。由式（6）、式（7）、图7 可知，kappa0 值与剪切波速VS20、VS30存在对数关系，拟合优度值约为0.4，拟合效果不理想，结合朱百慧（2016)的研究成果，笔者认为kappa0 值随着VS20、VS30的增大逐渐减至0（Ktenidou 等，2013）。

图7 kappa0 值与剪切波速VS20、VS30 拟合关系曲线Fig. 7 The fitting relationship between kappa0 and shear-wave（VS20、VS30）velocity

3 结论

本文利用Anderson 等（1984）提出的经典方法，通过计算得出以下结论：

（1）鄂尔多斯块体西缘地区强震台站kappa 值与震中距呈线性相关，且拟合优度值较高；

（2）kappa0 值为0.015 45～0.065 60 s；

（3）kappa0 值与VS20、VS30呈指数相关关系，离散性较大，拟合优度值约为0.4；

（4）随着VS20、VS30的增大，kappa0 值逐渐减至0。（5）本研究采用的地震数据较少，分析得到的kappa 值与震中距线性拟合受远震的影响较大，因此计算kappa 值时应尽量考虑震中距500 km 以内的地震数据。

致谢 本文在撰写过程中得到中国地震局工程力学研究所研究员陶正如、新疆地震局工程师李文倩、广东地震局高级工程师吴华灯、云南地震局研究员崔建文、防灾学院教授万永革、宁夏地震局高级工程师吕俊强及工程师李鸿廷的悉心帮助，在此向他们表示衷心感谢，同时感谢评审专家提出的宝贵意见。