云南是一个多震区域,进行地下水动态变化观测的台站也较多。大多数台站的地下水动态观测有可能会监测到地震前兆异常信息,但常有较多的环境干扰因素混淆其中。因此,如何及时排除干扰识别地震异常信息,对于准确判断把握地震前兆异常信息是至关重要的(孙小龙等,2016)。一般在识别地下水动态的干扰分析中,主要把握“四个相关性”,即成因、空间、时间及强度的相关性(车用太等,2011)。通过这些相关性的分析,可以定性定量地认识地壳浅层观测井水位、水温等物理观测量受环境等外界因素的影响。

近年来,许多学者结合水物理和水化学方法,广泛开展了地下水异常等分析研究工作(张磊等,2014; 张涛等,2016; 方震等,2016; 胡小静等,2018),取得较多研究成果。本文利用云南腾冲地震台的地下水动态观测资料,应用水物理方法和水化学方法相结合的实验分析方法,对多年的观测资料变化进行了分析研究,探讨其异常变化的原因。

腾冲地震台(以下简称“腾冲台”)流体观测井位于腾冲县城东南郊腾越镇董库村地震台内,距离城区约3 km。井孔位于NE向大盈江断裂与SN向怒江断裂交汇处(图1a),处于腾冲盆地东侧近山麓地带,盆地属熔岩盆地堆积地貌类型,东山区为近代火山岩分布区,属火山穹丘地貌类型,火山呈近南北向线状排列。该井成井于1999年,井深120 m,非自流井,属静水位观测。井壁套管下至井下60 m,直径为160 mm。包网滤水管在井下60～117 m,直径为108 mm,水位埋深为5.8 m左右,水的矿化度为255.44 mg/L,观测水层为第四系英安岩、安山岩及熔结凝灰岩沙层,地下水类型属裂隙水(图1b),水化学类型为HCO^-₃-Ca²⁺,Mg²⁺型。

距观测井的直线距离约100 m处(图2)有天然出露的2个低温温泉——大、小闫家塘,表面积分别为7 990 m²和5 099 m²,储水量分别为8 500 m³和 6 200 m³。每年1—4月,大、小闫家塘均会由村民进行放水捕鱼、蓄水供游泳等。水位观测日志记录显示闫家塘的放水、蓄水均引起了腾冲流体观测井井水位的大幅度下降、上升变化。

图1 腾冲台井构造位置分布(a)及井孔柱状图(b)
Fig.1 Tectonic map(a)and columnar section of well(b)at Tengchong seismic station

图2 腾冲台流体观测井与闫家塘相对位置示意图
Fig.2 Map showing the Tengchong seismic station fluid observation wells with Yanjiatang relative position

一般情况下腾冲台井水位能记录到清晰的固体潮变化,日变形态为双峰单谷型,潮差不超过0.02 m。年变规律主要受降雨干扰,呈夏秋高冬春低的变化形态,动态基本稳定。每年1—4月井水位都会出现一次以上超过正常日变幅度10倍的大幅度快速下降又快速上升的变化过程(表1)。尤其2018年几乎每个月井水位都出现大幅快速升降异常变化(图3)。

腾冲台井水位大幅快速升降异常变化有以下特点:水位下降历时3～6天,水位上升历时5～15天,水位下降速率高于上升速率。 下降幅度在0.246～0.775 m内波动变化,最大降幅为2015年的0.755 m,最小降幅为2017年的0.246 m。上升幅度在0.256～0.693 m波动变化,最大升幅为2011年的0.693 m,最小升幅为2008年的0.256 m。下降幅度和上升幅度均达到或超过正常水位日变动态幅度的10～20倍。

当水位处于趋势性下降时段时,水位下降幅度略大于上升幅度。腾冲台井水位在时间尺度上的下降与上升的动态变化过程大致存在“快速下降—平稳—快速上升”和“快速下降—快速上升”2种类型。

影响地下水动态的因素主要包括大气降水、水文、地震活动和人为活动(抽水与注水)等。其中,大气降水和地下水开采对井水位的影响最为普遍(车用太等,2011)。由图3a可知,腾冲台井水位大幅快速升降异常变化大都出现在每年的1—4月,这段时间属于腾冲的旱季,降雨相对较少,无雨或少量的降雨不足以引起水位产生如此大的变化。查阅腾冲台前兆观测日志和工作日志重大事项记载发现,每次腾冲台水位出现大幅

升降变化时,都对应闫家塘放水、蓄水记录。结合腾冲台井水位变化情况,笔者认为腾冲台流体

表1 腾冲台水位大幅快速升降异常变化记录表
Tab.1 Records of abnormal changes of water level at the Tengchong seismic station in a large and rapid manner

图3 2007—2018年(a)、2018年(b)腾冲台水位与降雨量对比图
Fig.3 Comparison of water level and rainfall at the Tengchong seismic station

观测井距闫家塘较近,井水位每次大幅快速升降异常变化的时间与闫家塘蓄水、放水的时间基本同步,井水位升降变化的幅度与闫家塘的放水量密切相关。因此,初步判断腾冲台井水位大幅快速升降异常变化是由闫家塘放水、蓄水引起的,井水位变化受闫家塘水位升降变化干扰的可能性较大。

由图3可见,腾冲台井水位异常映震的概率较高,在2008年8月21日盈江5.9级地震、2011年3月10日盈江5.8级地震、2011年6月20日及8月9日腾冲5.2级双震和2014年5月30日盈江6.1级地震前都出现了高水位前兆现象,震中距最近为20 km,最远为97 km。腾冲台井水位在近震前出现了前兆异常,表现为异常升高,异常重复性较好。由上文分析可知,腾冲台流体观测井水来自大气降水,水位的升高主要是由大气降水增多引起的。赵洪声(2006)通过对1966年邢台7.2级大震和1970年通海7.8级大震水-震过程的分析以及40多年来云南涝旱剧变与M≥6.2强震的相关研究得出,地表水和地下水的大面积、大幅度涨落,很可能是孕震区应变非线性增长并驱使孕震体进入不稳定状态乃至失稳的重要原因。腾冲过量降雨的加载对该地区构造活动具有调制作用,降雨量超出正常变化水平后,腾冲台井水位会出现高水位异常变化,这种异常变化能够提供一定的地震前兆异常信息,反映台站附近地区的部分地下构造活动。

本文对大、小闫家塘放水、蓄水过程与腾冲台井水位的变化的实验数据进行分析,证明两者有很好的相关性。闫家塘放水水位下降或蓄水水位上升,都对腾冲台井水位的变化起到主要影响作用。

水化学离子浓度分析结果与氢氧同位素的分析结果表明,腾冲台流体观测井井水和小闫家塘水主要来源于大气降水,大闫家塘水来源除大气降水外,还有其它方式。

腾冲台井水位在近震前出现水位异常,表明该井水位变化能够提供一定的地震前兆异常信息,对台站附近地区的地下构造活动有反映。

实验表明大、小闫家塘的蓄放水过程对腾冲台水位井周围的地表也产生了一定的加卸载效应。闫家塘与腾冲井流体观测井直线距离约100 m,在地质构造上具备较复杂的区域地质构造应力关系。大闫家塘储水量约为8 500 m³,小闫家塘储水量约为6 200 m³,按普通水的密度计算,大闫家塘水质量约位8 500 t,小闫家塘水质量约6 200 t,具有加卸载作用。

当闫家塘放水时,腾冲台井水位同步下降。但当放水结束后,闫家塘水体对腾冲井孔周围地表载荷降低,井孔周围的岩体孔隙有变化,井水位有持续下降的滞后过程。当周围环境达到基本平衡时,腾冲台井水位停止变化。当闫家塘蓄水时,产生与放水时相反的变化过程,腾冲台井水位同步上升,并有持续上升滞后并达到平衡的过程。

根据实验数据及变化情况分析,可以得出腾冲台井水位受环境的影响是明显的,产生大幅快速升降变化是闫家塘放水、蓄水过程所产生的。因此,在使用腾冲井水位观测资料时要充分考虑受环境影响的具体原因,客观正确地分析观测资料的变化。

感谢中国地震台网中心田雷,云南省地震局预报研究中心高文斐、张翔在样品采集和测试过程中提供的帮助与技术指导; 感谢审稿专家对本文提出的建设性意见。