四川长宁地区地震活动空间分布特征研究.pdf

书书书第 卷第 期（）年 月中 国 地 震 扈桂让，张广伟，梁姗姗 四川长宁地区地震活动空间分布特征研究 中国地震，（）：四川长宁地区地震活动空间分布特征研究扈桂让）张广伟）梁姗姗）山西省地震局，太原 ）应急管理部国家自然灾害防治研究院，北京 ）中国地震台网中心，北京 摘要近年来，四川长宁地区地震活动频繁，特别是 年 月 日发生的长宁 地震造成长宁及周边地区较为严重的人员伤亡和财产损失，引起社会和科研工作者的广泛关注。为探究长宁地区地震的发震构造和机理，利用双差定位方法对长宁 地震前后一年时间段内的地震进行重定位，并采用波形拟合方法获得了 次 地震的震源机制解，综合分析了长宁地区地震活动的空间分布特征，获得以下主要认识：年长宁 地震的发震构造为长宁背斜 向的高倾角断层，地震序列的震源深度西深东浅，该地震是位于滑脱面之上先存断裂活化的构造地震事件；南部建武向斜地区地震活动频度高，且表现出丛集性的分布特征，推测与零散分布的页岩气开采活动密切相关；长宁地区南北 个区域地震活动的差异性受控于局部应力场方向的变化，北侧历史盐矿开采、南侧页岩气开采以及岩石介质物性的横向差异直接影响了地震的空间特征。关键词：四川长宁地震双差定位震源机制丛集分布 文章编号 （）中图分类号 文献标识码 收稿日期 修定日期 项目类别中国地震局地质研究所基本科研业务专项（）、应急管理部国家自然灾害防治研究院基本科研业务专项（）共同资助 作者简介 扈桂让，男，年生，高级工程师，主要从事地震活动研究。：张广伟，通讯作者，男，年生，副研究员，主要从事地震定位研究。：引言长宁地区位于四川盆地的东南缘，是小震活动较为频繁的区域，在其周边发生过 年汶川 、年芦山 以及 年鲁甸 地震等（图 （）。近年来，长宁地区地震频度明显增强，其中 年连续发生 次 地震，而 年 月 日长宁 地震造成长宁县及其周边地区较为严重的人员伤亡和财产损失（易桂喜等，；，；梁姗姗等，），引起社会和科研工作者的广泛关注。这些中等地震并未造成明显的地表破裂，且均远离已知的区域活动断裂（图 （），其成因机理受到广泛关注。另外，长宁地区历史上具有长时间的盐矿开采活动（，），且近十几年来区域内开展了密集的页岩气开采（，；，）（图 （），因此详细分析该区域地震空间分布特征对于探讨发震构造尤为重要。中国地震 卷注：（）研究区周边断层及近年来中强地震分布；（）长宁地区 地震及主要构造分布，其中盐井位置据 等（），页岩气井位置据 等（）；（）横切建武向斜和长宁背斜的人工地震剖面，其中红线代表滑脱面，虚线代表断层（据何登发等（）。图 研究区域构造背景关于长宁地区地震的发震机制，不同学者开展了大量工作。其中，针对长宁 地震序列的精定位结果显示，发震破裂面呈 走向，断层面较为陡立，且震源深度表现出西深东浅的趋势（易桂喜等，；徐志国等，；郭志等，）。震源机制研究结果表明长宁 主震为逆断型兼左旋走滑性质（梁姗姗等，），震源破裂位置位于长宁背斜上的高倾角发震断层（图 （）、（）。区域三维速度模型揭示出长宁地区高波速比和低波速比相间存在的特征，暗示地震触发机制与区域结构和流体注入密切相关（，；，；，）。此外，等（）利用高分辨率的地震反射剖面揭示出长宁地区断层深部三维几何结构，提出基底先存断层系统为地震的孕育提供了构造条件。上述关于长宁地区发震机理的研究对认识四川盆地边缘地震的孕震环境及地震活动特征具有启示意义。然而，精细的震前 震后地震序列空间分布和震源破裂机制对探讨区域地震活动水平，理解发震断层破裂延展范围具有重要的作用。为此，本研究收集了 年长宁 地震发生前后近一年时间段内的区域台网震相数据，通过双差地震定位方法（，）获取地震的空间位置；进一步利用宽频带地震波形数据，采用波形拟合方法（，）反演震源机制解，分析长宁地区地震空间分布的差异性，并讨论不同构造位置地震的发震机理，为探究长宁资源开采区地震发震机理提供新的认识。期扈桂让等：四川长宁地区地震活动空间分布特征研究 资料与方法 地震重定位本研究使用 年 月 日 年 月 日四川地震台网的观测报告震相到时资料，选择震中距 范围内、至少 个到时记录的地震事件（图 （），最终符合条件的地震共 个，包含 条 波和 条 波走时数据（图 （）。采用双差定位法（，）开展重新定位研究，该方法利用台站相对走时残差来修定地震初始位置，基于相邻 个事件传播到台站的射线路径几乎相同，有效降低了速度模型不准确性对定位结果的 影 响。在 地 震重 定位 过 程 中，将 事 件 对 之 间 的 最 大 距 离 设 为 ，构 建出 条 波和 条 波走时差数据。初始一维模型参考张广伟等（）和 等（）的研究结果（图 （）。注：（）本研究所使用的地震台站（蓝色三角形）和地震（白色圆圈）分布；（）波和 波观测走时曲线；（）地震定位和震源机制反演采用的一维速度模型。图 地震定位和震源机制反演资料 震源机制反演震源机制反演采用近震波形拟合方法 （，），该方法将地震波形分为 和 波（或面波），并对两部分赋予不同权重，计算理论和实际波形的拟合误差函数。反演前，首先对原始波形去除仪器响应、去倾斜、去平均值及旋转到大圆弧路径；在反演过程中，和 波滤波范围分别为 和 ，走向、倾角和滑动角的搜索间隔为 ，深度搜索间隔为 。格林函数采用频率 波数法（法）计算（，），采样间隔设为 ，采样点为 个。中国地震 卷 结果与讨论 地震空间分布本研究获得四川长宁地区 个地震事件的重定位结果，其 、和垂直方向的平均相对误差分别为 、和 ，较大提高了地震的定位精度。为直观展示震中位置随时间的空间变化，以 年 月 日为初始时间起点，给出地震的平面分布图，如图 所示。由图可以看出，地震活动的初始位置主要集中在建武向斜，该区域正是页岩气开采的主要位置；随着时间的推移，地震活动向南部扩展，至 年 月在长宁背斜发生 向的条带状地震活动。图 长宁地区地震事件随时间变化的空间分布进一步给出重定位后地震平面及深度分布图，如图 所示，由图可见地震在空间上呈现丛集性分布，主要分为 个丛集（图 （）。其中，丛集 为长宁 地震的发震区域，该区位于盐矿开采区，地震序列表现为 向条带状分布，与长宁背斜走向方向相一致，长度范围约 ，呈现出北西深、南东浅的趋势，横切剖面显示断层面较为陡立（图 （）。丛集 位于建武向斜，其地震分布的优势方向性并不明显，年 月 日 和 年 月 日 地震均位于该丛集，重定位结果也显示出 次地震的深度位于人工反射剖面揭示的滑脱面上（何登发等，），表明向斜构造对地震的发生起到控制作用（图（）。丛集 地震显示出倾斜的柱状分布特征，并且震源深度较浅，该位置处于页岩气开采区，推测该区地震活动与开采活动相关（图 （）。丛集 地震序列展现出明显的 向条带分布，震源深度由北向南逐渐加深，横切剖面清楚地刻画出陡立的断层面，表明该处存在一条近 向的隐伏断层（图 （）。丛集 显示出一条 向的断层，长度约 ，且断层面较为陡立，与其他丛集不同，该丛集震源深度相对较深，达到 （图 （）。丛集 地震序列在平面上表现为窝状分布，与丛集 类似，其发震机制可能与页岩气开采活 期扈桂让等：四川长宁地区地震活动空间分布特征研究动相关（图 （）。注：表示 个地震丛集，不同丛集给出沿矩形框长轴和横切剖面图；图（）中 剖面的灰色线表示滑脱面；每个横切剖面的虚线代表可能的断层面，星号代表 地震；剖面纵横坐标比例相等。图 长宁地区重定位地震事件平面和纵剖面分布 震源机制分布特征采用近震波形反演 方法，获得了 次 地震的震源机制解（表 ）。图 展示出长宁 地震震源机制反演的理论波形和实际波形拟合图，整体上波形拟合较好，波形相关系数较高，表明反演结果可靠；通过不同深度反演，在矩心深度 处所得到的震源机制最优。最终得到 地震的双力偶机制解节面走向 ，倾角 ，滑动角 ；节面走向 ，倾角 ，滑动角 ；矩震级为 。断层面几何参数显示此次地震 中国地震 卷表 地震震源机制解事件号发震时刻（年 月 日 时：分：秒）矩心深度 矩震级（）节面 （）节面 （）走向倾角滑动角走向倾角滑动角 ：注：波形下方数字为各段理论地震波形相对实际观测波形的移动时间，以及波形拟合的相关系数（）；波形图左侧字母为台站代码，其下方数字分别为台站震中距（）和理论值与观测值拟合偏移（）。图 长宁 地震震源机制解的波形拟合（）及不同深度震源机制拟合残差（）期扈桂让等：四川长宁地区地震活动空间分布特征研究为一次逆冲为主、兼少量左旋走滑分量的地震事件，断层面较为陡立，且倾向 。结合震源区构造走向（易桂喜等，；何登发等，）和余震空间分布特征，推断发震断层面与 向的长宁背斜主轴走向一致，是位于背斜构造上的高倾角断层（，）。图 （）显示发生在建武向斜的地震事件 和 表现为正断和走滑性质，这与震源发生的位置处在滑脱层的特殊构造上密切相关（，）。地震事件 为走滑机制，倾角较为陡立，与精定位地震序列显示的近 向断层条带具有较好的一致性，进一步表明该处存在一条隐伏断层。地震事件 位于 向的断层上，与地震空间展布方向也较为一致，印证了震源机制解的可靠性。共 个地震事件均位于长宁 地震序列的长宁背斜上，表明区域地质构造控制了地震展布的方向，褶皱及伴生断层为发震断层（图 （）；另外，震源破裂机制主要以逆冲挤压型为主，断层面的走向自 至 向存在少许差异，震源机制节面走向的差异暗示此次地震序列受到区域构造应力场的控制，但是存在小尺度的应力方向变化（胡幸平等，）；同时，三维速度结构也显示出在长宁背斜的 向存在一个高速异常体（，），断层的 方向可能受到该高速异常体的阻挡，这种结构的横向不均匀性也可能导致应力场方向产生了局部变化，从而造成震源破裂机制的复杂性。注：图（）长宁地区 地震震源机制解，地震事件号与表 一致，其中震源球上灰色区域代表压缩区，白色代表拉张区，震源球采用下半球投影，震源球上符号“”表示反演使用台站投影。图 震源机制（）与区域地质构造简图（）反演获得的震源矩心深度均较浅，主要集中在 之间（表 ），这暗示地震的发生可能与开采活动密切相关。在研究区的长宁背斜，历史上具有长时间段的盐矿开采（阮祥等，；，），而在建武向斜区，近十几年来页岩气开采频度显著增加（，）。大量的研究结果揭示出开采注水活动与地震发生存在密切关系（，、），本研究获得较浅的地震矩心深度也表明高压注水可能触发了区域浅部先存断层的活动。同时，矩心深度较浅的地震对地表建筑物破坏性也较大，这从另外一个角度解释了长宁 地震造成地表灾害严重的重要原因。中国地震 卷综上，基于精细的地震空间位置和震源破裂机制，并结合人工地震剖面结果（何登发等，；，），给出长宁地区的发震构造示意图（图 ）。整体上，震源区主要受到长宁背斜和建武向斜的控制，长时间的开采注水活动触发了不同尺度、不同几何产状的先存断层活动。北部地震活动主要与盐矿开采有关，在背斜核部存在高角度的 向断层，因此地震的 向条带也较为明显；南部地震活动主要与页岩气开采活动有关，开采井分布较为零散，因此地震的分布也呈现出多个丛集式的展布。胡幸平等（）指出长宁地区地壳应力场的最大主应力轴在整个区域内基本均处于近水平状态，然而北侧和南侧存在显著差异，这种局部改变是长宁地区复杂地震活动的必要力学基础。本研究也认为介质岩性的横向差异和高压流体注入直接影响了长宁地区地震活动的空间分布格局。图 长宁地区地震发震构造示意图 结论本研究利用双差定位方法对 年长宁 地震发生前后一年时间段内的地震进行重定位，并采用波形拟合方法获得了 次 地震的震源机制解，综合分析了长宁地区地震活动的空间分布特征，获得以下主要认识：（）年长宁 地震的发震构造为长宁背斜 向高倾角左旋逆冲兼走滑断层，地震序列的震源深度呈现西深东浅，是 向先存断裂活化的构造地震事件。（）南部建武向斜地区地震频度远远高于北部，并且表现出丛集性的分布特征，推测其与页岩气开采活动密切相关。（）长宁地区南、北 个区域地震活动的差异性受控于局部应力场方向的变化，而介质岩性的横向差异和高压流体注入直接影响了地震活动的空间分布特征。致谢：感谢 审 稿 专 家 提 出 的 宝 贵 修 改 意 见；感 谢 中 国 地 震 台 网 中 心 国 家 地 震 科 学 数 据 中 心（：）为本研究提供数据。期扈桂让等：四川长宁地区地震活动空间分布特征研究参考文献郭志，高星，路珍 年 月 日四川长宁地震重定位及震源机制研究 地震学报，（）：何登发，鲁人齐，黄涵宇，等 长宁页岩气开发区地震的构造地质背景 石油勘探与开发，（）：胡幸平，崔效锋，张广伟，等 长宁地区复杂地震活动的力学成因分析 地球物理学报，（）：梁姗姗，徐志国，盛书中，等 年四川长宁 级地震主震及中强余震（）的震源机制及其应力场 地震地质，（）：阮祥，程万正，张永久，等 四川长宁盐矿井注水诱发地震研究 中国地震，（）：徐志国，梁姗姗，盛书中，等 年四川长宁 地震序列重定位和震源特征分析 地震学报，（）：易桂喜，龙峰，梁明剑，等 年 月 日四川长宁 地震序列震源机制解与发震构造分析 地球物理学报，（）：张广伟，雷建设，梁姗姗，等 年 月 日云南鲁甸 地震序列重定位与震源机制研究 地球物理学报，（）：，：，（）：，（）：，（）：，（）：，（）：，：，：，（）：，：，（）：，（）：中国地震 卷 ，），），），），），：，：；