水文作业10月5号交 达西定律的理论推导水文作业10月5号交 达西定律的理论推导.doc

达西定律的理论推导 　　达西定律可以从多孔介质中层流运动所遭遇的阻力关系推导出来。 　　图1为沿流线方向s取得单元微分体, 长为ds, 断面积为dA; ?? 图1 渗透水体的受力 作用在单元柱体上的力有: 两端的孔隙水压力, 孔隙水流的自重及水流受到颗粒孔隙道的摩阻力F。沿土柱方向写渗流的三力平衡式(略去水流的惯性力) 　　pndA - ( p+ dp ) ndA - γndsdA sinθ- F = 0 　　因为 dz/ds= sinθ, h =p/γ+ z , dp= ?γ( dh - dz ) 　　代入上式则得 　　γndA dh + F = 0 ( 1)引用司托克斯对于一个颗粒上的层流阻力的公式D=3πμdν , 式中D 常被称为拖引力; d 为颗粒直径; v为颗粒周围沿渗流方向的局部平均流速; !为水的动力粘滞性; 为一个系数,决定于邻近颗粒的影响(对于无限水体中的圆球。若土柱中土颗粒数为N, 并引用一个球体系数β(圆球时β=π/6), 则总阻力应为 　　 F = DN [( 1- n ) dA ds]/[βd3] ?μdv? ( 2) 将( 2)式代入( 1)并考虑到断面上平均流速v= nv?及渗流坡降J= -dh/ds, 则得v=[βn2]/[令C=[βn2]/再令达西渗透系数k= Cd2 γ/μ, 即得达西定律表达式v= kJ。　渗透系数 k 是一个代表土的渗透性强弱的定量指标，也是渗流计算时必须用到的一个基本参数。不同种类的土，k 值差别很大。因此，准确的测定土的渗透系数是一项十分重要的工作。 渗透系数的测定方法 　　渗透系数的测定方法主要分“实验室测定”和“野外现场测定“两大类。 　　1.实验室测定法 　　目前在实验室中测定渗透系数 k 的仪器种类和试验方法很多，但从试验原理上大体可分为”常水头法“和变水头法两种。 　　常水头试验法 就是在整个试验过程中保持水头为一常数，从而水头差也为常数。 如图： ?? 常水头法测渗透系数k 　　试验时，在透明塑料筒中装填截面为A，长度为L的饱和试样，打开水阀，使水自上而下流经试样，并自出水口处排出。待水头差△h和渗出流量Q稳定后，量测经过一定时间 t 内流经试样的水量V，则 　　V = Q*t = ν*A*t 　　根据达西定律，v = k*i，则 　　V = k*（△h/L）*A*t 　　从而得出 　　k = V*L / A*△h*t 　　常水头试验适用于测定透水性大的沙性土的渗透参数。粘性土由于渗透系数很小，渗透水量很少，用这种试验不易准确测定，须改用变水头试验。 　　变水头试验法就是试验过程中水头差一直随时间而变化，其装置如图： ?? 变水头法测渗透系数 水从一根直立的带有刻度的玻璃管和U形管自下而上流经土样。试验时，将玻璃管充水至需要高度后，开动秒表，测记起始水头差△h1，经时间 t 后，再测记终了水头差△h2，通过建立瞬时达西定律，即可推出渗透系数 k 的表达式。 　　设试验过程中任意时刻 t 作用于两段的水头差为△h，经过时间dt后，管中水位下降dh，则dt时间内流入试样的水量为 　　dVe = －a dh 　　式中 a 为玻璃管断面积；右端的负号表示水量随△h的减少而增加。 　　根据达西定律，dt时间内流出试样的渗流量为： 　　dVo = k*i*A*dt = k*（△h/L）*A*dt 　　式中，A——试样断面积；L——试样长度。 　　根据水流连续原理， 应有dVe = dVo，即得到 　　k = （a*L/A*t）㏑（△h1/△h2） 　　或用常用对数表示，则上式可写为 　　k = 2.3*（a*L/A*t）㏒（△h1/△h2）