论文部分内容阅读
摘要：对地震活动强烈地区进行水文地球化学监测可以很好地获取与地震相关的深部流体地球化学信息。通过研究龙门山断裂带及其周围13个温泉及观测井的水文地球化学特征，建立该断裂带温泉水文循环模型，揭示其水化学变化与地震活动的关系。结果表明：①龙门山断裂带温泉水主要来自其周围0.8～3.2 km高山的大气降水;②温泉水化学类型沿龙门山断裂从北到南、从西到东依次为重碳酸型、硫酸型、氯化物型，温泉水的循环深度、水岩反应程度及微量元素富集因子不断增加;③在汶川MS8.0地震、芦山MS7.0地震发生后，震中距200 km以内的温泉的离子组分呈下降趋势，这可能是在震后愈合过程中，龙门山断裂内部渗透性减弱、水岩反应程度衰减所致。
关键词：温泉;水文地球化学;地震;龙门山断裂带
中图分类号：P314.1 文献标识码：A 文章编号：1000-0666（2021）02-0170-15
0 引言
温泉多沿断裂带出露，是沟通地壳深部与浅部的窗口，和深部地下水系统有着密切的关系（Zhang et al，2015;Zhou et al，2017）。研究活动断层的温泉流体地球化学特征有助于深入了解该断层的活动方式、活动强度以及切割深度等断层运动特性，对探索断层带深部、特别是震源区附近介质的物理化学环境、研究各类地球化学组分的地震前兆机制、寻找特定地区的前兆观测组分有重要意义（Lambrakis，Kallergis，2005;Liu et al，2009）。因此，在地震易发区对温泉的水文地球化学特征进行监测，可捕捉来自地壳深部活动的信息（蒲小武等，2013;王云等，2014;Walraevens et al，2018;Xu et al，2019）。
近年来，国内外许多专家对温泉动态与地震活动的关系进行了研究（Tsunogai，Wakita，1995;Shi et al，2020）。Skelton等（2014）对冰岛2012年10月以及2013年4月2次5级地震之间以及地震前后4～6个月的地下水中稳定同位素比率和钠、钙等溶质的浓度变化规律进行研究，获得了该地区地震前后地下水水文地球化学变化的规律;Okuyama等（2016）研究发现1965—1968年日本松代地震群前后水中某些溶解组分随井深呈线性增长;Pérez 等（2008）研究1995年西班牙加利西亚4.6级地震以及1997年5.3级地震发现，在震中距为90 km的地下水观测井中Cl-浓度明显增加。另外，Shi等（2020）的研究显示，氢氧同位素以及微量元素的变化对地震的响应也较为明显，2018年通海MS5.0地震后，江川观测井的水化学以及氢氧同位素明显上升，微量元素（Li+、Pb2+、As2+等元素）的浓度呈明显下降。
龙门山断裂带是青藏高原东缘山脉与扬子板块的结合部位，地质构造复杂，地震活动剧烈（Lin et al，2018）。自2008年汶川MS8.0地震后，我国在该地区建立了大量的流动观测网、温泉和热水井固定观测点，丰富了该地区的流体地球化学背景（程万正等，2013）。目前，关于龙门山断裂带地下流体的研究主要集中在地震前后温泉水温、水化学组分、水氡以及温泉逸出气异常等方面（Chen et al，2015;晏锐等，2015;Zhou et al，2015;杨耀等，2019）。Chen等（2015）对2008—2010年川西地区32个温泉的水化学变化进行了研究，发现了汶川地震后温泉的离子浓度呈下降变化;Zhou等（2015）对该地区32个温泉气进行了3次采样，结果表明汶川MS8.0地震后，四川西部地区有更多地幔流體进入地壳，地幔成分的流体上涌可能与大地震的发生有关系;杨耀等（2019）研究了川西地区地下流体观测井的地下水的水化学特征，确定了四川西部地区主要观测井的水文地球化学背景。
前人对龙门山断裂带温泉水文地球化学特征的研究主要集中于川西地区附近，没有对整个龙门山断裂带展开更为细致的分析，鉴于此，本文对龙门山断裂带13个温泉及观测井的水文地球化学特征进行分析，针对该断裂带不同断裂计算温泉水的热储温度与循环深度，确定温泉水的补给来源和补给高程，判断其水文地球化学的控制因素和影响因素，建立龙门山断裂带温泉水循环模型。通过探讨2008—2016年龙门山断裂带温泉水文地球化学变化与汶川、芦山地震活动的关系，为该地区震情跟踪、异常落实和地震预测提供基础研究。
1 研究区地质概况
龙门山断裂带是青藏高原东缘的边界山脉，北起广元、南至天全，长约500 km、宽约30～50 km，呈NE—SW向展布，形成于印支晚期，受多期复杂构造运动影响（邓起东等，1994;Clark，Royden，2000;Godard et al，2009;Kouketsu et al，2017）。该构造带是由一系列大致平行的逆冲、逆冲兼走滑的断裂构成，自西向东发育有茂汶—汶川断裂（F1）、北川—映秀断裂（F2）、安县—灌县断裂（F3）和广元—大邑隐伏断裂带（F4）以及一系列龙门山山前断裂（图1）（王二七等，2001;徐锡伟等，2008;付碧宏等，2008）。
龙门山地区位于高原气候区与季风气候区交界处，受大气降水和高山冰雪融水的共同补给，水系较为发育（李奋生等，2015），从西南到东北依次发育青衣江、岷山、沱江、涪江和嘉陵江等一系列斜交于龙门山断裂带的河流（王志等，2010）。龙门山断裂带由多期构造变动形成，水文地质条件复杂，温泉多沿断裂带发育且呈带状分布。研究区主要的含水层包括第四系松散堆积层孔隙含水层、侏罗系—白垩系碎屑岩裂隙孔隙—基岩裂隙含水层、震旦系—三叠系碳酸盐岩岩溶裂隙—基岩裂隙水以及前震旦系基底含水层4大含水层，且在褶皱和断裂构造共同作用下，各含水层之间相互联系（陶广斌，2019）。前人已收集到的地震数据显示龙门山地区最早记录的地震为1161年北川MS4.45地震，根据中国地震局监测预报司统计，自公元1161年至2020年7月，本区共发生MS>43/4地震196次，其中MS5.0～5.9地震112次，MS6.0～6.9地震21次，MS>7.0地震5次，因此，该区具有重要的研究价值（张岳桥等，2008）。 2 研究方法
笔者对2008—2016年龙门山断裂带13个温泉（表1）进行了9次调查，共59个温泉水样，每个温泉水样采集4瓶用于测定其主量元素、微量元素、氢氧同位素及SiO2含量。野外采样时，采样容器为250 mL的聚乙烯塑料瓶，将聚乙烯瓶置于水中润洗2～3次后装样，取样时将瓶内空气排净，避免水样受空气影响，并现场测定温泉水的水温、pH、电导率，并向测定微量元素的样品中加入1.00 mL的浓硝酸（pH<2），在室温条件下保存。
温泉水样品中常量元素的测定在中国地震局地震预测重点实验室完成，使用仪器为Dionex ICS-900离子色谱仪及AS40自动进样器，检测限为0.01 mg/L。温泉水中的HCO-3和CO2-3浓度使用ZDJ-3D型电位滴定仪通过标准滴定程序进行测量，滴定所用盐酸浓度为0.05 mol/L。为了校准色谱，在测量每批水样品之前需测量标准样品，测量偏差在±2‰以内。离子平衡（ib）计算公式如下：
ib=∑阳离子-∑阴离子（∑阳离子+∑阴离子）×0.5×100%（1）
式中：阴阳离子的测量误差均小于5%。
微量元素的含量分析在中国地震局地壳应力研究所采用Element型电感耦合等离子体质谱仪ICP-MS完成，样品的分析精度优于5%。氢氧同位素分析采用气体同位素质谱仪MAT253，以国际标准即标准平均海洋水（SMOW）为标准，分别记为δD和δ18O，氢氧同位素参考标准物质为国家标准物质（GBW04401，04402），样品的分析精度为δD<±1‰和δ18O<±0.2‰。利用电感耦合等离子体发射光谱仪Optima-5300DV检测SiO2含量。
3 温泉水样分析结果
对所采集的温泉水样进行分析，结果见表2。研究区温泉水样水温变化范围为10.8℃～64.0 ℃;pH变化范围为6.86～9.50，平均值为7.96，均为典型的弱碱性水;电导率变化范围为154.00～132 400.00 μs/cm，TDS的变化范围为0.08～201.37 g/L，其中WJJ为卤水，LFS、ZGS为盐水，HSW、TC、PL为咸水，GYG、SZZ、JYG为淡水，LX、WC、ZBS、ZLJ为微咸水温泉。从表2还可看出，水样中的主量元素相对丰度：阳离子主要为Na+、Mg2+、Ca2+;阴离子主要为SO2-4、HCO-3、Cl-;水样中δD和δ18O测值的分布范围分别为-6.18‰～-17.80‰和-51.50‰～-119.20‰。微量元素含量较低，其丰度结果见表3，SiO2含量的分布范围为0.81～29.70 mg/L。
4 讨论
4.1 温泉水来源
Craig（1961）通过总结全球各地区大气降水中的δD、δ18O数值变化范围及二者之间的关系，最先提出了用于判断地下水现代补给来源的全球大气降水线（GMLW）。在此之后，不断有学者提出了一系列降水线方程，本文选取了Xu等（2015）对龙门山地区建立的大气降水线方程（LMLW）：δD=7.6δ18O+11.3，绘制龙门山断裂带温泉δD-δ18O关系图。从图2可观察到，龙门山断裂带的温泉水同位素数据主要分布在大气降水线两侧，表明研究区温泉水及观测井水主要来源于大气降水。
研究区水样δD和δ18O测值的分布范围分别为-6.18‰～-17.80‰和-51.50‰～-119.20‰，由于龍门山地区海拔变化较大，氢氧同位素会受到高程效应的影响，因此可以利用同位素的高程效应估算温泉水的补给高程（Clark，Royden，2000;于津生，1997;Zhou et al，2010;刘成龙等，2020）。
我国西部降水的δD值与高程H的关系为（于津生，1997）：
δD=-0.026H-30.2（2）
中国大气降水δD的高程效应可以表示为（Zhou et al，2010）：
δD=-0.03HALT -27（3）
式中：HALT为海拔高程值。
根据这些公式计算研究区温泉补给高程再取平均值（表4）。由表4可见，研究区温泉的补给高程大致为0.8～3.2 km;断裂带南西段温泉的补给高程为1.4～3.2 km，补给高程较高;北东段温泉的补给高程为0.8～1.9 km，相对较低，这与龙门山断裂带总体呈现由南向北、由西向东倾斜的趋势一致。
4.2 水化学成分特征
4.2.1 常量元素
本文使用了舒卡列夫分类法（毫克当量百分比分别大于25%的阴阳离子参与命名）将龙门山断裂带的温泉水划分为Na-HCO3、Ca·Na-HCO3、Na-HCO3·SO4、Ca-HCO3·SO4、Ca·Mg-SO4、Na-HCO3·Cl、Na-SO4·Cl、Na-Cl八种水化学类型，利用origin软件绘制Piper三线图，观测离子分配规律（图3），从空间上看，该地区水化学类型复杂，从时间上看，该地区温泉水化学类型在采样期间没有明显改变。
位于茂县—汶川断裂带的GRG、JYG温泉水化学类型为Na-HCO3、Na-HCO3·SO4，矿化度较低，其含水层为分别为侏罗系莲花口组钙质砂岩、粉砂岩以及震旦系灯影组白云岩、页岩。温泉中Na+主要来源于钠长石、钾长石、钙长石等矿物发生的溶解作用和阳离子吸附的交替作用。HCO-3主要来源于灰岩、白云岩等碳酸盐岩的溶解，一般来说，随着地下水入渗深度的增加，溶滤作用增强，温泉水中的SO2-4含量也会增加。WC、LX温泉的水化学类型为Ca·Mg-SO4，出露于志留系-泥盆系泥质灰岩夹黄铁矿岩层中，含水介质中黄铁矿的氧化作用可加快灰岩及白云岩的溶解，从而使WC、LX温泉富集Ca2+、Mg2+及SO2-4。 位于龙门山地区北川—映秀断裂、安县—灌县断裂及广元—大邑隐伏断裂带南段附近的WJJ、ZGS、HSW温泉，靠近四川盆地，矿化度高，为卤水和咸水，水化学类型为Na-Cl型及Ma-SO2-4型，而LFS、ZBS、ZLJ温泉水化学类型为Na-HCO3·Cl型。这可能是由于四川盆地从震旦系到白垩系各个地层均赋存有地下卤水（周训等，2015），而在构造应力相对集中带、断裂破碎带和裂隙发育带是卤水的主要富集地带（林耀庭，2001;杨耀等，2019），因此WJJ、ZGS、HSW温泉的形成可能起源于古海水，在水体循环过程中以浓缩作用为主，最终形成高矿化度的深部卤水（曹琴等，2015）。而LFS、ZBS、ZLJ温泉可能是地下水体循环过程中经过含盐地层，水体溶滤大量氯化物岩类（周训，1993），从而导致地下水的矿化度增高。研究区中段SZZ温泉水化学类型为Ca·Na-HCO3型，而北段PL、TC温泉的水化学类型为Ca-SO4·HCO3，都出露于三叠系嘉陵江组灰岩、白云岩。而北段温泉中SO2-4含量增加是由该温泉含水介质中三叠系嘉陵江组岩层中夹厚层硬石膏所致。
4.2.2 微量元素
深循环地下热水中微量元素的含量特征与区域地球化学背景值、元素性质、地下水的物理化学条件、水-岩相互作用等密切相关（Zhang et al，2016）。研究区温泉围岩以白云岩、泥质灰岩、碎屑岩及砂岩为主，温泉水循环于碳酸盐岩中，微量元素富集因子的大小可以定性判断温泉的水岩反应程度（Soto-Jiménez，Páez-Osuna，2001），其计算方法为：
EFi=（Ci/CR）w /（Ci/CR）r（4）
式中：CR为选定的参比元素含量;Ci为样品中元素含量;w为水样中元素浓度;r为岩石中元素浓度。因此将龙门山地区方解石及石英岩脉的微量元素含量（薛钧月，2009）作为标准，使用Ti浓度对微量元素数据进行归一化，计算研究区微量元素富集因子EF（图4）。从图4中可以明显看到微量元素的富集状况，结合微量元素含量，对研究区微量元素地球化学特征进行分析，发现断裂带南西段温泉微量元素相对富集，且微量元素富集因子比北东段大。
4.3 水岩反应程度
Na-K-Mg三角图能够反应温泉的水岩平衡状态和平衡温度（Giggenbach，1988）。从研究区的水岩反应图（图5）可以看出，WJJ温泉落在完全平衡区。林耀庭等（1997）和Zhou等（2015）根据四川盆地地下卤水的δD和δ18O值将其划分为4种类型：大气水淋滤型卤水、海相沉积型卤水、海相沉积-大气降水叠加型卤水以及海相沉积-岩浆水叠加型卤水。WJJ温泉矿化度较高，其δD和δ18O落于大气降水线上部，表明其成因与大气降水有关，为大气淋滤型卤水，且水-岩反应达到离子平衡状态。HSW06、ZLJ06、LFA06、ZBS06、ZGS06、ZGS07为部分平衡水，其它温泉为未成熟水在图中Mg1/2的一侧，表明水岩反应程度低，并受到一些浅层冷水混合作用。
4.4 热储温度与循环深度
目前常用的地热温标主要分为两大类，一类是阳离子温标，另一类为SiO2温标。阳离子温标方法一般用于估算部分平衡水及平衡水的热储温度（Fournier，Rowe，1966），阳离子温标法中K-Mg温标适用于中低温地热水。研究区WJJ平衡水温泉及HSW、ZLJ、LFA、ZBS、ZGS温泉等部分平衡水为中低温地热水系统，因此采用K-Mg温标计算热储温度更适合（王逸凌等，2020）。而SiO2相对于其他矿物而言，稳定性较高，能够指示未成熟水的热储温度（庞忠和等，1990），因此该地区未成熟水温泉采用SiO2温标。通过计算得到龙门山地区热储温度范围为20.1 ℃～123.4 ℃（表5）。
研究区温泉及观测井为大气降水入渗补给通过深循环获得深部热流加热之后在升流至地表形成，其热水循环深度计算公式如下：
H=（T-T0）/g+h（5）
式中：H为循环深度（km）;T为地热水的热储温度（℃）;T0为研究区的平均气温（℃），研究区内平均气温取17 ℃;g为地温梯度，为24.5 ℃/km;h为常温带深度，为30 m（崔希林，2014）。
由表5得到龙门山断裂带内热水循环深度范围为0.1～4.4 km（表5）。其北东段TC及PL温泉循环深度为1.1～1.5 km，中部GRG等温泉的循环深度为0.7～2.5 km，其中以Cl-为主的LFS、ZBS、ZLJ温泉，其循环深度为0.9～2.8 km，北部HSW、WJJ、ZGS温泉循环深度为4.0～4.4 km。结合温泉常量元素及微量元素特征可以发现，以HCO-3为主的温泉，循环深度较浅，SO2-4次之，而以Cl-为主要阴离子的温泉循环深度较深。结合Na-K-Mg三角图（图4）、水化学组分特征以及循环深度综合分析，龙门山断裂带温泉的水化学类型沿茂县—汶川断裂带、北川—映秀断裂、安县—灌县断裂及广元—大邑隐伏断裂带从北到南、从西到东依次为重碳酸型、硫酸型、氯化物型温泉水，其水岩反应程度、水化学组分含量、循环深度沿断裂分布不断加深。
4.5 温泉水化学变化与地震的响应
构造活动和地震活动会打破地壳中流体原有的水-岩平衡状态，使温泉循环深度，热储温度，深部热水循环程度，水-岩相互作用改变，造成离子组分，氢氧同位素等发生变化（Skelton et al，2014）。通过龙门山地区温泉的各离子浓度绘制主要离子相对变化率（（测量值-平均值）/平均值×100%）随时间变化的折线图（图6），根据采集样品离子含量的方差占平均值的百分比来判断研究区异常值的临界。
通过对比距离汶川MS8.0地震震中108 km以内的GRG、JYG、SZZ、LX、WC温泉在震后1个月及5个月采集的两次温泉水样（表6），发现随着时间的变化，温泉中的主要离子组分（Na+、Cl-、SO2-4）及微量元素含量（Li+、Pb2+、As2+）呈下降趋势，同时2008年6月至2010年6月采集水样的δD和δ18O也回归至大氣降水线附近。从δD-δ18O关系图（图2）发现，2008年6月和10月采集的GRG01、GRG02、JYG01、LX01、SZZ01、WC02的氧同位素偏离于大气降水线的右下方。这可能是由于区域构造活动的改变使断层封闭性改变，增加了温泉水与矿物的接触，使其水岩反应增强，水化学离子含量增加。 2013年芦山MS7.0地震后，在2013年4月及2014年8月采集了震中距为65 km以内HSW和WJJ温泉以及震中距为195 km以内的SZZ、WC的水样，其化学组分部分呈下降趋势，但变化不明显，也有少量元素呈增加的趋势（表6）。这可能与芦山地震后研究区多次发生5级以上余震，且在2014年11月发生了康定MS6.5地震（见图6中红色框），其深部热水供应变化不大有关。从δD-δ18O关系图可以发现，2013年5月采集的JYG06、LFS06、PL06、SZZ06、WC06、WJJ06、ZBS06、ZGS06、ZLJ06氢氧同位素偏离于大气降水线的左上方。因此，推测2013年芦山地震时各温泉离子浓度发生变化的原因为地震引起的各含水层间渗透性增强，地下水混合所致。
地震活动还会改变应力场变化，断裂带裂隙发生变化，引起渗透率发生变化，从而导致地震震中区周围的流体行为改变（Zhou et al，2017）。周永胜等（2008）、宋娟和周永胜（2013）研究发现川滇地区的龙门山断裂以及红河断裂带发生内部流体孔隙压力周期性变化以及断层带脆-塑性转化、裂缝张开与愈合等现象。强震发生后，在岩石塑性变形、压溶、静态重结晶、水岩反应等作用下，断层内部由地震形成的断层和微裂缝逐渐被愈合，把断层深部的流体密封。Luo等（2017）研究发现地下热水的氧同位素向右偏移的原因为地下热水在深循环过程中与富集δ18O围岩发生水岩作用，伴随氧同位素的交换，导致地下热水δ18O的升高，这与2008年氢氧同位素的变化表现一致。2013年5月，氢氧同位素偏离大气降水线的左上方，其原因为温泉水与深部来源的CO2相互溶解的过程中发生了δ18O同位素交换（Benavente et al，2016），这与2013年氢氧同位素的变化表现一致。
地震的孕育和发生过程中也会改变含水层及其周围岩土环境，引起渗透率变化或水体混合，从而引起水文地球化学变化。地震后，断层愈合，各离子浓度又大幅度回落，这可能是由地震后断裂内部渗透性减弱、深部流体供应减少、水岩反应程度衰减以及当地降雨量的增多、大量浅层水混入所致。
4.6 温泉成因模式
结合已有的地球物理资料和本文的研究成果建立模型，龙门山断裂带周围3.2 km的高山为温泉的补给区，大气降水入渗沿断裂带汇聚补给地下水，地下水沿断层破碎带或裂隙带经历深循环，获得来自深部热源加热，此时水体的温度增加至123 ℃，循环深度增大至地下4 km左右，与周围的岩石发生水岩反应，并在水头差驱动下上升，在地形较低的山谷和河谷处出露地表（图7）。而龙门山断裂带两侧中上地壳存在低速高导层，为龙门山中上地壳内浅源地震的孕育提供了深部构造条件。地震的孕育及发生可以改变地下应力应变状态，触发断裂带内部含水层破碎，导致深部流体以及浅层冷水混入从而影响地下热水的水岩反应程度，并改变研究区温泉的水化学组分及同位素特征。本文通过以上研究建立了龙门山地区温泉水化学模型，结合研究区温泉水化学变化对地震的响应，为该地区的震情跟踪、异常落实和地震预测提供基础研究。
5 结论
本文对龙门山断裂带13个温泉及观测井的水文地球化学特征进行分析研究，得到以下结论：
（1）龙门山断裂带温泉水的主要来源为周围3.2 km高山的大气降水。
（2）根据Na-K-Mg三角图、水化学组分特征以及循环深度综合分析，龙门山断裂带温泉的水化学类型沿龙门山断裂带从北到南、从西到东依次为重碳酸型、硫酸型、氯化物型，水岩反应程度、水化学组分含量、循环深度沿断裂分布不断增加。其中茂县—汶川断裂带温泉为未成熟水，水岩反应程度较弱，循环深度达到0.7～2.5 km;北川—映秀断裂、安县—灌县断裂及广元—大邑隐伏断裂带北部为未成熟水，水岩反应程度较弱，循环深度达到1.1～1.5 km;中部为部分平衡水，水岩反应程度较强，循环深度达到0.9～2.8 km;南部为平衡水及部分平衡水，水岩反应程度强，循环深度达到4.0～4.4 km。
（3）对比汶川MS8.0地震、芦山MS7.0地震震后几个月内温泉水化学组分（Na+、Cl-、SO42-、Li、Sr、As、δ18O及δD）的变化，发现震中距200 km以内的温泉水化学组分呈下降趋势，结合氢氧同位素的变化推测这可能是由震后断层愈合、断层渗透系数减弱、水岩反应程度衰减所致。以上研究可为该地区震情跟踪、异常落实和地震预测提供基础。
感谢中国科学院地质与地球物理研究所兰州油气资源研究中心与中国地震局地壳应力研究所的相关工作人员在样品分析過程中给予的帮助，在此表示衷心的感谢！
参考文献：
曹琴，周训，张欢，等.2015.四川盆地卧龙河储卤构造地下卤水的水化学特征及成因[J].地质通报，34（5）：990-997.
陈志.2014.汶川地震后川西地区温泉水地球化学研究[D].合肥：中国科学技术大学.
程万正，官致君，李军.2013.对汶川8.0级地震前四川地区地下流体观测异常的研究[J].四川地震，（2）：1-8.
崔希林.2014.四川盆地中西部温泉成因模式研究[D].成都：成都理工大学.
邓起东，陈社发，赵小麟.1994.龙门山及其邻区的构造和地震活动及动力学[J].地震地质，16（4）：389-403.
付碧宏，时丕龙，张之武.2008.四川汶川MS8.0大地震地表破裂带的遥感影像解析[J].地质学报，82（12）：1679-1687.
柯斌.2014.四川茂县吉鱼温泉成因模式及开发远景分析[D].成都：成都理工大学.
李奋生，赵国华，李勇，等.2015.龙门山地区水系发育特征及其对青藏高原东缘隆升的指示[J].地质论评，61（2）：345-355. 李生紅.2013.川东褶皱带中、下三叠统水文地质特征研究[D].成都：成都理工大学.
林耀庭，唐庆，熊淑君，等.1997.四川盆地卤水的氢、氧同位素地球化学特征及其成因分类研究[J].地质地球化学，（4）：20-26.
林耀庭.2001.四川盆地三叠纪卤水成藏条件[J].化工矿产地质，23（1）：19-24.
刘成龙，王广才，史浙明，等.2020.云南硫磺洞温泉水文地球化学特征和成因分析[J].地震研究，43（2）：278-286.
庞忠和，汪集旸，樊志成.1990.利用SiO2混合模型计算漳州地热田热储温度[J].科学通报，（1）：57-59.
蒲小武，武银，狄国荣，等.2013.甘肃清水温泉井与临夏井水温在汶川地震前后的异常变化特征[J].地震研究，36（3）：269-274.
宋娟，周永胜.2013.断层带流体对断层强度和强震孕育的影响[J].国际地震动态，420（12）：5-16.
陶广斌.2019.川西淡矿化温泉地球化学特征及成因研究[D].成都：成都理工大学.
王二七，孟庆任，陈智樑，等.2001.龙门山断裂带印支期左旋走滑运动及其大地构造成因[J].地学前缘，8（2）：165-174.
王逸凌，穆文清，史浙明.2020.鲜水河断裂南段温泉水文地球化学特征[J].地震研究，43（2）：287-295.
王云，赵慈平，刘峰，等.2014.小江断裂带及邻近地区温泉地球化学特征与地震活动关系研究[J].地震研究，37（2）：228-243.
王志，赵琳娜，张国平，等.2010.汶川地震灾区堰塞湖流域面雨量计算方法研究[J].气象，36（6）：7-12.
徐锡伟，闻学泽，叶建青，等.2008.汶川MS8.0地震地表破裂带及其发震构造[J].地震地质，30（3）：597-629.
薛钧月.2009.龙门山构造带中—北段构造流体地球化学特征及其与成藏关系的探讨[D].成都：成都理工大学.
晏锐，官致君，刘耀炜.2015.川西温泉水温观测及其在芦山MS7.0地震前的异常现象[J].地震学报，37（2）：347-356.
杨耀，周晓成，官致君，等.2019.川西地下流体观测井水文地球化学特征[J].矿物岩石地球化学通报，38（5）：966-976.
于津生.1997.中国同位素地球化学研究[M].北京：科学出版社.
张岳桥，杨农，施炜，等.2008.青藏高原东缘新构造及其对汶川地震的控制作用[J].地质学报，82（12）：1668-1678.
周晓成，王万丽，李立武，等.2020.金沙江—红河断裂带温泉气体地球化学特征[J].岩石学报，36（7）：2197-2214.
周训，曹琴，尹菲，等.2015.四川盆地东部高褶带三叠系地层卤水和温泉的地球化学特征及成因[J].地质学报，89（11）：1908-1920.
周训.1993.四川盆地龙女寺储卤构造深层地下卤水的水文地球化学特征及成因[J].现代地质，7（1）：83-92.
周永胜，何昌荣，杨晓松.2008.中地壳韧性剪切带中的水与变形机制[J].中国科学：地球科学，38（7）：819-832.
Benavente O，Tassi F，Reich M，et al.2016.Chemical and isotopic features of cold and thermal fluids discharged in the Southern Volcanic Zone between 32.5°S and 36°S：Insights into the physical and chemical processes controlling fluid geochemistry in geothermal systems of Central Chile[J].Chemical Geology，420（1）：97-113.
Chen Z，Zhou X C，Du J G，et al.2015.Hydrochemical characteristics of hot spring waters in the Kangding district related to the Lushan MS= 7.0 earthquake in Sichuan，China[J].Natural Hazards & Earth System Ences Discussions，2（6）：1149-1156.
Clark M K，Royden L H.2000.Topographic ooze：Building the eastern margin of Tibet by lower crustal flow[J].Geology，28（8）：703-706.
Craig H.1961.Isotopic Variations in Meteoric Waters[J].Science，133（3465）：1702-1703.
Fournier R O，Rowe J J.1966.The deposition of silica in hot springs[J].29（1）：585-587.
Giggenbach W F.1988.Geothermal solute equilibria.Derivation of Na-K-Mg-Ca geoindicators[J].Geochimica Et Cosmochimica Acta，52（12）：2749-2765. Godard V，Pik R，Cattin R，et al.2009.Late Cenozoic evolution of the central Longmen Shan，eastern Tibet：Insight from（U-Th）/He thermochronometry[J].Tectonics，28（5）：TC5009.
Kouketsu Y，Shimizu I，Wang Y，et al.2017.Raman spectra of carbonaceous materials in a fault zone in the Longmenshan thrust belt，China;comparisons with those of sedimentary and metamorphic rocks[J].Tectonophysics，699（1）：129-145.
Lambrakis N，Kallergis G.2005.Contribution to the study of Greek thermal springs：hydrogeological and hydrochemical characteristics and origin of thermal waters[J].Hydrogeology Journal，13（3）：506-521.
Lin X，Dreger D，Ge H，et al.2018.Spatial and Temporal Variations in the Moment Tensor Solutions of the 2008 Wenchuan Earthquake Aftershocks and Their Tectonic Implications[J].Tectonics，37（3）：989-1005.
Liu Y，Ren H，Wang B.2009.Application of environmental isotopes and tracer techniques to seismic subsurface fluids[J].Earth Science Frontiers，16（1）：369-377.
Luo J，Pang Z H，Kong Y K，et al.2017.Geothermal potential evaluation and development prioritization based on geochemistry of geothermal waters from Kangding area，western Sichuan，China[J].Environmental Earth Ences，76（9）：343.
Mokadem N，Demdoum A，Hamed Y，et al.2016.Hydrogeochemical and stable isotope data of groundwater of a multi-aquifer system：Northern Gafsa basin- Central Tunisia[J].Journal of African Earthences，114（1）：174-191.
Pérez N M，Hernández P A，Igarashi G，et al.2008.Searching and detecting earthquake geochemical precursors in CO2-rich groundwaters from Galicia，Spain[J].Geochemical Journal，42（1）：75-83.
Okuyama Y，Funatsu T，Fujii T，et al.2016.Mid-crustal fluid related to the Matsushiro earthquake swarm（1965-1967）in northern Central Japan：geochemical reproduction，Tectonophysics，679（3）：61-72.
Shi Z，Zhang H，Wang G.2020.Groundwater trace elements change induced by MS5.0 earthquake in Yunnan[J].Journal of Hydrology，581：124424.
Skelton A，Andren M，Kristmannsdottir H，et al.2014.Changes in groundwater chemistry before two consecutive earthquakes in Iceland[J].Nature Geoscience，7（10）：752-756.
Soto-Jiménez M F，Páez-Osuna F.2001.Distribution and Normalization of Heavy Metal Concentrations in Mangrove and Lagoonal Sediments from Mazatlán Harbor（SE Gulf of California）[J].Estuarine Coastal & Shelf Ence，53（3）：259-274.
Tsunogai U，Wakita H.1995.Precursory Chemical Changes in Ground Water：Kobe Earthquake，Japan[J].Science，5220（269）：61-63.
Walraevens K，Bakundukize C，Mtoni Y E，et al.2018.Understanding the hydrogeochemical evolution of groundwater in Precambrian basement aquifers：A case study of Bugesera region in Burundi[J].Journal of Geochemical Exploration，188（3）：24-42. Xu Q，Hoke G D，Jing L Z，et al.2015.Stable isotopes of surface water across the Longmenshan margin of the eastern Tibetan Plateau[J].Geochemistry Geophysics Geosystems，15（8）：3416-3429.
Xu Z，Zhen F，Li SX，et al.2019.The 2018 MS5.9 Mojiang Earthquake：Source model and intensity based on near-field seismic recordings[J].Earth and Planetary Physics，3（3）：88-101.
Zhang Y F，Tan H B，Zhang W J，et al.2015.A New Geochemical Perspective on Hydrochemical Evolution of The Tibetan Geothermal System[J].Geochemistry International，53（12）：1090-1106.
Zhang Y F，Tan H B，Zhang W J，et al.2016.Geochemical constraint on origin and evolution of solutes in geothermal springs in western Yunnan，China[J].Chemie der Erde-Geochemistry，76（1）：63-75.
Zhou X，Jiang C，Zhao J，et al.2015.Occurrence and resource evaluation of the subsurface high-K brines in the Pingluoba brine-bearing structure in western Sichuan Basin[J].Environmental Earth Sciences，73（12）：8565-8574.
Zhou X C，Li C，Ju X M，et al.2010.Origin of Subsurface Brines in the Sichuan Basin[J].Groundwater，35（1）：53-58.
Zhou X C，Liu L，Chen Z，et al.2017.Gas geochemistry of the hot spring in the Litang fault zone，Southeast Tibetan Plateau[J].Applied Geochemistry，79（1）：17-26.
Hydrogeochemical Characteristics of the Hot Springsin the Longmenshan Fault Zone
YAN Yucong1，LIU Fengli1，GUO Lishuang2，ZHOU Xiaocheng1，OUYANG Shupei1，LI Jingchao1，WANG Wanli1
（1.Institute of Earthquake Forecasting，China Earthquake Administration，Beijing 100036，China）
（2.National Institute of Natural Hazards，Ministry of Emergency Management，Key Laboratory ofCrustal Dynamics，China Earthquake Administration，Beijing 100085，China）
Abstract
Hydrogeochemical Monitoring in the regions with strong earthquake activities would facilitate the obtaining of the geochemical information about deep fluids correlated to seismic anomalies.In this paper the hydrogeochemical characteristics of 13 hot springs（wells）in the Longmenshan Fault Zone were investigated，so as to establish the hydrogeochemical cycle model of the hot springs and reveal the relationship between hydrogeochemical variations and seismic activity in the Fault Zone.The results show that：①The hot springs in the Longmenshan Fault Zone were recharged by the meteoric water at elevations ranging from 0.8 to 3.2 km.②The chemical types of the hot springs from north to south，and from west to east along the Longmenshan Fault Zone were respectively carbonic acid，sulfuric acid and chloride.Likewise，the circulation depth of spring water deepened;the water-rock reaction intensified;and enrichment factors of the trace element of the hot springs continuously increased.③After the Wenchuan MS8.0 earthquake and the Lushan MS7.0 earthquake，the dissolved hydrochemical composition got declining in hot springs within the epicenter of 200 kilometers，which could be the result of the permeability decrease and the attenuation of water-rock reaction inside the fault in the healing process of the Longmenshan Fault after earthquake，
Keywords：hot spring;hydrogeochemistry;earthquake;the Longmenshan Fault Zone
收稿日期：2021-01-11.
基金项目：国家重点研发计划（2017YFC1500501）、国家自然科学基金面上项目（41673106）和中国地震局地震预测研究所基本科研业务费（2017IES010205，2016IES010304，2018IEF010104，2020IEF0604，2020IEF0703）联合资助.
第一作者简介：颜玉聪（1996-），女，在读硕士研究生，主要从事地震地球化学研究.E-mail：[email protected].
通訊作者简介：郭丽爽（1983-），高级工程师，主要从事地下流体地球化学方面的研究.E-mail：[email protected].
关键词：温泉;水文地球化学;地震;龙门山断裂带
中图分类号：P314.1 文献标识码：A 文章编号：1000-0666（2021）02-0170-15
0 引言
温泉多沿断裂带出露，是沟通地壳深部与浅部的窗口，和深部地下水系统有着密切的关系（Zhang et al，2015;Zhou et al，2017）。研究活动断层的温泉流体地球化学特征有助于深入了解该断层的活动方式、活动强度以及切割深度等断层运动特性，对探索断层带深部、特别是震源区附近介质的物理化学环境、研究各类地球化学组分的地震前兆机制、寻找特定地区的前兆观测组分有重要意义（Lambrakis，Kallergis，2005;Liu et al，2009）。因此，在地震易发区对温泉的水文地球化学特征进行监测，可捕捉来自地壳深部活动的信息（蒲小武等，2013;王云等，2014;Walraevens et al，2018;Xu et al，2019）。
近年来，国内外许多专家对温泉动态与地震活动的关系进行了研究（Tsunogai，Wakita，1995;Shi et al，2020）。Skelton等（2014）对冰岛2012年10月以及2013年4月2次5级地震之间以及地震前后4～6个月的地下水中稳定同位素比率和钠、钙等溶质的浓度变化规律进行研究，获得了该地区地震前后地下水水文地球化学变化的规律;Okuyama等（2016）研究发现1965—1968年日本松代地震群前后水中某些溶解组分随井深呈线性增长;Pérez 等（2008）研究1995年西班牙加利西亚4.6级地震以及1997年5.3级地震发现，在震中距为90 km的地下水观测井中Cl-浓度明显增加。另外，Shi等（2020）的研究显示，氢氧同位素以及微量元素的变化对地震的响应也较为明显，2018年通海MS5.0地震后，江川观测井的水化学以及氢氧同位素明显上升，微量元素（Li+、Pb2+、As2+等元素）的浓度呈明显下降。
龙门山断裂带是青藏高原东缘山脉与扬子板块的结合部位，地质构造复杂，地震活动剧烈（Lin et al，2018）。自2008年汶川MS8.0地震后，我国在该地区建立了大量的流动观测网、温泉和热水井固定观测点，丰富了该地区的流体地球化学背景（程万正等，2013）。目前，关于龙门山断裂带地下流体的研究主要集中在地震前后温泉水温、水化学组分、水氡以及温泉逸出气异常等方面（Chen et al，2015;晏锐等，2015;Zhou et al，2015;杨耀等，2019）。Chen等（2015）对2008—2010年川西地区32个温泉的水化学变化进行了研究，发现了汶川地震后温泉的离子浓度呈下降变化;Zhou等（2015）对该地区32个温泉气进行了3次采样，结果表明汶川MS8.0地震后，四川西部地区有更多地幔流體进入地壳，地幔成分的流体上涌可能与大地震的发生有关系;杨耀等（2019）研究了川西地区地下流体观测井的地下水的水化学特征，确定了四川西部地区主要观测井的水文地球化学背景。
前人对龙门山断裂带温泉水文地球化学特征的研究主要集中于川西地区附近，没有对整个龙门山断裂带展开更为细致的分析，鉴于此，本文对龙门山断裂带13个温泉及观测井的水文地球化学特征进行分析，针对该断裂带不同断裂计算温泉水的热储温度与循环深度，确定温泉水的补给来源和补给高程，判断其水文地球化学的控制因素和影响因素，建立龙门山断裂带温泉水循环模型。通过探讨2008—2016年龙门山断裂带温泉水文地球化学变化与汶川、芦山地震活动的关系，为该地区震情跟踪、异常落实和地震预测提供基础研究。
1 研究区地质概况
龙门山断裂带是青藏高原东缘的边界山脉，北起广元、南至天全，长约500 km、宽约30～50 km，呈NE—SW向展布，形成于印支晚期，受多期复杂构造运动影响（邓起东等，1994;Clark，Royden，2000;Godard et al，2009;Kouketsu et al，2017）。该构造带是由一系列大致平行的逆冲、逆冲兼走滑的断裂构成，自西向东发育有茂汶—汶川断裂（F1）、北川—映秀断裂（F2）、安县—灌县断裂（F3）和广元—大邑隐伏断裂带（F4）以及一系列龙门山山前断裂（图1）（王二七等，2001;徐锡伟等，2008;付碧宏等，2008）。
龙门山地区位于高原气候区与季风气候区交界处，受大气降水和高山冰雪融水的共同补给，水系较为发育（李奋生等，2015），从西南到东北依次发育青衣江、岷山、沱江、涪江和嘉陵江等一系列斜交于龙门山断裂带的河流（王志等，2010）。龙门山断裂带由多期构造变动形成，水文地质条件复杂，温泉多沿断裂带发育且呈带状分布。研究区主要的含水层包括第四系松散堆积层孔隙含水层、侏罗系—白垩系碎屑岩裂隙孔隙—基岩裂隙含水层、震旦系—三叠系碳酸盐岩岩溶裂隙—基岩裂隙水以及前震旦系基底含水层4大含水层，且在褶皱和断裂构造共同作用下，各含水层之间相互联系（陶广斌，2019）。前人已收集到的地震数据显示龙门山地区最早记录的地震为1161年北川MS4.45地震，根据中国地震局监测预报司统计，自公元1161年至2020年7月，本区共发生MS>43/4地震196次，其中MS5.0～5.9地震112次，MS6.0～6.9地震21次，MS>7.0地震5次，因此，该区具有重要的研究价值（张岳桥等，2008）。 2 研究方法
笔者对2008—2016年龙门山断裂带13个温泉（表1）进行了9次调查，共59个温泉水样，每个温泉水样采集4瓶用于测定其主量元素、微量元素、氢氧同位素及SiO2含量。野外采样时，采样容器为250 mL的聚乙烯塑料瓶，将聚乙烯瓶置于水中润洗2～3次后装样，取样时将瓶内空气排净，避免水样受空气影响，并现场测定温泉水的水温、pH、电导率，并向测定微量元素的样品中加入1.00 mL的浓硝酸（pH<2），在室温条件下保存。
温泉水样品中常量元素的测定在中国地震局地震预测重点实验室完成，使用仪器为Dionex ICS-900离子色谱仪及AS40自动进样器，检测限为0.01 mg/L。温泉水中的HCO-3和CO2-3浓度使用ZDJ-3D型电位滴定仪通过标准滴定程序进行测量，滴定所用盐酸浓度为0.05 mol/L。为了校准色谱，在测量每批水样品之前需测量标准样品，测量偏差在±2‰以内。离子平衡（ib）计算公式如下：
ib=∑阳离子-∑阴离子（∑阳离子+∑阴离子）×0.5×100%（1）
式中：阴阳离子的测量误差均小于5%。
微量元素的含量分析在中国地震局地壳应力研究所采用Element型电感耦合等离子体质谱仪ICP-MS完成，样品的分析精度优于5%。氢氧同位素分析采用气体同位素质谱仪MAT253，以国际标准即标准平均海洋水（SMOW）为标准，分别记为δD和δ18O，氢氧同位素参考标准物质为国家标准物质（GBW04401，04402），样品的分析精度为δD<±1‰和δ18O<±0.2‰。利用电感耦合等离子体发射光谱仪Optima-5300DV检测SiO2含量。
3 温泉水样分析结果
对所采集的温泉水样进行分析，结果见表2。研究区温泉水样水温变化范围为10.8℃～64.0 ℃;pH变化范围为6.86～9.50，平均值为7.96，均为典型的弱碱性水;电导率变化范围为154.00～132 400.00 μs/cm，TDS的变化范围为0.08～201.37 g/L，其中WJJ为卤水，LFS、ZGS为盐水，HSW、TC、PL为咸水，GYG、SZZ、JYG为淡水，LX、WC、ZBS、ZLJ为微咸水温泉。从表2还可看出，水样中的主量元素相对丰度：阳离子主要为Na+、Mg2+、Ca2+;阴离子主要为SO2-4、HCO-3、Cl-;水样中δD和δ18O测值的分布范围分别为-6.18‰～-17.80‰和-51.50‰～-119.20‰。微量元素含量较低，其丰度结果见表3，SiO2含量的分布范围为0.81～29.70 mg/L。
4 讨论
4.1 温泉水来源
Craig（1961）通过总结全球各地区大气降水中的δD、δ18O数值变化范围及二者之间的关系，最先提出了用于判断地下水现代补给来源的全球大气降水线（GMLW）。在此之后，不断有学者提出了一系列降水线方程，本文选取了Xu等（2015）对龙门山地区建立的大气降水线方程（LMLW）：δD=7.6δ18O+11.3，绘制龙门山断裂带温泉δD-δ18O关系图。从图2可观察到，龙门山断裂带的温泉水同位素数据主要分布在大气降水线两侧，表明研究区温泉水及观测井水主要来源于大气降水。
研究区水样δD和δ18O测值的分布范围分别为-6.18‰～-17.80‰和-51.50‰～-119.20‰，由于龍门山地区海拔变化较大，氢氧同位素会受到高程效应的影响，因此可以利用同位素的高程效应估算温泉水的补给高程（Clark，Royden，2000;于津生，1997;Zhou et al，2010;刘成龙等，2020）。
我国西部降水的δD值与高程H的关系为（于津生，1997）：
δD=-0.026H-30.2（2）
中国大气降水δD的高程效应可以表示为（Zhou et al，2010）：
δD=-0.03HALT -27（3）
式中：HALT为海拔高程值。
根据这些公式计算研究区温泉补给高程再取平均值（表4）。由表4可见，研究区温泉的补给高程大致为0.8～3.2 km;断裂带南西段温泉的补给高程为1.4～3.2 km，补给高程较高;北东段温泉的补给高程为0.8～1.9 km，相对较低，这与龙门山断裂带总体呈现由南向北、由西向东倾斜的趋势一致。
4.2 水化学成分特征
4.2.1 常量元素
本文使用了舒卡列夫分类法（毫克当量百分比分别大于25%的阴阳离子参与命名）将龙门山断裂带的温泉水划分为Na-HCO3、Ca·Na-HCO3、Na-HCO3·SO4、Ca-HCO3·SO4、Ca·Mg-SO4、Na-HCO3·Cl、Na-SO4·Cl、Na-Cl八种水化学类型，利用origin软件绘制Piper三线图，观测离子分配规律（图3），从空间上看，该地区水化学类型复杂，从时间上看，该地区温泉水化学类型在采样期间没有明显改变。
位于茂县—汶川断裂带的GRG、JYG温泉水化学类型为Na-HCO3、Na-HCO3·SO4，矿化度较低，其含水层为分别为侏罗系莲花口组钙质砂岩、粉砂岩以及震旦系灯影组白云岩、页岩。温泉中Na+主要来源于钠长石、钾长石、钙长石等矿物发生的溶解作用和阳离子吸附的交替作用。HCO-3主要来源于灰岩、白云岩等碳酸盐岩的溶解，一般来说，随着地下水入渗深度的增加，溶滤作用增强，温泉水中的SO2-4含量也会增加。WC、LX温泉的水化学类型为Ca·Mg-SO4，出露于志留系-泥盆系泥质灰岩夹黄铁矿岩层中，含水介质中黄铁矿的氧化作用可加快灰岩及白云岩的溶解，从而使WC、LX温泉富集Ca2+、Mg2+及SO2-4。 位于龙门山地区北川—映秀断裂、安县—灌县断裂及广元—大邑隐伏断裂带南段附近的WJJ、ZGS、HSW温泉，靠近四川盆地，矿化度高，为卤水和咸水，水化学类型为Na-Cl型及Ma-SO2-4型，而LFS、ZBS、ZLJ温泉水化学类型为Na-HCO3·Cl型。这可能是由于四川盆地从震旦系到白垩系各个地层均赋存有地下卤水（周训等，2015），而在构造应力相对集中带、断裂破碎带和裂隙发育带是卤水的主要富集地带（林耀庭，2001;杨耀等，2019），因此WJJ、ZGS、HSW温泉的形成可能起源于古海水，在水体循环过程中以浓缩作用为主，最终形成高矿化度的深部卤水（曹琴等，2015）。而LFS、ZBS、ZLJ温泉可能是地下水体循环过程中经过含盐地层，水体溶滤大量氯化物岩类（周训，1993），从而导致地下水的矿化度增高。研究区中段SZZ温泉水化学类型为Ca·Na-HCO3型，而北段PL、TC温泉的水化学类型为Ca-SO4·HCO3，都出露于三叠系嘉陵江组灰岩、白云岩。而北段温泉中SO2-4含量增加是由该温泉含水介质中三叠系嘉陵江组岩层中夹厚层硬石膏所致。
4.2.2 微量元素
深循环地下热水中微量元素的含量特征与区域地球化学背景值、元素性质、地下水的物理化学条件、水-岩相互作用等密切相关（Zhang et al，2016）。研究区温泉围岩以白云岩、泥质灰岩、碎屑岩及砂岩为主，温泉水循环于碳酸盐岩中，微量元素富集因子的大小可以定性判断温泉的水岩反应程度（Soto-Jiménez，Páez-Osuna，2001），其计算方法为：
EFi=（Ci/CR）w /（Ci/CR）r（4）
式中：CR为选定的参比元素含量;Ci为样品中元素含量;w为水样中元素浓度;r为岩石中元素浓度。因此将龙门山地区方解石及石英岩脉的微量元素含量（薛钧月，2009）作为标准，使用Ti浓度对微量元素数据进行归一化，计算研究区微量元素富集因子EF（图4）。从图4中可以明显看到微量元素的富集状况，结合微量元素含量，对研究区微量元素地球化学特征进行分析，发现断裂带南西段温泉微量元素相对富集，且微量元素富集因子比北东段大。
4.3 水岩反应程度
Na-K-Mg三角图能够反应温泉的水岩平衡状态和平衡温度（Giggenbach，1988）。从研究区的水岩反应图（图5）可以看出，WJJ温泉落在完全平衡区。林耀庭等（1997）和Zhou等（2015）根据四川盆地地下卤水的δD和δ18O值将其划分为4种类型：大气水淋滤型卤水、海相沉积型卤水、海相沉积-大气降水叠加型卤水以及海相沉积-岩浆水叠加型卤水。WJJ温泉矿化度较高，其δD和δ18O落于大气降水线上部，表明其成因与大气降水有关，为大气淋滤型卤水，且水-岩反应达到离子平衡状态。HSW06、ZLJ06、LFA06、ZBS06、ZGS06、ZGS07为部分平衡水，其它温泉为未成熟水在图中Mg1/2的一侧，表明水岩反应程度低，并受到一些浅层冷水混合作用。
4.4 热储温度与循环深度
目前常用的地热温标主要分为两大类，一类是阳离子温标，另一类为SiO2温标。阳离子温标方法一般用于估算部分平衡水及平衡水的热储温度（Fournier，Rowe，1966），阳离子温标法中K-Mg温标适用于中低温地热水。研究区WJJ平衡水温泉及HSW、ZLJ、LFA、ZBS、ZGS温泉等部分平衡水为中低温地热水系统，因此采用K-Mg温标计算热储温度更适合（王逸凌等，2020）。而SiO2相对于其他矿物而言，稳定性较高，能够指示未成熟水的热储温度（庞忠和等，1990），因此该地区未成熟水温泉采用SiO2温标。通过计算得到龙门山地区热储温度范围为20.1 ℃～123.4 ℃（表5）。
研究区温泉及观测井为大气降水入渗补给通过深循环获得深部热流加热之后在升流至地表形成，其热水循环深度计算公式如下：
H=（T-T0）/g+h（5）
式中：H为循环深度（km）;T为地热水的热储温度（℃）;T0为研究区的平均气温（℃），研究区内平均气温取17 ℃;g为地温梯度，为24.5 ℃/km;h为常温带深度，为30 m（崔希林，2014）。
由表5得到龙门山断裂带内热水循环深度范围为0.1～4.4 km（表5）。其北东段TC及PL温泉循环深度为1.1～1.5 km，中部GRG等温泉的循环深度为0.7～2.5 km，其中以Cl-为主的LFS、ZBS、ZLJ温泉，其循环深度为0.9～2.8 km，北部HSW、WJJ、ZGS温泉循环深度为4.0～4.4 km。结合温泉常量元素及微量元素特征可以发现，以HCO-3为主的温泉，循环深度较浅，SO2-4次之，而以Cl-为主要阴离子的温泉循环深度较深。结合Na-K-Mg三角图（图4）、水化学组分特征以及循环深度综合分析，龙门山断裂带温泉的水化学类型沿茂县—汶川断裂带、北川—映秀断裂、安县—灌县断裂及广元—大邑隐伏断裂带从北到南、从西到东依次为重碳酸型、硫酸型、氯化物型温泉水，其水岩反应程度、水化学组分含量、循环深度沿断裂分布不断加深。
4.5 温泉水化学变化与地震的响应
构造活动和地震活动会打破地壳中流体原有的水-岩平衡状态，使温泉循环深度，热储温度，深部热水循环程度，水-岩相互作用改变，造成离子组分，氢氧同位素等发生变化（Skelton et al，2014）。通过龙门山地区温泉的各离子浓度绘制主要离子相对变化率（（测量值-平均值）/平均值×100%）随时间变化的折线图（图6），根据采集样品离子含量的方差占平均值的百分比来判断研究区异常值的临界。
通过对比距离汶川MS8.0地震震中108 km以内的GRG、JYG、SZZ、LX、WC温泉在震后1个月及5个月采集的两次温泉水样（表6），发现随着时间的变化，温泉中的主要离子组分（Na+、Cl-、SO2-4）及微量元素含量（Li+、Pb2+、As2+）呈下降趋势，同时2008年6月至2010年6月采集水样的δD和δ18O也回归至大氣降水线附近。从δD-δ18O关系图（图2）发现，2008年6月和10月采集的GRG01、GRG02、JYG01、LX01、SZZ01、WC02的氧同位素偏离于大气降水线的右下方。这可能是由于区域构造活动的改变使断层封闭性改变，增加了温泉水与矿物的接触，使其水岩反应增强，水化学离子含量增加。 2013年芦山MS7.0地震后，在2013年4月及2014年8月采集了震中距为65 km以内HSW和WJJ温泉以及震中距为195 km以内的SZZ、WC的水样，其化学组分部分呈下降趋势，但变化不明显，也有少量元素呈增加的趋势（表6）。这可能与芦山地震后研究区多次发生5级以上余震，且在2014年11月发生了康定MS6.5地震（见图6中红色框），其深部热水供应变化不大有关。从δD-δ18O关系图可以发现，2013年5月采集的JYG06、LFS06、PL06、SZZ06、WC06、WJJ06、ZBS06、ZGS06、ZLJ06氢氧同位素偏离于大气降水线的左上方。因此，推测2013年芦山地震时各温泉离子浓度发生变化的原因为地震引起的各含水层间渗透性增强，地下水混合所致。
地震活动还会改变应力场变化，断裂带裂隙发生变化，引起渗透率发生变化，从而导致地震震中区周围的流体行为改变（Zhou et al，2017）。周永胜等（2008）、宋娟和周永胜（2013）研究发现川滇地区的龙门山断裂以及红河断裂带发生内部流体孔隙压力周期性变化以及断层带脆-塑性转化、裂缝张开与愈合等现象。强震发生后，在岩石塑性变形、压溶、静态重结晶、水岩反应等作用下，断层内部由地震形成的断层和微裂缝逐渐被愈合，把断层深部的流体密封。Luo等（2017）研究发现地下热水的氧同位素向右偏移的原因为地下热水在深循环过程中与富集δ18O围岩发生水岩作用，伴随氧同位素的交换，导致地下热水δ18O的升高，这与2008年氢氧同位素的变化表现一致。2013年5月，氢氧同位素偏离大气降水线的左上方，其原因为温泉水与深部来源的CO2相互溶解的过程中发生了δ18O同位素交换（Benavente et al，2016），这与2013年氢氧同位素的变化表现一致。
地震的孕育和发生过程中也会改变含水层及其周围岩土环境，引起渗透率变化或水体混合，从而引起水文地球化学变化。地震后，断层愈合，各离子浓度又大幅度回落，这可能是由地震后断裂内部渗透性减弱、深部流体供应减少、水岩反应程度衰减以及当地降雨量的增多、大量浅层水混入所致。
4.6 温泉成因模式
结合已有的地球物理资料和本文的研究成果建立模型，龙门山断裂带周围3.2 km的高山为温泉的补给区，大气降水入渗沿断裂带汇聚补给地下水，地下水沿断层破碎带或裂隙带经历深循环，获得来自深部热源加热，此时水体的温度增加至123 ℃，循环深度增大至地下4 km左右，与周围的岩石发生水岩反应，并在水头差驱动下上升，在地形较低的山谷和河谷处出露地表（图7）。而龙门山断裂带两侧中上地壳存在低速高导层，为龙门山中上地壳内浅源地震的孕育提供了深部构造条件。地震的孕育及发生可以改变地下应力应变状态，触发断裂带内部含水层破碎，导致深部流体以及浅层冷水混入从而影响地下热水的水岩反应程度，并改变研究区温泉的水化学组分及同位素特征。本文通过以上研究建立了龙门山地区温泉水化学模型，结合研究区温泉水化学变化对地震的响应，为该地区的震情跟踪、异常落实和地震预测提供基础研究。
5 结论
本文对龙门山断裂带13个温泉及观测井的水文地球化学特征进行分析研究，得到以下结论：
（1）龙门山断裂带温泉水的主要来源为周围3.2 km高山的大气降水。
（2）根据Na-K-Mg三角图、水化学组分特征以及循环深度综合分析，龙门山断裂带温泉的水化学类型沿龙门山断裂带从北到南、从西到东依次为重碳酸型、硫酸型、氯化物型，水岩反应程度、水化学组分含量、循环深度沿断裂分布不断增加。其中茂县—汶川断裂带温泉为未成熟水，水岩反应程度较弱，循环深度达到0.7～2.5 km;北川—映秀断裂、安县—灌县断裂及广元—大邑隐伏断裂带北部为未成熟水，水岩反应程度较弱，循环深度达到1.1～1.5 km;中部为部分平衡水，水岩反应程度较强，循环深度达到0.9～2.8 km;南部为平衡水及部分平衡水，水岩反应程度强，循环深度达到4.0～4.4 km。
（3）对比汶川MS8.0地震、芦山MS7.0地震震后几个月内温泉水化学组分（Na+、Cl-、SO42-、Li、Sr、As、δ18O及δD）的变化，发现震中距200 km以内的温泉水化学组分呈下降趋势，结合氢氧同位素的变化推测这可能是由震后断层愈合、断层渗透系数减弱、水岩反应程度衰减所致。以上研究可为该地区震情跟踪、异常落实和地震预测提供基础。
感谢中国科学院地质与地球物理研究所兰州油气资源研究中心与中国地震局地壳应力研究所的相关工作人员在样品分析過程中给予的帮助，在此表示衷心的感谢！
参考文献：
曹琴，周训，张欢，等.2015.四川盆地卧龙河储卤构造地下卤水的水化学特征及成因[J].地质通报，34（5）：990-997.
陈志.2014.汶川地震后川西地区温泉水地球化学研究[D].合肥：中国科学技术大学.
程万正，官致君，李军.2013.对汶川8.0级地震前四川地区地下流体观测异常的研究[J].四川地震，（2）：1-8.
崔希林.2014.四川盆地中西部温泉成因模式研究[D].成都：成都理工大学.
邓起东，陈社发，赵小麟.1994.龙门山及其邻区的构造和地震活动及动力学[J].地震地质，16（4）：389-403.
付碧宏，时丕龙，张之武.2008.四川汶川MS8.0大地震地表破裂带的遥感影像解析[J].地质学报，82（12）：1679-1687.
柯斌.2014.四川茂县吉鱼温泉成因模式及开发远景分析[D].成都：成都理工大学.
李奋生，赵国华，李勇，等.2015.龙门山地区水系发育特征及其对青藏高原东缘隆升的指示[J].地质论评，61（2）：345-355. 李生紅.2013.川东褶皱带中、下三叠统水文地质特征研究[D].成都：成都理工大学.
林耀庭，唐庆，熊淑君，等.1997.四川盆地卤水的氢、氧同位素地球化学特征及其成因分类研究[J].地质地球化学，（4）：20-26.
林耀庭.2001.四川盆地三叠纪卤水成藏条件[J].化工矿产地质，23（1）：19-24.
刘成龙，王广才，史浙明，等.2020.云南硫磺洞温泉水文地球化学特征和成因分析[J].地震研究，43（2）：278-286.
庞忠和，汪集旸，樊志成.1990.利用SiO2混合模型计算漳州地热田热储温度[J].科学通报，（1）：57-59.
蒲小武，武银，狄国荣，等.2013.甘肃清水温泉井与临夏井水温在汶川地震前后的异常变化特征[J].地震研究，36（3）：269-274.
宋娟，周永胜.2013.断层带流体对断层强度和强震孕育的影响[J].国际地震动态，420（12）：5-16.
陶广斌.2019.川西淡矿化温泉地球化学特征及成因研究[D].成都：成都理工大学.
王二七，孟庆任，陈智樑，等.2001.龙门山断裂带印支期左旋走滑运动及其大地构造成因[J].地学前缘，8（2）：165-174.
王逸凌，穆文清，史浙明.2020.鲜水河断裂南段温泉水文地球化学特征[J].地震研究，43（2）：287-295.
王云，赵慈平，刘峰，等.2014.小江断裂带及邻近地区温泉地球化学特征与地震活动关系研究[J].地震研究，37（2）：228-243.
王志，赵琳娜，张国平，等.2010.汶川地震灾区堰塞湖流域面雨量计算方法研究[J].气象，36（6）：7-12.
徐锡伟，闻学泽，叶建青，等.2008.汶川MS8.0地震地表破裂带及其发震构造[J].地震地质，30（3）：597-629.
薛钧月.2009.龙门山构造带中—北段构造流体地球化学特征及其与成藏关系的探讨[D].成都：成都理工大学.
晏锐，官致君，刘耀炜.2015.川西温泉水温观测及其在芦山MS7.0地震前的异常现象[J].地震学报，37（2）：347-356.
杨耀，周晓成，官致君，等.2019.川西地下流体观测井水文地球化学特征[J].矿物岩石地球化学通报，38（5）：966-976.
于津生.1997.中国同位素地球化学研究[M].北京：科学出版社.
张岳桥，杨农，施炜，等.2008.青藏高原东缘新构造及其对汶川地震的控制作用[J].地质学报，82（12）：1668-1678.
周晓成，王万丽，李立武，等.2020.金沙江—红河断裂带温泉气体地球化学特征[J].岩石学报，36（7）：2197-2214.
周训，曹琴，尹菲，等.2015.四川盆地东部高褶带三叠系地层卤水和温泉的地球化学特征及成因[J].地质学报，89（11）：1908-1920.
周训.1993.四川盆地龙女寺储卤构造深层地下卤水的水文地球化学特征及成因[J].现代地质，7（1）：83-92.
周永胜，何昌荣，杨晓松.2008.中地壳韧性剪切带中的水与变形机制[J].中国科学：地球科学，38（7）：819-832.
Benavente O，Tassi F，Reich M，et al.2016.Chemical and isotopic features of cold and thermal fluids discharged in the Southern Volcanic Zone between 32.5°S and 36°S：Insights into the physical and chemical processes controlling fluid geochemistry in geothermal systems of Central Chile[J].Chemical Geology，420（1）：97-113.
Chen Z，Zhou X C，Du J G，et al.2015.Hydrochemical characteristics of hot spring waters in the Kangding district related to the Lushan MS= 7.0 earthquake in Sichuan，China[J].Natural Hazards & Earth System Ences Discussions，2（6）：1149-1156.
Clark M K，Royden L H.2000.Topographic ooze：Building the eastern margin of Tibet by lower crustal flow[J].Geology，28（8）：703-706.
Craig H.1961.Isotopic Variations in Meteoric Waters[J].Science，133（3465）：1702-1703.
Fournier R O，Rowe J J.1966.The deposition of silica in hot springs[J].29（1）：585-587.
Giggenbach W F.1988.Geothermal solute equilibria.Derivation of Na-K-Mg-Ca geoindicators[J].Geochimica Et Cosmochimica Acta，52（12）：2749-2765. Godard V，Pik R，Cattin R，et al.2009.Late Cenozoic evolution of the central Longmen Shan，eastern Tibet：Insight from（U-Th）/He thermochronometry[J].Tectonics，28（5）：TC5009.
Kouketsu Y，Shimizu I，Wang Y，et al.2017.Raman spectra of carbonaceous materials in a fault zone in the Longmenshan thrust belt，China;comparisons with those of sedimentary and metamorphic rocks[J].Tectonophysics，699（1）：129-145.
Lambrakis N，Kallergis G.2005.Contribution to the study of Greek thermal springs：hydrogeological and hydrochemical characteristics and origin of thermal waters[J].Hydrogeology Journal，13（3）：506-521.
Lin X，Dreger D，Ge H，et al.2018.Spatial and Temporal Variations in the Moment Tensor Solutions of the 2008 Wenchuan Earthquake Aftershocks and Their Tectonic Implications[J].Tectonics，37（3）：989-1005.
Liu Y，Ren H，Wang B.2009.Application of environmental isotopes and tracer techniques to seismic subsurface fluids[J].Earth Science Frontiers，16（1）：369-377.
Luo J，Pang Z H，Kong Y K，et al.2017.Geothermal potential evaluation and development prioritization based on geochemistry of geothermal waters from Kangding area，western Sichuan，China[J].Environmental Earth Ences，76（9）：343.
Mokadem N，Demdoum A，Hamed Y，et al.2016.Hydrogeochemical and stable isotope data of groundwater of a multi-aquifer system：Northern Gafsa basin- Central Tunisia[J].Journal of African Earthences，114（1）：174-191.
Pérez N M，Hernández P A，Igarashi G，et al.2008.Searching and detecting earthquake geochemical precursors in CO2-rich groundwaters from Galicia，Spain[J].Geochemical Journal，42（1）：75-83.
Okuyama Y，Funatsu T，Fujii T，et al.2016.Mid-crustal fluid related to the Matsushiro earthquake swarm（1965-1967）in northern Central Japan：geochemical reproduction，Tectonophysics，679（3）：61-72.
Shi Z，Zhang H，Wang G.2020.Groundwater trace elements change induced by MS5.0 earthquake in Yunnan[J].Journal of Hydrology，581：124424.
Skelton A，Andren M，Kristmannsdottir H，et al.2014.Changes in groundwater chemistry before two consecutive earthquakes in Iceland[J].Nature Geoscience，7（10）：752-756.
Soto-Jiménez M F，Páez-Osuna F.2001.Distribution and Normalization of Heavy Metal Concentrations in Mangrove and Lagoonal Sediments from Mazatlán Harbor（SE Gulf of California）[J].Estuarine Coastal & Shelf Ence，53（3）：259-274.
Tsunogai U，Wakita H.1995.Precursory Chemical Changes in Ground Water：Kobe Earthquake，Japan[J].Science，5220（269）：61-63.
Walraevens K，Bakundukize C，Mtoni Y E，et al.2018.Understanding the hydrogeochemical evolution of groundwater in Precambrian basement aquifers：A case study of Bugesera region in Burundi[J].Journal of Geochemical Exploration，188（3）：24-42. Xu Q，Hoke G D，Jing L Z，et al.2015.Stable isotopes of surface water across the Longmenshan margin of the eastern Tibetan Plateau[J].Geochemistry Geophysics Geosystems，15（8）：3416-3429.
Xu Z，Zhen F，Li SX，et al.2019.The 2018 MS5.9 Mojiang Earthquake：Source model and intensity based on near-field seismic recordings[J].Earth and Planetary Physics，3（3）：88-101.
Zhang Y F，Tan H B，Zhang W J，et al.2015.A New Geochemical Perspective on Hydrochemical Evolution of The Tibetan Geothermal System[J].Geochemistry International，53（12）：1090-1106.
Zhang Y F，Tan H B，Zhang W J，et al.2016.Geochemical constraint on origin and evolution of solutes in geothermal springs in western Yunnan，China[J].Chemie der Erde-Geochemistry，76（1）：63-75.
Zhou X，Jiang C，Zhao J，et al.2015.Occurrence and resource evaluation of the subsurface high-K brines in the Pingluoba brine-bearing structure in western Sichuan Basin[J].Environmental Earth Sciences，73（12）：8565-8574.
Zhou X C，Li C，Ju X M，et al.2010.Origin of Subsurface Brines in the Sichuan Basin[J].Groundwater，35（1）：53-58.
Zhou X C，Liu L，Chen Z，et al.2017.Gas geochemistry of the hot spring in the Litang fault zone，Southeast Tibetan Plateau[J].Applied Geochemistry，79（1）：17-26.
Hydrogeochemical Characteristics of the Hot Springsin the Longmenshan Fault Zone
YAN Yucong1，LIU Fengli1，GUO Lishuang2，ZHOU Xiaocheng1，OUYANG Shupei1，LI Jingchao1，WANG Wanli1
（1.Institute of Earthquake Forecasting，China Earthquake Administration，Beijing 100036，China）
（2.National Institute of Natural Hazards，Ministry of Emergency Management，Key Laboratory ofCrustal Dynamics，China Earthquake Administration，Beijing 100085，China）
Abstract
Hydrogeochemical Monitoring in the regions with strong earthquake activities would facilitate the obtaining of the geochemical information about deep fluids correlated to seismic anomalies.In this paper the hydrogeochemical characteristics of 13 hot springs（wells）in the Longmenshan Fault Zone were investigated，so as to establish the hydrogeochemical cycle model of the hot springs and reveal the relationship between hydrogeochemical variations and seismic activity in the Fault Zone.The results show that：①The hot springs in the Longmenshan Fault Zone were recharged by the meteoric water at elevations ranging from 0.8 to 3.2 km.②The chemical types of the hot springs from north to south，and from west to east along the Longmenshan Fault Zone were respectively carbonic acid，sulfuric acid and chloride.Likewise，the circulation depth of spring water deepened;the water-rock reaction intensified;and enrichment factors of the trace element of the hot springs continuously increased.③After the Wenchuan MS8.0 earthquake and the Lushan MS7.0 earthquake，the dissolved hydrochemical composition got declining in hot springs within the epicenter of 200 kilometers，which could be the result of the permeability decrease and the attenuation of water-rock reaction inside the fault in the healing process of the Longmenshan Fault after earthquake，
Keywords：hot spring;hydrogeochemistry;earthquake;the Longmenshan Fault Zone
收稿日期：2021-01-11.
基金项目：国家重点研发计划（2017YFC1500501）、国家自然科学基金面上项目（41673106）和中国地震局地震预测研究所基本科研业务费（2017IES010205，2016IES010304，2018IEF010104，2020IEF0604，2020IEF0703）联合资助.
第一作者简介：颜玉聪（1996-），女，在读硕士研究生，主要从事地震地球化学研究.E-mail：[email protected].
通訊作者简介：郭丽爽（1983-），高级工程师，主要从事地下流体地球化学方面的研究.E-mail：[email protected].