在地壳动力作用信息(又称为地下水微动态)的观测与研究中,观测井地下水的承压性是至关重要的条件,因此观测井建设的相关规范与标准都明确要求地震地下水观测井必须是承压井。然而,由于我国地震地下水观测井多为由地质、石油、水利、矿产等行业提供的勘探井经筛选而成,且在成井时也仅是定性的从地下水种类、成因、补给区与分布区的关系和动态特征等方面进行区分,观测井地下水类型主要有潜水、承压水以及二者的混合型(混合水或半承压水); 另一方面,由于这些井已有四五十年的历史,井孔条件发生了很大变化,如井孔坍塌、井管变形与破裂、滤水管渗水孔被堵塞、井底泥沙沉积等,井孔原始的承压性已经改变,这种状况不仅影响了井水位对地壳动力作用的响应能力,还严重影响观测数据分析与微动态特征的识别,特别是地震前兆信息的识别与提取,因此,对观测井地下水埋藏类型的研究与评价迫在眉睫。

对于严格按照相关规范与标准新建的观测井,其地下水承压性的判定与研究在钻井施工时即可进行,在成井报告中也会给出明确的说明,但对于观测了几十年的观测井而言,这就成了值得研究的问题。本文提出了基于井水位微动态数据分析及井水化学特性分析来判定观测井现今承压性的4种方法,详述其判定原理与具体操作过程,以期对地下水观测井映震能力的评估、观测井质量的评价与观测条件的改造与改善、提升地震地下水监测的有效性和科学性、认识和理解地下水运移和动态变化提供一定地理论基础和技术支撑。

进行地下水埋藏类型综合评价主要通过以下步骤完成(图1):①采集观测井至少48 h的水位、整理气压和理论固体潮的整点值数据。②利用卷积回归法中水位对气压的滞后时间与阶跃响应函数的关系式分析观测井的地下水埋藏类型。③将水位整点值数据进行傅立叶变换得到水位的频谱分析结果,利用各潮汐波主要分量的频率特性分析观测井的地下水埋藏类型。④利用维尼迪柯夫调和分析方法得到相位滞后参数,利用潮汐波群相位超前或滞后方法分析观测井的地下水埋藏类型。⑤采集观测井的水样(Na⁺,K⁺,Mg²⁺,Ca²⁺,Cl^-,SO^2-₄,CO^2-₃,HCO^-₃),进行水化学成分分析,通过Na-K-Mg三角图,依据未成熟水、部分平衡水或完全平衡水3种类型,分析待测观测井的水-岩平衡状态。

图1 地下水埋藏类型综合评价示意图
Fig.1 Comprehensive evaluation of the type of the embeded groundwater

将观测井的水位和气压数据整理成数据文本文件,按年逐月检查、整理,要求每天必须有24个数据,缺数的地方按照要求用“99999”进行标记,并改变为统一格式备用。水位单位为m,气压单位为hPa。

结合3次样条插值和一般多项式分段拟合值对数据进行替换。并对数据进行粗差检查,目的是检查整个数据段内的基线值是否统一,数据是否有明显的错误甚至非法。

准备一个各井点的excel文件,包括:日期、水位、气压和理论固体潮,且是等间隔的整点值。水位须由埋深值换算成水头高度值,单位为m; 气压单位为hPa; 理论固体潮单位为nm/s²。

井水位对气压的阶跃响应函数方法能区分2种类型的地下水,即承压水和非承压水(或半承压水); 频谱分析方法可区分3种类型,即承压水、半承压水和非承压水; 潮汐波群相位超前或滞后方法也只能区分两种类型,即承压水和非承压水或半承压水; 水-岩平衡状态分析方法能区分3种类型,即未成熟水、部分平衡水和完全平衡水。因此,可通过上述4种方法的组合来综合判定观测井的地下水埋藏类型,其依据是:①若使用3种定量方法——井水位对气压的阶跃响应函数、井水位潮汐波的频谱分析、井水位潮汐波群的相位超前或滞后都判定为承压水,且不管水-岩平衡状态方法判定结果如何,则综合判定为承压水。②当使用水-岩平衡状态方法判定为完全平衡水时,直接判定为承压水。③当使用3种定量方法——井水位对气压的阶跃响应函数、井水位潮汐波的频谱分析、井水位潮汐波群的相位超前或滞后,其中有1～2种方法都判定为承压水,且水-岩平衡状态方法判定结果为未成熟水或部分平衡水,则综合判定为混合水。④当使用井水位对气压的阶跃响应函数方法和井水位潮汐波群的相位超前或滞后方法判定为非承压水或半承压水、使用频谱分析方法判定为非承压水,且水-岩平衡状态方法表明是部分平衡水,则该井综合判定为混合水。⑤当使用井水位对气压的阶跃响应函数方法和井水位潮汐波群的相位超前或滞后方法判定为非承压水或半承压水,同时,使用频谱分析方法判定为非承压水,且水-岩平衡状态方法表明是未成熟水,则综合判定为非承压水。

以上可见,除了上述①、②和⑤条外,其余都可将观测井水综合判定为混合水。

以川滇地区的泸州13井、南溪井、川32井、红河开远井、腾冲井、丽江井、通海高大井、德阳井、姚安井、永胜井、保山井和德宏法帕井等12口井为例,利用上述4种方法对地下水埋藏类型进行了综合判定(表2),结果表明:①通海高大井、永胜井、姚安井基于井水位对气压的阶跃响应函数方法、频谱分析方法、潮汐波群相位超前或滞后方法诊断和判别的地下水埋藏类型都为承压水,且与3口井依据成井时的贮存埋藏条件定性判别的结果(裂隙承压水和孔隙承压水)一致。从水-岩平衡状态来看,虽然3口井都具有承压性,但受控于区域地形地貌、构造特征以及裂隙发育程度等,使得井孔中的地下水与外界有着一定的水力联系,受大气降水的地表径流、侧向径流和越层补给影响较大。②德阳川08井基于3种方法虽然没有诊断和判别为承压水,但水-岩平衡状态显示,该井水-岩反应较为充分,补给来源主要为深层地下水。③其余8口井不同程度地表现出承压水的特征,但与外界水力联系关系密切,受大气降水的地表径流、侧向径流和越层补给影响较显著。因此,综合判定通海高大井、永胜井、姚安井和德阳川08井为承压水; 腾冲井、德宏法帕井、川32井、开远井、保山井、丽江井、泸州川13井、南溪川12井为混合水。

表2 多方法综合判定川滇地区12口井地下水埋藏类型结果
Tab.2 Results from comprehensive determination of groundwater types of 12 wells in Sichuan-Yunnan region by multiple methods

运用本文提出的多种方法综合判定的地震观测井的地下水埋藏类型,与依据成井时的贮存埋藏条件定性判别的地下水埋藏类型有很多不一致。究其原因是:①有些井点虽然在承压区(具有承压性)内,但由于比较靠近补给区或排泄区,受控于区域地形地貌、构造特征以及裂隙发育程度等,使得井孔中的地下水与外界有着一定的水力联系。②有些井点虽然在承压区内,但井管年久失修,存在腐蚀破损现象,地表水和浅层水很容易进入井管内,与外界的水力联系也很密切。③远离补给区或排泄区,且井况良好的井点受地表水和浅层水的干扰较少,水-岩反应相对较弱和充分,补给来源主要为深层地下水,Na-K-Mg三角图反映的水-岩平衡状态能提供很好的佐证。

地震观测井水位观测作为地震流体观测最有效的手段之一,过去的几十年里在地震监测预报方面发挥了重要的作用。但受观测环境干扰日益加剧、井点年久失修等因素的综合影响,其产出资料的质量和信度不尽如人意。本文基于井水位微动态数据及井水化学特性分析,参考井水位对气压的阶跃响应函数方法、频谱分析方法、潮汐波群相位超前或滞后方法、水-岩平衡状态分析方法,提出了地震观测井地下水埋藏类型判定方法及综合判定依据,并以川滇地区的12口地震观测井为例,分析了多方法综合判定结果与依据成井时的贮存埋藏条件定性判别的地下水埋藏类型结果之间有很多不一致的原因。

因此,应用本文提出的4种方法对地震观测井地下水埋藏类型判定进行研究,是井孔观测质量和可信度评价、水位异常的甄别、观测井(点)的勘选等的一种有效途径。