贵州西部多流域水库区地壳三维速度结构.pdf

书书书第 卷第 期（）年 月中 国 地 震 杜兴忠，左可桢，段龙飞，等 贵州西部多流域水库区地壳三维速度结构 中国地震，（）：贵州西部多流域水库区地壳三维速度结构杜兴忠）左可桢）段龙飞）赵翠萍）中国电建集团贵阳勘测设计研究院有限公司，贵阳 ）中国地震局地震预测研究所，北京 摘要使用基于机器学习构建的贵州西部多流域水库区地震目录和震相报告，采用波速比一致性约束的双差层析成像方法联合反演得到了该地区的地震位置和三维、结构。结果显示，研究区内的地壳速度结构具有明显的不均匀性，不同大地构造变形分区展现出不同的速度结构特征。在 以上深度，横穿黔西中部 向的威宁构造变形区显示出显著的低速异常条带，揭示了威宁水城深大断裂带的影响深度和范围。受库区岩性和流体渗透的影响，深度的速度结构显示普遍的低波速和高波速比特征，包括夹岩、平寨、光照和马马崖水库区域。重定位后，地震的空间分布勾勒出大量隐伏断层的几何展布特征，结合三维速度结构，推测该区域的地震活动与水库周围的断层活化有关。关键词：贵州西部水库地震地震空间分布三维速度结构 文章编号 （）中图分类号 文献标识码 收稿日期 修定日期 项目类别 中国电力建设股份有限公司科技项目（）资助 作者简介 杜兴忠，男，年生，正高级工程师，主要从事地球物理勘探及地震监测相关技术研究工作。：左可桢，通讯作者，男，年生，助理研究员，主要从事震源参数和地震层析成像研究。：引言贵州地区在地质构造上位于西部特提斯构造区域与东部濒太平洋构造的交接处，历史上构造地震活动较弱。其中，贵州西部地震活动相对较强，但是地震监测台网建立以来发生的构造地震较少，因此一直以来对该地区的地震学研究，包括地壳三维速度结构等基础性研究较为缺乏。自 年以来，随着贵州西南地区乌江、北盘江等几大流域大型梯级水电站的建设，水库地震监测专用台网的监测结果表明黔西南地区水库区域地震活动呈现出相对活跃的状态，这为利用地震数据开展水库地震等相关研究提供了条件。等（）利用黔西南 个水库台网 个台站记录的连续观测数据，采用基于人工智能技术的 流程（，），得到 年 月 年 月期间大约 倍于人工处理的地震目录和震相数据，地震目录的完备震级从 下降至 。本文使用该研究得到的地震目录和震相数据，进一步反演黔西南地区的三维速度和波速比结构，为认识贵州西南地区的地壳结构特征，分析地震活动的深部环境，以及其他地学研究提供速度结构模型。期杜兴忠等：贵州西部多流域水库区地壳三维速度结构 贵州西南区域地质背景和地震活动情况受青藏高原隆起的影响，贵州高原山岭纵横，地表崎岖。从中元古宙至今的地史时期，贵州地壳经历了多次构造运动，造成了较为复杂的地质构造格局。岩石分布以沉积岩分布最广，发育最佳。本文研究区贵州西南地区为黔北台隆，属于二级构造。黔北台隆内部的威宁水城深断裂（属于垭都紫云断裂的组成部分）将黔北台隆划分为六盘水断陷和遵义断拱 个三级构造单元。其中，六盘水断陷由威宁 向构造变形区和普安旋钮构造变形区组成，以 向构造为主，走滑、逆冲断层均比较发育，特别是其南、北盘江流域三叠系陆源碎屑岩区，褶皱紧密，变形较强。遵义断拱则包括毕节 向构造变形区、凤岗 向构造变形区和贵阳复杂构造变形区（徐伟钧等，）。除 向的威宁水城逆冲断裂外，研究区内的主要活动断裂还包括位于北部乌江流域附近、近 向的开阳纳雍断裂和马场断裂，以及南部北盘江流域附近的 向晴隆、盘县断裂。其中，开阳纳雍断裂位于乌江以南，属于走滑断层，其切断了地台基底和盖层构造；马场断裂位于乌江以北，具有逆冲和走滑性质；北盘江南侧 走向的晴隆断裂活动性质为正断，其控制了一系列第四系盆地的发育，历史上曾发生多次 级以上地震，该断裂最新活动时代为中更新世晚期，晚更新世以来不活动。研究区内沿北盘江流域、乌江干流三岔河、六冲河分别建设了黔西北、黔中、黔西南水利工程，其中沿北盘江流域建设有三座大型水库，分别为光照水库、马马崖水库和董箐水库。光照水库坝高 ，库容 亿，年开始蓄水；马马崖水库坝高 ，库容 亿，年开始蓄水；董箐水库坝高 ，库容 亿，年开始蓄水。北盘江流域的水库库区基岩以碳酸盐岩为主，岩溶、地下暗河发育。已有研究表明，岩溶发育的地区在水库蓄水后容易诱发地震（胡毓良等，），光照水库和董箐水库在蓄水后地震活动性明显增强（王尚彦，）。三岔河中游是黔中水利枢纽，平寨水库位于三岔河中上游六枝与织金交界河段，最大坝高 ，库容 亿。平寨水库大坝左岸岩溶十分发育，受岩性和构造影响共发育五套岩溶系统（孙浩淼，），这也是河流以北地震相对较多且集中在上游几个支流附近的原因。乌江干流六冲河上的洪家渡水库坝高 ，库容 亿，库区附近岩溶发育、地质构造复杂、裂隙节理丰富，处于相对较高的应力状态区（姜朝松等，）。该水库于 年正式开闸蓄水，在蓄水后诱发了有影响的地震（王尚彦，）。位于其上游的夹岩水库坝高 ，库容 亿，于 年 月 日开始正式蓄水。上述 个流域水库地震监测台网的监测数据在联合处理后，可极大地提升整个黔中区域的地震监测能力。与西侧的川滇地区相比，研究区内（，）历史地震活动相对较弱。公元前 年至 年共记录到 级以上地震 次；年至 年共记录到 级以上地震 次，级以上地震 次，其中 次位于大型水库附近（图 （）；年至 年期间共记录到研究区内地震超过 次，大量地震丛集位于北盘江流域的董箐水库和光照水库库区附近（图 （）。地震活动频次、强度与水库蓄水具有一定的相关性，其中光照库区蓄水后发生的最大地震震级为 ，董箐库区发生的最大地震震级为 。中国地震 卷注：（）研究区位置及历史地震 震 中 分 布；（）年 月 至 年 月 地 震 震 中 分 布；（）利 用 得到的 年 月至 年 月地震震中分布；其中蓝色曲线代表河流，棕色曲线代表已知断裂，黑色矩形代表水库大坝，三角形代表台站，其颜色表示不同的台网；图（）中白色五角星代表公元前 年至 年 级以上地震事件，红色五角星代表 年至 年 级以上地震事件；图（）、（）中白色圆代表地震震中，大小代表震级。图 研究区地震及台站分布 数据与方法 数据流域型梯级水库的水库地震监测台站沿江布设在 范围内（图 （），每个台网的台站均沿河分布且台网数据单独进行定位处理，线状分布的台站布局对于库区内外的地震监测和定位存在缺陷，也未能充分发挥对整个区域地震监测的作用。为了使得这些监测数据得到更为有效的利用，获取贵州西南区域更加完整和准确的地震目录，对这些台网记录的 年 月至 年 月的连续波形数据进行合并，利用 自动定位流程开展基于机器学习的库区地震识别和定位研究，获取了该区域 年 月 日至 年 月 日的地震目录（图 （）和震相观测报告，并基于获得的高精度定位结果分析该阶段黔西南水库库区的地震活动特征（，）。本文进一步采用波速比一致性约束的双差层析成像方法（，），利用获取的定位结果和震相观测报告进行地震位置和三维速度结构联合反演，研究黔西地区的地壳结构特征。为保证结果的可靠性，从中挑选至少被 个台站记录到的地震事件用于震源位置和速度结构的联合反演。此外，要求所使用的震相数据的震中距不超过 ，且根据走时曲线（图 ）删除偏差较大的震相数据后，参与反演的每个事件至少包括 个震相到时。经过筛选，共有 个地震事件和 个台站用于反演计算，其中 波震相 个，波震相 个，所使用的地震和台站分布如图 所示。方法流程本研究所使用的波速比一致性约束的双差层析成像方法是在传统双差层析成像方法（，）的基础上发展而来的。该方法可以同时使用 波、波以及 波和 波的走时差（）等 种数据的绝对走时和相对走时来联合反演研究区的地震震源位置和三 期杜兴忠等：贵州西部多流域水库区地壳三维速度结构图 波（）和 波（）震相走时曲线图 研究区内数据及反演网格节点分布维、和 结构，并在反演过程中对 模型除以 模型得到的波速比模型和由 数据直接反演得到的波速比模型添加一致性约束，使二者之间差异最小，从而得到更为精确的地下介质结构和震源位置（，；，）。在地震对匹配时，设置每个地震对之间的最大距离为 ，每个地震最多可以和 个地震组成地震对。基于数据分布情况，经测试将研究区内水平网格间距设为 （图 ），深度方向上的网格节点分别位于 、和 深度处。反演使用的最优阻尼因子和光滑因子通过 曲线方法（，）获得，分别为 和 （图 ）。反演使用的一维速度模型与定位时使用的速度模型相同。中国地震 卷图 确定最优阻尼因子（）和光滑因子（）的 曲线 结果可靠性分析采用检测板测试方法（，）评估反演得到的研究区三维速度模型的可靠性。首先，在联合反演所使用的一维速度模型的基础上，按照实际反演时采用的网格节点划分方式，对 波速度模型相邻的网格节点添加 的扰动，对波速比模型添加 的扰动。由于初始 波速度模型是由 波速度模型和波速比计算得到的，因此 波速度模型相邻网格节点的扰动量约为 和 。然后，基于该理论检测板模型和实际反演时使用的地震数据，采用伪弯曲射线追踪方法（，）生成一个理论的走时数据集。最后，采用与实际反演时相同的反演策略对合成的理论数据集进行反演，分析反演得到的模型与初始检测板模型的一致性。最终的检测板测试结果如图 所示，从图中可以看出，在 深度，对于研究区内地震和台站分布比较集中的区域，其、和 模型均能得到较好的恢复（红色曲线所包围的区域）。因此，下文中主要讨论这部分区域的三维速度结构。结果与讨论经过联合反演，最终得到 个地震的重定位结果，重定位前后的地震震中分布如图 所示，重定位前，地震在研究区内的分布比较离散（图 （），重定位后的地震分布明显收敛（图 （）。研究区内地震活动在北盘江的董箐水库大坝、光照大坝及其上游库区、晴隆断裂与北盘江交汇区域形成了地震丛集，在三岔河流域两侧分布着大量地震，在六冲河两侧远离河流的地区存在密集的地震分布。结合反演得到的三维、和 结构对该地区的地震活动进行进一步分析。不同深度水平方向上的、和 模型反演结果如图 所示。通常情况下，射线分布密集、检测板模型恢复较好的区域其速度结构反演结果比较可靠。根据研究区内射线分布（图 ）和检测板测试结果（图 ），本文仅展示速度结构反演结果较为可靠的区域，即图 中红色曲线所包围的区域。反演结果显示，研究区内的地震波速度结构具有明显的横向不均匀性，不同大地构造变形分区展现出不同的速度结构特征。总体来讲，在 、和 深度层，研究区中部六盘水断陷与遵义断拱分界带的 向威宁构造变形区（图 ）期杜兴忠等：贵州西部多流域水库区地壳三维速度结构注：红色曲线所包围的区域检测板模型恢复较好。图 （）、（）和 （）模型检测板测试结果注：黑线代表断层；红色实线为图 中剖面位置，红色和紫色虚线分别为二级、三级构造单元边界。图 重定位前（）、重定位后（）研究区内的地震震中分布 中国地震 卷 期杜兴忠等：贵州西部多流域水库区地壳三维速度结构注：白色圆圈为地震在距离最近的剖面上的垂直投影；蓝线代表河流；红线代表断层，：马场断裂，：开阳纳雍断裂，：威宁水城断裂，：盘县断裂，：晴隆断裂。图 研究区不同深度的（）、（）和 （）模型显示出显著而持续一致的低、特征（图 （）、（）。其西南侧的六盘水断陷在 和 深度层以低、为主，在 和 深度层显示出显著的高、特征。在其北东侧的遵义断拱在 和 深度层显示出低、特征，在 和 深度层显示低、高分布。在 深度层，波和 波速度均显示出较好的一致性和均匀性。在浅层地表附近（深度），夹岩、平寨和马马崖水库大坝所在位置表现为 波和 波低速异常，光照水库大坝及其上游区域也为低速异常区（图 （）、（）。对于浅层的波速比结构，夹岩、平寨、光照和马马崖水库大坝所在区域为高波速比区（图 （）。浅层的低速特征与横穿研究区的 条河流及碳酸盐岩为主的岩性特征有关。地质研究表明，研究区内碳酸盐岩广泛分布，六冲河干流上溶洞发育，夹岩水库库区伏流溶洞与明流河道呈相间分布（杨周等，）。溶洞、地下暗河等岩溶地貌是水体渗透扩散的有利通道（易立新等，；纪星星等，），流体沿通道向深部渗透，能够使地下介质的波速发生变化（曹颖等，）。平寨水库坝址区出露地层为三叠系夜郎组、永宁镇组、关岭组以及第四系堆积物，库区岩溶和裂隙十分发育，构成了复杂的岩溶 裂隙网络（简文星等，）。光照水库库区构造以短轴褶皱为主，断裂和岩溶比较发育（陈本金等，）。马马崖水库区的河谷边坡岩性为软岩、或上软下硬、或夹层发育，受构造切割等因素的影响，岸边多处分布有第四系覆盖层（张秀宏，）。因此，夹岩、平寨、光照和马马崖水库区域的低速、高波速比异常反映了该 中国地震 卷区域浅层的地质环境和流体渗透的影响。在 深度处，波低速异常主要分布在毕节大方纳雍一带和六盘水平寨晴隆一带，贞丰县西北侧和赫章县北侧则为显著的 波高速异常。波低速异常主要分布在六盘水 平寨晴隆一带，该区域同时具有较高的波速比。晴隆断裂两侧的速度结构表现出明显的差异。将研究区内不同深度的地震投影在距离最近的水平切片上（图 ），发现地震主要分布在赫章至纳雍之间的高低波速区的过渡带、平寨西北侧的高波速区和晴隆县附近低波速体的边缘。注：黑色圆圈为剖面两侧 范围内地震的垂直投影。图 不同深度剖面上、和 结构在 深度处，研究区显著的 波低速异常主要分布在六冲河两侧，并具有较低的波速比，而北盘江区域尤其是晴隆贞丰一带则为 波高速异常区和高波速比区。北盘江两侧的 波速度结构具有明显差异，其西南侧为 波高速区，东北侧为 波低速区。在 深度，晴隆贞丰一带仍然为显著的高、低 和高波速比区。整体上，研究区内的地震丛集主要位于高、低速异常体的边缘等速度结构变化剧烈的区域，对应于高波速比区或高、低波速比区的过渡带。为进一步分析研究区内的三维介质结构及地震活动，给出了穿过主要地震丛集区的 条深度剖面（剖面位置见图 （），其中 剖面（图 （）显示，夹岩水库西侧、六冲河南岸 向的地震丛集（图 （）其震源深度位于 ，呈近直立状分布，地震丛集南侧为低速区，北侧波速相对较高，地震丛集下方 以下深度存在明显的低波速比异常体。夹岩水库北侧的地震丛集由于受台站覆盖的影响，在深度上可靠性较差，分布较为离散。剖面（图 （）显示，北盘江和三岔河附近的 期杜兴忠等：贵州西部多流域水库区地壳三维速度结构地震震源深度同样较浅，主要集中在 以浅的深度。北盘江西南侧的地震丛集位于相对高速异常体内，在深度上具有倾向 的条带状特征。沿 剖面（图 （）地震主要位于北盘江光照水库、三岔河平寨水库和开阳纳雍断裂附近。其中，平寨水库位于研究区的中部，周围台站包围较好，重定位后地震收敛较好，在深度方向上呈现具有不同倾向和倾角的多条地震条带，反映了该区域的地震活动为多条小规模断层发生活动所致。平寨水库所在区域速度结构较为复杂，其南、北两侧速度结构和波速比具有明显差异。沿 剖面（图 （）地震主要集中在马马崖大坝附近的晴隆断裂上，位于高低速异常的交界处和低波速比区域。董箐大坝附近的地震活动受台站覆盖的影响，地震在深度上的分布比较离散。结论由于贵州地区的构造地震活动水平较弱，基于地震观测资料开展的地壳三维速度结构研究缺乏。本文使用基于机器学习方法构建的地震目录和观测报告，反演得到了该地区水平分辨率 、垂直分辨率 的三维、和 结构，并分析了地震的空间分布特征。研究区内分布有多个大型水库，受库区岩性和流体渗透的影响，浅层的速度结构普遍具有低波速和高波速比特征。上地壳速度结构呈现出与大地构造变形分区相关的横向不均匀性。重定位后的地震丛集勾勒出了大量隐伏断层的几何展布特征，其主要位于高、低速异常体的边缘等速度结构变化剧烈的区域，对应于高波速比区或高、低波速比区的过渡带上，推测该区域的地震活动与水库周围的断层活化有关。致谢：中国科学技术大学张海江教授提供 软件，文中图件使用 软件绘制，在此一并表示感谢。参考文献曹颖，黄江培，付虹 云南小湾水库蓄水后 波速度结构的双差地震层析成像研究 中国地震，（）：陈本金，欧品智，刘军 北盘江光照水库地震活动特征研究 科技创新导报，（）：胡毓良，陈献程 我国的水库地震及有关成因问题的讨论 地震地质，（）：纪星星，黄世源，罗浩，等 年 月 日贵州沿河 地震震害特征及发震构造分析 中国地震，（）：简文星，唐辉明，刘佑荣，等 平寨坝址左岸岩溶发育规律及其对建坝条件的影响 工程地质学报，（）：姜朝松，彭万里 贵州洪家渡电站水库诱发地震前期研究 地震研究，（）：孙浩淼 黔中水利枢纽平寨水库左岸岩溶溶蚀分析 水利与建筑工程学报，（）：王尚彦 贵州水库地震研究 贵州科学，（）：徐伟钧，黄波，范中学 贵州省大地构造特征分析 有色金属（矿山部分），（）：杨周，刘冬梅 岩溶地区水库设计洪水分析 以夹岩水利枢纽为例 水利科技与经济，（）：易立新，车用太 水库诱发地震及其水文地质条件和诱震机理 中国地质灾害与防治学报，（）：张秀宏 马马崖水电站堆积体对坝址选择的影响 贵州水力发电，（）：，（）：，（）：中国地震 卷 ，：，（）：，（）：，：，（）：，：，（）：，（）：），），），），），：；