朱鹏涛,沈 宁,李自芮,马禾青

（1.宁夏回族自治区地震局,银川 750001；2.宁夏回族自治区地震局银川基准台,银川 750001）

0 引言

地下水的基本特性是上地壳中广泛存在,具有流动性,对上地壳中发生的各种地壳动力作用的响应具有灵敏性,其动态对地震活动与构造活动的响应十分灵敏,同时受多种自然环境的影响也较大。因此,提高利用地下水动态预测地震的效能,必须不断强化地下水异常调查与落实的工作力度,力争事前能够识别和排除各种干扰异常,并对其前兆异常性质做出肯定的确认[1]。首先,确定地下水的主要干扰因素,一般可分为两类,一类是通过含水层水量变化引起的水位变化（其主要因素是降雨和开采）；另一类是含水层应力应变状态改变引起的水位变化,例如气压、固体潮、地面荷载、地震波等附加应力对含水层的影响[2]。井水位观测一直被看作是用于地震预测、预报的有效手段之一[3],在排除干扰异常之后,研究震例是确定地震前兆异常的重要手段,通过大量的震例表明在许多地震前可以观测到井水位的异常变化。因此,本文仔细分析研究西吉王民井水位的变化影响因素,期望研究成果对今后的地震预报有着重要的意义。

1 井孔基础资料

1.1 井周边水文地质条件

西吉王民井观测含水层为第三系向斜承压水,岩性为砂岩。由月亮山区的大气降水和基岩裂隙水补给。地下水补给源远,径流途径长,受地表水干扰小,沟谷地带一般为富水性地带。在观测井东侧有一隐伏断裂从区域穿过,且断错了中新统（N1）泥质砂岩与砂质泥岩和渐新统（E3）的泥岩等。

1.2 井孔基本情况

西吉王民井是宁夏地质局第二水文队于1982年9月16日在水文地质普查工作中建成的普查井。井层底深357.62 m,现有井深351.20 m,水位埋深10.16 m,套管直径127 mm,套管深 度326.44 m,水 温12.0 ℃。1984年 成 为 地震水位观测专用井,架设SW40-1型水位自记仪观测,2007年6月对该观测仪器进行数字化改造,使用LN-3A水位仪和SZW-1A水温仪进行观测,配有气温、气压和降雨等辅助观测项目。水位探头投放深度为16.00 m；水温探头投放深度为180 m。同时安装了公用数字化通信设备、供电、避雷等系统。水位观测段位于井下225.50～324.80 m间,过滤管位于井下225.50～324.80 m间。周围10 km无其它水井抽水干扰,观测条件比较理想。

2 水位影响因素分析

2.1 水位对气压的响应

2.1.1 频谱分析

将时域信号变换至频域加以分析的方法称为频谱分析。频谱分析的目的是把复杂的时间历程波形,经过傅里叶变换分解为若干单一的谐波分量来研究,以获得信号的频率结构以及各谐波和相位信息。选取2008年数字化正式观测以来,王民井2个不同时段（2012年1月至3月、2017年1月至3月）的水位整点值数据,得到该井水位2时段的频谱结果（图1）,通过频谱分析发现,引起王民水位变化的主要潮汐分波有M2波、S2波、O1波、K1波。其中,S2波和K1波由气压引起。

2.1.2 回归分析

（1） 线性回归分析

线性回归是利用称为线性回归方程的最小平方函数对一个或多个自变量和因变量之间关系进行建模的一种回归分析。选枯水期两个月（2018年11月—12月）水位整点值作为因变量与气压做线性回归计算,从置信水平95%条件下的回归结果看,水位和气压存在明显的负相关,并得到相关系数R=-0.813（图2）,从而得到气压对水位有一定的影响。

图1 王民井2时段内,原始水位的主要潮汐分波振幅谱Fig.1 The main tidal partial wave amplitude spectrum of the original water level in the 2nd period of Wangmin well

图2 2018年11—12月剔除气压（a）、拟合结果（b）及王民水位（c）对比曲线Fig.2 Contrast curves of rejecting air pressure (a), fitting result (b) and Wangmin water level (c) from November to December in 2018

（2）卷积回归分析

考虑水位与气压间的非线性关系及滞后时间的卷积回归法,已成为近年来气压校正的主要方法之一[4-11]。该方法利用相对应的水位、气压数据来拟合阶跃响应函数的最佳值,再由该最佳阶跃响应函数对水位进行校正,最终计算出校正后的水位。

井水位对气压的阶跃响应函数可表示为：

式中,A(i)为井水位对气压的阶跃响应函数,可由气压单位脉冲响应函数α(i)累加求和得到。

校正后的水位可表示为：

选取2018年12月的水位整点值,通过两种方法（线性回归和卷积回归）进行气压校正对比,两种方法剔除气压之后的固体潮清晰,曲线变化形态一致（图3）,再次证实西吉王民井水位受气压影响。

图3 西吉王民两种方法剔除气压之后的水位对比图Fig.3 Comparison of the water level of Wangmin wella fter removing the air pressure by the two methods(dotted line is convolution regression, solid line is linear regression)

2.2 水位对降雨的响应

西吉王民井水位数字化以来由于更换芯片,电瓶试验以及更换仪器产生多个台阶,选取2013年以来变化相对平稳的水位整点值进行分析。由于该井点的降雨辅助观测资料不全,故向宁夏气象局申请西吉县降雨量月值（表1）,分析西吉王民水位与降雨量之间的相关性。考虑到其他因素的影响,将原始水位进行去气压、固体潮及趋势,用剩余水位与降雨量做对比分析,结果显示该井水位与区域降雨有关,降雨量增大时水位也上升（图4）。进一步利用2013年以来西吉王民水位月均值与降雨量月值做线性回归分析,得到相关系数（见表2）,从表2可以得到在2013年、2017年和2018年雨季降雨量比往年同期降雨量偏高的情况下,相关系数也相对较高。证实西吉王民水位的变化受降雨影响。

表1 西吉县历年降雨月值Table 1 Monthly rainfall value of Xiji County (according to Ningxia Meteorological Bureau)

表2 西吉王民水位月均值与降雨量月值之间相关性统计Table 2 Correlation statistics between the monthly mean value of the water level of Wangmin well and the monthly value of rainfall

图4 2013年以来王民水位（a）、降雨量（b）对比曲线Fig.4 Comparison of Wangmin water level (a) and rainfall (b) since 2013

2.3 李马河水库对井水位的可能影响分析

2.3.1 水库基本情况

李马河水库于2015年下半年开建,2016年下半年完工,位于西吉王民测点上游1.2 km处（东经106°31′39.9”,北纬35°55′03.6”）,距西吉县城50 km,距王民乡东4.5 km。建设过程中开始蓄水,总库容628万m³,目前蓄水量30万m³,主要用于农田灌溉。集水面积54 km2,河道总长60.6 km,河道平均比降30.7‰,流域跨越甘肃、宁夏两个省区,涉及甘肃省静宁县的原安乡和宁夏西吉县的王民、平峰乡,共计8个行政村（原安乡的陈义村、关音村,王民乡的王民村、三岔村、学杨村、平峰乡的金塘村、西坡村、张武村）。库区属黄土丘陵地貌单元。由于地壳运动及河流长期切割作用,使河谷形成“U”字型,两岸并形成相称的Ⅰ、Ⅱ级阶地,并有数条冲沟发育,岸坡较陡。第三系基岩强风化厚2 m,中等风化3 m,库区地形简单,临水面多形成陡坎。库区地处黄土丘陵地貌单元,第四系堆积物较厚,两岸山体高耸且宽厚,壤土属弱透水层。

2.3.2 区域地形及高程变化情况

王民井处于典型的黄土沟壑区域,冲积沟发育,井水位变化受季节降雨影响较大。在观测井SW方向的冲积沟内有2016年上半年完工的李马河水库（图5中的AA剖面）,主要用于农田灌溉。值得注意的是,李马河水库所在的冲积沟（AA剖面）越接近观测井,其海拔高度越低（也即王民井处于该冲积沟的下游）,更有利于降雨及地表水的汇集。位于观测井北侧稍远的BB剖面,其海拔高度变化不明显,且整体要比西南侧的海拔高度低（图6）。

图5 王民井周边地形及AA和BB剖面示意图Fig.5 Schematic diagram of the surrounding terrain and AA and BB sections of Wangmin well

图6 AA和BB剖面海拔高度变化情况Fig.6 Altitude changes of AA and BB profiles

结合区域水文地质条件,西吉王民井观测含水层为第三系向斜承压水,岩性为砂岩,由月亮山区的大气降水和基岩裂隙水补给。该井上游有利的地形条件和李马河水库的蓄水对井水位动态变化可能产生很大影响。也即王民井水位自2017年年初出现的下降速率变缓现象与李马河水库的蓄水可能有很大关系。

3 震例分析

西吉王民井自观测以来,水位对发生在周边300 km范围内的地震有一定映震能力。1990年10月20日和2000年6月6日先后在甘肃景泰地区发生的MS5.9级和MS6.2级地震前,均出现了异常变化。两次地震前水位均出现明显的上升,在上升速率减缓转为下降的过程中发生地震,且异常持续时间较长,而此次该井水位自2017年年初至今出现下降速率减缓,并未出现明显的上升变化,与之前异常形态不一致（图7、表3）,是地震前兆异常的可能性不大。

图7 1985—2018年西吉王民水位月均值曲线Fig.7 Curve of monthly average water level of Wangmin well in Xiji from 1985 to 2018

表3 西吉王民井水位历史震例异常变化总结Table 3 Summary of abnormal changes in the historical earthquake case of Wangmin well

4 结语

西吉王民井水位2013年以来处于下降变化,2017年年初至今该井水位出现下降速率减缓的现象。本文利用频谱分析得到引起王民水位变化的潮汐分波主要有M2波、S2波、O1波、K1波。其中,S2波和K1波由气压引起。采用两种方法（卷积回归法和线性回归法）对水位进行气压校正,校正后的曲线变化形态一致。对水位和降雨量做线性回归分析,从结果可以得到,降雨量增多时相关系数增大。并且经过走访调查发现水位的变化与新建李马河水库蓄水有很大关系,李马河水库建成后区域地下水可能达到一种新的动态平衡,但由于水库建成时间不长,蓄水资料积累较少,后续需密切关注李马河水库蓄水以及区域降雨对该井水位动态变化的影响。通过对西吉王民井水位的变化影响因素分析,得出以下结论：

（1）气压和降雨是西吉王民井水位动态变化的主要影响因素。

（2）西吉王民井水位2017年年初至今出现下降速率减缓的现象可能是李马河水库建成后区域地下水达到了一种新的动态平衡。