2021, Vol. 37
2. 中国地震局地壳动力学重点实验室, 北京 100085;
3. 北京大学地球与空间科学学院, 北京 100871;
4. 四川省地震局, 成都 610041;
5. 中国地震局地球物理研究所, 北京 100081
2. Key Laboratory of Crustal Dynamics, China Earthquake Administration, Beijing 100085, China;
3. School of Earth and Space Sciences, Peking University, Beijing 100871, China;
4. Sichuan Earthquake Agency, Chengdu 610041, China;
5. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
地壳表层水循环系统中的不同类型水体，如大气降水、地表水和地下水具有不同的同位素组成特征(Clark and Fritz, 1997)。2H和18O作为水分子的组成部分，它们的稳定同位素是研究水的来源和水文循环过程的最有效的示踪工具(Craig，1961；Gat，1996)。温泉作为一种天然环保地热资源，是链接大气降水、地表水、地下水和地下深部流体的水热载体，是地球水循环过程中的重要组成部分，因此，对温泉开展氢氧稳定同位素研究，能够有效识别温泉的流体来源、水岩作用与水文循环过程等(周训等，2015；Li et al., 2018)，为温泉资源的开发与利用提供重要依据。
四川西部高原地区分布着大量隆起山地型地热资源，主要受断裂带分布控制(王贵玲，2018)。尤其是位于川西青藏高原东缘大型活动断裂带-鲜水河和安宁河断裂带上分布着大量中高温温泉。李晓等(2018)通过水化学和氢氧稳定同位素手段研究了鲜水河断裂带康定至道孚段的水热活动区热水的形成和联系，提出了不同水源补给的热水系统形成模式。Guo et al.(2017a)分析了康定高温地热水的水化学组成、气体组分和同位素组成，识别了影响热水组成的水文地球化学过程(包括混合和脱气作用)以及储层温度。Luo et al.(2017)利用地热温标、氢氧稳定同位素和水文地球化学模拟等手段评价了康定南北部地热开发潜力和优先顺序。也有学者通过对青藏高原东缘温泉的氦同位素、重力、磁和地震数据，获得了巴塘到康定区域的地下热结构和热液活动与地热动力学的关系(Tang et al., 2017)。所有这些研究为探讨青藏高原东缘温泉的热源和形成模式提供了依据。
鲜水河-安宁河断裂带及其周边大地震频发(Wen et al., 2008)，如发生在鲜水河断裂带上的1955年M7.5、1973年M7.6、1981年M6.9和2014年M6.3级地震。地震与构造活动可引起泉水物理和水化学性质发生变化，反之断裂带地下水尤其是温泉水中蕴含丰富的地下构造活动信息(赵永红等，2011；周晓成等，2020)，对其深入研究有助于评价构造活动对地下水资源的影响并可能获得可靠的地震信息(Skelton et al., 2014, 2019)。Chen et al.(2014)和Li et al.(2019)等发现康定区域温泉在大地震前后出现明显的氢氧稳定同位素变化，揭示了构造活动引起了断裂带含水层渗透性的增加或减小。氢氧稳定同位素技术为研究断裂带水循环过程及其与构造活动之间的关系提供了证据。
前人对青藏高原东缘鲜水河-安宁河断裂带温泉的氢氧稳定同位素研究多集中在鲜水河断裂带上康定区域，对鲜水河-安宁河断裂带上分布的众多温泉及地表水的氢氧稳定同位素的系统对比研究还较缺乏。在解决鲜水河-安宁河断裂带不同位置温泉的氢氧稳定同位素分布差异原因、区域地表水对地热水的补给等方面仍需深入开展相关工作，这对于深入了解活动断裂带条带分布地下水的循环特征有重要意义。为此本文以鲜水河-安宁河断裂带上的温泉和地表水作为一个系统研究对象，开展氢氧稳定同位素对比分析，同时结合温泉水中化学组分特征，初步揭示活动断裂带上地热水的同位素地球化学及水文循环特征。这对于研究青藏高原东缘断裂带上地热水的水文循环过程及其与构造活动的关系具有重要的意义。
1 研究区域概况研究区域(27.7°~31.2°N、101.1°~102.4°E，)主要位于四川西部道孚-康定-石棉-西昌一带，地处青藏高原东缘和四川盆地西缘(图 1)。气候类型主要为高原型季风气候和山地气候，降雨集中在每年5月至9月。鲜水河断裂带位于四川甘孜至石棉一带，为北西向的左旋走滑断裂带。安宁河断裂带位于四川石棉至西昌一带，以左旋走滑为主，近南北走向，长约160km。
|
图 1 研究区域位置与采样点分布示意图 Fig. 1 The study area and the sampling locations |
强烈的构造活动以及高山深切峡谷地形有利于温泉的形成，研究区域大多数温泉沿断裂呈带状分布(罗来麟，1994)。鲜水河断裂带上的温泉主要集中在道孚、中谷、康定和泸定磨西等4个热水区(罗来麟，1994)，多为高中温温泉，属于川西高原高中温地热区，热储岩性主要为花岗岩、石灰岩、大理岩、变质砂板岩以及第四系松散层胶结等(王贵玲，2018)。鲜水河断裂带地下冷水主要是基岩裂隙水，其补给为大气降水和高山冰雪融水(赵庆生，1984)。安宁河断裂带温泉主要分布在喜德和西昌周边，属于川西南中低温地热区，热储岩性一般为碳酸盐岩(罗来麟，1994；王贵玲，2018)。
2 样品采集与测试本研究于2019年5月共采集了鲜水河-安宁河断裂带上38个样品，包括温泉水样品24个，温泉附近河水12个，冷泉水样品1个，雪线附近积雪融水1个，采样位置见图 1和表 1。采样点命名中，S代表温泉水，S后数字代表采样点编号，R代表温泉附近的河流，R后数字为相应位置的温泉编号，CS1为冷泉水，M1为积雪融水(图 1、表 1)。在温泉主出水口处采集水样，共采集了中高温温泉(60℃≤T＜90℃)4处、中温温泉(40℃≤T＜60℃)16处和低温温泉(25℃≤T＜40℃)4处。采集温泉附近的流动河水样品，同时采集了道孚七美一处山顶雪线处积雪融水样品，康定一处冷泉样品。采样现场水样使用0.45μm过滤膜过滤，水样装满后密封保存在高密度聚乙烯塑料瓶内。在采样现场，使用AZ8821型号数字温度计测量水温。
|
表 1 采样点位置和氢氧稳定同位素组成 Table 1 The sampling locations and analyzed δ2H and δ18O data |
水样测试在中国地震局地壳动力学重点实验室完成。氢氧稳定同位素使用液态水同位素分析仪(LGR LWIA 912-0008)测量，结果以相对于维也纳标准平均海水(Vienna standard mean ocean water，V-SMOW)的千分差表示(δ2H和δ18O，‰)。同位素校准使用国家标准物质GBW04458、GBW04459和GBW04460。三种同位素标准物质的值分别为δ18OV-SMOW=-0.15‰和δ2HV-SMOW=-1.7‰；δ18OV-SMOW=-8.61‰和δ2HV-SMOW=-63.4‰；δ18OV-SMOW=-19.13‰和δ2HV-SMOW=-144.0‰。样品测试精度δ18O＜0.1‰，δ2H＜0.3‰。温泉和冷泉水离子质量浓度(Ca2+、Mg2+、Na+、K+、SO42-、Cl-、NO3-和F-)使用离子色谱仪(Dionex ICS-2100)测量，使用水中离子成分分析标准物质(GBW(E)080519~GBW(E)080529)，分别配置成阴离子(SO42-、Cl-、NO3-和F-)和阳离子(Ca2+、Mg2+、Na+、K+)的混合标准溶液，绘制标准曲线。HCO3-和CO32-采用滴定法测得。水中阴阳离子平衡相对误差(E)结果控制在±5%以内，为了保持数据的整体性，对于部分E超出±5%范围的离子数据亦列入表 2作为分析。微量元素(Li和Sr)使用电感耦合等离子体质谱(ICP-MS)(Thermo X Series II)测量，仪器测量误差＜5%。
|
|
表 2 温泉和冷泉水化学组成结果(mg/L) Table 2 The analyzed chemical compositions and related parameters (mg/L) |
研究区域温泉的海拔高度为1175~3946m，水温变化范围为37.29~74.81℃，水温平均值为52.02℃(表 1)。温泉最高水温为康定灌顶温泉(74.81℃)，最低水温为西昌河西温泉(37.29℃)。研究区域河水的海拔高度为1734~4197m，河水水温变化范围为6.48~23.82℃，水温平均值为12.66℃(表 1)，呈现出海拔高水温低、海拔低水温高的特征，这与采样点处的气温有关系。
研究区温泉、河水和积雪融水的δ18O和δ2H显示出差异性的分布特征(图 2)。温泉δ18O变化范围为-19.04‰~-12.71‰，平均值为-16.42‰；δ2H变化范围为-144.07‰~-88.63‰，平均值为-122.37‰(n=24)。在康定采集的1处冷泉(T=19.88℃)，δ18O和δ2H分别为-13.66‰和-106.74‰。研究区域河水的δ18O变化范围为-15.90‰~-10.85‰，平均值为-13.86‰；δ2H变化范围为-118.21‰~-71.12‰，平均值为-98.99‰(n=12)。雪线处积雪融水的δ18O和δ2H分别为-10.27‰和-65.41‰，相比温泉和河水样品富集重同位素。
|
图 2 水样δ18O-δ2H分布及其与全球和区域大气降水线的关系 Fig. 2 Plot of δ18O vs. δ2H and their correlations with GMWL and LMWL |
通过水样品的δ18O和δ2H与全球和区域大气降水线对比，能够了解水的来源与补给、水岩作用以及地下水混合作用等水文循环特征(Pang et al., 2017)。全球大气降水线方程(Global Meteoric Water Line，GMWL)为δ2H=8δ18O+10(Craig，1961)。研究区域位于青藏高原东缘，采用青藏高原大气降水线为区域大气降水线(Local Meteoric Water Line，LMWL)，降水线方程为δ2H=8.41δ18O+16.72(Kong et al., 2019)。δ18O和δ2H对比结果(图 2)显示出研究区域的大部分温泉和河水的δ18O和δ2H靠近GMWL和LMWL，说明补给源均为大气降水成因。但是不同类型采样点上同位素组成具有显著的差异性，如24个温泉样品具有4个明显的分布区域，河水R1~R17与R18~R23同样具有显著的同位素组成差异。上述氢氧稳定同位素的离散分布特征，表明了大气降水补给源的差别，这可能与不同位置大气降水同位素分布差异因素有关。冷泉水同位素呈现出向右侧偏离大气降水线的特征(图 2)。冷泉水的氢氧稳定同位素代表当地大气降水的平均同位素组成(王基华等，2000)，其向右偏离大气降水线的分布特征表明了补给过程中的蒸发作用引起了同位素的富集。
温泉与其附近的河流同位素值存在较大差异，温泉比其附近河流δ18O和δ2H值低(图 2)。这种特征同样出现在西藏地热田中，如羊八井地热田(郑淑蕙等，1982)。上述分析结果表明了近距离的河流和周边大气降水不是温泉的主要补给源，温泉的补给源为远距离的大气降水下渗经地下裂隙和管道运移补给。但是S17温泉与附近河流(R17)同位素值接近(图 2)，表明了该温泉的补给来源于附近的大气降水。这与该温泉为贡嘎山环绕的小区域气候环境有直接关系。
3.2 同位素高程效应研究区域位于青藏高原东缘的不同位置，采样点之间存在着明显的海拔高度差。前人研究表明，高程效应是影响大气降水同位素的重要因素，表现出随着海拔升高同位素值下降的特征(Clark and Fritz, 1997)。高程效应能够区分不同海拔处降水对含水层的补给(Blasch and Bryson, 2007)，例如Song et al. (2006)通过怀沙河流域内泉水的氧同位素高程效应确定了地下水的补给范围和水流路径。温泉的δ18O和δ2H与泉水出露高程之间存在着显著的负相关关系(图 3a, b)。虽然温泉均是上升泉，但是其补给为区域内大气降水(图 2)，而且沿着青藏高原东缘鲜水河-安宁河断裂带分布的研究区域温泉高程变化明显，因此温泉的同位素高程效应可能与大气降水的同位素高程效应有关。
|
图 3 温泉(a、b)和河水(c、d)采样点海拔与δ18O和δ2H的关系 Fig. 3 Plot of altitude vs. δ18O and δ2H for hot springs (a, b) and for rivers (c, d) |
研究区域温泉同位素分布特征受高程效应控制，温泉海拔与同位素值具有强负相关性(图 3)，同位素高程效应拟合公式为：
|
(1) |
|
(2) |
式中，h为温泉泉口出露高程(m)。因此得到氧同位素高程效应为-0.23‰/100m，氢同位素高程效应为-1.95‰/100m，即海拔每增加100m，δ18O值减少0.23‰，δ2H值减少1.95‰。Yu et al.(1984)通过对贵州-重庆-四川-西藏一线的水同位素研究，提出了大气降水氧同位素高程效应为-0.31‰/100m，氢同位素高程效应为-2.6‰/100m，大于本文研究区域同位素效应，这与其研究区域海拔高度差(4980m)明显大于本研究(2770m)有直接关系。本研究的同位素高程效应更适用于鲜水河-安宁河断裂带上的四川道孚-康定-石棉-西昌区域。温泉水同位素与温泉高程的强相关性(图 3)，表明了引起研究区域温泉水同位素离散分布的主要因素为高程变化。河水的δ18O与海拔之间的R2值表明了断裂带河水存在着高程效应(图 3c, d)，但是河水的高程效应R2值明显小于温泉，这可能与温泉的补给和水汽来源地稳定有关，而河水的补给存在着以大气降水为主的多种来源。
3.3 温泉氧同位素漂移地热水表现出两种典型的氧同位素分布特征，一类是氧同位素漂移并表现出水岩作用特征。由于补给热储的大气降水与岩石之间同位素的不平衡，会发生水与岩石矿物之间氧同位素的交换，导致地热水重同位素富集，氧同位素向大气降水线右侧漂移(顾慰祖等，2011)。氧同位素漂移可作为地热系统深部温度的判定指标(卫克勤等，1983)。如西藏羊八井地热水氧同位素偏离大气降水线2‰，云南热海地热田的偏移1‰~3‰，属于与火山和岩浆有关的地热田类型(卫克勤等，1983；Guo et al., 2017b)。另一种氧同位素靠近大气降水线，未发生明显的氧同位素漂移。本文报道的温泉氢氧稳定同位素靠近大气降水线的分布特征，反映出温泉的补给与循环速度快，并未发生明显的水和岩石氧同位素交换作用(除S15康定龙头沟温泉)。S15温泉泉口出露温度72.17℃，地下水循环过程中与岩石发生了氧同位素交换作用，显示出高温作用下的地热水的一定程度水岩作用的特征。
3.4 水化学特征及水岩作用温泉中离子质量浓度变化范围较大，如Ca2+为0.86~245.1mg/L，Mg2+为0.03~61.97mg/L，Na+为46.77~521.5mg/L，K+为0.11~81.19mg/L，HCO3-为42.70~1940mg/L，SO42-为7.33~ 256.7mg/L，Cl-为2.65~228.3mg/L(表 2)。温泉中微量元素Li质量浓度变化范围为0.03~4.81mg/L，平均值为1.52mg/L，最高Li质量浓度为道孚七美1#温泉，最低为西昌河西温泉。温泉中微量元素Sr质量浓度变化范围为0.04~2.80mg/L，平均值为1.05mg/L，最高Li质量浓度为道孚八美牦牛沟1#温泉，最低为西昌河西温泉(表 2)。离子和微量元素质量浓度的分布特征表明了地下水化学成分形成过程的差异性。
地下水的化学组成与水岩反应过程中发生的组分溶解/沉淀、水岩之间离子交换、氧化还原作用等化学过程密切相关，因此水化学组成被用来识别水岩相互作用(Redwan et al., 2016)。相关研究发现应用(Na++K+-Cl-)/SO42-和(Ca2++Mg2+)/HCO3-等水化学组分比值能够用来指示补给水的深浅部来源以及新鲜程度(Soltan，1998；Redwan et al., 2016)，Na-K-Mg三角图能够用来评价地下水为完成平衡水、部分平衡水或未成熟水(Giggenbach，1988)，以此指示水岩反应程度。大气成因指数(MGI)能够用来确定地下水来源，其计算方法为(Na++K+-Cl-)/SO42-，其中各离子浓度单位为meq/L。当计算结果MGI＜1时，指示地下水来源于深层大气降水渗滤，当MGI>1时，指示地下水来源于浅层大气水渗滤(Soltan，1998)。MGI结果显示大部分温泉和冷泉为浅层大气水渗滤类型。温泉和冷泉(Ca2++Mg2+)/HCO3-(各离子浓度单位为meq/L)值均＜1表明了泉水补给来源于新鲜的水或者大气降水来源对地下水的补给(表 2)。上述水化学离子指标说明温泉样品的补给为新鲜的大气水补给，进一步说明其水岩作用并不强烈。
在Na-K-Mg三角图上(图 4)，大部分温泉样品和冷泉样品落在Mg端元附近，表明了地下水的水岩作用尚未达到离子平衡状态，为循环相对较快的浅层冷水补给。S9、S21和S24温泉处在部分平衡水区域，但是靠近Mg端元，显示出了一定程度的水岩作用特征。地热水中的微量元素与深部流体活动密切相关，为深入研究影响流体组成的各种地质过程提供了依据(Kaasalainen et al., 2015)。S12和S13温泉与CS1冷泉同处康定同一区域，相隔距离近，温泉Li和Sr含量略大于冷泉含量(表 2)，显示出温泉水岩作用弱的特征。来自氢氧稳定同位素和水化学组成的结果，说明鲜水河-安宁河断裂带上的大部分温泉以大气降水补给为主且循环相对较快，小部分温泉存在着一定程度的水岩反应。
|
图 4 泉水样品Na-K-Mg三角图 Fig. 4 Na-K-Mg triangular diagram for springs |
同处康定区域的温泉S9~S16与冷泉CS1有着明显的同位素分布差异(图 2)，显示出温泉同位素值较冷泉值低的特征。这种现象与我国东部郯庐断裂带和胶辽断块区温泉与冷泉的同位素分布特征相同，即地热水的δ18O和δ2H值低于当地冷泉的相应同位素值，同位素差异与地下水深循环过程产生的同位素物理分馏有关(Shangguan et al., 1998)。Shangguan et al.(1998)发现地热水δ2H和δ18O的变化幅度与循环深度之间存在着强相关性，并能用来评价地热水的循环深度，计算方法如下：
|
(3) |
|
(4) |
式中，d为地热水最大循环深度(km)，Δδ18O为地热水与冷泉δ18O的差，Δδ2H为地热水与冷泉δ2H的差。公式(3)和(4)计算康定地区8口温泉的d(δ18O)和d(δ2H)平均值估算得到温泉循环深度为3.9~10.2km。使用上述方法估算的温泉循环深度与δ18O和δ2H具有强相关性(图 5)，表明随着温泉的循环深度增加，同位素值下降。由于研究区域多为高山峡谷区，温泉存在着强的高程效应(图 3)，温泉同位素与循环深度存在强相关性(图 5)，分析温泉与冷泉同位素的差异性可能与补给区高程以及地下水深浅循环过程中的同位素分馏有关。
|
图 5 温泉循环深度与δ18O和δ2H的对比关系 Fig. 5 Plot of circulation depths vs. δ18O and δ2H for hot springs |
除上述氢氧稳定同位素方法计算循环深度，还可以应用地热温标计算温泉的循环深度。常用的SiO2地热温标和阳离子地热温标等方法利用经验公式估算热储温度(王莹等，2007)。王逸凌等(2020)对鲜水河断裂带的温泉采用了SiO2地热温标、阳离子地热温标以及多矿物平衡温标综合得出了鲜水河断裂带不同区域的热储温度，与本文研究康定区域类似范围内的热储温度为180~ 250℃。如果温泉的温度受其循环所处的地热控制，可通过下式估算热水循环深度(李娟等，2007)：
|
(5) |
式中Z为循环深度(m)，G为地热增温率，取20.4m/℃(张云辉，2018)，TZ为热储温度(℃)，T0为补给区年平均气温，取7℃，Z0为年常温带深度，取30m。根据此式推算出康定区域温泉的循环深度在3.6~5.0km。利用氢氧稳定同位素计算的循环深度较利用热储温度方法计算的值偏大，可能与同位素的高程效应导致同位素值偏低有关。S9温泉的Cl-和Li和Si含量较低(表 2)，也说明其循环深度可能比氢氧稳定同位素计算的10.2km要小。该温泉的海拔明显较其他康定区域温泉高，可能同位素高程效应引起该温泉与浅循环的冷泉氢氧同位素差值较大，进而得到的循环深度偏高。对同一小范围内的温泉S12、S13和冷泉CS1进行计算得到的循环深度(图 5)也较利用公式(5)计算的偏高，可能与冷泉的蒸发效应引起的同位素值偏大，引起Δδ18O和Δδ2H变化范围偏大有关。因此，在利用δ2H和δ18O的变化幅度评价高山区域地热水的循环深度时，应考虑同位素高程效应以及冷泉的蒸发效应的影响。考虑到由不同经验公式估算循环深度可能存在不确定性，而且不同种方法之间存在着误差，康定区域温泉的循环深度可能由公式(5)估算的更加适用。
3.6 温泉成因与同位素监测意义鲜水河断裂带为下切到地幔的深大活动断裂带，温泉逸出气体中氦同位素和碳同位素结果显示出幔源流体作用(Zhou et al., 2015)，但是氢氧稳定同位素未显示出明显的向左或向右偏离大气降水线的水岩作用和二氧化碳气体作用特征(图 2)，而且水化学揭示的大部分温泉水岩作用特征不强烈。这可能与补给份额以大气降水补给为主有关，地下水循环速度快、滞留时间短，而且深部补给份额少不足以引起水的同位素变化。温泉和地表水氢氧稳定同位素的分布特征，表明了温泉的形成过程，即远距离补给的大气降水和地表水通过地下裂隙运移补给地下水，地下水在地下热储层加热后通过断裂通道上升到地表形成温泉，这符合王贵玲和蔺文静(2020)提出的川西温泉主要为经典的对流型水热系统成因类型特征。Tang et al.(2017)和Zhang et al.(2017)对川西巴塘-理塘-康定高温水热活动区进行了地球物理(地震、航磁和重力)和氦同位素的观测，结果显示青藏高原东缘地热系统的成因为大气降水和地表水沿断裂带渗入和深循环，被地壳热源加热后经循环上升到地表形成温泉，热源来自于逃逸的幔源热流、花岗岩产生的放射性热和断层摩擦热。Shi et al.(2017)对川西温泉的水化学聚类分析研究证明了康定区域温泉水化学成因来源于浅层水混和。氢氧稳定同位素结果与地球物理和水化学统计分析结果是一致的，而且氢氧稳定同位素更易识别补给源的远近及其与地表水的补给关系。
活动断裂带温泉氢氧稳定同位素监测是研究构造活动特征的重要手段之一。首先，氢氧稳定同位素揭示出了温泉的补给、循环与水岩作用特征，能够评价是否适宜进行地震观测(Zhang et al., 2020)。其次，深循环地下水的同位素变化是稳定的(Clark and Fritz, 1997)，近距离地表水的干扰较小，这在本研究温泉和河流同位素区别分布上得到证实。该性质对于排除小区域大气降水和河流补给等环境干扰对温泉氢氧稳定同位素的影响十分重要，有利于识别与构造活动有关的异常。再次，活动地震断裂带上温泉的氢氧稳定同位素值在强震活动前后会发生明显的变化(Skelton et al., 2014；Zhang et al., 2020)。已有的报道显示出2008年汶川地震后，鲜水河-安宁河断裂带上部分温泉出现明显的同位素变化，比如位于康定区域的二道桥温泉(S14)、龙头沟温泉(S15)和灌顶温泉(S16)发生了明显的氧同位素偏移(Chen et al., 2014；Li et al., 2019)。同位素变化与地壳应力变化引起含水层性质变化有关，如含水层渗透性发生变化，致使不同类型的地下水混合或者补给比率发生变化，或者前震与同震促使含水层发生破裂，水与新鲜裂隙面接触增加了水岩反应程度(Thomas，1988；Claesson et al., 2007)。那么，定期开展活动断裂带上多口温泉的氢氧稳定同位素背景变化观测，对于了解构造活动强烈的鲜水河-安宁河断裂带的地下水循环特征、断层/含水层破裂与愈合以及大地震的孕育过程有重要意义。
4 结论(1) 鲜水河-安宁河断裂带温泉氢氧稳定同位素空间变化明显，其δ18O变化范围为-19.04‰~-12.71‰，平均值为-16.42‰，δ2H变化范围为-144.07‰~-88.63‰，平均值为-122.37‰。氢氧稳定同位素主要分布在全球和区域大气降水线上，表明了温泉的大气降水成因。温泉水缺少明显的氧同位素漂移现象，表明了大部分温泉弱的水岩作用特征，并在相关离子比值、Na-K-Mg三角图、Li和Sr元素等指标上得到证实。
(2) 研究区域温泉同位素值具有明显的同位素高程效应，温泉同位素与海拔的拟合关系式为δ18O=-0.0023h-10.011，δ2H=-0.0195h-67.813。氧同位素高程效应为-0.23‰/100m，氢同位素高程效应为-1.95‰/100m。
(3) 温泉与其附近地表水的δ18O和δ2H具有明显的差异性，表明了温泉的补给源为远距离的大气降水。氢氧稳定同位素特征、水岩作用特征和循环深度表明温泉的成因为远距离大气降水运移补给地下水，地下水在地下热储层加热后通过断裂通道上升到地表形成温泉。
致谢 衷心感谢两位审稿专家对本文提出的宝贵意见！
谨以此文祝贺沈其韩院士百岁华诞，祝沈先生健康长寿！
Blasch KW and Bryson JR. 2007. Distinguishing sources of ground water recharge by using δ2H and δ18O. Groundwater, 45(3): 294-308 DOI:10.1111/j.1745-6584.2006.00289.x
|
Chen Z, Du JG, Zhou XC, Yi L, Liu L, Xie C, Cui YJ and Li Y. 2014. Hydrochemistry of the hot springs in western Sichuan Province related to the Wenchuan MS 8. 0 earthquake. The Scientific World Journal, 2014: 901432
|
Claesson L, Skelton A, Graham C and Mörth CM. 2007. The timescale and mechanisms of fault sealing and water-rock interaction after an earthquake. Geofluids, 7(4): 427-440 DOI:10.1111/j.1468-8123.2007.00197.x
|
Clark ID and Fritz P. 1997. Environmental Isotopes in Hydrogeology. Boca Raton: CRC Press
|
Craig H. 1961. Isotopic variations in meteoric waters. Science, 133(3465): 1702-1703 DOI:10.1126/science.133.3465.1702
|
Gat JR. 1996. Oxygen and hydrogen isotopes in the hydrologic cycle. Annual Review of Earth and Planetary Sciences, 24: 225-262 DOI:10.1146/annurev.earth.24.1.225
|
Giggenbach WF. 1988. Geothermal solute equilibria. Derivation of Na-K-Mg-Ca geoindicators. Geochimica et Cosmochimica Acta, 52(12): 2749-2765 DOI:10.1016/0016-7037(88)90143-3
|
Gu WZ, Pang ZH, Wang QJ and Song XF. 2011. Isotope Hydrology. Beijing: Science Press (in Chinese)
|
Guo Q, Pang ZH, Wang YC and Tian J. 2017a. Fluid geochemistry and geothermometry applications of the Kangding high-temperature geothermal system in eastern Himalayas. Applied Geochemistry, 81: 63-75 DOI:10.1016/j.apgeochem.2017.03.007
|
Guo QH, Liu ML, Li JX, Zhang XB, Guo W and Wang YX. 2017b. Fluid geochemical constraints on the heat source and reservoir temperature of the Banglazhang hydrothermal system, Yunnan-Tibet Geothermal Province, China. Journal of Geochemical Exploration, 172: 109-119 DOI:10.1016/j.gexplo.2016.10.012
|
Kaasalainen H, Stefánsson A, Giroud N and Arnórsson S. 2015. The geochemistry of trace elements in geothermal fluids, Iceland. Applied Geochemistry, 62: 207-223 DOI:10.1016/j.apgeochem.2015.02.003
|
Kong YL, Wang K, Li J and Pang ZH. 2019. Stable isotopes of precipitation in China: A consideration of moisture sources. Water, 11(6): 1239 DOI:10.3390/w11061239
|
Li B, Shi ZM, Wang GC and Liu CL. 2019. Earthquake-related hydrochemical changes in thermal springs in the Xianshuihe Fault zone, western China. Journal of Hydrology, 579: 124175 DOI:10.1016/j.jhydrol.2019.124175
|
Li J, Zhou X, Fang B and Yang YX. 2007. Formation of hot springs at Wenquanpu, Qinhuangdao City, Hebei, China, and suggestions of their exploitation and utilization. Geological Bulletin of China, 26(3): 344-349 (in Chinese with English abstract)
|
Li JX, Yang G, Sagoe G and Li YL. 2018. Major hydrogeochemical processes controlling the composition of geothermal waters in the Kangding geothermal field, western Sichuan Province. Geothermics, 75: 154-163 DOI:10.1016/j.geothermics.2018.04.008
|
Li X, Wang JJ, Huang X and Yu ZY. 2018. Chemical and isotopic characteristics of hot water in the Kangding-Daofu section of Xianshuihe fault zone, Sichuan, China. Journal of Chengdu University of Technology (Science & Technology Edition), 45(6): 733-745 (in Chinese with English abstract)
|
Luo J, Pang ZH, Kong YK and Wang YC. 2017. Geothermal potential evaluation and development prioritization based on geochemistry of geothermal waters from Kangding area, western Sichuan, China. Environmental Earth Sciences, 76(9): 343 DOI:10.1007/s12665-017-6659-9
|
Luo LL. 1994. Inquisition of the distribution and cause of the hot springs in western Sichuan. Journal of Chongqing Teachers College (Natural Science Edition), 11(2): 39-47 (in Chinese with English abstract)
|
Pang ZH, Kong YL, Li J and Tian J. 2017. An isotopic geoindicator in the hydrological cycle. Procedia Earth and Planetary Science, 17: 534-537 DOI:10.1016/j.proeps.2016.12.135
|
Redwan M, Abdel Moneim AA and Amra MA. 2016. Effect of water-rock interaction processes on the hydrogeochemistry of groundwater west of Sohag area, Egypt. Arabian Journal of Geosciences, 9(2): 111 DOI:10.1007/s12517-015-2042-x
|
Shangguan ZG, Du JK, Zang W, Wang JH, Kong LC and Gao SS. 1998. Modern hot spring geochemistry at the Tanlu fault and Jiaoliao block in eastern China. Science in China (Series D), 41(1): 87-94 DOI:10.1007/BF02932426
|
Shi ZM, Liao F, Wang GC, Xu QY, Mu WQ and Sun XY. 2017. Hydrogeochemical characteristics and evolution of hot springs in eastern Tibetan Plateau geothermal belt, western China: Insight from multivariate statistical analysis. Geofluids, 2017: 6546014
|
Skelton A, Andrén M, Kristmannsdóttir H, Stockmann G, Mörth CM, Sveinbjörnsdóttir Á, Jónsson S, Sturkell E, Guorúnardóttir HR, Hjartarson H, Siegmund H and Kockum I. 2014. Changes in groundwater chemistry before two consecutive earthquakes in Iceland. Nature Geoscience, 7(10): 752-756 DOI:10.1038/ngeo2250
|
Skelton A, Liljedahl-Claesson L, Wästeby N, Andrén M, Stockmann G, Sturkell E, Mörth CM, Stefansson A, Tollefsen E, Siegmund H, Keller N, Kjartansdóttir R, Hjartarson H and Kockum I. 2019. Hydrochemical changes before and after earthquakes based on long-term measurements of multiple parameters at two sites in northern Iceland: A review. Journal of Geophysical Research: Solid Earth, 124(3): 2702-2720 DOI:10.1029/2018JB016757
|
Soltan ME. 1998. Characterisation, classification, and evaluation of some ground water samples in Upper Egypt. Chemosphere, 37(4): 735-745 DOI:10.1016/S0045-6535(98)00079-4
|
Song XF, Liu XC, Xia J, Yu JJ and Tang CY. 2006. A study of interaction between surface water and groundwater using environmental isotope in Huaisha River basin. Science in China (Series D), 49(12): 1299-1310 DOI:10.1007/s11430-006-1299-z
|
Tang XC, Zhang J, Pang ZH, Hu SB, Tian J and Bao SJ. 2017. The eastern Tibetan Plateau geothermal belt, western China: Geology, geophysics, genesis, and hydrothermal system. Tectonophysics, 717: 433-448 DOI:10.1016/j.tecto.2017.08.035
|
Thomas D. 1988. Geochemical precursors to seismic activity. Pure and Applied Geophysics, 126(2-4): 241-266 DOI:10.1007/BF00878998
|
Wang GL. 2018. Geothermy of China. Beijing: Science Press (in Chinese)
|
Wang GL and Lin WJ. 2020. Main hydro-geothermal systems and their genetic models in China. Acta Geologica Sinica, 94(7): 1923-1937 (in Chinese with English abstract)
|
Wang JH, Lin YW, Liu CL and Liu WZ. 2000. Evidence to the hydrogen-oxygen stable isotopes and gas composition formed from hot springs in southern Zhangjiakou. Hydrogeology and Engineering Geology, 27(4): 30-33 (in Chinese with English abstract)
|
Wang Y, Zhou X, Yu Y, Liu CH and Zhou HY. 2007. Application of geothermometers to calculation of temperature of geothermal reservoirs. Geoscience, 21(4): 605-612 (in Chinese with English abstract)
|
Wang YL, Mu WQ and Shi ZM. 2020. Hydrogeochemical characteristics of hot springs in south segment of the Xianshuihe fault zone. Journal of Seismological Research, 43(2): 287-295 (in Chinese with English abstract)
|
Wei KQ, Lin RF and Wang ZX. 1983. Hydrogen and oxygen stable isotopic composition and tritium content of waters from Yangbajain geothermal area, Xizang, China. Geochimica, (4): 338-346 (in Chinese with English abstract)
|
Wen XZ, Ma SL, Xu XW and He YN. 2008. Historical pattern and behavior of earthquake ruptures along the eastern boundary of the Sichuan-Yunnan faulted-block, southwestern China. Physics of the Earth and Planetary Interiors, 168(1-2): 16-36 DOI:10.1016/j.pepi.2008.04.013
|
Yu JS, Zhang HB, Yu FJ and Liu DP. 1984. Oxygen and hydrogen isotopic compositions of meteoric waters in the eastern part of Xizang. Geochemistry, 3(2): 93-101
|
Zhang J, Li WY, Tang XC, Tian J, Wang YC, Guo Q and Pang ZH. 2017. Geothermal data analysis at the high-temperature hydrothermal area in Western Sichuan. Science China (Earth Sciences), 60(8): 1507-1521 DOI:10.1007/s11430-016-9053-2
|
Zhang L, Guo LS, Wang Y, Liu DY, Liu YW and Li J. 2020. Continuous monitoring of hydrogen and oxygen stable isotopes in a hot spring: Significance for distant earthquakes. Applied Geochemistry, 112: 104488 DOI:10.1016/j.apgeochem.2019.104488
|
Zhang YH. 2018. Research on genesis and development of the geothermal system in the Kangding-Moxi segment of the Xianshuihe fault. Ph. D. Dissertation. Chengdu: Chengdu University of Technology (in Chinese with English summary)
|
Zhao QS. 1984. The hydrogeochemical characteristics and forming model of hot water in the Xianshuihe fracture zone. Journal of Chengdu University of Science and Technology, (2): 77-88 (in Chinese with English abstract)
|
Zhao YH, Bai JT, Li XF, Jia K and Chen H. 2011. Correlation between hydrogen isotope in underground water near active fault and earthquakes. Acta Petrologica Sinica, 27(6): 1909-1915 (in Chinese with English abstract)
|
Zheng SH, Zhang ZF, Ni BL, Hou FG and Shen MZ. 1982. Hydrogen and oxygen isotopic studies of thermal waters in Xizang. Acta Scientiarum Naturalium Universitatis Pekinensis, (1): 99-106 (in Chinese with English abstract)
|
Zhou X, Cao Q, Yin F, Guo J, Wang XC, Zhang YS, Wang LD and Shen Y. 2015. Characteristics of the brines and hot springs in the Triassic carbonates in the high and steep fold zone of the eastern Sichuan basin. Acta Geologica Sinica, 89(11): 1908-1920 (in Chinese with English abstract)
|
Zhou XC, Wang WC, Chen Z, Yi L, Liu L, Xie C, Cui YJ, Du JG, Cheng JW and Yang LM. 2015. Hot spring gas geochemistry in western Sichuan Province, China after the Wenchuan Ms 8. 0 earthquake. Terrestrial, Atmospheric and Oceanic Sciences, 26(4): 361-373 DOI:10.3319/TAO.2015.01.05.01(TT)
|
Zhou XC, Wang WL, Li LW, Hou JM, Xing LT, Li ZP, Shi HY and Yan YC. 2020. Geochemical features of hot spring gases in the Jinshajiang-Red River fault zone, Southeast Tibetan Plateau. Acta Petrologica Sinica, 36(7): 2197-2214 (in Chinese with English abstract) DOI:10.18654/1000-0569/2020.07.18
|
顾慰祖, 庞忠和, 王全九, 宋献方. 2011. 同位素水文学. 北京: 科学出版社.
|
李娟, 周训, 方斌, 杨燕雄. 2007. 河北秦皇岛市温泉堡温泉的形成与开发利用建议. 地质通报, 26(3): 344-349. DOI:10.3969/j.issn.1671-2552.2007.03.012 |
李晓, 王金金, 黄珣, 余中友. 2018. 鲜水河断裂带康定至道孚段热水化学与同位素特征. 成都理工大学学报(自然科学版), 45(6): 733-745. DOI:10.3969/j.issn.1671-9727.2018.06.09 |
罗来麟. 1994. 四川西部温泉分布及成因初探. 重庆师范学院学报(自然科学版), 11(2): 39-47. |
王贵玲. 2018. 中国地热志. 北京: 科学出版社.
|
王贵玲, 蔺文静. 2020. 我国主要水热型地热系统形成机制与成因模式. 地质学报, 94(7): 1923-1937. DOI:10.3969/j.issn.0001-5717.2020.07.002 |
王基华, 林元武, 刘成龙, 刘五洲. 2000. 张家口南部地区温泉形成的氢氧稳定同位素及气体组成证据. 水文地质工程地质, 27(4): 30-33. DOI:10.3969/j.issn.1000-3665.2000.04.009 |
王莹, 周训, 于湲, 柳春晖, 周海燕. 2007. 应用地热温标估算地下热储温度. 现代地质, 21(4): 605-612. DOI:10.3969/j.issn.1000-8527.2007.04.003 |
王逸凌, 穆文清, 史浙明. 2020. 鲜水河断裂南段温泉水文地球化学特征. 地震研究, 43(2): 287-295. DOI:10.3969/j.issn.1000-0666.2020.02.010 |
卫克勤, 林瑞芬, 王志祥. 1983. 西藏羊八井地热水的氢、氧稳定同位素组成及氚含量. 地球化学, (4): 338-346. DOI:10.3321/j.issn:0379-1726.1983.04.002 |
张云辉. 2018. 鲜水河断裂康定-磨西段地热系统成因及开发利用研究. 博士学位论文. 成都: 成都理工大学
|
赵庆生. 1984. 鲜水河断裂带热水水文地球化学特征及形成模式. 成都科技大学学报, (2): 77-88. |
赵永红, 白竣天, 李小凡, 贾科, 陈辉. 2011. 活动断裂带附近地下水中的氢同位素变化与地震关系研究. 岩石学报, 27(6): 1909-1915. |
郑淑蕙, 张知非, 倪葆龄, 侯发高, 沈敏子. 1982. 西藏地热水的氢氧稳定同位素研究. 北京大学学报, (1): 99-106. |
周训, 曹琴, 尹菲, 郭娟, 王晓翠, 张永帅, 王黎栋, 沈晔. 2015. 四川盆地东部高褶带三叠系地层卤水和温泉的地球化学特征及成因. 地质学报, 89(11): 1908-1920. |
周晓成, 王万丽, 李立武, 侯建民, 邢蓝田, 李中平, 石宏宇, 颜玉聪. 2020. 金沙江-红河断裂带温泉气体地球化学特征. 岩石学报, 36(7): 2197-2214. |