潜水—--微承压含水层给水参数确定方法及研究.doc

.页眉..页脚潜水——微承压含水层给水参数确定方法的研究陈庆秋(华南理工大学南方水政策研究中心)摘要:基于考虑潜水层重力释水滞后作用条件下潜水——微承压含水层完整单井非稳定流的解析解,讨论了布尔顿给水强度公式中给水度的物理内涵,指出布尔顿给水强度公式中的给水度是“动态稳定给水度”;探讨了布尔顿潜水井流模型中延迟指数的物理意义,得出了如下猜想性的定义:布尔顿延迟指数1/?表征潜水含水层在完整单井定流量抽水时,重力疏干迟后性的一个水文地质参数;对于确定的含水层,该参数是抽水流量和抽水时间函数,当抽水流量一定时,在抽水过程中的某一时刻t的1/?值等于潜水面从埋深为h(t’)外下降一个单位深度后,在埋深为h(t’)-1[h(t)坐标向下为正]的单位面积上获得给水度?大小的水量所需。该文在探讨了布尔顿给水强度公式中的给水度及延迟指数的物理意义后,还提出了一种考虑潜水层重力释水滞后作用条件下潜水——微承压含水层完整单井非稳定流参数的确定方法。考虑到所提出的潜水——微承压含水层完整单井非稳定流参数的确定方法借鉴布尔顿的第二潜水井模型的参数确定方法,该文还讨论了布尔顿潜水井流模型适应性。关键词:潜水——微承压含水层;含水层参数;方法1考虑潜水层重力释水滞后作用条件下潜水————微承压含水层完整单井非稳定流的数学描述潜水层:??????????????????????0201010)(111*11210zztttsssdestsbssK???????.页眉..页脚微承压层:上式中::2s??????????????????????????????????TQrsrsstsTbssTKrsrrsrrt??200122022221122221K为潜水层的渗透系数;为潜水层的水位降深;1s为潜水层的水位降深;为微承压含水层的水头降深;?1为潜水层的延迟指数;b为潜水层的初始厚度;*1?为潜水层的弹性释水系数;1?为潜水层的给水度;2?为微承压含水层的储水系数;T 为微承压含水层的导水系数; Q 为抽水流量。??????????????ttdestsrsrrsT0)(21222222)1(???????Qzr1,?zKb*2,?TO图 1 潜水——..页脚式中:潜水——微承压含水层的微承压层水头降深的解析解:式中:x为积分变量。;111bKbK??????;121??????xdBrxJshuutxchuexTQsu)(}]2)1([1{24'0222202???????????;1211??????????;2)1(21xtu????;24)1(22222xxtu??????;1'??TB?;22uueechu???;22uueeshu???)(0J为第一类零阶贝塞尔函数。.页眉..:计算出上式积分项,即可得到潜水层水位降深的解析解:式中井函数为:——“动态稳定给水度”笔者在《关于布尔顿给水强度公式中给水度的新概念》一文中提出了一个新概念——“动态稳定给水度”。他认为,当地下水水位以一不为零的速度匀速下降时,给水度不会达到完全给水度,而是在地下水水位下降一定时间后,稳定于小于完全给水度的某一常数,该常数即是“动态稳定给水度”。??????desKbsstt)(021121???????),,(4)1('1212111BrWKTQbsKbs??????????)(2]4)1(11[1{2),,('2'12220'??????????????????????uueexxxexBrWuutt.)(}])(24)1(11['0'2'12222xdBrxJuueexxxtuu????????????????)1(212'1??xu?????222'24)1(21xxu???.页眉..页脚“动态稳定给水度”具有如下性质:不同的地下水水位匀速下降速度,有不同的“动态稳定给水度”值;当下降速度越小时,“动态稳定给水度”越大;当下降速度趋于零时,“动态稳定给水度”趋于完全给水度。图2为地下水水位速度匀速下降时土壤含水量分布变化图。从该图可以看出:当地下水水位下降到一定深度后,土壤含水率剖面将随水位下降而平行下移,水位下降单位深度时,排出的水量为常数,含水层给水度在水位动态下移过程中达到相对稳定状态。水位下降速度越小,土壤含水率剖面平行下移收敛处的含水率越小,下降单位深度释放的水量越大,即“动态稳定给水度”越大。