地下水向河渠间的运动.ppt
第三章 地下水向河渠的运动,均质含水层中地下水向河渠的运动,承压水向河渠一维稳定运动,无入渗潜水向河渠二维稳定运动,隔水底板水平,隔水底板倾斜,无入渗潜水向河渠三维稳定运动,平面流线呈辐射状,渗流断面复杂变化,均匀入渗潜水向河渠二维稳定运动,承压水向河渠一维不稳定运动,非均质含水层地下水向河渠的运动,1.定义为地下水运动要素是否随时间发生变化,变化为非稳定流,不变为稳定流;强调流场内所有点的运动要素都随时间变化。,2.产生稳定流的条件 流入 流出,必要条件,首先必须保持补给区和排泄区边界的水头 保持不变。,两者缺一不可。,稳定流与非稳定流计算公式不同,对地下水资源评价意义重大。,3.1 均质含水层中地下水向河渠的运动,一、稳定与不稳定流,充分条件:要求所研究的渗流区段内补给量排泄量。,二、承压水向河渠一维稳定运动-物理模型,1、物理模型(水文地质模型描述)条件:均质、等厚、承压含水层,两条平行河流完整切割含水层。两河水位分别为H1,H2,当两河水位稳定时,地下水可形成稳定流动,地下水可形成稳定流动。这时,流网显示地下水流线是一条平行的直线。,二、承压水向河渠一维稳定运动数学模型与求解(1),二、承压水向河渠一维稳定运动-数学模型与求解(2),二、承压水向河渠一维稳定运动-数学模型与求解(2),三、无入渗潜水向河渠二维稳定运动-(一)隔水底板水平,此问题属于剖面二维流动(vz0),潜水面是流线,由于其水力坡度不仅沿流线变化,而且过水断面也发生变化。,引入裘布依假定 把二维流(x,z)问题降为一维流(x)问题处理。,由于无垂向补排,故q沿0l不变,积分从断面1 至断面2,对比两式,若令z=0,即取基准面与底板一致,水头线方程,改变积分限(0 x),此水头线的特点:它是以x轴为对称轴的抛物线(上半支的一部分);,2.它与渗透系数K值的大小无关,(解法一),水头线方程,数学模型,(解法二),三、无入渗潜水向河渠二维稳定运动-(二)隔水底板倾斜,沿水平方向取x轴,它和底板夹角为;H轴和井轴一致。基准面可取在底板以下任意高度水平(00)。当 20,渗流长度可以用以水平孔距l来近似表示,水力坡度。即引入裘布依假设。,o,流量方程和水头线方程推导,水头线方程讨论,12渗流段的流量公式,1x渗流段的流量公式,水均衡原理,四、无入渗潜水向河渠三维稳定运动(一)平面流线辐射状,底板水平时,渗流宽度沿流向呈线性变化,水流在x、y、z三个方向都有分流速,根据裘布依假设,忽略垂向分速度,则可将水流简化为平面二维流。,四、无入渗潜水向河渠三维稳定运动(一)平面流线辐射状,四、无入渗潜水向河渠三维稳定运动(一)平面流线辐射状,流量公式,水头线方程,四、无入渗潜水向河渠三维稳定运动-(二)渗流断面复杂变化,潜水含水层隔水底板倾斜且不平整,呈三维流动,若允许忽略垂向分流速,则可利用裘布依微分方程,分离变量并积分得流量公式为,水头线方程,若任意,五、均匀入渗潜水向河渠二维稳定运动,1、入渗强度(W):单位时间内入渗补给地下水的水量。,2、问题描述:两条完整切割潜水含水层的平行河流,潜水含水层隔水地板水平,入渗强度W分布均匀,Wconst,流场为剖面二维流。,3、取单位渗流宽度的河间地块为研究对象,流网如图。以潜水面接受入渗补给为流线起点,由于两河水位不等,存在分水岭。,(一)流量方程推导,若x断面在分水岭的左侧,即xa,则,若x端面取在分水岭的右侧,则有,可知无论x在何处,均可得相同均衡式,取坐标系:规定流向与x方向一致q为正,入渗W0,蒸发W0,(一)流量方程推导,引入裘布依假定,分离变量,由断面1至断面x积分,当x=l时,h=h2,单宽流量方程:断面1 断面2,任意断面,(一)流量方程推导,引入裘布依假定,分离变量积分,单宽流量方程:断面1 断面2,流量方程的讨论,当 该式为无入渗补给潜水剖面二维稳定流动,此时河间地段呈单向流动。,2.当 向两侧河流的排泄量相等,各为补给量的一半。,(据河1断面流量q 方程),1,流量方程的讨论,说明:(1)在分水岭处水流不满足裘布依假定(2)在地下水排入河流的河床壁面,在河水位之上存在“出渗面”,也不满足裘布依假定。(3)只有离河边界和分水岭边界,水平距离l1.52.0M的垂直面才视为等水头面。,(二)水头线(浸润曲线)方程,讨论:1.当W0时,水头线是椭圆曲线的上半支 当W0时,水头线是双曲线方程 当W=0时,水头线是抛物线方程,由断面1至断面x积分得,2.有无入渗水头线方程相比较,前者多一项,当W 0,说明同一断面处有入渗条件比无入渗 条件的水位高。当 即河间地块中间断面水位抬高最大。,3.水头线与K有关,K值小,由于入渗引起的水位抬高值越大。,4.水头线方程可用于排水渠的设计。,当h1=h2,两渠水位相等时:处h为极大值,用h a表示,若两渠(沟)的水位已定,可以根据当地土质情况以不发生盐渍,(三)地下水分水岭位置的确定,分水岭公式的应用:,判断水库是否发生渗漏,分水岭公式的应用:,库水位的极限高度hmax,由图可见到水库蓄水过程,分水岭不断向水库方向移动,而当a=0时,是库水位得极限高度值。,指导野外调查工作 分析影响渗漏的因素(a0)K愈大,愈易渗漏。调查时水库要避开喀斯特发育带、构造破碎带或古河道发育带。渗流途径l小,即两河之间距离越短越易渗漏。调查时要避免将库址选在分水岭过于狭窄的地带。入渗补给量W愈小,愈易渗漏。在干旱地选址时,要避开存在渗透性差的覆盖层。邻河水位愈低(h2愈小),愈易渗漏。选址时应注意选在邻河水位高的地段。,分水岭公式的应用:,(四)入渗强度(W)的计算,若已知河间地段任意断面的水流值h和岩层的渗透系数K,就可以利用上式计算入渗强度W.,若未知K,则可值W/K,,可代W/K入分水岭公式,以判断水库是否发生渗漏,六、承压含水层中地下水向河渠一维不稳定运动,(一)定流量沟(渠)流,矿山水平排水廊边、河渠等。如果有闸门控制排水量保持其不变,即是定流量沟流问题,条件:承压含水层,均质、各向同性、等厚且半无限,一侧为定流量抽水(注水),(一)定流量沟(渠)流,erf(u)误差函数,见,(二)定降深沟渠流,被淹没的矿山、巷道(由水闸门关闭)的放水;,渠边的输退水、河道上建闸蓄水或开闸放水;,洪水期水位突然上涨。,假定渗流区含水层均质、各向同性、含水层厚度不变、侧向无限延伸;,初始水平面水平,即初始水力坡度为零,河渠水位突然下降。,然后保持不变,在它的影响下渗流区的地下水为一维不稳定运动。,解:,所以任意x断面的流量:,对于沟渠一侧的单宽流量,只要上式令,就可获得:,3.2 非均质含水层中地下水向河渠的运动-、分段法,一、分段法(一)水平层状非均质含水层中地下水运动稳定运动问题。三个均质等厚的水平岩层组成承压含水系统、其平面及剖面上流线互相平行,属于一维流动。由于按流面划分可将总水流划分成三个互不干扰的均质岩层地下水流,采取不同方法将非均质岩层转换成等效均质岩层中的地下水流问题来解决,常用的有分段法、等效厚度法、吉林斯基势函数法。,(一)分段法求解水平层状非均质问题,所以根据每一个单层计算单宽的公式有:,因为,流线在各层平行,在剖面上等水头线与铅垂线一致,故有:,显然,若存在几个含水层,有,取一等效渗透系数,厚度为,则有,(二)分段法求解透水性沿流向突变的非均质含水层中的地下水稳定运动问题。,条件:河流阶地附近潜水含水层中的地下水运动。隔水底板水平,阶地两侧岩性截然不同,但分别为均质岩层接触面近似垂直,潜水面十分平缓,满足裘布依假定。,根据潜水单层q公式:,(三)分段法小结,1、分段法:将一个复杂的渗流分解成几个简单的分渗流段而使问题得到解答的方法。2、两条要求:(1)各分渗流段的渗流状况,即运动要素或流网,与总渗流相应部分应保持一致。即分段之后,不能“走样”,否则各分渗流段之和不等于原渗流。(2)每一分渗流段应有现成的解答(即流量。水头线方程已知)或解答容易求得,否则分段法就没有优越性了。3、实现方法(1)分段法必须从分析流网开始。流线(面)隔水边界 等势线 等水头边界 所以分段界面应取流面或等水头面。,(2)分段总数应满足“每个分渗流段有现成解”的前提下、而且越少越好。,4、应用(1)承压无压流动:通常按有已知解的承压流和无压流两段求解。(2)复杂的三维或剖面二维流动,若存在一水平或接近水平的流面,将其作为分段界面。一个有隔水底板的分渗流段,一个有隔水顶板的分渗流段。,(3)复杂渗流边界(水工建筑物),5、总流量方程等于分段流量的并联或串联。,二、等效厚度法,因此,断面1和断面2的含水层假想厚度分别为:,二、等效厚度法,任意断面x处的含水层厚度(等效厚度)为:,任意断面x处的含水层厚度(等效厚度)为:,h为x断面处上含水层厚度。,三 吉林斯基势函数方法(1946),(一)原理,研究透水性在垂线上渐变的含水层的地下水流动问题。假定隔水底板水平,基准面取在隔水底板上z0,渗透系数沿垂直方向变化,沿水平方向不变,在含水层任一铅垂线上,取微分厚度dZ,对应KK(z),通过dZ断面的微分单宽流量:,b:含水层的顶面高度,对于承压:b=M,对于潜水:bh。吉林斯基将渗透系数沿垂向变化的含水层中的流量势定义为:,H为水头值,从隔水底板算起。,将吉林斯基势函数用于层状非均质含水层:,四、直接积分法,(一)渗透系数呈线性变化的含水层中的地下水运动,假定潜水含水层隔水底板水平,K呈线性变化,即:,我们若改变积分限,即自断面1至断面x,则得到水头线方程:,(二)渗透系数变化复杂的含水层中的地下水运动,引入裘布依假定,则:,实际上,上式右端的积分就是导水系数T的定义,即,分离变量进行积分:,(二)渗透系数变化复杂的含水层中的地下水运动,Tm是断面1、2之间的平均导水系数。如果导水系数由断面1的T1单调地变化到断面2的T2,则可近似取:,其中,T1和T2分别是断面1和断面2的平均渗透系数,