第七章 水文地质参数计算课件

1、第七章 水文地质参数计算 水文地质参数是表征含水层水文地质特征的重要参数，其数值大小是含水层各种性能的综合反映，是地下水资源评价的重要基础资料。水文地质参数的计算，必须在分析勘察区水文地质条件的基础上，合理地选用公式（选用的公式应注明出处）。 主要包括：含水介质的渗透系数和导水系数、承压含水层的释水系数、潜水含水层的重力给水度、弱透水层的越流系数及水动力弥散系数等， 水文参数表征与岩土性质、水文气象等因素有关的参数，如：降雨入渗系数、潜水蒸发强度、灌溉入渗补给系数等。 主要方法有：野外试验、实验室测试及根据地下水动态观测资料采用理论公式求取，或采取数值法反演求参等7.1渗透系数（K） 渗透系数

2、的大小与介质的结构（颗粒大小、排列、空隙充填等）和水的物理性质（液体的粘滞性、容重等）有关。 表征岩石空隙渗透性的指标是渗透率，它取决于空隙的大小、形状、空隙度等因素。与渗透系数的关系为：kK 渗透系数渗透率液体重率液体的粘滞系数由于水的物理性质在常温下变化不大，可以忽略。利用裘布依公式求渗透系数时的注意事项：1）含水层的井壁边界条件2）影响半径 含水层很少能满足该条件。理论上，在抽水后的实际降落漏斗范围内，只有当观测孔与抽水井的距离（r）小于0.178倍的R时，水位降深s与r才属对数关系；当r大于0.178倍后就变为贝塞尔函数关系，而贝塞尔函数斜率小于对数函数，这也就是观测孔越远计算的k值越

3、大的原因。3）天然水力坡度的影响 利用上游观测孔求得的k值偏小，下游求得的k值偏大，主要是地下水是在天然水力坡度作用下运动的。4）抽水降深大小的影响 在实际计算中应选择抽水降深较大，同时抽水井周围的三维流或紊流影响较小，且s与r保持对数关系的观测孔资料计算渗透系数。利用泰斯公式求水文地质参数时应注意：（1）对于承压完整井抽水，当井内流速达到一定程度（1m/s）时，在井附近会产生三维流区，此时，利用主孔资料或布置在三维流区内的观测孔求解时，将产生三维流影响的水头损失，所以，应对实测降深值进行修正。（2）由于地下水运动天然水力坡度，在地下水上下游所计算的参数数值差异较大。解决的方法是在形成的降落漏

4、斗范围内布置较多的观测孔，求水文地质参数的平均值。（3）注意边界条件的影响。三、利用数值法求水文地质参数 数值法按其求解方法可分为：试估校正法和优化计算法，一般多采用试估校正法。 这种方法利用水文地质工作者对水文地质条件的认识，给出参数初值及其变化范围，用正演计算求解水头函数，将计算结果和实测值进行拟合比较，通过不断调整水文地质参数和反复多次的正演计算，使计算曲线和实测曲线符合拟合要求，此时的水文地质参数即为所求。 然而，求参结果的可靠性和花费时间的多少，除了取决于原始资料精度外，还取决于调参者的经验和技巧。-1001020304050607080G322孔200620042002200019

5、98199619941992计 算 水 位实 测 水 位地 下 水 位 (米 )时 间 (年 )地下水动态拟合曲线图 -1001020304050607080梅 花 寨 孔20062004200220001998199619941992计 算 水 位实 测 水 位地 下 水 位 (米 )时 间 (年 )地下水动态拟合曲线图 7.2给水度 给水度是表征潜水含水层给水能力或储水能力的一个指标，给水度和饱水带的岩性有关，随排水时间、潜水埋深、水位变化幅度及水质的变化而变化。给水度的确定方法 确定的方法有非稳定流抽水试验、实验室法、野外现场试验法、根据潜水位动态观测资料用有限差分法确定、经验数值法等。

6、给水度强弱相对比较（分选性较好时，典型推荐值）砾石0.12 粗砂 0.1 中砂 0.08 细砂 0.06 亚粘土 0.01 粘土0.002岩 性给水度岩 性给水度粘 土0.020.035细 砂0.080.11亚粘土0.030.045中细砂0.0850.12亚砂土0.0350.06中 砂0.090.13黄土状亚粘土0.020.05中粗砂0.100.15黄土状亚砂土0.030.06粗 砂0.110.15粉 砂0.060.08粘土胶结的砂岩0.020.03粉细砂0.070.010裂隙灰岩0.0080.10漏斗疏干法 在潜水面平缓、天然地下径流量很小的地区，抽水井所抽出的水主要来自降落漏斗疏干的水量，

7、随着抽水时间的延长，降落漏斗在不断扩展，只要将某一时刻以前抽出的水量，除以该时段的降落漏斗体积，即可得到给水度，计算公式如下： （6-5）式中抽水井的流量（m3/d） t抽水开始以后的时间（d） V降落漏斗的体积（m3） 降落漏斗的形状一般为不规则的漏斗形。当观测孔较多时，可根据实测的降落漏斗形状和漏斗内降深等值线，算出漏斗体积V。VtQ 7.3降水渗入系数的确定（一）基本概念 降水渗入系数是指降水渗入量与降水总量的比值。 值的大小取决于地表土层的岩性和土层结构、地形坡度、植被覆盖以及降水量的大小和降雨形式等，一般情况下地表土层的岩性对值的影响最显著。 一般情况下，地表土层的岩性对值的影响最显

8、著。降水入渗系数可分为次降水入渗补给系数、年降水入渗补给系数、多年平均降水入渗补给系数，它随着时间和空间的变化而变化。岩性P年（mm）粘土亚粘土亚砂土粉细砂砂卵砾石5000.020.010.050.020.070.050.110.080.121000.010.030.020.060.040.090.070.130.100.152000.030.050.040.100.070.130.100.170.150.214000.050.110.080.150.120.200.150.230.220.306000.080.140.110.200.150.240.200.290.260.368000.09

9、0.150.130.230.170.260.220.310.280.3810000.080.150.140.230.180.260.220.310.280.3812000.070.140.130.210.170.250.210.290.270.3715000.060.120.110.180.150.2218000.050.100.090.150.130.19不同岩性和降水量的平均年降水入渗补给系数值 注 东北黄土与表中亚粘土相近，陕北黄土含有裂隙，其与表中亚砂土相近（引自水利电力部水文局中国地下水资源）。 （二）降水入渗补给系数的确定方法（二）降水入渗补给系数的确定方法1.利用地下水位动态资料

10、计算降水入渗补给系数 这种方法适用于地下水位埋藏深度较小的平原区。我国北方平原区地形平缓，地下径流微弱，地下水从降水获得补给，消耗于蒸发和开采。在一次降雨的短时间内，水平排泄和蒸发消耗都很小，可以忽略不计。 根据降水过程前后的地下水位观测资料计算潜水含水层的一次降水入渗系数，可采用下式近似计算降水入渗补给系数：Xthhh/max式中 一次降水入渗系数； 降水后观测孔中的最大水柱高度（m）； 降水前观测孔中的水柱高度（m）； 临近降水前，地下水水位的天然平均降（升）速（m/d）； t 观测孔水柱高度从 变到 的时间（d）； Xt日内降水总量（m）。注：这种方法的适用条件是几乎没有水平排泄的潜水。

11、在水力坡度大、地下径流强烈的地区，降水入渗补给量不完全反映在潜水面的上升中，而有一部分水从水平方向排泄掉了，则会导致计算的降水入渗系数值偏小。如果是承压水，水位的上升不是由于当地水量的增加，而是由于压力的变化。以上情况本方法不适用。maxhmaxhhhh径流条件强地区2.回归分析法确定降水入渗补给系数 根据降水量和地下水的排泄量的系列资料，可用逐步回归分析法估算降水入渗补给量，从而可以确定降水入渗补给系数。 地下水的排泄量（主要指水平排泄，不包括潜水蒸发等垂直排泄）对于降水量存在时间滞后，用地下水的排泄量来代表同期的降水入渗补给量是错误的。但只要有地下水排泄量（河流基流量、矿坑排水量、泉流量等

12、）的系列资料和降水量的系列资料，用逐步回归分析法求出各时段降水量对于地下水排泄量的贡献，可以大致估算出降水入渗补给量。以当月的降水量和前几个月的降水量为自变量，共有n+1个自变量 。自变量可适当取多一些。地下水排泄量为因变量y。根据引入变量的数目和样本的大小，算出自由度f1和f2。给定显著水平，查表得到临界的F检验值F，然后进行逐步回归计算。如果算得的F值大于F，则认为这一自变量的影响是重要的，引入回归方程。否则予以剔除。最后得到的方程只包含对因变量y影响显著的 个自变量 即对地下水排泄量影响不大的月份的降水量不引入回归方程。最后得到的回归方程如下：式中 对地下水排泄量有显著影响的个月的降水量

13、； 回归系数lloooxbxbxbxbay2211nxxxx210,lxxxx210,lxxxx210,lbbbba2100,1l 考虑到入渗的机制，只有当降水量大于临界降水量后，才能产生入渗。可对作为自变量的降水量减去某一常数后，再进行回归计算。如再考虑到雨型的影响，可对七、八月份等有暴雨月份的自变量再乘以某一系数。 7.4潜水蒸发系数潜水蒸发系数 潜水蒸发是指潜水在土壤水势作用下运移至包气带并蒸发成为水汽的现象。在潜水埋深较小的地区，潜水蒸发是潜水的主要排泄途径，直接影响到潜水位的消退。单位时间的潜水蒸发量成为潜水蒸发强度，潜水蒸发强度的变化既受潜水埋深的制约，又受气象、土壤、植被等因素的

14、影响。（一）潜水蒸发强度确定1.经验公式法我国水利电利部水文局1982年11月编制的地下水资源调查和评价工作技术细则推荐下列三个经验公式计算潜水蒸发强度。（1）阿维里扬诺夫公式（1965年）： 式中 时段内日平均水面蒸发量（mm/d）； 时段内地下水平均埋深（m）； 潜水停止蒸发时的地下水埋深，又称潜水蒸发极限埋深（m）；n与包气带土质、气候有关的指数，一般取13。（2）沈立昌双曲线型公式（1979年）：式中 标志土质、植被、水文地质条件及其它因素的综合指数；a，b指数（无因次）；n)/1 (0000bak)1/()(0k（3）叶水庭指数型公式（1977年）：式中 衰减指数；其它符号意义同前。

15、 值的确定： 潜水蒸发极限埋深一方面与土质有关，其次也与作物等因素有关。例如，根据安徽省水利科学研究院五道沟试验场资料，在无作物的条件下，亚粘土的潜水蒸发极限埋深2.32.5m，亚砂土为4.0m；在有作物条件下，亚粘土的潜水蒸发极限埋深3.23.5m，亚砂土为5m左右。a100a02. 地中渗透仪法地中渗透仪是测定潜水蒸发、降水入渗、凝结水量的装置。 地中渗透仪示意图潜水蒸发系数潜水蒸发系数C泰森多边形法泰森多边形法60( )( )QFC ht蒸发7.5灌溉入渗补给系数 当引外水灌溉时，灌溉水经过渠系引入田间，灌溉水入渗补给地下水称为灌溉入渗补给，分为渠系的渗漏补给（条带状下渗）与田间灌溉入渗

16、补给（面状下渗）两类。 有的地区利用当地的水源（如抽取地下水）进行灌溉，灌溉水入渗后，地下水得到的补给应称之为灌溉回渗，它是当地水资源的重复量，不能作为地下水补给量。所以，灌溉入渗补给与灌溉回归应区别开来。 渠系渗漏系数m、田间灌溉入渗补给系数以及井灌回归系数的计算方法如下。（一）渠系渗漏补给系数的确定方法渠系渗漏补给系数m为渠系渗漏补给地下水的水量与渠首引水量的比值，即：令： 则：为简化起见，对 乘以折减系数，以代替右端应减去 ，写成如下形式：式中 Q引渠首引水量，可用实测的水文资料和调查资料；Q净经由渠系输送到田间的净灌水量； Q损渠系输水过程中的损失水量，包括水面蒸发损失、湿润渠底、两侧

17、土层的水量损失及退水填底损失等总和； 渠系有效利用系数； 修正系数，反应渠道在输水过程中消耗于湿润土层、浸润带蒸发损失水量的参数。引损净引QQQQm/ )(引净QQ /引损）（QQm/1）（1）（引损QQ /）（1mm分区衬砌情况渠床下岩性地下水埋深（m）渠系有效利用系数修正系数渠系渗漏补给系数长江以南地区和内陆河流域灌溉农业区未衬砌亚粘土亚砂土0.30.60.550.90.220.6部分衬砌亚粘土亚砂土40.40.70.300.800.180.45衬砌亚粘土亚砂土40.40.80.350.800.160.45北方半干旱半湿润区未衬砌亚粘土40.550.320.144亚砂土0.40.50.37

18、0.500.180.3亚粘土亚砂土0.40.550.320.140.3部分衬砌亚粘土40.550.700.300.090.135亚砂土40.520.730.350.100.17衬砌亚粘土亚砂土40.550.730.320.400.090.17亚粘土40.650.880.320.040.112亚砂土0.570.730.370.100.6 渠系渗漏补给系数值表 * 引自水利电力部水文局，中国地下水资源 （二）灌溉入渗补给系数的确定方法灌溉入渗补给系数是指某一时段田间灌溉入渗补给量与灌溉水量的比值，即式中 灌溉入渗补给系数； 灌溉入渗补给量（mm）； 灌溉水量（mm）。 灌溉入渗补给系数可采用试验方

19、法加以测定。试验时，选取面积为F的田地，在田地上布设专用观测井。测定灌水前的潜水位，然后让灌溉水均匀地灌入田间，测定灌水流量，并观测潜水位变化（包括区外水位）。经过时段后，测得试验区地下水位平均升幅，用下列公式计算：灌hhr/rh灌h式中 给水度； 计算时段（s）； 计算时段内试验区地下水位平均升幅（m）； 计算时段内流入试验区的灌水流量（m3/s）； 试验区面积（m2）。灌溉入渗补给系数主要的影响因素是岩性、地下水位埋深和灌溉定额 tQFhhhr灌/thQF地下水埋深灌水定额（m3/亩）岩性亚粘土亚砂土粉细砂黄土黄土状亚粘土1000.150.270.200.350.300.400.200.3

20、04m8m40700.080.140.120.18701000.100.220.150.250.200.300.150.200.150.201000.100.250.150.300.250.350.200.250.200.258m40700.050.060.050.10701000.050.150.060.150.050.200.060.100.050.131000.050.180.060.200.100.250.060.130.050.15田间灌溉入渗补给系数 值表* 1.引自水利电力部水文局，中国地下水资源。 2.表中井灌回归补给系数 取用表中相应栏的下限值 。 井井井（三）井灌回归系数的

21、确定方法（三）井灌回归系数的确定方法 在抽取当地地下水灌溉的井灌区，灌溉水的一部分下渗返回补给地下水，这种现象称为地下水灌溉回归。 井灌回归系数 是指灌溉水回归量与灌水量的比值，其测定方法与灌溉入渗补给系数相同。值得注意的是，试验时地下水处于开采过程中，则地下水位变幅中包括开采造成的变幅值，应予以考虑。井灌回归系数一般取值范围0.10.3。井7.6贮水率和贮水系数 贮水率和贮水系数是含水层中的重要水文地质参数，它们表明含水层中弹性贮存水量的变化和承压水头（潜水含水层中为潜水水头）相应变化之间的关系。 贮水率表示当含水层水头变化一个单位时，从单位体积含水层中，应水体积膨胀（或压缩）以及介质骨架的

22、压缩（或伸长）而释放（或贮存）的弹性水量，用 表示，它是描述地下水三维非稳定流或剖面二维流中的水文地质参数。s 贮水系数表示当含水层水头变化一个单位时，从底面积为一个单位、高等于含水层厚度的柱体中所释放（或贮存）的水量，用S表示。潜水层水层的贮水系数等于贮水率与含水层的厚度之积再加上给水度，潜水贮水系数所释放（贮存）的水量包括两部分，一部分是含水层由于压力变化所释放（贮存）的弹性水量，二是水头变化一个单位时所疏干（贮存）含水层的重力水量，这一部分水量正好等于含水层的给水度，由于潜水含水层的弹性变形很小，近似可用给水度代替贮水系数。承压含水层的贮水系数等于其贮水率与含水层厚度之积，它所释放（或贮存）的水量完全是弹性水量，承压含水层的贮水系数也