2022年四川泸定M_S6...的自动构建与地震活动性分析_刘泽民.pdf

第 卷 第期 年月地球物理学报 ，刘泽民，李俊，苏金蓉等 年四川泸定 地震余震序列的自动构建与地震活动性分析地球物理学报，（）：，：，.（），（）：，：年四川泸定 地震余震序列的自动构建与地震活动性分析刘泽民，李俊，苏金蓉，于子叶，黄春梅，祁玉萍，陈依伲，王伟涛，中国地震局地球物理研究所，北京 四川省地震局，成都 浙江省地震局，杭州 中国地震局震源物理重点实验室，北京 摘要北京时间 年月日 时 分，四川甘孜州泸定县发生 地震（以下简称泸定地震）本文基于神经网络方法从泸定地震周边 之内 个台站 年月 日的连续波形资料中拾取 震相到时，对由优选的震相识别阈值筛选后的震相到时进行震相关联和定位分析，经上述处理流程获得自动检测目录（）包含 个地震事件，为同期台网人工目录（即中国地震台网中心目录，）数量的 倍，完备性震级由 下降至 基于波形互相关的双差定位方法获得的 个地震的时空分布特征表明，年泸定 地震的主震及余震主要发生在鲜水河断裂磨西段，沿着该断裂在 方向分布并表现出明显分段特征，北段总体深度较浅，深度集中在 以上，中段和南段相对较深，深度分布在 之间 此外，余震序列在鲜水河断裂上展现出高倾角的特征，而在主震 方向的海螺沟断裂附近与鲜水河断裂呈共轭展布，推测在中段海螺沟附近有一条走向、倾向 的低倾角分支断裂 研究结果表明，利用人工智能方法快速实现中强地震余震序列的自动化检测是获取地震发震构造、破裂方向等信息的有效方法关键词泸定地震；神经网络；余震序列；地震精定位 ：中图分类号 收稿日期 ，收修定稿基金项目国家自然科学基金项目（），中国地震局地球物理研究所基本业务专项（）和地震科技星火计划项目（和 ）资助第一作者简介刘泽民，男，在读博士，主要从事地震学方面的研究 ：通讯作者李俊，男，高级工程师，主要从事地震监测和地震预测方面的工作 ：，期刘泽民等：年四川泸定 地震余震序列的自动构建与地震活动性分析 ，（），（），（），（）（），；引言中强地震往往会伴随有大量的余震，余震的快速、准确测定对于分析地震发震构造和机理、评估地震灾害分布等具有重要意义 震源机制结合余震分布可有效判别发震主断裂，并且高精度的余震序列可用 于 更 好 地 刻 画 发 震 断 层 的 精 细 结 构 相 比 等基于地表信息的研究，余震序列可提供深部的信息，余震的时空分布特征也对后续地震的分析研判具有重要参考价值（杨智娴等，；房立华等，；张广伟，）因此，开展快速的余震序列分析对研究中强地震的发震构造和孕育环境以及分析地震危险性意义重大中强地震引发的余震在主震发生之后密集发生，随时间推移余震发生频次逐渐降低 不同震级余震数量的统计关系符合 关系，即小震级余震数量较多 而 大 震 级 余 震 数 量 较 少（，）增强小震监测能力，获得较为完备的余震序列目录有助于刻画断层空间形态和分析强余震危险性地震台站数量的增加可提升地震监测能力，但也导致余震序列的人工处理分析难以兼顾时效性和完备性 近年来基于机器学习的震相识别和地震检测方法得到了快速发展（，；于子叶等，；，；，；，；，；，；，），在多次大地震的余震目录构建中表现出高效、快速的优势 相关方法应用于 年美国 和 地震、年青海玛多 地震、云南漾濞 地震和双柏 地震的余震序列分析（，；，；廖诗荣等，；苏金波等，；郭慧丽等，），也可用于中强地震余震序列的回溯性研究（，），甚至在监测非天然地震活动方面也具有一定的优势（，）对相应技术进行更多应用，可进一步细化流程，提升技术成熟度据中国地震台网中心（）测定，北京时间 年月日 时 分，在四川省甘孜州泸定县（，）发生 地震（图）泸定地震发生在川滇菱形地块东边界的鲜水河断裂带南段，该区域构造十分复杂，是鲜水河断裂带南段、安宁河断裂带北段和龙门山断裂带南段交汇区，俗称三岔口地区鲜水河断裂带在历史上曾多次发生级以上地震，其中最大为 年四川康定南的级地震（距离本次主震约 ），最近一次为地 球 物 理 学 报（）卷 年 月 日康定 级地震 本次地震主震震中就位于鲜水河断裂带磨西段附近泸定地震最高烈度为度，等震线长轴呈 走向，长轴 ，短轴 ，不仅引起了当地建筑的损坏（，），同时还诱发了严重的地质灾害，尤其是滑坡分布广泛，导致地表调查十分困难（，）余震序列的时空分布成为分析地震构造特征，快速评判地震破坏范围的重要方式 根据中国地震台网目录，截至 年月 日 时，泸定地震序列共记录到 余震 次，最大余震为月日（）地震人工定位结果显示主震深度约 ，破裂主要向浅部和 方向传播，余震总体沿 向展布，长度约 ，深度在 ，初步显示断层形态较为复杂完整的地震余震序列目录有利于对发震断层进行细致的分析 为完善自动化构建地震目录的流程，本文分析了地震序列构建中震相识别阈值的选取策略，利用神经网络算法拾取了距泸定地震主震 范围内 个地震台站（图）年月日至 日连续记录中的震相到时，通过震相关联和地震定位，构建了泸定地震高精度的余震序列目录 高精度的余震目录刻画了泸定地震余震的时空演化特征，较好地揭示了 向主断层、向共轭断层以及北部雁列状断裂的深部构造特征，有助于深入认识泸定地震的孕震环境和发震机制数据泸定地震区域为川滇菱形地块、巴颜喀拉地块和华南地块的交汇区，西临四川省内最高峰贡嘎山本研究选取了 主 震 震 中 范围内的 个测图研究区域和本文使用的地震台站分布图红色五角星是泸定地震主震震中位置；蓝色三角形是泸定地震周边 范围内的地震台站（个）；灰色圆点表示 年以来该区域发生的 级以上的地震事件；震源机制解来自于 ：在子图中，青藏高原地块；，滇缅地块；，华南地块；，西域地块；，东北亚地块；，华北地块 ：，期刘泽民等：年四川泸定 地震余震序列的自动构建与地震活动性分析震台站，其中包括四川数字地震台网的 个固定台站和 个重庆台网和云 南 台网 架 设 的 流 动 台 站（、）在地震检测中，近场台站可提取可靠的 震相，良好的台站方位覆盖有利于震相关联和定位 如图所示，在距离主震 范围内有 个地震台站，范围内有 个地震台站 从台站分布上，除在主震东南 范围台站较为稀疏外，已有台站的分布相对震源区有较好的距离和方位角覆盖，台站的整体分布利于余震序列的自动化构建本文使用了 个台站 年月 日的三分量连续数据来构建泸定地震余震目录 为对网络拾取阈值进行定量分析，利用了其中 个长期运行的高数据连续率台站 年月的数据进行了拾取测试，以优化拾取结果方法余震目录的自动构建主要包含三个步骤：（）震相拾取；（）震相关联；（）地震重定位 详述如下：基于神经网络的波和波震相拾取及置信度阈值优选基于神经网络的震相拾取算法提高了震相到时的拾取精度和准确度，相比经典的长短时窗法，神经网络拾取的震相中含有震相的类型信息（如 ，），为后续的震相关联提供了便利（，；，；，）目前的震相拾取神经网络多利用监督学习方式训练，基于特定研究区数据样本训练得到的网络在应用于研究区之外的数据时，往往会引起震相拾取精度和准确度的降低（，）加入目标研究区数据样本进行迁移学习以及采用覆盖研究区的较大数据集进行训练是提升网络性能的可行方式（赵明等，；，）等（）利用中国地震台网 年的震相数据，结合 地震波形对几种震相拾取网络进行了训练和评估，结果表明感 受 野（）较 大 的 以 及 网络的拾取效能较优 本文中拾取和分析震相的神经网络是在 网络的基础上改进的 网络，该网络去除了 网络的 层，省去了其地震事件判断的概率输出，利用 函数直接输出和震相的概率，以提高模型的训练效率 在模型构建中，使用了中国大陆地区 年共 年级以上地震的人工标注信息以及其对应的震相波形作为数据集，数据采样率选择了通用的 模型以 年的数据作为训练，年的数据进行测试 等（）结果表明，改进的 模型结构相比于 有较大程度的精简，利于训练收敛，且在精度上略有提升神经网络震相拾取算法利用网络得到数据点置信度来进行震相拾取，当置信度超过特定阈值时，输出对应的震相时间和类型 置信度阈值的选取会影响最终地震目录的构建 当选择较高阈值时，拾取震相和后续关联的可信度增加，同时也可能会降低较小震级地震的检测能力；反之降低阈值，可提高震相识别数量，但也可能会引入部分虚假震相影响后续目录构建 同时，受网络结构、训练数据、加权策略等的影响，不同的网络对同一数据的置信度表现也有差异 多数研究往往会使用经验性的阈值来开展余震检测研究（，；苏金波等，；郭慧丽等，）根据固定台网长期运行以及四川多震的特点，本文利用研究区历史地震记录寻找最优的置信度阈值，以进一步优化震相拾取流程 本文利用图台站中数据连续率高的 个固定台站 年月至 年月的连续波形，基于历史地震在对应台站上的震相编目信息，对不同震相拾取阈值进行测试 根据人工目录中对已有地震事件进行定位所使用的台站上的 和 震相信息，选取对应的波形，利用神经网络识别波形中的 和 震相，通过统计震相识别结果的准确率、召回率以及调和平均来实现拾取阈值的优选 ，（），（），（）其中 为正确拾取的震相数量，在本文中定义为神经网络拾取震相到时与人工标注到时相差在 之内的震相；是神经网络识别为但按上述标准无对应人工标注的震相数量高准确率可实现更加可信的震相识别，而高召回率则可保证不遗漏有效 震 相，可 以作为二者的 均衡评估指 标（，；，）通过相应评估，可判断所使用的网络需要多大的置信度阈值才能够较完整地识别人工能够判别的 震相，从而指导震相到时的智能拾取地 球 物 理 学 报（）卷利用 年月 个 的地震在 个台站上评估结果如图所示，随着置信度的提高，波拾取准确率上升（图），但是 波的准确率在置信度 之后下降（图），可能由于波波形复杂，网络判断波的最高置信度偏低引起曲线是准确率和召回率的调和平均，可作为网络整体性能的表征 从图中可以看出，网络整体性能均衡点所对应的 波置信度阈值约为，波为 基于人工智能自动化构建地震目录的目标是获得更加完整的地震目录，选用较低的阈值可提升震相查全率，有利于检测更多地震 本文通过将震相拾取的置信度阈值降低至 来提升震相的查全率召回率，从而减少遗漏震相 图 也分别显示了对应 波的拾取精度，从中可以看出 波到时差的均值为 ，标准差为 ；波到时差的均值为 ，标准差为 该拾取阈值可较好的满足地震定位的震相拾取精度需求图神经网络震相拾取效能与置信度阈值选取（，）和 震相到时拾取的准确率（）、召回率（）和调和平均与其置信度的关系图；（，）置信度阈值为 时 震相拾取偏差的分布 左上角注明了召回率、准确率、调和平均、平均值及标准差 （），；（）（）；（，），（），（），（），（）利用优选的阈值，本文对 个台站 年月 日的三分量数据进行了处理分析数据采样率为 ，仅做去均值和去趋势的预处理图显示了 台站上 窗口内的拾取结果从中可以看出，高信噪比震相置信度可超过 ，而低信噪比震相的置信度偏低 选取 为阈值，可以对信噪比较低的震相进行准确的到时拾取和类别判定，同时不引入过低置信度的信号（图，附近）基于网格搜索的震相关联利用基于网格搜索的 程序来对拾取到期刘泽民等：年四川泸定 地震余震序列的自动构建与地震活动性分析图神经网络震相到时拾取结果 的震相进行震相关联和初步定位（，），该算法依据一维速度获取理论走时曲线，通过格点搜索将震相数量以及发震时刻的标准差作为指标，确认对应地震的发震时刻和概略位置 文中采用的格点搜索策略为：以台站为中心，水平向搜索范围为 ，搜索步长为；深度方向的搜索范围为 ，搜索步长为 成功关联的震相数量策略为：震相数量不少于个，震相数量不少于个，震相不少于个 为保持关联稳定性，要求同时具有和震相识别的台站数目不少于个 在网格搜索中，震相关联最大范围为 图显示了最终关联结果以及所使用的一维速度模型，该模型是以康定速度模型（易桂喜等，）为初始模型，通过多阶段定位法（，）对泸定地震的余震序列进行定位反演获得 通过震相关联，共从拾取的震相中获取具有时空关联的 条波震相和 条波震相（图）地震重定位以震相关联获得地震的概略发震时间和位置为初始值，本文调用了 系统中的 定位方法对关联到的地震进行了绝对定位，最终获得的自动检测目录（）包含了 个余震事件（图）之后，基于图 所示的一维速度模型，利用双差定位方法（，图震相关联结果的走时曲线图（）与震相关联和地震重定位所使用的一维波（红）和波（蓝）速度模型（）（）；（），）对这些余震进行了精定位分析 在双差定位中，选择震源距离小于 的事件对计算走时差，震中和台站之间的最大距离小于 ，并只对有个及以上相邻事件的地震事件进行重定位，定位过程中采用阻尼最小二乘算法进行四次迭代，利用震相走时差获得的地震目录（）中包含 个地震（图）之后，加入基于波形互相关的震相走时差数据来提高定位精度，在计算震相波形的互相关时，选取相对于震相到时，地 球 物 理 学 报（）卷图不同目录下泸定地震主震和余震空间分布图（）中国地震台网中心（）目录；（）自动检测（）目录；（）基于震相走时差的双差定位（）目录；（）基于波形互相关的双差定位（）目录 （）（）；（）（）；（）（）；（）（）的震相波形并对其进行 的带通滤波对于互相关的计算结果，本文中仅保留互相关系数不小于、信噪比大于的震相走时差，最终获取了高精度的精定位目录（），该目录中含有 个地震事件（图）由图可知，相比人工目录，自动生成的目录中期刘泽民等：年四川泸定 地震余震序列的自动构建与地震活动性分析地震数量有大幅的增加 相对完备的余震目录易于勾勒断层结构，双差定位之后地震分布较好地显示了地震的丛集性分布特征 基于直接震相走时差以及波形互相关的重定位目录中地震数目变化不大，整体分布类似，表明本研究自动拾取的震相到时具有较高的精度结果与讨论 定位结果对比为验证自动余震检测结果的可靠性，本文将获得的自动 检 测 目 录 与台网人工目录 进 行 了 对比 当两个目录中的地震发震时刻差在以内，水平定位偏差在 之内时，将其视为同一地震 在图所示的区域内，目录在 年月 日 内 包 含 个 地 震 事 件（图 ），目录识别了其中的 个，地震检测率达到 由图可知，在检测到的地震中，发震时刻偏差为 （），水平定位偏差主要在 之间，深度偏差主要在 以内，震级误差为 （），展示了自动检测结果的可靠性图震源参数偏差分布图 目录和 目录共同事件的（）发震时刻差、（）水平定位偏差、（）震级差以及（）深度偏差 （）；（）；（）（）本文所得到的自动检测余震目录在数量上为人工目录的倍，地震总数量是人工地震目录的 倍（图）震级频度关系（关系）可表达为 （），（）式中，表示地震震级，（）为震级不小于的地震数量，和为常数 为研究余震序列的 关系，本文统计了 ，）震级范围内的地震数量（），其中震级在 之间取值，取值间隔为 将地震数量最大时对应的作为该目录的完备震级 基于公式（），采用最小二乘算法拟合了与 （）之间的关系（图）结果表明，自动检测余震目录较好地满足 关地 球 物 理 学 报（）卷图 年月 日 目录（黑色）与 目录（灰色）地震数量对比（）每小时频次图，灰色实线表示每天 目录与 目录地震数量之比；（）累计数量图 （）（），（）（）（）图 目录（）和 目录（）的关系图 （）（）系，其完备震级由人工目录的 下降至，说明在检测人工难以完整识别的小震级余震方面，基于人工智能的自动检测方法具备较大优势 相对于人工目录，自动检测目录在地震数量及最小震级方面均有提升，也有利于计算更为精确的值利用其构建的 关系外推估计最大余震的震级为，与实际最大余震 符合更好 泸定地震余震分布与发震构造本次泸定 地震发生在鲜水河断裂带磨西段附近，地震所处区域有多条断层发育，构造复杂 鲜水河断裂带是青藏高原东缘一条 走向的左旋走滑全新世活动断裂带（，），断裂带长约 ，以惠元寺拉分盆地为界，呈现出北西段和南东段两个结构不同的段落（闻学泽等，），北西段断层结构特征比较单一，自北西向南东由炉霍段、道孚段和乾宁段呈左阶羽列而成，晚第四纪水平滑动速率达 （闻学泽等，；唐荣昌等，）康定以南是鲜水河断裂带磨西段（图），长度达 ，断层带宽约 ，水平 滑动速率在 左右，走向 ，期刘泽民等：年四川泸定 地震余震序列的自动构建与地震活动性分析图 目录中余震震源位置与断层分布，鲜水河断裂磨西段；，大渡河断裂；，安宁河断裂；，燕子沟断裂；，磨子沟断裂；，海螺沟断裂；，石棉；，得妥；，泸定；，西店子；，加郡；，冷碛镇断层数据来源于邓起东等（）和屈春燕（）；红色短棒示意断层在剖面上的位置 ，；，；，；，；，；，；，；，；，；，；，；，（）（）；倾向，呈高倾角左旋走滑性质（闻学泽，）泸定地震震中位置比较靠近 走向的龙门山断裂带大渡河断裂 以及 走向的安宁河断裂 与鲜水河断裂三者交汇处，易于应力积累 此外，该区域还发育有 走向的三组近似平行的断裂，靠近震中的海螺沟断裂 走向 ，在东地 球 物 理 学 报（）卷段偏转为 ，该断裂为压扭性，以右旋水平错动为主，断层破碎带近乎直立，略向南倾，南盘相对北盘上升；靠北部的磨子沟断裂 和燕子沟断裂 则为高倾角右旋走滑断层（冯静等，）图显示了本文中自动构建的泸定地震的余震序列基于波形互相关数据进行双差定位后的分布，其空间位置清晰地刻画了断层的深部位置和形态余震整体沿着鲜水河断裂 向展布，余震区长约 ，并具有明显的分段特征 主震发生在余震序列的中段，该区域余震分布较为分散，而北段、南段地震活动相对集中 定位结果显示，沿着鲜水河断裂的 剖面北段的余震总体深度较浅，深度集中在 内，中段和南段相对较深，优势分布范围在 之间，在深度空间上北段和南段呈现出一个台阶（图）余震在三个区段展现出不同的特征：在北西段，余震序列展现出了一组 走向呈雁列状的细小断裂（剖面），长度约（剖面），深度在 以上（图，）在中段西南侧，存在一个余震密集区，从精定位的结果看，沿着 剖面存在一条 走向的次级断层，倾角为 左右（图，）剖面结果显示该断裂倾向，与震源区周边发育的 走向的次级断层，如海螺沟断层、燕子沟断层相符 在南东段，地震分布显示安宁河断裂高倾角的特点（图）从余震频次和随时间演化过程来看，主震发生后的短时间内余震沿南段发生（图 ，），后开始向北段迁移（图 ，），后才在震中附近西南侧发震，之后序列开始在整个余震区遍布 余震区北段地震的频次、强度均弱于南段，主要以 以下小震活动为主（图 ，）空间和时间演化过程表明此次地震为东南向单侧破裂，主震破裂区余震分布长度约为 ，北西段和中段西南侧为主震破裂触发区域应力释放产生的地震活动，相应余震演化与主震的破裂模型相符（图）泸定地震主震的震源机制解的两个节面分别为 和 （：），其中节面 走向与鲜水河断裂基本一致，地震呈现左旋走滑兼逆冲性质，与高倾角的磨西断裂具有较好的对应关系；节面的走向与海螺沟断裂近似平行，倾角稍陡，右旋走滑的性质也与海螺沟断裂的压扭性右旋水平错动的特性一致，而海螺沟断裂的地表迹线显示其穿过鲜水河断裂，并且东侧的走向相较于西侧顺时针偏转 （冯静等，；李娴，），与鲜水河断裂的左旋走滑特性相匹配 由于鲜水河断裂磨西段水平滑移速率 左右（闻学泽，；李乐等，），而海螺沟断裂附近的余震活动未跨过鲜水河断裂带往东（图），结合地震破裂模型的结果，本文认为泸定地震的发震断层为节面 表征的鲜水河断裂磨西段 地震分布与震源区速度特征刘影等利用体波和面波联合反演构建了高精度的川滇公共速度模型 （：），其 中 包 含 了 波和波速度，该模型在上地壳的分辨率可以达到 ，中下地壳的分辨率可达到 本文对比了在 及 剖面上精定位的余震和速度及波速比的空间分布关系（图）高程较高区域对应为贡嘎山，沿 剖面，高程从北西到南东下降幅度较大；相应地，波速度结构在 剖面北段 深处有个明显的高低速界面，其向南延伸至海螺沟附近，之后阶跃式下降至深度 左右，形成了一个波高低速不连续面（图 ）结合图所显示的定位结果，地震深度分布与速度界面的阶跃变化吻合较好，总体上余震深度在北段位于 以上，中段在 之间，南段则在 间都有分布 我们推断 波速度界面的不连续性与贡嘎山抬升有关，随着贡嘎山的抬升，导致了北段总体向上，在中部形成一个跃阶，同时产生了一些可能的次级断裂（；图 ）波速比与地下介质的岩性、岩石损伤、部分熔融以及含水饱和度等相关，是对波速的有益补充（，；，）在地壳中，波速比与断裂带的破碎程度以及流体有关（，；，）图 所显示的波速比在 剖面的 左右深度存在一个低波速比条带，条带北段较窄，中南段较宽 泸定地震的余震分布基本在条带浅部边界区域，推断可能与断裂带内部岩石的损伤和破坏程度有关，也有可能与浅层岩性或者流体状态有关 从 的波速比剖面还可以看出沿西南方向的断裂分布与低波速比带分布形态相似（图 ），推测可能为贡嘎山抬升过程中产生的裂隙性断裂所致 受相关速度模型分辨率的限制，具体的影响因素尚需要开展更为精细的研究予以确定结论本文利用深度神经网络方法对 年月期刘泽民等：年四川泸定 地震余震序列的自动构建与地震活动性分析图 不同时间段 目录中余震震源位置和断层滑移量分布（，）主震后 ；（，）主震后 ；（，）主震后 （，）；（，）（，）（，），地 球 物 理 学 报（）卷图 目录中余震震源位置与波速度和波速比分布，剖面位置见图 ，日四川泸定 地震余震序列进行了检测和定位，获得了震后 天内的余震序列目录 研究结果表明，自动构建的余震目录增加了小震数量，降低了完备震级 利用波形互相关方法获得的 个地震的双差定位结果表明，本次泸定地震余震主要沿鲜水河断裂带磨西段，在 方向展布，余震深度总体呈北浅南深，与地壳内速度界面形态相符，推测与贡嘎山抬升对应的构造背景有关 余震分布在南段和北段相对集中，中段分布较为松散，部分余震沿方向与鲜水河断裂呈共轭展布，可能有一条低倾角分支断层存在 本文中所获得的余震序列可为分析泸定地震的发震机制和构造提供有益参考致谢感谢四川省地震局杜方、易桂喜研究员以及四川地震台监测预报人员提供的宝贵建议，感谢审稿专家的专业建议，感谢中国地震局地球物理研究所郭祥云副研究员等提供的震源机制解结果以及中国地震局地球物理研究所张旭副研究员等提供的震源破 裂 过 程（：）文中图件使用 软件（，）绘 制，计 算 波 形 互 相 关 的 程 序 为 （：）川滇公 共 速 度 模 型 可 由 ：下载，本文使用的震相拾取模型可由 ：下载 ：，（）：，：，（），（）：，（）：，：期刘泽民等：年四川泸定 地震余震序列的自动构建与地震活动性分析 ，（）：，：，（），（）：，：，（），（）：，：，（），（）：，：，（）：，：，（）：，：，：，（）：，：，（），（）：，：（）：，（），（）：，：，（）：，：，（）：，：，（）：，：，（）：，：，（）：，：，？：，（）：，：，：，（）：，：（），（）：，（）：，：，：，（）：，：，（）：，：，（）：，：，（），（）：，：，（），（）：，：，（）：，：，：，（）：，：，（），（）：，地 球 物 理 学 报（）卷 （），（）：，（）：，：，（），（）：，（），（）：，：，（），（）：，：，：，（）：，：，（）：，：（），（）：，（）：，：，：，（）：，：，（），（）：，：，：，（），（）：，：，（）：，：，：，（）：，：附中文参考文献邓起东，张培震，冉勇康等 中国活动构造与地震活动地学前缘，（）：房立华，吴建平，王未来等 四川芦山 级地震及其余震序列重定位科学通报，（）：房立华，吴建平，苏金蓉等 四川九寨沟 地震主震及其余震序列精定位科学通报，（）：冯静，孔军，康宏等 年月泸定地震序列重定位和震源机制研究地球物理学进展，（）：，：郭慧丽，常利军，鲁来玉等 基于深度学习震相拾取和密集台阵数据构建青海玛多 地震震源区高分辨率地震目录地球物理学报，（）：，：李乐，陈棋福，钮凤林等 鲜水河断裂带南段深部变形的重复地震研 究地 球 物 理 学 报，（）：，：李娴 贡嘎山地区旅游地学特征及开发模式研究 博士论文成都：成都理工大学廖诗荣，杨婷，张红才等 年云南双柏 地震序列的快速检测与精定位研究地震研究，（）：屈春燕 最新 万中国活动构造空间数据库的建立地震地质，（）：苏金波，刘敏，张云鹏等 基于深度学习构建 年月 日云南漾濞 地震序列高分辨率地震目录地球物理学报，（）：，：唐荣昌，张耀国，黄祖智等 四川石棉西昌地区地震区划研究地震研究，（）：闻学泽，罗灼礼等 鲜水河全新世断裂带的分段性、几何特征及其地震构造意义地震学报，（）：闻学泽 四川西部鲜水河安宁河则木河断裂带的地震破裂分段特征地震地质，（）：杨智娴，于湘伟，郑月军等 中国中西部地区地震的重新定位和三维地壳速度结构地震学报，（）：易桂喜，龙锋，闻学泽等 年 月 日康定 级地震序列发震构造分析地球物理学报，（）：，：于子叶，储日升，盛敏汉 深度神经网络拾取地震和波到时地球物理学报，（）：，：张广伟 云南地区地震的重新定位及值研究中国地震，（）：赵明，唐淋，陈石等 基于深度学习到时拾取自动构建长宁地震前震目录地球物理学报，（）：，：（本文编辑何燕）