2024年4月23日发(作者：联想m6100t配置)

洱源温泉水氡预测效能检验及其水化学特征

作者：钟骏 王博 晏锐 李庆

来源：《地震研究》2021年第04期

摘要：利用Molchan图表法对洱源温泉水氡观测资料进行预测效能检验，并开展了温泉水

化学组分、氢氧同位素和水循环特征研究。结果表明：①DMolehan图表法检验结果显示，洱

源水氡出现高值时，对观测点附近中强以上地震反应较灵敏，优势对应地震时间段为90d以

内;②洱源温泉井的水化学类型为HCO3·SO4-Na，与周边地表水体及浅层地下水无直接水力联

系;③洱源温泉井水属于“部分平衡水”，具有深部来源特征。

关键词：洱源水氡;预测效能;水化学特征

中图分类号：P315.72 文献标识码：A文章编号：1000-0666（2021）04-0564-08

0引言

氡作为地球内部放射性元素铀、镭衰变的产物，是地震科学研究中观测历史最久、研究程

度最深的地下流体前兆手段之一（刘伯礼等，1975;刘耀炜，2006）。前人观测实践表明，氡

的浓度变化是地球深部信息的重要表现，特别是地震孕育及发生阶段，水中氡的浓度会出现不

同程度的异常变化（付虹等，1997;Ge et al ，2014;孙小龙等，2020）。然而，氡浓度的变化与

地震活动的关系较为复杂，变化形态往往也不一致。如龙陵邦腊掌温泉水氡在多次中强地震前

出现低值波动异常（李加稳，孙文波，2004），但在汶川M8.0地震前，位于龙门山断裂带附

近的多个水氡观测项目的自相关系数则出现了明显趋势性上升的异常现象（晏锐等，2011）。

地下水溶解氡浓度在地震前上升或者下降的异常变化，在一定程度上反映了应力作用条件下所

引起的水动力条件和局部水文地球化学环境的改变（Muir-Wood，King，1993）。

针对氡的异常成因，前人提出了诸如膨胀扩散、破裂混合、热液爆沸、氡团混入、超声振

动、岩浆冲熔及水动力作用等机理和模式（Birchard，Libby ，1978;郭增建等，1979;冯玮等，

1981;李宣瑚，1981;Thomas ，1988;蔣凤亮等，1989;车用太等，1997）。在此基础上，张慧等

（2005）提出了孕震体介质微破裂生长的观点，进一步揭示了水氡前兆场的时空变化特征及其

物理机制。温泉作为连接地球深部与浅部的重要通道，与构造活动和地震活动都有较好的联

系，能反映较多的地球深部信息，而温泉水氡观测是十分理想的前兆观测手段（苏鹤军等，

2020）。因此，深入分析温泉水氡的水化学特征及其与地震活动的关系，对使用温泉水氡进行

地震预测实践具有重要意义。

日常震情跟踪发现，洱源温泉水氡与周边地区一定范围内的中强地震存在一定的关联性，

而针对该泉点资料的归纳梳理和系统研究较少。因此，本文以洱源温泉为研究对象，使用

Molchan图表法对洱源温泉水氡进行映震效能检验，从水化学组成和氢氧同位素的角度分析洱

源温泉水化学特征和氡的来源问题，并对其前兆机理进行了讨论。

1观测点概况

洱源温泉水氡（以下简称洱源水氡）观测点位于洱源断陷盆地，构造上属于扬子准地台、

松潘—甘孜褶皱带与三江褶皱带的交界部位（房艳国等，2015），主要出露地层岩性为二叠

系、三叠系灰岩（李晓露，2017）。区内活动断裂发育，红河断裂北段、鹤庆—洱源断裂、龙

蟠—乔后断裂和通甸—巍山断裂等多条断裂交汇于此（图1），构造活动较为强烈，是滇西地

区中强地震的密集活动区域（林爱文，1997）。

洱源水氡自2006年开始在兴兴温泉取样，水温为62.0℃左右，流量为（0.75±0.03）L/s

（段美芳等，2020），观测仪器为FD-105K型测氡仪，测量灵敏度为（0.37±0.15）Bq/

（L·格·min），采样率为日值，无明显年变规律，受气象因素影响较小，浓度变化范围为

38.2～50.1Bq/L。2006—2007年洱源水氡因影响观测稳定因素的出现，观测数据出现大幅波

动，之后趋于稳定并且连续可靠。在2008年8月21日盈江M6.0、2009年7月9日姚安

M6.0、2012年6月24日宁蒗M5.7、2013年8月31日德钦M5.9、2015年10月31日昌宁

M5.1等多次地震前，洱源水氡都出现明显的异常变化，表现为水氡浓度值快速上升或波动幅

度增大，显示出较强的映震能力，说明其异常可能与地震活动存在一定联系（图2）。

2地震预测效能检验

虽然洱源水氡在一些地震前出现过明显异常变化，但异常与地震并非一一对应，很可能存

在其他影响因素，导致地震前兆异常被掩盖或非地震前兆信息被放大。因此，需要使用有效的

方法对观测资料进行预处理，从而提取有效异常指标，并依据预测效能评价方法——Molchan

图表法（Molchan，1990）来检验这种异常提取方法的映震水平，获取最优预测指标，为更好

地利用观测资料进行地震监测预测提供方法依据。

根据模糊数学原理，引入从属函数方法中的斜率K和相关系数R的乘积Y来反映观测数

据变化斜率随时间的变化程度（冯德益等，1984）：

式中：K为观测数据的斜率，反映数据的变化速率;R1为观测数据的自相关系数，反映数

据内在质量的好坏。K与R可以利用指定窗长内连续观测数据及其相应的时间序列作一元线性

回归分析求得，即：

式中：n为滑动窗长，即一元线性回归分析时样本的个数，通常取n=5;Y为观测数据;T为

与观测数据相应的时间序列。

当|R不变时，Y随|K的增大而增大，表明观测值变化速率越大，异常的可能性越高;当K

不变时，R越大，则数据的相关性越好，说明|K作为判别异常的信度越大。

本文选取20082018年洱源水氡的观测数据进行分析，并依据震级与震中距一般对应关

系：M5.0～5.9地震的震中距小于250km、M6.0～6.9地震的震中距小于300km（孙小龙等，

2018），筛选出该时间段内洱源水氡观测点250km范围内发生的M3≥5.5地震共21次，其

中，M5.5～5.9地震17次、M6.0～6.9地震4次（图1），地震参数详见表1。在实际计算过程

中，为减小数据中的系统误差和偶然误差，本文采用45天滑动平均对洱源水氡的原始数据进

行趋势提取，如图3所示，可以看出水氡浓度呈现明显的上升-下降变化，预处理后曲线相对

于原始曲线的震前异常特征更为突出、更易于识别。

Molchan图表法通过对预测结果与观测目标地震的差异度进行检验，直观反映观测数据的

预报效能，得到最佳阈值所对应的异常识别指标（蒋长胜等，2011;孙小龙等，2018）。每一

个阈值对应一组漏报率r和异常时空占有率v，r-v值连线与横、纵坐标围成的面积代表检验的

预测效果，且与面积大小呈反比。此外，报准率h所对应的显著性水平α越低、概率增益Gain

越大时，预测效果越好。

本文利用Molchan图表法对洱源水氡預处理后的斜率数据（图3c）与周边地震关系进行

检验，得到该测项的整体预测效能、最佳预测指标（即异常的最佳判别指标）及优势对应时

间。图3c中空心五角星为选取的21个地震，其中红色代表有效预测期内报准的地震，蓝色代

表有效预测期外漏报的地震，红色虚线为最佳阈值线（0.0175）;图4a为洱源水氡Y值整体预

报效能与概率增益关系图;图4b为洱

源水氡Y值漏报率v与时间占有率r关系曲线及不同显著性水平α等值线。图中红色粗阶

梯线为预测阈值从大到小滑动得到的r-v曲线，蓝色圆圈代表最接近原点（0，0）的r-v值，

其阈值可视为最佳阈值。图中右下角（1，0）点表示地震全部报准，但其时间占有率也最大，

覆盖数据整个时间段;左上角（0，1）点表示地震全部漏报，其时间占有率最小，相当于未做

出预测。

由图4a可见，洱源水氡在不同预测时长的概率增益均大于1，预测效能基本在0.6以上，

具有显著的统计和检验意义。尤其是短期（90d内），概率增益均大于2，预测效能均在0.7

以上，表明洱源水氡短期预测效果更好。因此选取预测时长80d进行计算，如图4b所示，在

80d所示位置（蓝色圆圈），21个地震中，报准15个、漏报6个，报准率约71%，时间占有

率约0.33，置信水平较高，说明预测效果较好。因此，洱源水氡的优势对应地震时间段为90d

以内，在预测时长为80d时，最佳阈值为0.0175，可作为相应时间段的“警报”阈值。

3水化学特征及温泉来源分析

水化学特征分析能够探明地下流体的来源与补给过程，是区别地下流体是否受构造活动影

响的重要技术方法之一（蒋凤亮等，1989）。地壳深部构造活动或浅层地表水补给会影响地下

流体的流量、温度和水化学组分的变化。为了判断氡的来源及其与深部构造活动的关系，笔者

于2018年5月8日对洱源水氡（SD）、洱源老街（LJ）和温泉疗养院（LYY）3处温泉点，

以及滇20井（D20）和附近河流（HL）等浅表水，共计5个水体进行了水样采集，采样点空

间分布如图1所示。样品检测由中国地震局地壳动力学实验室完成，使用Dionex ICS -2100离

子色谱仪进行离子组分检测，包括Ca2+、Mg2+、Na+、K+、SO、Cl-和NO3，浓度检测限为

0.01mg/L;HCO3采用酸碱中和滴定法求得;氢氧同位素δ0与8D使用LGR LWIA -24-EP激光液

态水稳定同位素分析仪测定，精度分别高于0.1%和0.4%。水化学测试结果见表2。

由表2及图5可以看出，洱源水氡和洱源老街优势阳离子成分为Na+、优势阴离子成分为

HCO3和SO4，两处水体的Na+、HCO3和SO4毫克当量百分比分别为92.5%、67.7%、29.4%

以及90.5%、67.2%、29.8%，水化学类型均为HCO3·SO4-Na;温泉疗养院优势阳离子成分为

Ca2+和Na+、优势阴离子成分为HCO3和SO，其各离子毫克当量百分比分别为45.3%、

40.8%、57.5%和41.3%，水化学类型为HCO3·SO4-Ca·Na;滇20井优势阳离子成分为Na+、优

势阴离子成分为HCO3，其离子毫克当量百分比分别为85.9%和93.8%，水化学类型为HCO3-

Na;河流优势阳离子成分为Ca2+和Mg2+、优势阴离子成分为HCO3，其离子毫克当量百分比

分别为60.6%、30.5%和89.8%，水化学类型为HCO3-Ca·Mg。由上可见，温泉水的大多数离

子含量明显高于浅表水，表明其来源较深、离子反应更强烈，从而导致不同水体的水化学类型

差别较大，呈现出明显分区特征。

从图5可见，洱源水氡与洱源老街在菱形区域中的投影点几乎重合，表明其来源较为一致;

而温泉疗养院、滇20井和河流落在不同的位置，说明其来源和形成过程不同，与洱源水氡与

洱源老街没有直接的水力联系。分析认为，洱源地区温泉水水化学组分的形成过程是地下热水

与围岩发生溶滤作用和离子交换作用的结果。由于该地区发育大量二叠系、三叠系灰岩，三叠

系和白垩系的碎屑岩中含有石膏（张春山等，2003），当地下热水流经灰岩层及膏盐层时与围

岩发生溶滤作用，使得Ca2+、HCO3和SO大量进入水中。随着循环深度和循环时间的增加，

CO2和O2等气体耗尽，溶滤作用减弱，离子交换作用开始增强，Ca2+不断被Na+替换，逐渐

形成HCO3·SO4-Ca·Na或HCO3·SO4-Na型水（王大纯，1986;吴桂武等，2015）。

Na-K-Mg三角图通过水-岩平衡状态，可以判断地下水循环深度及混合程度。水-岩平衡程

度越高，地下水循环越深或循环周期越长（Giggen-bach，1988）。由图6可见，所有采样点均

靠近右下角Mg端元，未达到饱和状态。其中，洱源水氡属于“部分平衡水”，洱源老街接近

“部分平衡水”，两者接近重合，说明其水-岩平衡状态相似，来源于深部的可能性更大。且相

比于后者，洱源水氡的水-岩反应程度更高，其接受大气降水补给或浅部冷水混合更少。结合

表2及水化学结果分析，可以清楚看到温泉疗养院具有深部来源特征（SO含量较高），但温

泉疗养院与浅表水体（滇20井、河流）同属“未成熟水”，均表现出大气降水补给的特征，说

明在循环过程中明显受到了地表水或浅部冷水混入影响（Ca2含量较高）。研究区位于多条断

裂交汇处，构造裂隙十分发育，特别是红河断裂北段和鹤庆—洱源断裂具有活动性强、切割深

度大的特点，有利于地下水的深循环和深部热物质的对流（颜鹍等，1997;王云等，2019）。

大气降水经洱源盆地周边山区入渗补给，沿构造裂隙向地下运移，经过深循环或深部热物质对

流后，在水头压力差或阻水断裂的阻挡作用下，沿压力较小的断裂破碎带向地表运移并与浅层

冷水发生不同程度的混合，进而表现出不同的水-岩平衡状态。

Molchan图表法通过对预测结果与观测目标地震的差异度进行检验，直观反映观测数据的

预报效能，得到最佳阈值所对应的异常识别指标（蒋长胜等，2011;孙小龙等，2018）。每一

个阈值对应一组漏报率r和异常时空占有率v，r-v值连线与横、纵坐标围成的面积代表检验的

预测效果，且与面积大小呈反比。此外，报准率h所对应的显著性水平α越低、概率增益Gain

越大时，预测效果越好。

本文利用Molchan图表法对洱源水氡预处理后的斜率数据（图3c）与周边地震关系进行

检验，得到该测项的整体预测效能、最佳预测指标（即异常的最佳判别指标）及优势对应时

间。图3c中空心五角星为选取的21个地震，其中红色代表有效预测期内报准的地震，蓝色代

表有效预测期外漏报的地震，红色虚线为最佳阈值线（0.0175）;图4a为洱源水氡Y值整体预

报效能与概率增益关系图;图4b为洱

源水氡Y值漏报率v与时间占有率r关系曲线及不同显著性水平α等值线。图中红色粗阶

梯线为预测阈值从大到小滑动得到的r-v曲线，蓝色圆圈代表最接近原点（0，0）的r-v值，

其阈值可视为最佳阈值。图中右下角（1，0）点表示地震全部报准，但其时间占有率也最大，

覆盖数据整个时间段;左上角（0，1）点表示地震全部漏报，其时间占有率最小，相当于未做

出预测。

由图4a可见，洱源水氡在不同预测时长的概率增益均大于1，预测效能基本在0.6以上，

具有显著的统计和检验意义。尤其是短期（90d内），概率增益均大于2，预测效能均在0.7

以上，表明洱源水氡短期预测效果更好。因此选取预测时长80d进行计算，如图4b所示，在

80d所示位置（蓝色圆圈），21个地震中，报准15个、漏报6个，报准率约71%，时间占有

率约0.33，置信水平较高，说明预测效果较好。因此，洱源水氡的优势对应地震时间段为90d

以内，在預测时长为80d时，最佳阈值为0.0175，可作为相应时间段的“警报”阈值。

3水化学特征及温泉来源分析

水化学特征分析能够探明地下流体的来源与补给过程，是区别地下流体是否受构造活动影

响的重要技术方法之一（蒋凤亮等，1989）。地壳深部构造活动或浅层地表水补给会影响地下

流体的流量、温度和水化学组分的变化。为了判断氡的来源及其与深部构造活动的关系，笔者

于2018年5月8日对洱源水氡（SD）、洱源老街（LJ）和温泉疗养院（LYY）3处温泉点，

以及滇20井（D20）和附近河流（HL）等浅表水，共计5个水体进行了水样采集，采样点空

间分布如图1所示。样品检测由中国地震局地壳动力学实验室完成，使用Dionex ICS -2100离

子色谱仪进行离子组分检测，包括Ca2+、Mg2+、Na+、K+、SO、Cl-和NO3，浓度检测限为

0.01mg/L;HCO3采用酸碱中和滴定法求得;氢氧同位素δ0与8D使用LGR LWIA -24-EP激光液

态水稳定同位素分析仪测定，精度分别高于0.1%和0.4%。水化学测试结果见表2。

由表2及图5可以看出，洱源水氡和洱源老街优势阳离子成分为Na+、优势阴离子成分为

HCO3和SO4，两处水体的Na+、HCO3和SO4毫克当量百分比分别为92.5%、67.7%、29.4%

以及90.5%、67.2%、29.8%，水化学类型均为HCO3·SO4-Na;温泉疗养院优势阳离子成分为

Ca2+和Na+、优势阴离子成分为HCO3和SO，其各离子毫克当量百分比分别为45.3%、

40.8%、57.5%和41.3%，水化学类型为HCO3·SO4-Ca·Na;滇20井优势阳离子成分为Na+、优

势阴离子成分为HCO3，其离子毫克当量百分比分别为85.9%和93.8%，水化学类型为HCO3-

Na;河流优势阳离子成分为Ca2+和Mg2+、优势阴离子成分为HCO3，其离子毫克当量百分比

分别为60.6%、30.5%和89.8%，水化学类型为HCO3-Ca·Mg。由上可见，温泉水的大多数离

子含量明显高于浅表水，表明其来源较深、离子反应更强烈，从而导致不同水体的水化学类型

差别较大，呈现出明显分区特征。

从图5可见，洱源水氡与洱源老街在菱形区域中的投影点几乎重合，表明其来源较为一致;

而温泉疗养院、滇20井和河流落在不同的位置，说明其来源和形成过程不同，与洱源水氡与

洱源老街没有直接的水力联系。分析认为，洱源地区温泉水水化学组分的形成过程是地下热水

与围岩发生溶滤作用和离子交换作用的结果。由于该地区发育大量二叠系、三叠系灰岩，三叠

系和白垩系的碎屑岩中含有石膏（张春山等，2003），当地下热水流经灰岩层及膏盐层时与围

岩发生溶滤作用，使得Ca2+、HCO3和SO大量进入水中。随着循环深度和循环时间的增加，

CO2和O2等气体耗尽，溶滤作用减弱，离子交换作用开始增强，Ca2+不断被Na+替换，逐渐

形成HCO3·SO4-Ca·Na或HCO3·SO4-Na型水（王大纯，1986;吴桂武等，2015）。

Na-K-Mg三角图通过水-岩平衡状态，可以判断地下水循环深度及混合程度。水-岩平衡程

度越高，地下水循环越深或循环周期越长（Giggen-bach，1988）。由图6可见，所有采样点均

靠近右下角Mg端元，未达到饱和状态。其中，洱源水氡属于“部分平衡水”，洱源老街接近

“部分平衡水”，两者接近重合，说明其水-岩平衡状态相似，来源于深部的可能性更大。且相

比于后者，洱源水氡的水-岩反应程度更高，其接受大气降水补给或浅部冷水混合更少。结合

表2及水化学结果分析，可以清楚看到温泉疗养院具有深部来源特征（SO含量较高），但温

泉疗养院与浅表水体（滇20井、河流）同属“未成熟水”，均表现出大气降水补给的特征，说

明在循环过程中明显受到了地表水或浅部冷水混入影响（Ca2含量较高）。研究区位于多条断

裂交汇处，构造裂隙十分发育，特别是红河断裂北段和鹤庆—洱源断裂具有活动性强、切割深

度大的特点，有利于地下水的深循环和深部热物质的对流（颜鹍等，1997;王云等，2019）。

大气降水经洱源盆地周边山区入渗补给，沿构造裂隙向地下运移，经过深循环或深部热物质对

流后，在水头压力差或阻水断裂的阻挡作用下，沿压力较小的断裂破碎带向地表运移并与浅层

冷水发生不同程度的混合，进而表现出不同的水-岩平衡状态。

Molchan图表法通过对预测结果与观测目标地震的差异度进行检验，直观反映观测数据的

预报效能，得到最佳阈值所对应的异常识别指标（蒋长胜等，2011;孙小龙等，2018）。每一

个阈值对应一组漏报率r和异常时空占有率v，r-v值连线与横、纵坐标围成的面积代表检验的

预测效果，且与面积大小呈反比。此外，报准率h所对应的显著性水平α越低、概率增益Gain

越大时，预测效果越好。

本文利用Molchan图表法对洱源水氡预处理后的斜率数据（图3c）与周边地震关系进行

检验，得到该测项的整体预测效能、最佳预测指标（即异常的最佳判别指标）及优势对应时

间。图3c中空心五角星为选取的21个地震，其中红色代表有效预测期内报准的地震，蓝色代

表有效预测期外漏报的地震，红色虚线为最佳阈值线（0.0175）;图4a为洱源水氡Y值整体预

报效能与概率增益关系图;图4b为洱

源水氡Y值漏报率v与时间占有率r关系曲线及不同显著性水平α等值线。图中红色粗阶

梯线为预测阈值从大到小滑动得到的r-v曲线，蓝色圆圈代表最接近原点（0，0）的r-v值，

其阈值可视为最佳阈值。图中右下角（1，0）点表示地震全部报准，但其时间占有率也最大，

覆盖数据整个时间段;左上角（0，1）点表示地震全部漏报，其时间占有率最小，相当于未做

出预测。

由图4a可见，洱源水氡在不同预测时长的概率增益均大于1，预测效能基本在0.6以上，

具有显著的统计和检验意义。尤其是短期（90d内），概率增益均大于2，预测效能均在0.7

以上，表明洱源水氡短期预测效果更好。因此选取预测时长80d进行计算，如图4b所示，在

80d所示位置（蓝色圆圈），21个地震中，报准15个、漏报6个，报准率约71%，时间占有

率约0.33，置信水平较高，说明预测效果较好。因此，洱源水氡的优势对应地震时间段为90d

以内，在预测时长为80d时，最佳阈值为0.0175，可作为相应时间段的“警报”阈值。

3水化学特征及温泉来源分析

水化学特征分析能够探明地下流体的来源与补给过程，是区别地下流体是否受构造活动影

响的重要技术方法之一（蒋凤亮等，1989）。地壳深部构造活动或浅层地表水补给会影响地下

流体的流量、温度和水化学组分的变化。为了判断氡的来源及其与深部构造活动的关系，笔者

于2018年5月8日对洱源水氡（SD）、洱源老街（LJ）和温泉疗养院（LYY）3处温泉点，

以及滇20井（D20）和附近河流（HL）等浅表水，共计5个水体进行了水样采集，采样点空

间分布如图1所示。样品检测由中国地震局地壳动力学实验室完成，使用Dionex ICS -2100离

子色谱仪进行离子组分检测，包括Ca2+、Mg2+、Na+、K+、SO、Cl-和NO3，浓度检测限为

0.01mg/L;HCO3采用酸碱中和滴定法求得;氢氧同位素δ0与8D使用LGR LWIA -24-EP激光液

态水稳定同位素分析仪测定，精度分别高于0.1%和0.4%。水化学测试结果见表2。

由表2及图5可以看出，洱源水氡和洱源老街优势阳离子成分为Na+、优势阴离子成分为

HCO3和SO4，两处水体的Na+、HCO3和SO4毫克当量百分比分别为92.5%、67.7%、29.4%

以及90.5%、67.2%、29.8%，水化学类型均为HCO3·SO4-Na;温泉疗养院优势阳离子成分为

Ca2+和Na+、优势阴离子成分为HCO3和SO，其各离子毫克当量百分比分别为45.3%、

40.8%、57.5%和41.3%，水化学类型为HCO3·SO4-Ca·Na;滇20井优势阳离子成分为Na+、优

势阴离子成分为HCO3，其离子毫克当量百分比分别为85.9%和93.8%，水化学类型为HCO3-

Na;河流优势阳离子成分为Ca2+和Mg2+、优势阴离子成分为HCO3，其离子毫克当量百分比

分别为60.6%、30.5%和89.8%，水化学类型为HCO3-Ca·Mg。由上可见，温泉水的大多数离

子含量明显高于浅表水，表明其来源较深、离子反应更强烈，从而导致不同水体的水化学类型

差别较大，呈现出明显分区特征。

从图5可见，洱源水氡与洱源老街在菱形区域中的投影点几乎重合，表明其来源较为一致;

而温泉疗养院、滇20井和河流落在不同的位置，说明其来源和形成过程不同，与洱源水氡与

洱源老街没有直接的水力联系。分析认为，洱源地区温泉水水化学组分的形成过程是地下热水

与围岩发生溶滤作用和离子交换作用的结果。由于该地区发育大量二叠系、三叠系灰岩，三叠

系和白垩系的碎屑岩中含有石膏（张春山等，2003），当地下热水流经灰岩层及膏盐层时与围

岩发生溶滤作用，使得Ca2+、HCO3和SO大量进入水中。随着循环深度和循环时间的增加，

CO2和O2等气体耗尽，溶滤作用减弱，离子交换作用开始增强，Ca2+不断被Na+替换，逐渐

形成HCO3·SO4-Ca·Na或HCO3·SO4-Na型水（王大纯，1986;吴桂武等，2015）。

Na-K-Mg三角图通过水-岩平衡状态，可以判断地下水循环深度及混合程度。水-岩平衡程

度越高，地下水循环越深或循环周期越长（Giggen-bach，1988）。由图6可见，所有采样点均

靠近右下角Mg端元，未达到饱和状态。其中，洱源水氡属于“部分平衡水”，洱源老街接近

“部分平衡水”，两者接近重合，说明其水-岩平衡状态相似，来源于深部的可能性更大。且相

比于后者，洱源水氡的水-岩反应程度更高，其接受大气降水补给或浅部冷水混合更少。结合

表2及水化学结果分析，可以清楚看到温泉疗养院具有深部来源特征（SO含量较高），但温

泉疗养院与浅表水体（滇20井、河流）同属“未成熟水”，均表现出大气降水补给的特征，说

明在循环过程中明显受到了地表水或浅部冷水混入影响（Ca2含量较高）。研究区位于多条断

裂交汇处，构造裂隙十分发育，特别是红河断裂北段和鹤庆—洱源断裂具有活动性强、切割深

度大的特点，有利于地下水的深循环和深部热物质的对流（颜鹍等，1997;王云等，2019）。

大气降水经洱源盆地周边山区入渗补给，沿构造裂隙向地下运移，经过深循环或深部热物质对

流后，在水头压力差或阻水断裂的阻挡作用下，沿压力较小的断裂破碎带向地表运移并与浅层

冷水发生不同程度的混合，进而表现出不同的水-岩平衡狀态。

Molchan图表法通过对预测结果与观测目标地震的差异度进行检验，直观反映观测数据的

预报效能，得到最佳阈值所对应的异常识别指标（蒋长胜等，2011;孙小龙等，2018）。每一

个阈值对应一组漏报率r和异常时空占有率v，r-v值连线与横、纵坐标围成的面积代表检验的

预测效果，且与面积大小呈反比。此外，报准率h所对应的显著性水平α越低、概率增益Gain

越大时，预测效果越好。

本文利用Molchan图表法对洱源水氡预处理后的斜率数据（图3c）与周边地震关系进行

检验，得到该测项的整体预测效能、最佳预测指标（即异常的最佳判别指标）及优势对应时

间。图3c中空心五角星为选取的21个地震，其中红色代表有效预测期内报准的地震，蓝色代

表有效预测期外漏报的地震，红色虚线为最佳阈值线（0.0175）;图4a为洱源水氡Y值整体预

报效能与概率增益关系图;图4b为洱

源水氡Y值漏报率v与时间占有率r关系曲线及不同显著性水平α等值线。图中红色粗阶

梯线为预测阈值从大到小滑动得到的r-v曲线，蓝色圆圈代表最接近原点（0，0）的r-v值，

其阈值可视为最佳阈值。图中右下角（1，0）点表示地震全部报准，但其时间占有率也最大，

覆盖数据整个时间段;左上角（0，1）点表示地震全部漏报，其时间占有率最小，相当于未做

出预测。

由图4a可见，洱源水氡在不同预测时长的概率增益均大于1，预测效能基本在0.6以上，

具有显著的统计和检验意义。尤其是短期（90d内），概率增益均大于2，预测效能均在0.7

以上，表明洱源水氡短期预测效果更好。因此选取预测时长80d进行计算，如图4b所示，在

80d所示位置（蓝色圆圈），21个地震中，报准15个、漏报6个，报准率约71%，时间占有

率约0.33，置信水平较高，说明预测效果较好。因此，洱源水氡的优势对应地震时间段为90d

以内，在预测时长为80d时，最佳阈值为0.0175，可作为相应时间段的“警报”阈值。

3水化学特征及温泉来源分析

水化学特征分析能够探明地下流体的来源与补给过程，是區别地下流体是否受构造活动影

响的重要技术方法之一（蒋凤亮等，1989）。地壳深部构造活动或浅层地表水补给会影响地下

流体的流量、温度和水化学组分的变化。为了判断氡的来源及其与深部构造活动的关系，笔者

于2018年5月8日对洱源水氡（SD）、洱源老街（LJ）和温泉疗养院（LYY）3处温泉点，

以及滇20井（D20）和附近河流（HL）等浅表水，共计5个水体进行了水样采集，采样点空

间分布如图1所示。样品检测由中国地震局地壳动力学实验室完成，使用Dionex ICS -2100离

子色谱仪进行离子组分检测，包括Ca2+、Mg2+、Na+、K+、SO、Cl-和NO3，浓度检测限为

0.01mg/L;HCO3采用酸碱中和滴定法求得;氢氧同位素δ0与8D使用LGR LWIA -24-EP激光液

态水稳定同位素分析仪测定，精度分别高于0.1%和0.4%。水化学测试结果见表2。

由表2及图5可以看出，洱源水氡和洱源老街优势阳离子成分为Na+、优势阴离子成分为

HCO3和SO4，两处水体的Na+、HCO3和SO4毫克当量百分比分别为92.5%、67.7%、29.4%

以及90.5%、67.2%、29.8%，水化学类型均为HCO3·SO4-Na;温泉疗养院优势阳离子成分为

Ca2+和Na+、优势阴离子成分为HCO3和SO，其各离子毫克当量百分比分别为45.3%、

40.8%、57.5%和41.3%，水化学类型为HCO3·SO4-Ca·Na;滇20井优势阳离子成分为Na+、优

势阴离子成分为HCO3，其离子毫克当量百分比分别为85.9%和93.8%，水化学类型为HCO3-

Na;河流优势阳离子成分为Ca2+和Mg2+、优势阴离子成分为HCO3，其离子毫克当量百分比

分别为60.6%、30.5%和89.8%，水化学类型为HCO3-Ca·Mg。由上可见，温泉水的大多数离

子含量明显高于浅表水，表明其来源较深、离子反应更强烈，从而导致不同水体的水化学类型

差别较大，呈现出明显分区特征。

从图5可见，洱源水氡与洱源老街在菱形区域中的投影点几乎重合，表明其来源较为一致;

而温泉疗养院、滇20井和河流落在不同的位置，说明其来源和形成过程不同，与洱源水氡与

洱源老街没有直接的水力联系。分析认为，洱源地区温泉水水化学组分的形成过程是地下热水

与围岩发生溶滤作用和离子交换作用的结果。由于该地区发育大量二叠系、三叠系灰岩，三叠

系和白垩系的碎屑岩中含有石膏（张春山等，2003），当地下热水流经灰岩层及膏盐层时与围

岩发生溶滤作用，使得Ca2+、HCO3和SO大量进入水中。随着循环深度和循环时间的增加，

CO2和O2等气体耗尽，溶滤作用减弱，离子交换作用开始增强，Ca2+不断被Na+替换，逐渐

形成HCO3·SO4-Ca·Na或HCO3·SO4-Na型水（王大纯，1986;吴桂武等，2015）。

Na-K-Mg三角图通过水-岩平衡状态，可以判断地下水循环深度及混合程度。水-岩平衡程

度越高，地下水循环越深或循环周期越长（Giggen-bach，1988）。由图6可见，所有采样点均

靠近右下角Mg端元，未达到饱和状态。其中，洱源水氡属于“部分平衡水”，洱源老街接近

“部分平衡水”，两者接近重合，说明其水-岩平衡状态相似，来源于深部的可能性更大。且相

比于后者，洱源水氡的水-岩反应程度更高，其接受大气降水补给或浅部冷水混合更少。结合

表2及水化学结果分析，可以清楚看到温泉疗养院具有深部来源特征（SO含量较高），但温

泉疗养院与浅表水体（滇20井、河流）同属“未成熟水”，均表现出大气降水补给的特征，说

明在循环过程中明显受到了地表水或浅部冷水混入影响（Ca2含量较高）。研究区位于多条断

裂交汇处，构造裂隙十分发育，特别是红河断裂北段和鹤庆—洱源断裂具有活动性强、切割深

度大的特点，有利于地下水的深循环和深部热物质的对流（颜鹍等，1997;王云等，2019）。

大气降水经洱源盆地周边山区入渗补给，沿构造裂隙向地下运移，经过深循环或深部热物质对

流后，在水头压力差或阻水断裂的阻挡作用下，沿压力较小的断裂破碎带向地表运移并与浅层

冷水发生不同程度的混合，进而表现出不同的水-岩平衡状态。