机器学习揭示的2010年Ms 7.1玉树地震带断层控制的水文地球化学空间演化

《Journal of Hydrology: Regional Studies》：Fault-controlled hydrogeochemical spatial evolution revealed by machine learning around the 2010 Ms 7.1 Yushu earthquake zone

　　

研究区域：甘孜-玉树断裂带（GYFZ）是青藏高原东部一条活动的走滑断层，分隔了巴颜喀拉地块和羌塘地块。研究焦点：活动走滑断层可同时充当地下水的屏障和通道，因此成为理解地震变形如何重新分配大气降水补给和深部挥发性组分的重要构造环境。在GYFZ中，沿断层走向的耦合状态差异、2010年玉树地震破裂以及非均质岩性，要求厘清断层构造如何控制地下水循环和水化学分异。区域水文地质新认识：来自31个水样品和9个气体采样点的水化学、水同位素、气体同位素以及无监督机器学习证据（主成分分析（PCA）、K-means聚类分析和均匀流形逼近与投影（UMAP））表明，GYFZ呈分隔特征。稳定同位素指示大气降水补给，但循环过程受力学状态控制。高度耦合的锁定段保留了稀释的大气降水补给簇2水体，并限制了与深部流体的垂直交换。破裂段和耦合转换域则为更深、更热的簇1流体提供了通道，产生了SiO₂、Na⁺、Cl^-和HCO₃^-的富集以及变化的挥发性组分输入（CO₂ = 1.02–99.04%；Rc = 0.10–0.51 Ra）。簇3由样品GQNY代表，记录了流经蒸发岩和碳酸盐岩地层的局部硫酸盐富集流动路径；其SO₄^2-富集反映的是矿物可用性和水文连通性，而非CO₂驱动的酸性条件。因此，断层耦合决定了通道-屏障行为，而岩性控制了反应指纹，从而解释了沿断层分布的局部挥发性异常。

基于断层带水文地球化学与机器学习耦合研究的论文解读

## 一、研究背景与科学问题

活动断裂带是结构非均质的地下水系统，变形作用、渗透率生成、挥发性组分运移和水-岩相互作用在广泛深度范围内相互耦合。在青藏高原这类大陆碰撞带中，深根走滑断层可能将浅层大气降水循环与更深部的热源和岩石圈挥发性储库连接起来。尽管已有研究认识到断层能够输运流体，但尚未解决的关键问题是：力学性质和岩性均非均质的断层带，如何将浅层补给、深部含挥发性组分的流体以及反应路径，划分为空间上不同的水化学模式。甘孜-玉树断裂带（GYFZ）因其沿走向的构造和力学差异明显而成为验证这一问题的理想场所。该断裂构成巴颜喀拉与羌塘地块的边界，被解释为切至约85 km深的莫霍面，并包含高度耦合的锁定段与低耦合域。2010年Ms 7.1玉树地震、孔隙流体压力升高的证据以及密集的破裂相关裂隙网络，表明锁定域、破裂影响区和耦合转换区的渗透率可能存在显著差异。此外，岩性非均质性意味着硫酸盐、碳酸盐和硅酸盐相关溶质不能仅作为单一区域来源解释。

本研究旨在回答三个具体问题：（i）GYFZ沿线的泉水能否被划分为具有水文意义的水化学组；（ii）这些组是否占据不同的构造域；（iii）δD-δ¹⁸O同位素、CO₂含量、He-Ne同位素特征和主要离子组合如何指示大气降水补给、深部挥发性组分输入和依赖岩性的水-岩相互作用。

## 二、研究数据与技术方法

研究人员于2024年8月沿甘孜-玉树断裂带开展泉流体地球化学调查。共采集13个水化学样品，其中9个点位同步采集冒泡气体样品，并整合了已发表数据集中18个热泉的水化学数据。所用主要方法包括：主成分分析（PCA），用于提取主要地球化学梯度并降低数据维度；K-means聚类算法，用于依据化学成分相似性划分水样组别；均匀流形逼近与投影（UMAP），用于非线性可视化并验证聚类结果的拓扑分离。研究通过Z-score标准化对水化学数据进行预处理，采用WCSS肘部法则和轮廓系数双重验证确定最优聚类数。水源样品测试分析在中国地震局地震预测研究所重点实验室完成，气体样品测试在中国科学院西北生态环境资源研究院油气资源研究中心完成。

## 三、主要研究结果

### 3.1 水化学和稳定同位素特征

研究区流体温度介于4.8–58.9 °C，pH值为6.3–8.3，电导率（EC）为343.9–3207.0 μS/cm。主要阳离子为Na⁺和Ca²⁺，阴离子以HCO₃^-为主。依据舒卡列夫分类，流体样品均为重碳酸盐型水，可划分为六种水化学类型。稳定同位素δ¹⁸O_V-SMOW介于-17.4‰至-13.1‰，δD_V-SMOW介于-132.7‰至-100.2‰，样品点沿全球大气降水线（GMWL）和当地大气降水线（LMWL）分布，指示大气降水为主要补给来源。

### 3.2 CO₂含量与He-Ne同位素组成

气体样品中CO₂含量变化极大，介于1.02%至99.04%之间。样品的⁴He/²⁰Ne比值在1.02–119.83之间，空气校正后的氦同位素比值（Rc）高度集中，介于0.10–0.51 Ra之间。

### 3.3 水化学簇的定量划分与演化

通过双验证策略确定了K = 3的最优聚类数。PCA结果显示前两个主成分共解释85.7%的方差。PC1（60.8%）作为水-岩相互作用强度和流体矿化程度的定量指标，PC2（24.9%）用于区分不同岩性风化路径。Cluster 2占据PC1负域，代表未演化的浅层大气降水补给水；Cluster 1跨越PC1正域，与SiO₂、Na⁺和Cl^-载荷向量密切相关，代表经历了较强水-岩相互作用的深循环较热流体；Cluster 3由样品GQNY单独构成，以极高的Ca²⁺和SO₄^2-浓度为特征，被解释为局部水化学极端，可能反映了流经石膏质蒸发岩透镜体和碳酸盐岩地层的局地循环。

### 3.4 非线性拓扑验证与地球化学指纹

UMAP投影展示了三个K-means簇在非线性拓扑空间中的清晰分离，验证了分类结果并非线性投影的数学假象。Z-score热图显示：Cluster 1富集SiO₂、碱金属和Cl^-，与碳酸增强的硅酸盐风化反应路径一致；Cluster 2在所有变量上均表现出相对亏损；Cluster 3以SO₄^2-（Z-score 4.19）、Ca²⁺（3.86）和Mg²⁺（3.10）的强正异常为特征，其硫酸盐富集主要归因于石膏等蒸发岩矿物的溶解。

### 3.5 深部流体上涌的构造控制

空间分析表明，Cluster 1深循环热流体和Cluster 3局部异常集中在2010年Ms 7.1主震震中附近，并沿北北西走向的线性轨迹延伸。这些异常与CO₂浓度超过80%以及区域氦同位素比值（Rc）为0.2–0.5 Ra的区域相吻合，指示深部地壳和岩石圈来源的挥发性组分上涌。Cluster 2样品则广泛分布于高度耦合的锁定段，这些段落作为低渗透率屏障限制了深部流体的垂直迁移。

## 四、讨论与结论

研究人员提出了一个概念模型：甘孜-玉树断裂带作为一个分隔的地下水系统，局部将深部含挥发性组分的流体与地表连接。高度耦合的锁定段充当低渗透率屏障，其水化学组成以未演化的浅层大气降水补给（Cluster 2）为主；而高-低耦合转换边界和地震破裂段则作为优势垂直通道，地震事件打破了锁定凹凸体的力学完整性，使次级裂隙网络扩张，形成类断层阀行为，促进深部地壳流体和岩石圈CO₂的聚焦上涌。上涌流体的化学演化路径取决于所遇岩性：穿经深部硅酸盐岩层时，在较高温度下发生更强的硅酸盐-水相互作用，形成Cluster 1的SiO₂、Na⁺和Cl^-富集特征；而当构造聚焦流路径穿经局部蒸发岩和碳酸盐岩透镜体时，则产生Cluster 3的Ca²⁺、Mg²⁺和SO₄^2-特征。

研究的主要结论为：GYFZ的水化学信号可被划分为三个具有物理意义的簇，其中前两个主成分解释了85.7%的方差；不同簇记录了不同的控制过程，Cluster 2受大气降水补给和稀释控制，Cluster 1受深循环、较热硅酸盐-水相互作用和沿渗透性断层通道上涌控制，Cluster 3反映局部流动路径穿经蒸发岩和碳酸盐岩地层；断层耦合和破裂构造提供了主要的本文控制作用，高度耦合锁定段为浅层大气降水循环主导的低渗透率域，而地震破裂域和耦合转换边界为深部含挥发性组分流体提供优势垂直通道。这种分隔构造解释了深部过程如何以局部地表水化学异常的形式在GYFZ中表现。该研究发表于《Journal of Hydrology: Regional Studies》。