流体在地球动力学过程中有不可或缺的作用，水文地球化学在地震监测预测方面得到广泛应用，地震常伴有流体各组分地球化学异常变化，地下流体异常对地震孕震区及其外围区域介质的变化有较为灵敏的响应，某些强震前会出现地下流体群体性异常变化^[1-3]。2007年普洱地震震前后温泉水中氟离子发生明显变化^[4]，2008年汶川大地震伴有地下水位和流体地球化学异常现象^[5]，2008年汶川8.0级地震前后的5月1日至7月4日，距震中约400 km的重庆北温泉群所有温泉的水温下降了1℃左右，震后海拔较低泉点K⁺、Na⁺、Ca⁺浓度大幅降低^[6]，2015年阿左旗5.8级地震，周围区域的水化学组成变化较大，对地震响应较为敏感^[7]。

2021年5月21日21时48分37秒（北京时间），在云南省大理白族自治州漾濞县（25.67°N，99.87° E）发生M6.4级地震，震源深度8 km。雷兴林等^[8]对漾濞地震特征和成因进行研究，认为漾濞地震受深部流体作用的影响明显，漾濞地震后的流体特征需要进一步研究。此次漾濞M_S6.4地震发生于滇西北地区，该区域历史上M≥6. 0地震频发，且大多数地震发生于主要断裂及其附近区域，2013年3月以来，滇西北地区M_L≥4. 0地震显著增强，且主要集中分布在维西—乔后断裂、程海断裂及其附近地区。刘自凤等^[9]研究了滇西北地区b值以及其应力张量方差的空间扫描，显示出程海断裂、龙蟠—乔后断裂和维西—乔后断裂南段及邻区低值变化特征相对显著，认为强震发生的可能性大。基于强震危险性的较大背景，从2019年开始，对大理地区主体断裂（维西—乔后断裂、程海断裂、红河断裂等）附近6个温泉每季度定期进行采样观测，主要测定断裂附近温泉水体的水氡、水汞、水体主要阴离子、固体CO₂、电导率、pH值、水温等地球化学量；此次漾濞6.4级地震后，也及时对区域内温泉开展了现场采样观测。共获得6个温泉震前震后8组连续流动观测数据。因此，本文通过分析漾濞6.4级地震前后大理地区温泉水化学组分，分析该地区温泉水文地球化学变化规律，探索温泉水地球化学特征变化与漾濞地震的关系。研究成果不仅对认识该区温泉水地球化学特征及其演化有重要意义，而且对探索和改进利用该区温泉水化学组分进行地震预测的方法也具有重要意义。

从构造上看，漾濞M6.4地震发生在川滇菱形块体西侧的印支块体上，印支地块自西向东分为三个次级构造单元：临沧花岗岩带、澜沧江断裂带和兰坪—思茅盆地。兰坪—思茅盆地内部构造复杂，发育了早期的逆冲断层、推覆构造与晚期的走滑构造断层^[10-11]。此次地震发生在兰坪—思茅盆地内，离东侧维西—乔后—巍山断裂约20 km （图1）。维西—乔后—巍山断裂带，北起雪龙山东廘白济汛一带，经维西、通甸、乔后，止于巍山盆地南端，走向北北西，倾向NE/NW，长约280 km。大致以玉狮场、平坡为界，可将断裂分为北段、中段和南段3段，是一条走滑兼正断的活动断裂^[12-13]。该区域是青藏高原物质向东南侧挤出的通道^[14]，也是特提斯—喜马拉雅构造域东南部弧形构造转折处，处于川滇菱形块体西南边界红河断裂带与近南北向丽江—大理断裂系统的交汇区及应力集中带^[15]。区域内强震频发，图1b给出了1925年有仪器记录以来的M≥6.0级地震M-T图，包括漾濞M6.4级地震在内，共发生了22次6级以上地震，占同时期云南地区地震的34%，最近的一次强震频发时段是1996—2003年，期间发生1次丽江7.0级地震，4次M≥6.0级地震。漾濞M6.4级地震与上一次地震相隔12年。大理地区位于红河断裂带北端，构造上表现为一系列近SN向及NNE向的活动断裂，如离此次漾濞地震约20 km的维西—乔后—巍山断裂，中新世以来的断陷盆地从南北两个方向逐渐向NW走向的红河断裂交会，构成断裂东北侧典型的滇西北伸展裂陷区^[16] 。

图1 区域构造与历史地震震中分布Fig.1 Regional tectonics and epicenters distribution of the historical earthquakes

此次定期观测的6个温泉为：洱源县牛街乡练渡村温泉（以下简称洱源温泉），永胜县片角镇热河村温泉（以下简称永胜温泉），弥渡县弥城镇高芹温泉（以下简称高芹温泉），弥渡县弥城镇白总旗温泉（以下简称白总旗温泉），云龙县漕涧下澡堂温泉（以下简称云龙温泉），剑川县老君山新和村温泉（以下简称剑川温泉），位置如图2所示，除剑川温泉离漾濞6.4级地震震中距为104 km，其余温泉震中距都在100 km内（表1）。洱源温泉位于红河断裂和丽江—剑川断裂交互处，该区内构造体系呈SN向展布，为地下热水提供循环通道;洱源盆地基底有大面积的多期岩浆侵入玄武岩，为地下热水提供热源;基底上覆中、新生界碎屑岩和冲湖积黏土层组成良好隔热保温层，三者构成良好的储热构造^[17]。弥渡高芹温泉、白总旗温泉近红河断裂南段，弥渡盆地及周边山区的碳酸盐岩破碎带地区是该区域地下水的补给区，大气降水沿岩溶裂隙及构造裂隙向深部运移获得热量形成深远高温热水，在盆地总体的由北至南的地下水流动方向上，运移至温泉出露区域，受致密岩体或者压性断裂带阻拦，热水沿构造上升出露成泉^[18] 。云龙温泉位于兰坪—永平断裂西侧，该区位于怒江及澜沧江沿岸地下热水密集区，属于中低温水热型地热类型，地壳地温梯度为4℃/100 m^[19]。区域内出露地层的岩性以紫红色砂岩、粉砂岩、粉砂质泥岩和泥岩为特征，统称红层^[20]。剑川温泉位于维西—乔后断裂附近，出露于黄色砂岩的裂缝中，处于兰坪—思茅盆地的边缘^[21]，位于我国著名的滇藏地热带内，地热活动强烈，构造复杂，不同构造体系的断裂互相重合、交接、穿插，为地下热水将深部热量带至地壳浅表、孕育高温水热系统、产生强烈水热活动提供了重要的构造条件^[22]。永胜温泉位于程海—宾川断裂红岩坡段，程海—宾川断裂可分为八段，为最南段，区域内活动断裂发育，程海断裂纵贯^[23]。

图2 温泉取样点分布图Fig.2 Distribution map of hot spring sampling points

表1 温泉基本信息表Table 1 Basic information sheet of hot spring

对6个温泉点分别于2020年5、7、9、11月， 2021年3、5、8、11月共8次观测采样，获得温泉水氡、水汞、固体CO₂ 、阴离子、pH、电导率的连续8次观测数据，以及温泉阳离子数据（无连续数据）。水氡、水汞分别采用RAD7测氡仪和RA-915W 测汞仪野外现场测定，其余取样实验室测定；温泉固体CO₂利用沉淀法获取，取50 mL温泉水样，加入过量饱和Ba（OH）₂溶液，将Ba CO₃沉淀洗净过滤、烘干，计算固体 CO₂ 值；用intergration离子色谱仪测定常量元素，仪器RSD%在5%以内，电导率为0. 1 nS/cm，HCO₃^-用滴定法测定。数据见表2、表3。

表2 温泉采样时间、水化学数据Table 2 Sampling time and water chemistry data of hot spring

表3 温泉水样主要阴阳离子组分（2021年8月） Table 3 Main anion and cation components of hot spring water samples（August 2021)

Piper图解和Na-K-Mg三角图解常被用来评价水—岩化学平衡状态、区分不同类型水样以及判断地下水循环深度^[24] 。从Piper图可见（图3），洱源温泉属于HCO₃·SO₄-Na型水，白总旗温泉属于HCO₃·Cl-Na型水，高芹、剑川、永胜、云龙温泉均属于HCO₃-Na型水。6个温泉水常量离子组分测试结果（表2 、表3）显示，6个温泉所含阴离子均以HCO₃^-占绝对优势，阳离子均以Na⁺占绝对优势，这可能和钠的强溶解性和强迁移性有关。洱源，高芹，白总旗，剑川温泉Na⁺占阳离子组分60%左右，且富含K⁺和Ca₂ ⁺，反映出这4个温泉可能受到岩浆岩和碳酸盐岩的共同影响，大理地区广泛分布石炭—二叠系海相沉积灰岩，较高的HCO₃^－、Ca₂ ⁺含量可能与灰岩的溶解有关，该区同时分布喜山期碱性岩浆岩，这或许导致了较高的Na⁺、K⁺；不同的是云龙和永胜温泉主体为Na⁺，但K⁺，Ca²⁺，Mg²⁺等阳离子基本很少，这可能和离子交换作用有关， K⁺， Ca²⁺， Mg²⁺吸附能力强于Na⁺，围岩中的钠离子置换钙镁离子。从样品的Na-K-Mg平衡图解（图4）可见，永胜温泉水属于部分成熟水，表明其发生部分水岩反应，温泉水与围岩的反应达到部分平衡状态；其余温泉为未成熟水，且靠近Mg^1/2端，表明其水岩反应较强，属于浅层水区域。6个温泉的pH为6.64～7.54之间，中性；电导率为0.22～1.34之间，6个温泉pH和电导率多次测值变化不大，数据较稳定。

图3 温泉水的Piper图解Fig.3 Piper diagram of hot spring water

图4 温泉水的Na-K-Mg图解Fig.4 Na-K-Mg ternary diagram of hot spring water

根据前文分析显示，漾濞6.4级地震前后区域内温泉水的化学参数主要表现在Cl^-和SO₄^2-浓度以及温泉水氡、水汞含量出现不同程度的变化，Cl^-和SO₄^2-变化最为显著，如剑川温泉Cl^-，永胜温泉SO₄^2-，相对震前采集的水样，下降幅度都达到98%以上，温泉的水温和固体CO₂变化不明显。总体显示出漾濞6.4级地震前大多数温泉Cl^-、SO₄^2-、水氡、水汞浓度震前出现明显上升，震后各组分含量显著下降。许多学者对地下水前兆异常的机制进行了研究，提出了许多解释地下水前兆异常机制的模型^[[27-30]，其中含水层破裂/流体混合模型（AB/FM模型）认为隔离含水层之间的水文屏障的前兆破裂，从而使不同含水不层发生混合，而这种前兆破裂是由地震前与应力积累相关的地壳膨胀引起的。大理地区6个温泉的水化学变化现象可以用AB/FM模型模型机制解释。

前人研究滇西北地区b值和应力张量方差的空间扫描结果显示程海断裂、龙蟠—乔后断裂和维西—乔后断裂南段及邻区低值变化特征相对显著，这可能是区域应力水平增强的表现^[9]，此次研究Cl^-和SO₄^2-震前升高，震后剧烈下降的剑川，永胜，白总旗，高芹等温泉刚好在这些断裂附近。剑川温泉位于维西—乔后断裂附近，永胜温泉位于程海—宾川断裂红岩坡段，白总旗和高芹温泉近红河断裂南段。既区域应力积累增加，孕震区的岩石的内部结构经历岩石的弹性压缩、原生裂隙的闭合、新微裂隙的产生、裂隙的弥散发展、裂隙的集结到岩石的宏观破裂的演化过程^[26]，演化过程中深部流体扩散上升、不同层位水体的混合以及孔隙流体的受挤压流体岩石间的相互作用导致地下流体水化组分的异常变化^[27-30]。洱源温泉属于 HCO₃ · SO₄-Na型水，白总旗温泉属于HCO₃·Cl-Na型水，高芹、剑川、永胜、云龙温泉均属于HCO₃-Na型水，虽然这些温泉有一定含量的Cl^-和SO₄^2-，却不是温泉的主量化学元素，而SO₄^2-及Cl^-除少量的地表来源外，主要为中、深部来源^[7]，这些温泉水较高的Cl^-、SO₄^2-浓度可能主要来源于深部高Cl^-、SO₄^2-浓度水的混入，地震孕育过程中，区域压应力增强导致高Cl^-、SO₄^2-浓度的深部水供给量的增加所致，地震后，随着区域应力的调整，深部供给的流体量减少，各温泉的离子浓度下降。

（1）大理地区主要活动断裂附近6个温泉主要为中高温温泉，洱源温泉温度最高，约74.5℃。温泉水的常量离子分析表明，阳离子主要以Na⁺为主，Ca²⁺、Mg²⁺含量不高且差异较大，位于程海—宾川断裂的永胜温泉和兰坪—永平断裂附近的云龙温泉Ca²⁺、 Mg²⁺含量极少，阴离子以HCO₃^-为主，较高的HCO₃^-和固体CO₂含量可能受大理地区广泛分布的石炭—二叠系海相沉积灰岩影响。洱源温泉属于HCO₃·SO₄-Na型水，白总旗温泉属于HCO₃·Cl-Na型水，高芹、剑川、永胜、漕涧温泉均属于HCO₃-Na型水。永胜温泉水属于部分成熟水，其余温泉为未成熟水。区域内温泉有较高水氡水汞背景值，不同断裂带及同一断裂带上不同区位的温泉水汞，水氡背景值各不相同，具有一定的区域性。

（2）漾濞6.4级地震前后大理地区主要断裂附近温泉水的Cl^-、SO₄^2-、水汞、水氡含量有明显变化，白总旗温泉Cl^-，永胜、洱源温泉SO₄^2-以及高芹温泉水汞含量在漾濞6.4级地震前明显上升，震后含量显著下降；可能是震前区域挤压应力增强过程高Cl^-和SO₄^2-，水体的混入，震后区域应力的调整，中深部流体量补给量减少，各化学组分多数含量下降；温泉水温和固体CO₂含量变化不明显。