1、第六章 地下水动态与均衡的研究1 地下水动态和均衡的概念地下水资源和其它矿产资源的最主要区别是，其量和质总是随着时间而不停地变化着。所谓地下水动态即是指表征地下水数量与质量的各种要素（如水位、泉流量、开采量、溶质成分与含量、温度及其它物理特征等）随时间而变化的规律。其变化规律可以是周期性的变化，也可以是趋势性的变化。变化的周期可以是昼夜的（如月球引力导致的固体潮），也可以是季节性的或者是多年的。其变化的速率，在天然状态下一般具较明显的周期性，或具极为缓慢的趋势性。在人为因素（开采或排除）的影响下，其变化率可大大加强。这种迅速的变化，可能对地下水本身和环境带来严重的后果。地下水的质与量之所以变化
2、，主要是由于水量和溶质成分在补充和消耗上的不平衡所造成的。所谓地下水均衡，就是指这种在一定范围、一定时间内，地下水水量、溶质含量及热量等的补充（流入）与消耗（流出）量之间的数量关系。当补充与消耗量相等时，地下水（量与质）处于均衡状态；当补充量小于消耗量时，地下水处于负均衡状态；当补充量大于消耗量时，地下水处于正均衡状态。地下水在天然条件下，一般多处于均衡状态；在人为活动影响下，则可能出现负均衡或正均衡状态。从上述概念可知，地下水动态与均衡之间存在着互为因果的紧密联系。地下水均衡是导致动态变化的实质，即导致动态变化的原因；而地下水动态则是地下水均衡的外部表现，即动态变化的方向与幅度是由均衡的性质
3、和数量所决定的。2 研究地下水动态与均衡的意义研究地下水动态与均衡，对于认识区域水文地质条件、水量和水质评价，以及水资源的合理开发与管理，都具有非常重要的意义。任何目的、任何勘查阶段的水文地质调查，都必须重视地下水动态与均衡的研究工作。由于对地下水动态规律的认识，往往要经过相当长时间的资料积累才能得出结论，因此在水文地质调查时，应尽早开展地下水动态与均衡研究。其研究意义具体表现在：（1）在天然条件下，地下水的动态是地下水埋藏条件和形成条件的综合反映。因此，可根据地下水的动态特征分析、认识地下水的埋藏条件、水量、水质形成条件和区分不同类型的含水层。（2）地下水动态是均衡的外部表现，故可利用地下水
4、动态资料去计算地下水的某些均衡要素。如根据次降水量、潜水位升幅和潜水含水层给水度计算大气降水的入渗系数；根据潜水位的升幅或降幅计算地下水的储存量及潜水的蒸发量等。（3）由于地下水的数量与质量均随着时间而变化，因此一切水量、水质的计算与评价，都必须有时间的概念。如对同一含水层来说，在雨季、旱季、丰水年、枯水年，其水资源数量与水质都可能大不一样。因此，地下水动态资料是地下水资源评价和预测时必不可少的依据。（4）用任何方法计算的地下水允许开采量，都必须能经受地下水均衡计算的检验；任何地下水开采方案，都必须受地下水均衡量的约束。为尽可能地减少开采地下水引起的负作用，开采量一般不能超过地下水的补给量，即
5、不应破坏地下水的均衡状态。（5）研究地下水的均衡状态，可预测地下水水量、水质及与地下水有关的环境地质作用的变化及总体发展趋势。因此，在各种目的的水文地质勘探中，都规定进行一定时期的地下水长期观测，以便进行地下水动态与均衡的研究。勘探阶段愈详细，长期观测工作量愈大，要求的精确度愈高。3 地下水动态和均衡研究的基本任务一、研究地下水动态的基本任务（1）正确布设地下水动态监测网点，对动态监测的频率、监测次数及监测时间作出科学的规定。地下水动态监测点的布置形式和位置，主要决定于水文地质调查的主要任务。动态监测成果要满足水文地质条件的论证，地下水水量、水质评价及水资源科学管理方案制定等方面的要求。对干不
6、同的勘查阶段，对以上要求各有侧重。为阐明区域水文地质条件服务的动态监测工作，其主要任务在于查明区域内地下水动态的成因类型和动态特征的变化规律。因此，监测点一般应布置成监测线形式。主要的监测线应穿过地下水不同动态成因类型的地段，沿着区域水文地质条件变化最大的方向布置。对不同成因类型的动态区，不同含水层，地下水的补给、径流和排泄区，均应有动态监测点控制。为地下水水量、水质计算与资源管理服务的动态监测工作，其主要任务是：为建立数学模型、水文地质参数分区及选择参数提供资料。鉴于地下水数值模型在地下水水量、水质评价与管理工作中的广泛应用，要求将相应的动态监测点布置成网状形式，以求能控制区内地下水流场及水
7、质变化。对流场中的地下分水岭、汇水槽谷、开采水位降落漏斗中心、计算区的边界、不同水文地质参数分区及有害的环境地质作用已发生和可能发生的地段，均应有动态监测点控制。地下水动态的监测点，除井、孔外，还应充分利用区内已有的地下水天然及人工水点。对有关的地表水体、各种污染源，以及有害的环境地质现象，亦应进行监测。科学规定地下水动态项目的监测频率、监测次数和时间，对于获得真实、完整的动态资料十分重要。对于不同的监测项目，监测的频率、次数和时间的具体要求虽有不同，但其总的原则是一致的，即要求按规定的监测频率、次数和时间所获得的地下水动态资料，应能最逼真地反映出年内地下水动态变化规律。以上问题的具体要求，可
8、参阅有关水文地质勘查和地下水动态观测规范。需强调的是，为了能从动态变化规律中分析出不同动态要素（监测项目）间的相互联系，对各监测项目的监测时间，在一年中至少要有几次是统一的。（2）根据所获得的地下水动态监测资料，分析地下水动态的年内及年际间的变化规律。动态变化规律，主要是指某种动态要素随时间的变化过程、变化形态及变幅大小等的水文地质意义，变化的周期性与趋势性，并通过不同监测项目动态特征的对比，确定它们之间的相关关系。（3）根据所获得的各种动态资料，考虑各种影响因素（水文、气象、开采或人工补给地下水等）的作用，确定区内地下水的成因类型。为认识区域地下水的埋藏条件，水质、水量的形成条件及有害环境地
9、质作用的产生和发展原因等，提供动态上的佐证。二、地下水均衡研究的基本任务（1）为进行均衡研究，首先要确定均衡区的范围及边界的位置与性质。当区域较大，各地段的地下水均衡要素组成又不相同时，应划分均衡亚区。为便于均衡计算，每个均衡区（或亚区）最好是一个相对独立的水文地质单元。均衡区的边界最好是性质比较明确、位置比较清楚的某一自然边界（或地质界线）。（2）确定均衡区内地下水均衡要素的组成及地下水水量或水质均衡方程的基本形式。在建立方程时，应考虑到，同一均衡区在不同的时段，其均衡要素的组成可能是不同的。因此，在均衡计算之前，还应划分出均衡计算的时段，即确定出均衡期。（3）通过直接（野外实测或室内测定）
10、或间接（参数计算）方法，确定出地下水各项均衡要素值，为地下水水量、水质的计算与预测提供基础数据。（4）通过区域水均衡计算，确定出区内地下水的均衡状态，预测某些水文地质条件的变化方向，为制定合理的地下水开发方案及科学管理措施提供基本依据。4 地下水动态与均衡的监测项目一、地下水动态监测项目对大多数水文地质勘查任务来讲，地下水动态监测的基本项目都应包括地下水水位、水温、水化学成分和井、泉流量等。对与地下水有水力联系的地表水水位与流量，以及矿山井巷和其它地下工程的出水点、排水量及水位标高也应进行监测。水质的监测，一般是以水质简分析项目作为基本监测项目，再加上某些选择性监测项目。选择性监测项目是指那些
11、在本地区地下水中已经出现或可能出现的特殊成分及污染物质，或被选定为水质模型模拟因子的化学指标。为掌握区内水文地球化学条件的基本趋势，可在每年或隔年对监测点的水质进行一次全分析。地下水动态资料，常常随着观测资料系列的延长而具有更大的使用价值，故监测点位置确定后，一般都不要轻易变动。二、地下水的均衡项目（或均衡要素）地下水的均衡包括水量均衡、水质均衡和热量均衡等不同性质的均衡。不同性质均衡方程的均衡项目（均衡要素），也就必然有所区别。在多数情况下，人们首先关注的还是水量问题，而水量均衡又是其它两种均衡的基础。因此，下面着重讨论水量均衡的组成项目。根据质量守恒定律，在任何地区，在任一时间段内，地下水
12、系统中地下水（或溶质或热）的流入量A（或补充量）与 流出量B（或消耗量）之差，恒等于该系统中水（溶质或热）储存量的变化量W。据此，我们可直接写出均衡区在某均衡期内的各类水量均衡方程。总水量均衡方程的一般形式为： BAW进一步写为（单位面积）： )()( 2211 RWZYRZYXPVh 式中：h潜水储存量的变化量，其中，为潜水位变动带内岩石的给水度或饱和差，h为均衡期内潜水位的变化值；V，P分别为 地表水体和包气带水储存量的变化量；X降水量；Y1，Y 2地表水的流入和流出量；Z1，Z 2凝结水量和蒸发量（包括地表水面、陆面和潜水的蒸发量）；W1，W 2地下径流的流入和流出量；R1，R 2人工引
13、入和排出的水量。潜水水量均衡方程的一般形式为： )()( 2211 RZWRZWYXhsff 式中：X f 降水入渗量；Z1，Z 2潜水的凝结补给量及蒸发量；Ws泉的流量；Yf地表水对潜水的补给量；R1，R 2人工注入量和排出量；其余符号同前式。承压水的水量均衡方程，比潜水为简，常见形式为： )()(21* kRWEh式中：*承压含水层的弹性给水度（贮水系数）；E1越流补给量；R2k承压水的开采量；其余符号同前式。对于不同条件的均衡区及同一均衡区的不同时间段，均衡方程的组成项可能增加或减少。如：当地下水位埋深很大时，Z 1和Z 2常常忽略不计。分析上述各水量均衡方程，可清楚地看到，一切水量均衡
14、方程均由三部分组成，即均衡期内水量的变化量(W)、地下水系统的 补给量(或流入量A) 和消耗量(或流出量B)。在补给量中，最重要的是降水入渗量（X f）、地表水入渗量（Y f）、地下径流的流入量（W 1）；在某些情况下，越流补给量(E 1)和人工注入量（R 1）也有较大意义；在 消耗量中，最重要的是潜水的蒸发量（Z 1）、地下径流的流出量（W 2）、地下水的人工排泄量（R 2和R 2k）；有时，泉水的溢出量（W s）和越流流出量（E 2）也很有意义。5 地下水动态的成因类型及主要特征地下水动态成因类型的划分，主要是根据地下水的水位动态过程曲线的特点，以及对地下水动态影响最大的自然和人为因素对地
15、下水动态成因类型进行划分。综合国内、外一些地下水动态成因类型分类方案，本书将地下水动态成因类型归纳为8种基本类型（见表61），而由基本类型又可组成多种混合成因类型。6 地下水均衡要素的测定方法一、潜水储存量变化量(h)的测定方法潜水储存量变化量由潜水位变化值h和水位变动带岩层的给水度（或饱和差）组成。h能通过水位观测孔实测获得。因此，确定潜水储存量变化量的关键在于值的测定。当潜水水位上升或下降时，值具有不同的物理意义。下降时，表征水位变动带地层的给水度；上升时则表征饱和不足量（或饱和差）。但当潜水面未上升到近地表的湿度变动带时，值仍可视为给水度。确定给水度的常用方法有：（1）室内参数测定法按要
16、求深度定期采取水位变动带内的岩土样，在室内测定饱和容水度（饱和含水量）、持水度、天然湿度（天然含水量），即可得出不同时段的值。此法取样繁琐，且难保证土样的天然结构不被破坏，而粘性土又难测出其持水度，故已很少使用。（2）根据抽水前后包气带上层天然湿度的变化来确定值据包气带中非饱和水流的运移和分带规律知，抽水前包气带内土层的天然湿度分布应如图 61 中的 oacd 线所示。抽水后，潜水面由 A 下降到 B（下降水头高度为己 h），故毛细水带将下移，由 aa段下移至 bb段，此时的土层天然湿度分布线则变为图中的oabd。对比抽水前后的两条湿度分布线可知，由于抽水水位下降，水位变动带将会给出一定量的水
17、。按水均衡原理，抽水前后包气带内湿度（含水量）之差（阴影面积），应等于潜水位下降 h 时包气带（主要是水位下降带）所给出之水量（h），即：hWZni ii121)(故给水度：hni ii121)(式中：Z i包气带湿度测定分段长度 (空间步长)；h抽水产生的潜水面下移深度（水位降深）；W1i，W 2i抽水前后Z i段内的土层天然湿度（含水量）；n取样数。土层的天然湿度，可采取原状土样在实验室测定，或利用中子水分计（中子仪）在钻孔中直接测定土层的含水量。（3）根据潜水水位动态观测资料用有限差分法确定值（卡明斯基有限差分法）如果潜水为单向流动（一维流），隔水层水平，含水层均质，可沿流向布置3个地下
18、水水位动态观测孔（图62），然后根据水位动态观测资料，按下式计算值： 2,2,3,12)(hthxtKttt 式中： szwh1,t，h 2,t，h 3,t1，2，3号观测孔 t时刻水位或含水层厚度；h 2t时段内2号孔水位变幅；垂向流入和流出量之和称综合补给强度；K渗透系数；x观测孔间距（空间步长）；t时间步长。如地下水流入和流出量以裘布衣公式表示，整理上图阴影部分水均衡式，也可得上式。如潜水为二维流，观测孔作方形网格布置（图63）时，仍可按上述方法得出下列计算式： 52,42,3,2,15 )( hthhxt ttttt 式中各符号意义同前式。在上式的参数中，值常未知，但可选择近乎常数的两
19、个时段（2个t），写出两个计算式，解出和值。此法优点在于，能确定较大范围内的值，可用于基岩和地下水深埋区，对不同边界适应性较强。由于未知，常取0。但此时平原区的h较小，计算的相对误差大。二、降水入渗补给量（X f）及蒸发量（Z 2）的确定1地中渗透仪（或蒸渗仪lysimeter）测定法这是较老但又是唯一可直接测到降水入渗补给量和潜水蒸发量的方法。此方法仪器的结构装置如图64所示。整个装置由左方的地中渗透计、右方的给水观测装置构成。地中渗透计的圆筒内装有均衡地段的标准土柱，土柱下方为砂砾和滤网组成的外滤层（图6一4中的2，3）。给水观测部分由供水（盛水）用的有刻度的马利奥特瓶（图中10）和控制地
20、中渗透计筒内水位高度的盛水漏斗（11）及量简（l4 ）组成。两部分以导水管连结，将两端构成统一的连通管。1土柱其工作原理如下：首先调整盛水漏斗的高度，使漏斗中的水面与渗透计中的设计地下水面（相当潜水埋深）保持在同一高度上。当渗透计中的土柱接受降水入渗和凝结水补给时，其补给水量将会通过连通管（4）和水管（13）流入量筒（14）内，可直接读出补给水量；当土柱内的水面产生蒸发时，便可由漏斗供给水量，再从马利奥特瓶读出供水水量（此即潜水蒸发消耗量）。在测定凝结补给量时，应在该渗透计上方加棚，以隔离降水。此法装置可用多个不同岩性和不同水位埋深的土柱，分别观测其降水补给和蒸发值。本方法缺陷是，很难如实模拟
21、天然的入渗补给条件，故其结果的可靠性有时值得商榷，而且此法只适用于松散岩层。2通量法包气带土壤水分运动所遵循的基本规律是非饱和水流运动的达西定律（Darcys law或Richards equation）和质量守恒原理（conservation of mass），在实际应用中，可以直接应用达西定律和质量守恒原理分析或解决水量均衡问题。下面介绍具有重要应用价值的土壤水分通量法来求降水入渗补给量及蒸发量。要研究降雨或灌溉对土壤及潜水的入渗补给，大气蒸发作用下土壤及潜水的消耗，从而分析四水（大气水、地表水、土壤水和地下水）的转化关系等，在田间监测土壤水分的分布和运动是十分必要的。田间土壤水分运动，可
22、近似视为一维垂向的流动。于是，连续方程可简化为上式由 z*至 z 积分，得式中，q（z*）和q（z ）分别表示高度为z* 和z处的土壤水分运动通量（单位时间，通过单位面积的量：cm3/cm2s）。当时间由t 1改变到 t2时，以Q （z* ）和Q（z）分别表示在此时段内通过z*和z处单位土壤断面面积上的水量(cm 3/cm2)，由上式积分或直接由质量守恒原理写出无源汇情况下的水量平衡方程为土壤含水率的分布（z，t）在田间可用中子测水仪（或其他方法）监测。若测得某一断面z*处的土壤 水分运动通量q（z*）或单位面积上的水量Q（z* ）（通量已知的断面），土壤中任一断面z 处的通量q（z）及单位面
23、积上的水量Q （ z）(Q（z）= q（z）t)（通量未知的断面）便可由上式计算出。确定某一断面z*处的通量，主要应用达西定律，其方法有零通量面法、表面通量法和定位通量法，统称为土壤水分运动通量法。（1）零通量面法土壤水的重力势 g由垂直坐标z决定，若用负压计测得土壤剖面各点的基质势 m，则可得到总水势 mz的分布，如图210土壤中任一点的土壤水分通量由达西定律Zz*zXq(z)*q(z)给出。当水势梯度 时，该处的通量 ，则称该处的水平zKqm)( 0z0q面为零通量面ZFP（zero flux plane ），位置记为z 0。土壤剖面中出现零通量面时，可根据水势的分布特点，区分为以下几种类
24、型：a单一聚合型零通量面。若在降雨（或灌溉）前，土壤长期处于蒸发状态，上层土壤的水分只由潜水补给，水分自下而上在土壤中运移，水势自潜水面向上是逐渐减小的。在连续降雨过程中，由于水分不断入渗，上部土壤的水势将随之增加。如果水分的入渗尚未对下部的土壤产生明显的影响，此时水势的分布则如图210中（a）所示。在此情况下，土壤中某处水势出现最小值。该处 ，即为零通量面ZFP，也可视其为入0z渗前锋面。由于这种情况下土壤水分由上下两侧向零通量面处迁移，故称为聚合型。b单一发散型零通量面。当降雨停止且入渗锋面已下移到潜水面以后，上部土壤开始蒸发，水分自下面上迁移。与此同时，下部土壤水分继续处于向下入惨状态。
25、这时，水势的分布如图210中（b）所示。在此情况下，土壤中某处水势达最大值，该处 ，0z即为零通量面ZFP。由于这种情况下土壤水分自零通量面处分别向上、向下运动，故称为发散型。c具有多个零通量面。如图 210 中（c）所示.，这发生在间隔降雨、入渗和蒸发交替出现的情况下，具有多个零通量面。 零通量面的位置不变时：i地表处蒸发量：当零通量面存在时，该断面即为通量已知的断面 z*。若由 t1 至 t2 这一时段内（t ），零通量面的位置不变，测得 t1 和 t2 时刻的土壤含水率 （ z，t 1）和 （z，t 2），如图211。利用式（267）可计算出 t 时段内任一断面处单位面积上所流过的土壤水的水量 Q（z）。若 以 Qs 表示地表处相应的水量，其值可由下式计算：
