第10章 地下水动态与均衡

10.1 地下水动态与均衡的概念

        地下水各种要素（水位、水量、化学组分、气体成分、温度、微生物等）随时间的变化，称为地下水动态（groundwater regime）。

        某一时段、某一范围内地下水量、盐量、热量等的收支状况，称为地下水均衡（groundwater budget)。 地下水动态与均衡的关系是：地下水动态是地下水均衡的外在表现，地下水均衡是地下水动态的内在原因。 地下水动态的研究包括：影响因素、类型及成果分析。

        地下水动态监测及成果分析，可以解决一系列理论与实际问题：

        ①检验并完善前期水文地质研究结论；

        ②监测地下水量、水位、水质等变化，为保护地下水及相关生态环境进行预警；

        ③为数学模拟提供依据；

        ④为拟定合理的地下水利用/防治方案及措施提供依据；

        ⑤检验实施中的利用/防治措施的合理性。

        地下水均衡的研究包括：均衡区和均衡期的确定，均衡方程式的确定，各收支项的求取，均衡计算结果的校核与分析。 地下水均衡研究，可以为拟定合理的地下水利用/防治方案及措施提供定量依据，检验并完善利用/防治方案及措施。 2015年7月，水利部水文局和中国地质环境监测院联合启动“国家地下水监测工程”三年计划，涉及我国七大流域和16个含水系统的重点地下水监测，监控面积约350×10* km2，计划建成国家级（自动）监测站点共计20401个（中国地质环境监测院，2017）。此外，还有各部门省市区根据各自需要布设的地下水监测站点。

10.2 地下水动态的影响因素

10.2.1 地下水动态的两类影响因素

        地下水动态是地下含水系统对外界激励（输入）产生的响应（输出）。 影响地下水动态的因素分为两类：一类是地下水诸要素（水量、盐量、热量、能量等）本身的收支变化，即外界激励（输入）因素；另一类是影响输入-输出关系转换的转换因素，即影响地下水动态曲线具体形态的因素。

        地下水动态是两类因素综合影响的结果。 以大气降水入渗补给抬升潜水位为例加以说明。一个降雨-地下水位抬升过程，可以看作一个脉冲转换为波形的过程。包气带的滤波作用将一次降雨脉冲转换为一个时间滞后和延迟（延续）的地下水位波峰。波峰与降雨相对应，波峰出现和延续的时间，以及波峰形态，取决于包气带岩性及地下水埋藏深度。

         图10.1表示包气带岩性和厚度不同的条件下，地下水位对一次降雨的响应。

        图10.1中波形变化的水位分别代表不同类型包气带：

        1为包气带渗透性良好（此时包气带厚度影响可以忽略）；

        2为包气带渗透性及厚度均为中等；

        3为包气带渗透性低且厚度大。

        三种条件下，地下水位抬升对降雨的响应分别为：

        1为时间滞后及时间延迟都短的尖峰；

        2为时间滞后及时间延迟都中等的波峰；

        3为时间滞后及时间延迟都很大的缓峰。

连续若干次降雨，可形成叠合波峰。

图10.1显示了不同条件下地下水位对降雨的响应，其中a1、a2、a3代表地下水位波峰对降雨响应的时间滞后；b1、b2、b3代表地下水位波峰对降雨响应的时间延迟。1代表岩溶发育区，2代表地下水埋藏深度中等的裂隙砂岩，3代表地下水埋藏深度大于100 m的黄土高原。

         由此可见，通过地下水动态分析，不仅能够获取地下含水系统收支的信息，还能获取地下含水系统结构的信息。

        通过收支影响地下水动态的因素随时间变动，包括：气象（气候）因素、水文因素、生物因素、地质营力因素、天文因素等。

        地下水动态的转换因素主要是地质结构及水文地质条件，如地质构造、含水层类型、岩性等。转换因素只在地质时间尺度内变化，对地下水动态而言，是不随时间变动的因素。

10.2.2气象(气候)因素影响下的地下水动态

        气象（气候）因素对地下水动态影响最为普遍。降水的数量及其时间分布决定地下水位与水量的时间变化。在干旱、半干旱平原和盆地中，地下水以蒸发和排泄为主，随着地下水位和水量变化，水质也随时间有规律地变化。 在气象（气候）的影响下，地下水动态呈现昼夜变化、季节变化、年际变化及多年变化。 地下水位的昼夜变化，可以由蒸发和植物蒸腾引起。白天植物生长旺盛，蒸腾强烈，地下水位下降；夜晚蒸腾停止，得到来自周围水的补充，地下水位上升，由此引起的昼夜变幅可达数厘米（Todd et al., 2005）。 我国大多数地区为季风气候，旱季和雨季分明，地下水位呈现明显季节变化。夏季多雨，地下水位抬升达到最高，雨季过后，地下水通过蒸散和（或）径流排泄，次年雨季以前，地下水位最低，全年地下水位呈单峰单谷形态（图10.2）。

        图10.2显示了季风气候影响下的北京市潜水位动态曲线。 气候的周期变化控制地下水动态的多年变化，其中周期约为11年的太阳黑子变化影响最为明显。太阳黑子平静期，降水丰盛，地下水位高，地下水储存量增加；太阳黑子活动期，降水稀少，地下水位低，地下水储存量减少（图10.3）。

图10.3显示了苏联卡明草原地下水位多年变化曲线，根据每年9月1日平水位资料编绘，实点代表实测水位，空心点代表可能水位。 

        大型地下水供水工程的设计，需要考虑干旱周期地下水供水能力能否满足需求。相反，大型地下水排水工程的设计则需要满足湿润周期的排水要求。如果长期地下水动态的观测资料不足，可以利用多年的气象、水文资料以及历史资料进行推求。

        大气压强可以通过井孔影响周边小范围地下水位。大气压强变大，井孔水位降低，大气压强变低，井孔水位抬升。由于大气压强变化引起的潜水井孔水位变化非常小，通常在1厘米左右；相反，大气压强变化引起的承压水井孔水位变化大，可超过10厘米（陈葆仁等，1988年）。因此，大气压强引起的井孔水位变化不能代表地下水位及水量的实际变化。 

10.2.3 水文因素影响下的地下水动态

        当河流补给地下水时，随着远离河流，地下水位抬升的时间滞后和延迟增大，波形趋于平缓（图10.4）。河岸及含水层的渗透性越强，地下水位响应的时间滞后和延迟越小；含水层给水度越大，则波形越平缓。通常情况下，河流对地下水质和温度的影响范围小于地下水位波动范围。

图10.4显示了莱茵河洪水对潜水位的影响，摘自卡明斯基 (1958)。1~5表示观测井潜水位，数字越大则距离河流越远，6表示莱茵河水位。

10.2.4 其他因素影响下的地下水动态

        地震、固体潮、潮汐、外部荷载等都可以引起地下水要素变化。地震孕震及发震阶段的地应力变化会引起地下水位、化学成分、气体成分等的变化。1975年2月4日，辽宁海城发生7.3级大地震，根据包括地下水动态在内的各种前兆，成功进行了地震预报，避免了大量伤亡。震前，附近地区地下水位上升，之后下降，震后大幅度回升（图10.5；王永林，1987；杨成双，1982）。1974年1月至地震发生时，地下水出现氡 (Rn) 异常（图10.6；张立海等，2007）。

        需要注意的是，地下水动态异常不一定是地震前兆。迄今为止，准确鉴别地下水异常是否是地震前兆仍然是一个有待解决的课题（车用太等，2005）。

        内陆地区的承压含水层可以观测到周期为12小时的测压水位波动。这是由于月亮和太阳对地球表面的引力相互作用而产生的固体潮引起的。当月亮运行到某地中天时，承压含水层荷载减少，轻微膨胀，测压水位出现厘米级波动。满月时，月亮达到中天位置，月亮和太阳对地球表面的引力合力最大，测压水位降低幅度最大（陈葆仁等，1988）。

图10.5显示了根据杨成双（1982）的资料，1975年海城7.3级地震前后丹东文斌井地下水位变化的情况。

图10.6显示了根据张立海等（2007）的资料，1975年海城7.3级地震前后地下水氡 (Rn) 异常变化的情况。

        海洋潮汐会增减承压含水层的荷载，从而使地下水位发生相应的升降。例如，湛江市海潮引起的地下水位升降幅度为0.1~1米，最大可达2.5米（章至洁等，1995年）。同理，火车停车及开动，会使附近承压含水层测压水位出现厘米级升降（陈葆仁等，1988年）。

        在上述各种情况下，外力通过改变含水层应力状态而引起地下水位变化，并不伴随地下水储存量的变化。

10.2.5人为活动影响下的地下水动态

        在天然条件下，气候因素在多年中趋于某种平均状态。因此，地下水的补给量、排泄量和储存量多年中保持平衡状态，地下水位围绕某一平均水位波动，水质稳定地趋向淡化或盐化。

        人类活动增加了新的补给来源或排泄去路，影响了地下水的天然均衡状态，从而改变了地下水动态，并影响了水质演变方向。

        井孔开采地下水是最常见的人工排泄。新增的人工排泄将减少甚至完全替代原有的天然排泄（如泉流量减少或枯竭、向河流泄流较少或停止、蒸发量减少）。有时候，还伴随着某些补给的增加（如地下水由补给河水而转变为接受河水补给，原先地下水埋藏深度过浅降水入渗受到限制的地段，因水位下降而入渗补给量增加）。

        如果新增的补给量和减少的天然排泄量之和等于人工排泄量，地下水将会达到新的均衡状态。地下水位将维持在较原先平均高程更低的位置，并以更大的幅度变动，但不会持续下降（图10.7）。

图10.7显示了据河北省第九地质大队观测资料编制的河北饶阳五公地下水位变化曲线。

        如果人工开采水量过大，新增的补给量及减少的天然排泄量之和不足以补偿人工排泄量，地下水位将持续下降（图10.8）。

图10.8 河北保定西部地下水位变化曲线(据河北省地矿局第四水文地质大队，1976)
Fig. 10.8 Changes of groundwater level，Baoding,Hebei
 

        修建水库、引用外来地表水灌溉等都会增加新的人工补给，抬高地下水位。例如，1974年初，河北冀州新庄的潜水埋藏深度大于4米，引用外来地表水灌溉后，到1977年雨季，潜水位接近地表，发生次生沼泽化，导致农业大幅度减产（图10.9）。

        地下水的人工补给或人工排泄都有可能打破原有的地下水均衡，形成新的地下水均衡，进而影响地下水所支撑的水文系统及生态环境系统。地下水人工排泄减少蒸发量，在干旱半干旱地区，可以减轻乃至消除原有的土壤盐渍化。地下水人工排泄大幅度降低地下水位，减少向河流泄流，则将损害水文系统以及相关生态环境系统的正常运行。地下水人工补给过多，促使地下水位抬升到离地面很近的位置，则将引起次生沼泽化。在干旱半干旱地区，还会导致土壤盐渍化。

        图10.9河北冀州新庄地下水位变化曲线
 

10.2.6影响地下水动态摇的转换因素

        含水系统的结构对输入信号进行滤波或增强，然后输出为观测得到的地下水动态。地质结构以及非变动性水文地质条件（如岩层渗透性、孔隙率、裂隙率等）是影响地下水动态的转换因素。

        大气降水入渗补给抬升潜水位时，包气带岩性及厚度对降水脉冲起滤波作用。饱水带岩性也会影响潜水位变幅大小。潜水储存量的变化，以给水度（饱和差）u与水位变幅Ah的乘积uAh表示。当入渗补给量相同时，给水度（饱和差）u愈大，潜水位抬升值△h便愈小。承压含水层获得补充水量或能量时，储存量的变化以弹性给水度ue(贮水系数S）与测压水位变幅△h的乘积uAh表示。由于弹性给水度比给水度小1~3个数量级，接受同量补给或增加同等应力时，承压水测压水位抬升幅度比潜水位大得多。

        潜水含水层水位的变化通过质量传输完成。承压含水层中测压水位的变化，则是压力传递的结果。压力传递速度远大于质量传输。例如，河水补给承压含水层时，测压水位的变化，滞后时间短，波及距离大。

        承压含水层的隔水顶板限制了承压水和大气及地表水的联系，只能在有限的范围接受补给。因此，承压水位动态变化通常小于潜水。构造越封闭，承压水的动态变化越不明显。

        地下水流系统的不同部位，地下水位的波动幅度不同。在区域系统的补给区，地下水位变幅最大，排泄区变幅最小。在局部系统的补给区，地下水位变幅较大，排泄区变幅较小。原因在于：排泄区附近获得补给时，受排泄区高程限制，水力梯度显著增大，径流明显加强，地下水位不可能明显抬升；补给区接受降水补给时，因远离排泄区，水力梯度无明显增加，径流增强也不大，水位得以累积抬升。随后，由排泄区向补给区，水力梯度溯源增大，补给区径流加强，水位逐渐下降。

        我国南方岩溶水区域水流系统的补给区，地下水位对降雨响应迅速，并且变动幅度很大，可以达到数十米。这是多种因素综合影响的结果。岩溶水的含水介质，具有空隙尺寸大、空隙率小的特点。空隙尺寸大，渗透性良好，有利于降雨大量快速入渗；雨季过后，良好的渗透性使得径流强烈，地下水位迅速降低。空隙率小（相当于给水度小），接受补给时，地下水位抬升幅度大；发生排泄时，地下水位下降幅度大。