2022 年青海门源MS6.9 地震后冷龙岭断裂未来强震的水平位错量评估1

吴 果 孙浩越 吕丽星 冉洪流 周 庆 周介元

（中国地震局地质研究所, 北京 100029）

引言

2022 年1 月8 日青海省海北州门源县发生MS6.9 地震（图1），震中位于37.77°N，101.26°E，震源深度10 km。发震断裂为青藏高原北东缘祁连-海原断裂中西段的冷龙岭断裂，震中位于断裂南侧3～4 km（Yuan 等，2013；李智敏等，2022）。本次地震产生了以左旋走滑为主的同震地表破裂带。破裂带分为南、北支，分别沿托莱山断裂东端和冷龙岭断裂西端展布，总长约27 km（潘家伟等，2022）。南支破裂带最大左行水平位错量约85 cm，北支破裂带最大左行水平位错量约3.7 m。

门源地震虽未造成人员伤亡，但使冷龙岭断裂西端产生破裂带，破裂带与兰新铁路相交（图1），对铁路大桥和隧道造成了严重破坏,其中大梁隧道被本次地震破裂带直接错断，受损最严重（图2）。在距地震破裂带稍远处，地震对隧道的破坏以外墙和拱顶开裂、剥落、掉块及变形为主（图2（a））。在地震破裂带与隧道交汇处，隧道被断裂左旋水平错动约2 m（图2（b）），这一水平错动使隧道在短期内难以修复，造成巨大的经济损失。拟定隧道修复方案时，需对冷龙岭断裂未来强震的水平位错量进行评估，为隧道抗震设计提供参考。

图1 门源MS6.9 地震构造背景Fig. 1 Tectonic background map of Menyuan MS6.9 earthquake

图2 大梁隧道破坏情况Fig. 2 Photos of Daliang tunnel damaged by Menyuan MS6.9 earthquake

活断层同震位错量是铁路、隧道、油气管线等长线路工程抗断设计中的关键性参数（曹毅渊等，2019）。目前，我国工程地震界普遍采用确定性方法评估活断层同震位错量，即通过活断层最大潜在地震震级和震级与最大位错量的经验关系估算对应的最大位错量。然而目前人类对地震的认识水平有限，发震断层震级分布、破裂长度、位错分布等均存在较大的不确定性。因此，地震学家尝试借鉴经典的概率地震危险性分析理论，采用概率性方法评价断层未来一段时间内可能的位错分布情况，即概率断层位错危险性分析（Probabilistic Fault Displacement Hazard Analysis, PFDHA）（Youngs 等，2003）。概率断层位错危险性分析理论和方法已被国内外学者广泛应用于断层同震位错量评估中（冉洪流等，2004a，2004b；Chen 等，2011；Petersen 等，2011；荆旭，2019；Valentini 等，2021；吴果等，2022）。

本文基于近年来冷龙岭断裂的最新地震活动性研究成果，选取适当的平均地震复发间隔、最大震级等参数，进而采用确定性方法和概率断层位错危险性分析方法对冷龙岭断裂未来强震的水平位错量进行评估。本研究可为大梁隧道修复工程提供参考，为其他长线路工程抗断设计提供技术支撑。

1 冷龙岭断裂地震活动参数

冷龙岭断裂在青藏高原东北缘的构造变形中起着重要的转换调节作用（Gaudemer 等，1995）。断裂西端与托莱山断裂呈左阶斜列分布，东端与天桥沟-黄羊川断裂及金强河断裂相接（李智敏等，2022）。断裂晚第四纪主要表现为左旋走滑运动，地震活动强烈（郭鹏等，2017）。

郭鹏（2019）通过探槽开挖和放射性碳定年技术揭示了冷龙岭断裂最近发生的6 次古地震事件和形成时代，由新到老的年龄分别为438-350 a BP 以来、2 951-1 155 a BP、4 016-3 609 a BP、5 325-4 476 a BP、7 284-6 690 a BP 和8 483-7 989 a BP。6 次古地震事件的平均复发间隔为（1 640±570） a，相应变异系数为0.34。

冷龙岭断裂在我国第五代地震动参数区划图中位于六盘山-祁连山地震带的冷龙岭潜在震源区内（潘华等，2013）。冷龙岭潜在震源区震级上限为M7.5，所在地震带的b值为0.75，本文采用以上参数开展相关研究。

2 水平位错量评估

2.1 确定性方法与评估结果

确定性方法是将断裂最大潜在地震震级代入震级与最大位错量的经验关系式，从而估算对应的最大位错量Dm。对于本研究，Dm特指最大水平位错量。已有学者基于不同数据集拟合了不同的经验关系式。考虑不确定性因素影响，本文同时引用邓起东等（1992）、冉洪流（2011）、Wells 等（1994）提出的关系式。

邓起东等（1992）建立的适用于我国青藏地区走滑断层的震级与最大位错量关系为：

式中，MW为矩震级。

由于Wells 等（1994）的关系式是基于矩震级，因此需将我国工程地震界常用的面波震级转换为矩震级（张力方等，2013；吴果等，2014；谢卓娟等，2020）。本文引用Cheng 等（2017）拟合的适用于我国大陆地区的震级转换关系，当MS≥7.0 时：

式中，MS为面波震级。

与冷龙岭断裂最大潜在地震震级对应的矩震级为Mw7.18，最大水平位错量计算结果如表1 所示。

表1 确定性方法给出的水平位错量评估结果Table 1 Assessment results of horizontal displacement by deterministic method

2.2 概率性方法与评估结果

2.2.1 基本原理

评估断裂同震位错量的过程中面临多个关键环节的不确定性，本文给出的震级与最大位错量的经验关系式均带有一定标准差，即存在不确定性。经验关系式的不确定性影响计算结果（吴果等，2022），而该影响在确定性方法中往往被忽视。此外，确定性方法未考虑断裂震级分布，默认采用最大潜在地震震级估算对应的最大位错量，这种方法可能会造成结果过于保守（孙建宝，2002）。事实上，大震产生的同震位错量大，复发间隔长；小震产生的同震位错量小，复发间隔小于大震。因此，断裂上大、小地震之间的频度关系（震级分布）与未来一定时段内可能发生的同震位错量直接相关。

冷龙岭断裂平均复发间隔为（1 640±570）a，存在较大的不确定性。此外，2022 年门源地震震级为MS6.9，小于M7.5，同样引起了显著的同震位错。这说明可能导致冷龙岭断裂错动的地震存在较大的震级范围。综上，冷龙岭断裂在复发间隔和发震震级上均存在较大的不确定性，且不符合经典的特征地震模型。因此，本研究在上述确定性方法的基础上，进行冷龙岭断裂概率断层位错危险性分析。

在概率断层位错危险性分析中，需同时考虑震级分布、震级与破裂尺度的经验关系、同震位错沿断层走向分布等的不确定性。概率断层位错危险性分析结果形式类似于概率地震危险性分析（Petersen 等，2011；潘华等，2017），以不同预设位错量对应的超越概率表示，不仅针对最大同震位错量。吴果等（2022）给出了概率断层位错危险性分析完整表达形式，但其涉及的环节较多、形式较复杂。

本文基于现有资料和认知水平，对概率断层位错危险性分析模型进行适当简化：①郭鹏（2019）基于错断地貌和古地震探槽研究给出了断裂平均复发间隔，隐含有平均每隔一定年限地震破裂穿过1 次该位置的含义。因此，本研究不考虑分段破裂、级联破裂等复杂过程，且不考虑这些过程涉及的破裂传播距离、破裂能否传播到目标点的问题。②实际震例中位错量是沿断层走向变化的，以2008 年汶川地震为例，地震在北川段产生的平均垂直位错量为3～4 m，而最大位错量为（6.5±0.5） m（Xu 等，2009）。然而，Wesnousky（2008）收集了全球多次地震的地表破裂数据，发现不同地震的位错分布规律差异巨大。基于上述原因，同时考虑到冷龙岭断裂运动性质是以左旋走滑为主，本研究暂不考虑位错量随地震破裂走向变化的复杂情形。③本研究假设同震位错全部分布在单一主干断裂上，暂不考虑同震位错中一部分被分配到主干断裂以外分支断裂的情形（Chen 等，2011）。

在此基础上，基于断层上不同震级档地震之间相互独立且年发生率极低的假设，可得到概率断层位错危险性分析模型的表达形式。对于某条断层上的任意目标点，断层上发生的地震在该点产生的位错量d大于等于预设位错量d0的年发生概率为：

式中，m0和mu分别为断层起算震级和最大潜在地震震级；n（m0）为断层起算震级m0以上地震年发生率；f（m）为断层震级概率密度函数；P（d≥d0|m）为断层上发生m级地震时，在目标点处产生的位错量大于等于预设位错量的概率，P（d≥d0|m）取决于震级与最大位错量的经验关系式及其不确定性。

在概率地震危险性分析中，当超越概率水平非常低时，地震动衰减关系的不确定性校正过程可能会产生失去物理意义的过大值（Zhang 等，2021）。在对震级与最大同震位错量的经验关系式进行不确定性校正时，同样需增加适当的约束条件（吴果等，2022）。本研究参考经典概率地震危险性分析方法，使用经验关系式时将结果限制在3 倍标准差以内。此外，历史上我国西部地区记录到的最大同震位错量为14 m，由1931 年新疆富蕴8 级地震产生（新疆维吾尔自治区地震局，1985）。因此本文将考虑标准差后的最大位错量限制在14 m 以内。

2.2.2 震级分布模型

潘家伟等（2022）指出在青藏高原地区通常矩震级Mw≥7 的地震才会产生同震地表破裂，7 级以下地震不产生明显同震地表破裂（顾功叙，1983），因此本研究将M7.0 作为起算震级，最大潜在地震震级为M7.5，M7.0～M7.5 之间的震级概率密度函数可通过地震带的b值约束。

假设起算震级以上地震年发生率服从泊松模型。泊松模型是现行概率地震危险性分析的基础模型，其假设时间上先后发生的地震之间是相互独立的。假设断层上地震年发生率为λ，则实际发生的地震数X为k次的概率为：

对于本研究而言，λ 即起算震级以上地震年发生率n（m0），其等于平均复发间隔的倒数，约为0.000 61。对该参数进行简单验证，采用我国第五代地震动参数区划图中六盘山-祁连山地震带的地震活动参数、冷龙岭潜在震源区的空间分布函数（潘华等，2013）、潜在震源区中冷龙岭断裂和托莱山断裂的长度比例，可推算出冷龙岭断裂M7.0 以上地震年发生率约为0.000 58，其与0.000 61 接近。

进一步可得到断裂在ΔT年内的发震概率为：

2.2.3 评估结果

基于概率断层位错危险性分析方法原理，平均复发间隔取1 640 a，分别计算50 年超越概率2%、100 年超越概率2%和100 年超越概率1%的水平位错量（表2），对应的重现周期分别为2 475 a、4 950 a 和9 950 a。冷龙岭断裂的平均复发间隔T为（1 640±570） a，不确定性较大。为考虑该不确定性的影响，将T取下限值1 070 a 后再次计算，对应的平均值*结果如表2 所示。

表2 概率断层位错危险性分析得到的水平位错量评估结果Table 2 Assessment results of horizontal displacement by probabilistic fault displacement hazard analysis

3 评估结果综合分析

（1）采用不同震级与最大位错量的经验关系式，确定性方法评估结果显示出较大差异，根据冉洪流（2011）给出的关系式计算得到的水平位错量最大，为4.36 m，根据Wells 等（1994）给出的关系式计算得到的水平位错量最小，为2.32 m。概率断层位错危险性分析结果受所用经验关系式的影响显著，可知采用确定性方法和概率性方法时均须同时采用多组经验关系式，以考虑结果的不确定性。

（2）当超越概率水平为50 年2%、100 年2%和100 年1%时，概率断层位错危险性分析结果均值分别为1.82 m、3.17 m、4.61 m，可知评估结果随着超越概率的降低而增大，这与概率地震危险性分析中地震动随超越概率的变化规律一致。概率断层位错危险性分析可根据工程设施的重要程度和设计使用寿命，提供不同超越概率水平的位错参数，较确定性方法更灵活。

（3）考虑平均复发间隔的不确定性，将其由1 640 a 改为下限值1 070 a 时，概率断层位错危险性分析结果有所增大，平均值的放大倍数随着超越概率的降低而减小，这说明断裂的平均复发间隔对高超越概率下的概率断层位错危险性分析结果影响更显著。

（4）总体上，确定性方法评估结果介于100 年超越概率2%和100 年超越概率1%的概率断层位错危险性分析结果之间，明显大于50 年超越概率2%的概率断层位错危险性分析结果。然而，如果断裂的平均复发间隔缩短至下限值1 070 a，100 年超越概率2%的概率断层位错危险性分析结果将略大于确定性方法评估结果。考虑到大部分工程设施的抗震设防要求均低于100 年超越概率1%的结果，此时确定性方法评估结果具备一定保守性。同时，确定性方法还具备简单易行的优点，这是其被广泛采用的原因之一。但对于跨越地震活跃断层的重要设施，开展概率断层位错危险性分析有助于认识评估结果的不确定性，且存在提供更保守的评估结果可能性。

（5）本文模型进行了一定简化，还有众多潜在的不确定性值得进一步探讨，如Guo 等（2019）的研究认为冷龙岭断裂上最大潜在地震震级可达Mw7.5～Mw7.8，大于本研究所用的MS7.5。此外，Guo 等（2019）认为历史上冷龙岭断裂的同震位错量是沿着断层走向变化的，在发生本次地震的断裂西端位置约为3 m。

在现有技术水平下，断裂未来强震的同震位错量评价仍面临诸多不确定性。因此，本文虽提供了不同方法的计算结果，但最终的取舍和决策仍需要专家判断。如果开展更深入的研究，可建立逻辑树，以减小认知不确定性的影响（Youngs 等，2003）。

4 结论

本文基于冷龙岭断裂最新的地震活动参数研究成果，同时采用确定性方法和概率断层位错危险性分析方法评估了该断裂未来强震的水平位错量，得出以下结论：

（1）不同研究者提供的震级与最大位错量的经验关系式之间差异较大，对于确定性方法评估结果和概率断层位错危险性分析结果均有显著影响，建议同时采用多组关系式，以考虑不确定性。

（2）通常确定性方法仅提供单一的位错量评估结果，而概率断层位错危险性分析结果随着超越概率的降低而增大，可供不同重要程度和设计使用寿命的工程设施选择相应的抗断参数。

（3）确定性方法评估结果介于概率断层位错危险性分析结果中100 年超越概率2%和100 年超越概率1%之间，显著大于50 年超越概率2%的结果。这说明确定性方法提供的同震位错量评估结果对于大部分工程设施而言是偏保守的。然而，随着断裂平均复发间隔的缩短，概率断层位错危险性分析结果会进一步增大。因此，对于地震活动十分活跃的断层，有必要开展概率断层位错危险性分析，以提供更保守的评估结果。

（4）本文基于现有资料和研究程度，对概率断层位错危险性分析中的破裂传播过程和位错分布模型进行了简化，后续需开展更深入的研究。在现有技术水平下，断裂未来强震的同震位错量评价过程中仍有诸多不确定性难以避免，不同方法评估结果的取舍仍需要专家的主观决策。在条件允许时可建立逻辑树，以减小认知不确定性的影响。

致谢 感谢中国地震局地质研究所郭鹏副研究员对本文提供地震地质数据和相应的解释，感谢审稿人给予的专业修改建议，感谢编辑老师对本文的润色。