《水资源评价方法》由会员分享，可在线阅读，更多相关《水资源评价方法（15页珍藏版）》请在金锄头文库上搜索。
1、水资源研究第 25卷 第 2 期（总第 91期）2004 年 6 月平原区地下水资源评价方法综述刘 予 伟 金 栋 梁长江水利委员会 水文局，湖北 武汉 430010）摘 要： 全面阐述了平原区地下水资源评价方法，包括水文地质参数的获得和选取以及地下水资源量的计算方法，并通过实例来评价其合理性和可靠性。关键词：平原区；地下水资源；评价方法；综述平原区包括一般沿江、沿湖、沿海平原和山间盆地平原两类。就长江流域而言，前者有 洞庭湖平原、江汉平原、鄱阳湖平原、太湖平原、长江中下游沿江平原和江苏、浙江沿海平 原。后者有成都盆地、汉中盆地和南阳盆地。地下水评价对象是与大气降水和地表水体有直 接联系的浅层
2、地下水，一般仅评价矿化度小于2 g/L 的多年平均淡水资源，以现状条件为评 价基础，以水均衡法为主评价出各项补给量和排泄量。1含水层参数的确定含水层参数是定量描述含水层物理特性的指标或系数，是评价含水层的主要依据。在计 算各项补给量和排泄量时都要根据准确的参数来计算。主要参数有:潜水变幅带给水度（M ），降水入渗补给系数（a ），潜水蒸发系数（c），渠系渗漏补给系数（m），灌溉 入渗补给系数（B ）和水稻田渗漏率（0 ）等。现将上述几种含水层参数的确定分述如下。1.1潜水变幅带给水度（“）给水度（p ）是指饱和岩土在重力作用下，自由排出重力水的体积与该饱和岩土相应 体积的比值，它是一个无因次大
3、于零而小于1 的数值。可通过简易测筒或地中渗透仪试验、 利用地下水动态观测资料分析、剖面含水率测量和抽水试验等方法求得。1.1.1 简易测筒和地中渗透仪法用一个金属圆筒，将被测给水度的原状土（即保持天然结构的土层）装入筒内，使土层充 水达到饱和状态，然后在上部加盖，但不密封，防止水分蒸发，筒的下部留有排水孔，在重 力作用下，筒中的水会自由地从排水孔中流出，测量排出水的体积。排水体积和筒内土体积 之比即为给水度。此种测筒，制作和操作都甚简便，曾在第一次全国水资源评价中广为使用。另一种类似于测筒的是地中渗透仪，图1（略）是地中渗透仪的示意图。利用潜水位控制， 可将左边测筒内土体积饱和到任意位置，然
4、后将连通管控制进水，测量由连通管自由流出水 的体积，使之与其土体体积相比，即得给水度。地中渗透仪虽造价较高，但由于它可进行多项参数的观测试验，故我国的黄淮海平原区 有多处此种实验装置。1.1.2 包气带剖面含水率法设有一均质土层，其颗粒组成较粗，颗粒之间的孔隙排水滞后作用时间短，假设在无蒸 发的条件下地下水位上升(或下降)值为 H。在水位变化前后分别测定水层剖面的土壤面含水率曲线如图 2；图中横座标代表土壤含水率，纵座标代表埋深，纵横座标所夹的面积 即代表含水量(以mm计)。由于是均质土，无土壤水蒸发，又不考虑滞后作用，因此在水位变化前后的土壤含水率 剖面线应是平行的，即 AB 平行 CD,
5、W值是田间持水率，W是饱和含水率。ArnB间C CD，D 形状完全一样，面积相等。即水量相等。因此，不难得出ABDC 的 面积等于 BDFE 。根据实测资料作出图2后，在图上便可量出 ABCD 的面积即含水量(以 mm计)，设其含水量为HA，则有给水度：ABCDH 一 H(1)卩=BDFE = ABCD(1)AHAH由上述可知，本方法原理十分简单，但在实际应用中可能会有困难，因为自然界非常复 杂的，需要在具体实践中经过适当处理求得满意结果。1.1.3 回归分析法根据水量平衡原理，在一定时期内，一定区域潜水或浅层地下水量的变化，应等于其收 入与支出量的差值、水量变化反映在水位变幅上。收入：降水入
6、渗补给量二入渗系数X降水量X面积二a PF侧向径流补给量=渗透系数X水力坡降X时段长X横断面=K1I1A1A t 支出：人工开采量=区域开采地下水体积=V开 开潜水蒸发量=给水度X地下水位降深X面积二p HF狈恫径流排泄量=渗透系数X水力坡降X横断面X时段长=呼2人24 t 故有：水量变化=收入-支出u HF=a PF+K I A t- V -p HF-K I A t(2)11 1开2 2 2式中U为区域地下水变幅带平均给水度； H为 t时段内，区域平均地下水变幅；F为计算区域面积；a为降水入渗系数； P为 t时段内区域平均降水量；KK2为入流断面和出流断面的渗透系数；AA2为入流和出流断面面
7、积；I2为 t时段内，入流和出流断面平均水力坡降；A t为计算时段长；V开为 F区域内在 t时 开段中地下水开采体积； h为 t时段内潜水蒸发引起的地下水位下降值。将式(2)除以( p F )得：aKA IA t 1VKAH = P + 廿（）-（十）-Ah -旷） 卩 F卩F卩a1KK令：a =a =-a = 1 a=212 3 4A I A tA I A tV又令： h=h =2 2,h=开入F出f开f将上述各式代入式(3)得：ah = a p - a h + a h - a h - a h(4)1 2 开 3 入 4 出式中h称为开采模数，h 、h分别为地下水流进流出的单位渗透模数，都以
8、mm计。开入出各变量 H、P、h h 、h 、卜h 之间为相关关系，故可将此式视为多元回归方程：开 入出ah = a p - a h + a h - a h ah + a(5)12 开 3 入 4 出0式中 a为常数项； a、a、a、a为待定回归系数，可以使用最小二乘法求得，利用01234a、a、a、a与p、a、K、K 的关系式，就可获得a、p、K、K 了。12341212在具体计算中，如果知道某些因素影响小可以忽略时，则计算工作大为简化。如地下径流微弱时，h入和h出可略而不计，则上式变为：入出ah + ah = a p - a h + a(6)1 2 开 0由于 h系数为1,故可移至左端，这
9、样式(5)成为二元线性回归方程了。 当埋深较大时，潜水蒸发可略而不计，则式(6)为：ah = a p - a h + a(7)12 开 0当时段内无降水、地下径流微弱、地下水埋深大，无潜水蒸发时，回归方程为：ah = -a h + a(8)2 开0式(8)成为一元线性回归方程。但必须注意的是开采量一项是不能缺少的，有了开采项，才能算得由此也可了解到通过抽水可以求得给水度 p 值。推求给水度的方法还有很多，如坑测法、入渗差值法、潜水位增幅法、优选法等。由于参数给水度在地下水资源评价中极其重要，它的精度直接影响资源估算的数量， 所以水文地质工作者投入试验研究的时间也是很多的，取得很多成果，根据淮委
10、的研究，现 在应用的给水度 p 是地下水变幅带的平均给水度，实际上给水度是随地下水埋深而变的， 在埋深0.2 m以内为最大值。大于1.0 m基本稳定不变。如图3所示。我国第一次水资源评 价时，全国各流域对各种岩性的给水度进行大量试验研究，经综合归纳后的给水度如表 1。(略)1.2降水入渗补给系数(a)降水入渗补给系数 (a ) 是地下水资源评价和系统管理中常用的重要参数，是地区水资 源主要补给来源，降水入渗系数选用是否准确合理对地下水资源的计算有着决定性的作用。降水入渗补给系数，为降水入渗补给地下水的量(P)与降水总量(P)的比值，即：ra = P / P = r ( a ) a h / P(
11、9)r式中a为降水入渗补给系数；P为时段降水入渗补给量，(mm)；rP 为时段降水总量， (mm)； h为时段降水入渗引起的地下水位升幅，(mm)；p ()为随埋深()而变的地下水位变幅带含水层的变给水度。 影响a值的因素很多：时段降水总量、降水强度、降水时间分布、地下水埋深、时段初包气 带含水量大小，土壤类别、结构、地表植被等。因此，a 值是随时间和空间变化的。但对某一特定地区，由于土壤岩性和气候条件变化不大。 影响a 值的主要因素是：降水总量、地下水埋深和时段初包气带含水量大小。a 值的确定方法主要有：地下水动态资料分析法、人工降雨模拟试验法和含水层参数率 定模型法等。含水层参数率定模型是
12、建立在地下水长观资料基础上的，对资料要求较高，一 般难以应用。小型人工降雨模拟试验，代表性不尽人意，一般只作验证性试验。目前确定a 值的主要方法还是地下水动态观测资料分析等法。兹分述如下：1.2.1 地下水动态资料分析法根据长观井的地下水动态资料，用如下的水均衡公式计算降水入渗补给系数：Q - Q - Q 士卩 AH Pd = 开 河 側（io）FP式中 a 为时段降水入渗补给系数；Q开为地下水开采量；Q河为河道渗漏补给量；Q側为侧向补给量；M 为给水度； H开河側为时段地下水位升幅； F 为流域面积； P 为时段降水总量。采用多元回归分析法（见给水度分析确定部分）可以求降水入渗系数。当流域内
13、无开采、无灌溉、无侧向补给、无河道渗漏时，仅根据地下水位的升幅及给水 度数据，便可计算出降水入渗系数 a 。1.2.2 补偿疏干法在开采条件下，在雨季所得到的补给量除满足当时开采外，并用以补偿地下水储存量， 因之，引起地下水位上升，其计算式为：V =V -V -V +u HF补开河側V =EV年补a =V /FP年 年 年式中V补为时段降水入渗补给量；补V 年为年降水入渗补给量；年V 河为时段河道渗漏补给量；河V 側为时段侧向补给量；側 H为时段地下水位升幅；F 为流域面积。1.2.3 岩溶区降水入渗补给量的推求据邵正介绍，选择岩溶区内枯季断流的泉，并确定其泉域（面积），并查清泉域内的厚度 较
14、大的由粘土亚粘土覆盖的非岩溶区面积。待雨季来临时，泉水涌流，在测定泉域内平均降 雨量外，还要测出泉的涌水量及泉域内人畜饮水、灌溉水量消耗、可算出降水入渗系数。例如：山东东平县中套泉，泉域面积11.875 km+2.泉域北部有洪坡积粘土和亚粘土覆 盖厚度达4m以上的非岩溶区。1983年5月至7月29日，泉水断流。7月27日及29日分别 降雨 53.1、148.5 mm。7月30日泉水开始外流。至1984年4月17日止，泉域总降水量471.5 mm。测得泉水溢出量为48.62万m3。在此期间，调查到引用泉水灌溉及人畜饮用共 耗水42.11万m3,因此泉水总溢出量为90.73万m3。由此计算得：90
15、.73 x 10 4=0.160.471 x11.875 x106覆盖土层厚度大于4 m的面积2.7 km2(根据当地土层确定a =0.13),则张夏灰岩区a 值如下：907300 0.471 x 2.7 x 10 6 x 0.13a = 0.170.471 x (11.875 2.7)x1061.2.4降水入渗系数的修正降水入渗系数由于受到多种因素的影响，某时段的a次值，几乎没有实用价值，所以 次一般仅采用其均值a 。为了消除影响因素时间上的变化，提高计算精度可作如下处理：P p = a - p(11)rP式中 Pr 为时段内降水入渗补给量， (mm)； P 为时段内总降水量， (mm)；P 为时段内多年平均总降水量，(mm)； a为时段内多年平均a 值。经过这样处理后的成果比较符
