基于随机振动模型的强地面运动数值模拟及其在广西灵山M63_4地震中的应用.pdf

书书书第 卷第 期（）年 月中 国 地 震 陈刚，李茂峰，李克华，等 基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用 中国地震，（）：基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用陈刚）李茂峰）李克华）唐勇）张清）马桂芳）申文豪，）姜文亮）中国南方电网有限责任公司超高压输电公司百色局，广西百色 ）广西壮族自治区地震局，南宁 ）应急管理部国家自然灾害防治研究院，北京 ）上海佘山地球物理国家野外科学观测研究站，上海 摘要 年 月 日广西灵山县平山镇东南罗阳山附近发生 地震，该地震是华南大陆自有地震记载以来发生的最大地震。本研究收集整理了灵山 地震的地质资料、活动断层探测最新成果等，选取适当的研究区域，利用随机振动有限断层模型计算区域内网格点的峰值加速度（）及峰值速度（）等参数，并且加入浅层横波速度结构 对地震动参数的影响，最终得到此次地震的地震动分布并分析了地震动特征。本研究将模拟结果与野外调查烈度数据和地震动衰减关系进行对比，结果显示模拟结果与调查烈度值和衰减关系在整体衰减特征、极震区的分布等方面均符合较好，模拟结果可为该地区未来地震危害性评估提供依据。本研究所使用的方法流程亦可应用于本地区地震烈度速报，为震后应急救援及决策指挥提供帮助。关键词：随机振动模型强地面运动模拟灵山 地震共轭破裂衰减关系 文章编号 （）中图分类号 文献标识码 收稿日期 修定日期 项目类别 中国地震 局“三 结 合”课 题（）、上 海 佘 山 地 球 物 理 国 家 野 外 科 学 观 测 研 究 站 开 放 基 金（）、国家自然科学基金（）共同资助 作者简介 陈刚，男，年生，高级工程师，主要从事强地震动模拟及电塔震害研究。：引言地震灾害的根源主要是断层突然运动产生地表位移，激发地震波造成地表强烈振动，即强地面运动。强地面运动可直接导致建、构筑物的倒塌、破坏，也可引起大量的崩塌、滑坡、滚石、地裂缝、砂土液化等，从而间接导致基础设施等遭受严重破坏。强地面运动研究的基础是强震记录数据，但是由于地震的随机性和不确定性，强震记录往往集中在某一些地区，地理环境和地质条件的不同使得许多数据的应用价值受到限制。为了在强震数据匮乏地区获得满足特定地理地质条件的地震数据，需要开展强地面运动模拟研究工作。广西是我国少震、弱震地区，其所在的华南地块是中国大陆一个相对稳定的构造单元，期陈刚等：基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用地震活动性较弱（黄强强等，）。虽然广西地震活动性水平较低，但是在历史上却发生过华南内陆地区自有地震记载以来最大的地震 年灵山 地震。年 月 日时 分，广西灵山县平山圩东罗阳山发生 强烈地震，震中位于灵山县东北约 的罗阳山西北麓一带，震中烈度达度或度强（李冰溯等，）。该地震造成 人死亡，人受伤，余间房屋倒塌破坏（陈恩民等，；李伟琦，；任镇寰等，）。对于缺乏仪器数据的历史地震，进行强地面运动模拟是评估该地区潜在地震危害性的重要手段。作为广西乃至华南内陆有记录以来最大的历史地震，开展灵山 地震的强地面运动模拟，进而补充本地区强震记录数据库，具有重要的科学价值和工程意义。但是由于历史久远和仪器记录匮乏等原因，自发震至 世纪初 余年的时间里，灵山地震的研究程度一直较低（李细光等，）。近年来，随着技术手段进步，相关学者对灵山地震开展了一系列研究，获得了一批突破性成果，已基本厘清此次地震的震级、震中位置、极震区等级范围和发震构造等关键信息（周本刚等，；何军等，；李细光等，、；李冰溯等，；李细光等，）。因此，针对灵山地震开展强地面运动模拟的条件已经成熟。本研究使用随机振动有限断层法对 年广西灵山 地震进行了强地面运动模拟，在收集地震地质资料的基础上，确定发震断层的走向、倾角、长度、宽度、震级、应力降等参数，选取适当的传播路径参数和场地效应参数，将上述参数输入随机有限断层模型进行强地面运动模拟，并将模拟结果与地震烈度野外调查数据和地震动衰减关系进行对比，验证模拟结果的可靠性。本研究最终获得 年广西灵山 地震区域及极震区峰值加速度（）、峰值速度（）及烈度的空间分布，相关结果可为工程抗震设计、政府决策和地震灾害风险评估提供依据。方法与原理 随机振动有限断层方法随机振动法是由 等（）和 （）于 世纪 年代提出的一种用随机振动来模拟高频地震波的方法。这种方法的本质是把远场地震动看作有限带宽、有限持时的高斯白噪声，并把影响地震动参数的一系列因素（包括震源、路径、场地）在频率域里表达为简单的乘积形式，然后通过窗函数调整噪声的傅里叶谱，最后通过傅里叶反变换得到地震动参数的时程曲线（，、）。对于一次大地震，合成地震动的有效方法是基于小震模拟。在随机有限断层模型中，断层破裂面被分成 个大小相等的矩形子源，每个子源即为一个点源，每个子源满足 平方衰减规律，通过设定震源模式和破裂传播速度，可以得到子源破裂的时间顺序，然后根据子源与场地的几何关系，计算每个子源对场地的影响。所有子源在观测点引起的地震动在时域中以适当的延迟时间叠加，便可获得场地的地震动时程（），即（）（）（）其中，、分别为断层沿长度和宽度划分的子断层数，和 分别代表沿断层走向和倾向方向的断层指标，为子源到场地的滞后时间，（）为子源利用随机点源方法得到的在观测点的地震动时程。针对有限断层模型采用部分子断层多次破裂的方式来保证地震矩守恒 中国地震 卷的不合理之处，等（）提出了“动力学拐角频率”的概念，认为子断层拐角频率随着破裂的传播而减小。每个子断层动力学拐角频率表达式为 （）（）（）（）其中，（）为第 个子断层的动力学拐角频率，为第 个子源被触发的时刻，（）为在时刻 已破裂子断层的累计数，为子断层的平均地震矩，为表示子断层辐射强度的一个常数。为了使断层总辐射能守恒，引入一个标度因子 。动力学拐角频率对子断层总辐射能的递减影响得到标度因子 的补偿，总辐射能保持不变，这样在保证辐射能守恒的前提下避免了子断层多次破裂这一非物理现象。至此，随机振动有限断层方法的理论框架得到最终完善。场地条件对地震动模拟结果的影响场地的浅层速度结构是引起地震场地效应的重要原因，直接影响地震灾害的严重程度。研究表明，基岩上的地震动参数值与考虑场地土层放大效应的参数值差别可达 倍（，），因此在强地面运动模拟中需要考虑浅层速度结构的影响。精细的浅层速度结构数据往往依赖于钻孔数据，而由于造价成本等因素，钻孔勘探数量很少，无法满足大规模区域性的强地面运动模拟的需求。为解决这一问题，等（）提出一种利用高精度地形数据计算地形坡度来求取浅层剪切波速 分布，进而得到近似场地校正放大因子的简单方法。这一方法已被成功应用于美国地质调查局（）所属的地震烈度速报系统（），并在实践中得到有效性的验证。的全球数据可以在 官方网站下载，不同学者则提出不同的计算场地放大因子的方法。（）计算了不同场地类型条件下短周期（）及中周期（）个基岩 和 输入档下的放大系数，这一结果无论在场地类型划分还是 或 分档方面均存在划分过粗的缺点。等（）在建立地震动衰减关系 （）时认为，对基岩 和 的放大作用可表达为 （）（）（）（）（）（）（）式中，和 为考虑了场地放大效应后的强地面运动参数。式（）避免了场地划分和地震动参数分档，适用于任何场地类型。通过式（）可将 的影响加入到强地面运动模拟的结果中，使得模拟结果能够充分反映浅层速度结构的影响。图 展示了本研究区内的地形及 分布情况，可以看到在平原、盆地等沉积层较厚的区域，数值较低，而在丘陵地带、山区及基岩裸露地区，数值较高。模型参数 发震构造发震断层的几何形态（包括走向、倾角、长度和宽度等）是随机振动有限断层模型的基本输入参数。灵山 地震所在的华南地区由于历史强震较少，人类改造活动频繁且降水较多，地震遗址保存非常困难，因此很长一段时间内关于此次地震的发震断层存在争议。陈国达（）最早研究灵山 地震时，认为灵山地震发震断层为 走向的灵山大容山断裂，而茹锦文等（）则认为发震断层应归属于 向的昆仑关断裂南延部分。前人对该地震发 期陈刚等：基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用图 研究区内地形（）及 分布（）震断层的研究大多从极震区震害和地裂缝的空间分布情况进行推测（陈国达，；茹锦文等，；李伟琦，；任镇寰等，），缺乏直接的地质学证据。年研究人员首次在灵山断裂北段发现灵山地震 条地表破裂带遗迹（李细光等，、），地表破裂带沿线开挖的探槽清晰揭露出断层对上覆土层的错动，这为灵山断裂北段晚更新世以来的活动提供了直接证据。探槽显示的最新一次地震事件距今约 年，与灵山地震发震时间吻合，证明了灵山断裂北段即为 年灵山 地震的发震断裂。另外，李伟琦（）和任镇寰等（）根据极震区的等震线形状和低烈度区的长轴方向，推测灵山地震可能为 向断裂与 向断裂共轭破裂的结果，其中 向构造为控震构造。李冰溯等（）综合极震区震害及地质学研究成果，认为灵山地震主发震构造为 走向的灵山断裂，罗阳山南麓的泗州断裂同时发生共轭破裂，造成局部烈度增强。因此，本文在开展模拟时将发震断层设定为 条，即灵山断裂北段和泗州断裂（图 ），条断裂的几何参数见表 。断层滑动模型断层滑动模型即断层面上的滑动量静态分布是进行强地面运动模拟的必要参数，一般通过震后地表形变数据或地震波数据反演获得（申文豪等，）。对于灵山 地震这类没有观测数据的历史地震，已很难通过反演方式获得滑动模型。除了反演模型之外，地震学者还根据地震波观测总结出断层滑动模型的空间分布规律，提出了通过正演方式生成断层滑动分布的方法（王海云，）。本文采用应用广泛的 平方模型正演获得断层面滑动分布。等（）提出的 平方模型二维波数谱为（，）珚 ()()()槡（，）（）中国地震 卷注：断层数据来自国家地震科学数据中心活断层数据分中心。图 灵山 地震震中区域地震构造简图其中，和 分别为断层破裂面的长度和宽度，为平均滑移，（，）为相位谱，和 分别为走向方向滑移分布和下倾方向滑移分布的空间拐角波数。等（）给出了震级 与空间拐角波数经验关系公式，即 （）（）（）本研究采取以下策略获取断层滑动模型：首先根据断层几何参数（表 ），由震级与断层破裂长度经验关系（等，）分别得到灵山断裂北段和泗州断裂的发震震级和平均滑动位移 ，再由公式（）得到 个断层滑动模型的空间拐角波数，最后根据公式（）分别生成灵山断裂北段和泗州断裂的滑动模型（图 ）。震源及路径传播参数根据现有文献及地质资料基本可以确定 年灵山 地震发震断层的破裂长度、走向、倾角、埋深、震中位置等几何参数（表 ）。除此之外，为描述断层破裂和传播特征，还需要确定地震矩、应力降、值、品质因子 模型、几何衰减模型等参数。龙政强等 ：期陈刚等：基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用表 模型计算输入参数参数名灵山断裂泗州断裂数据来源断层方位走向 倾角 走向 倾角 张沛全等（）矩震级 等（）应力降 龙政强等（）主断层尺度 张沛全等（），周本刚等（）地壳剪切波速度 梁飞等（）破裂速度 何欣娟等（）脉冲面积百分比 何欣娟等（）地壳平均密度 梁飞等（）品质因子 龙政强等（）值 黄强强等（）几何衰减模型 ，（），（），杨晓瑜等（）图 正演获得的灵山断裂（）和泗州断裂（）断层面滑动分布模型 中国地震 卷（）选取广西数字地震台网记录到的 年广西主要断裂带 地震事件，采用 方法反演得到介质品质因子 。根据 等（），值 可以由经验关系式确定，即 （）（）其中，为上地壳 深度内平均横波速度。根据广西中西部地区背景噪声成像结果，灵山地区 约为 （黄强强等，），因此 值计算结果为 。杨晓瑜等（）利用中国国家地震台网 个固定地震台站记录的远震波形资料，通过 波接收函数分析估算了中国华南地区的地壳厚度，其结果显示灵山地震震源区地壳厚度约为 ，根据三段式几何衰减模型的定义（等，），该地区几何衰减模型可以表示为（），（），（）其中，为震中距。龙政强等（）研究显示，广西地区 及其以上的高应力降地震主要分布在桂东南，且历史上发生过 级以上地震的震源区及其附近仍具有高应力背景。年 月 日发生的灵山 地震其应力降为 ，该地震位于 年灵山 地震和 年 月 日灵山 地震震中附近，这说明灵山地震区具有高应力背景，本研究将应力降的值取为 。结果分析与讨论灵山 地震虽然缺乏仪器记录，但具有丰富的烈度调查数据。本研究针对宏观区域（调查烈度覆盖区域及其扩展区域）和极震区（度和度区）分别划分不同的计算区域，并设置不同尺度的计算网格，利用随机振动有限断层模型计算每个网格点的 和 ，通过经验转换关系（丁宝荣等，）分别转化为烈度，最后取两种烈度的加权平均为最终烈度值，并将其与野外调查结果进行比较。丁宝荣等（）给出的 、与烈度 转换关系为 （）（）（）宏观烈度模拟宏观烈度模拟区域划定在 、范围内，区域网格大小为 ，经度及纬度方向网格精度分别为 和 。为更好地展示强地面运动参数衰减趋势，绘制了 （）和 （）空间分布图，如图 （）、（）所示。从图中可以看到模拟 的 强 地 震 动 参 数 具 有 以 下 特 征：最 大 值 在 左 右，最 大 值 达 ，转化为烈度均已超过度；和 均随距断层距离增加而迅速衰减，且 的衰减趋势更快；和 均受到场地条件的影响，在平原等沉积层较厚区域存在不同程度的放大效应。图 （）展示了模拟烈度与野外调查烈度（陈恩民等，）的对比，图中白色实线即为野外调查烈度等值线。可以看到，模拟获得的地震烈度在空间分布的大致走向、烈度等级范围与野外调查结果基本一致。模拟结果显示断层附近区域极震区的最大烈度已经超过 期陈刚等：基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用图 基于随机振动模型模拟得到的灵山 地震宏观区域强地震动参数分布度，其中灵山县城位于度区，南宁市区位于度区，玉林和钦州均位于度区，均与调查结果相同。另外需要注意的是，由于采样点的有限性，野外调查烈度图中的烈度等值线大多是经验性地圈划为圆形或者椭圆形，而本研究的数值模拟结果由于考虑了 的作用，烈度分布的形态明显受到场地条件的影响。其主要体现在一些沉积层较厚的冲积平原和盆地出现高烈度异常区，这些地区的 普遍较小，从而产生了局部场地放大效应，这可能是造成高烈度异常的主要原因。极震区模拟为更好地研究双断层共轭破裂造成的强地面运动特点，本研究针对极震区开展了网格尺度更精细的数值模拟，模拟区域划为 、，区域网格大小为 ，经度及纬度方向网格精度为 。极震区地震动模拟结果如图 所示，图中白色实线为野外调查度和度区等值线，黑色实线为模拟烈度度和度区等值线。国家地震局全国地震烈度区划编图组（）汇编的 中国地震等烈度线图集 中，灵山 地 中国地震 卷震的震中烈度最大为度，并且形状为一长轴方向呈 向的椭圆形。中国地震目录（顾功叙，）和 中国近代地震目录（中国地震局震害防御司，）均采用这一结果。此图 基于随机振动模型模拟得到的灵山 地震极震区强地震动参数分布后，陈恩民等（）、李伟琦（）以及任镇寰等（）对极震区烈度进行了重新调查和评定，结果均显示极震区的烈度分布有 和 两个优势方向，呈“”型分布。由图 可以看到，无论是 、分布，还是烈度分布，本研究的模拟结果较好地反映了极震区烈度分布的这种“”型特征，尤其是在度区。这一结果从模拟的角度验证了模型参数的合理性，也证实了灵山地震极震区分布是由灵山断裂及泗州断裂共轭破裂造成的。两断层发生共轭破裂的控制因素可能是罗阳山地区最大有效力矩准则约束下的以 向为主压应力方向的力学系统（张沛全等，）。图 （）还显示出与调查烈度不同的是，模拟 期陈刚等：基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用烈度度区在北部及西南部有两处明显的烈度异常区，而这两处烈度异常区均为河流冲积平原，北部为郁江平原区，西南部为钦江平原区。这两处平原区沉积层较厚，水平较低，从而导致明显的场地放大效应。注：（）和（）分别为宏观区域模拟 和 与 对比；（）和（）分别为极震区内模拟 和 与 对比。图 模拟得到的 、与衰减关系 对比 衰减关系对比除了与调查烈度进行空间分布对比之外，本研究还将宏观区域和极震区模拟得到的地震动参数与衰减关系 （等，）进行对比。图 （）、（）分别展示了宏观区域模拟 、与 的衰减趋势，图 （）、（）分别展示了极震区内模拟 、与 的衰减趋势，图中黑色实线为 衰减关系曲线，黑色虚线为 置信区间曲线。从图 可以看到，模拟结果与 结果的衰减趋势大体一致，在极震区内模拟结果与 曲线符合程度较高，但在震中距大于 时，地震动模拟结果整体略高于 预测值。图 展示了本研究模拟的烈度与王继等（）针对本地区烈度衰减公式的对比，其中蓝色实线为长轴衰减公式，蓝色虚线为短轴衰减公式。从图 可以看到，模拟烈度与理论公式在不同震中距尺度上具有相似的衰减特征。在震中距小于 的范围内，基于烈度衰减公式预测的最大烈度值为度强，显然低于模拟烈度和实际调查烈度。造成这一差异原因可能是王继等（）所使用的地震事件中超过 级的地震样本数较少，从而可能会低估大震级事件极震区内的烈度水平。在震中距大于 处，模拟烈度则低于基于烈度衰减公式 中国地震 卷预测的烈度水平。考虑到在远场区域地震动主要以面波形式出现，而品质因子 对其有着重要影响，在本研究过程中只采用了单一的 衰减模型。事实上，许多研究已经表明 值是依赖于震中距变化的（等，），因此，在今后工作中需要建立随距离变化的 值模型来改善模拟结果。图 模拟得到的烈度值与华南地区烈度衰减关系（王继等，）对比 结论与展望本研究收集整理了 年灵山 地震的地质资料、活动断层探测最新成果等，利用随机振动有限断层模型计算了宏观区域及极震区内网格点的峰值加速度及峰值速度等参数，并且加入了浅层横波速度结构对模拟结果的影响，最终得到此次地震的地震动分布，分析了地震动特征。本研究将模拟结果与历史地震烈度调查数据及地震动衰减关系进行对比，结果显示模拟结果与调查烈度值在整体特征、极震区的分布等方面均符合较好。由于广西地区中强地震少发，缺乏足够的强地面运动记录回归拟合得到本地区地震动衰减关系，而模拟得到的强震数据可以补充强震数据的不足，为建立广西地区地震动衰减关系曲线提供有益补充。另外，由于随机振动有限断层模型未考虑地震波在介质中的传播过程，不需要像复合震源模型、有限差分法或者谱元法等方法耗费大量时间用于地质体网格剖分和格林函数计算，具备时效性强的突出特点，这一点正是地震烈度速报所需要的。因此，未来随机振动有限断层模型可以应用于本地区地震烈度速报，为震后灾情快速评估和应急救援提供有效帮助。参考文献陈恩民，黄咏茵 华南十九次强震暨南海北部陆缘地震带概述 华南地震，（）：陈国达 广东灵山地震志 两广地质调查所特刊，（）：丁宝荣，孙景江，杜轲，等 地震烈度与峰值加速度、峰值速度相关性研究 地震工程与工程振动，（）：顾功叙 中国地震目录：公元前 公元 年 北京：科学出版社，国家地震局全国地震烈度区划编图组 中国地震等烈度线图集 北京：地震出版社，何军，刘怀庆，黎清华，等 广西防城灵山断裂带北东支灵山段活动性初探 华南地质与矿产，（）：何欣娟，潘华 年云南漾濞 地震的强地面运动模拟 地震地质，（）：期陈刚等：基于随机振动模型的强地面运动数值模拟及其在广西灵山 地震中的应用黄强 强，沈 旭 章，王 文 天，等 广 西 北 流 地 震 震 源 区 短 周 期 密 集 台 阵 背 景 噪 声 成 像地 球 物 理 学 报，（）：李冰溯，李细光，潘黎黎 年广西灵山 地震参数讨论 地震学报，（）：李伟琦 广西新构造分区特征及其与地震的关系 华南地震，（）：李细光，李冰溯，潘黎黎，等 广西灵山 年 级地震地表破裂带新发现 地震地质，（）：李细光，潘黎黎，李冰溯，等 广西灵山 年 级地震地表破裂类型与位错特征 地震地质，（）：李细光，潘黎黎，李冰溯，等 广西灵山断裂北段古地震事件初步研究 地学前缘，（）：梁飞，阎春恒，向巍，等 广西地区不同类型地震震源参数特征 华北地震科学，（增刊）：龙政强，谢夜玉，韦东，等 广西主要断裂带中小地震应力降特征分析 山西地震，（）：，任镇寰，杨廉法，邓业权 年广西灵山 级地震极震区震害和地震影响场的研究 中国地震，（）：茹锦文，黄瑞照 灵山地震发震机制及邻区区域稳定性评价 华南地震，（）：申文豪，李永生，焦其松，等 联合强震记录和 结果的四川九寨沟 级地震震源滑动分布反演及其地震学应用 地球物理学报，（）：王海云 近场强地震动预测的有限断层震源模型 博士学位论文 哈尔滨：中国地震局工程力学研究所王继，俞言祥 华中、华南中强地震区地震烈度衰减关系研究 震灾防御技术，（）：杨晓瑜，李永华 中国华南地区地壳厚度与波速比分布特征及其地质意义 地球物理学报，（）：张沛全，李冰溯 最大有效力矩准则约束下的广西罗阳山地震、地质与地貌效应 地质力学学报，（）：中国地震局震害防御司 中国近代地震目录（公元 年，）北京：中国科学技术出版社，周本刚，杨晓平，杜龙 广西防城灵山断裂带活动性分段与潜在震源区划分研究 震灾防御技术，（）：，（）：，（）：，（）：，：（），（）：，（）：，（）：，：，（）：，（）：，（）：，（）：中国地震 卷 ），），），），），），），），），），），：；