资源描述：
《非饱和土的渗透系数_孙大松.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、DOI:10.13243/j.cnki.slxb.2004.03.012水利学报2004年3月SHUILIXUEBAO第3期文章编号:0559-9350(2004)03-0071-05非饱和土的渗透系数1122孙大松,刘鹏,夏小和,王建华(1.江苏省交通规划设计院,江苏南京210000;2.上海交通大学建筑工程与力学学院,上海200030)摘要:非饱和土的渗透系数是分析水分和物质迁移的重要参数,直接测量非饱和土渗透系数的代价较高,且直接测量的精度较差,因此间接估算非饱和土渗透系数成为很好的选择。分形理论适合用来描述多孔介质的结构和透水性。本文建立了土体孔隙分布
2、的分形模型,导出用分维和进气值表示的水分特征曲线和渗透系数的理论表达式。与实验结果的比较表明,用分形模型计算得到的水分特征曲线和渗透系数与试验结果一致。关键词:渗透系数;非饱和土;分形;分维;水分特征曲线;孔隙分布;吸力中图分类号:P642.11文献标识码:A非饱和土的渗透系数是石油工程和环境工程中分析水、油和污染物迁移必不可少的重要参数。非饱和土渗透系数的试验测量结果表明,不同土样在不同的含水量时的渗透系数差别较大。即使在近饱[1]和状态下(吸力介于0～1kPa)吸力的微小变化也能引起渗透系数改变1～3个量级。由于非饱和土渗透系数对饱和度的依赖性大,造成测量
3、的困难,特别是在低饱和度时,非饱和土渗透系数测量很费时间。一般情况下,用实验方法测量非饱和土的渗透系数不太现实。经验公式和理论预测成为确定非饱和土渗透系数的不可回避的选择。土体中水分运移与其孔隙结构密切相关,土体的孔隙是各向等性[2,3]的,且可以视为均匀分布。分形理论很适合表示土体的结构性状。徐永福和孙德安用分形理论表[2]示了土体的孔隙分布,并建立了孔隙分布分维与非饱和土体导水系数间的相互关系。本文根据土体孔隙分布的分形模型,导出了非饱和土的水分特征曲线和非饱和土渗透系数的理论表达式,并将理论表达式与已有的非饱和土的水分特征曲线和非饱和土渗透系数进行比较,
4、验证本文提出的理论公式。1土体孔隙分布的分形模型[2]分形可以用来表示非饱和土的孔隙结构。分形孔隙结构可以通过以下方法生成:假定一个立方3块体的边长为单位长度1,被分成1/(r)个边长为r的次一级的小立方块体,从中拿走N个次一级的332小立方块;将剩下的1(r)-N个次一级小立方块再分成1(r)个边长为r的次二级立方块体,从中3拿走N[1(r)-N]个次二级的小立方块,如此重复下去,直到孔隙的孔径为无穷小,这样就形成一个分[3]形分布的孔隙结构。根据分形理论,孔隙的孔径与其对应的个数有下面的关系-DN=Cr(1)式中:C为常数;r为孔径;D为孔隙分布的分维。r
5、3孔隙的体积VP由VP=∫(4πr3)dN算出,即03-DVp=Ar(2)收稿日期:2002-11-22基金项目:国家自然科学基金资助项目(40201024)作者简介:孙大松(1967-),男,江苏人,硕士,高级工程师,主要从事土力学、公路路基处理和设计等工作。—71—式中:A=4πCD3(3-D)。土体的孔隙体积与孔径间相互关系可以用水银压入法测出来,如果孔隙体积与孔径在双对数坐标上回归线的斜率是κ,那么土体孔隙分布的斜率为D=3-κ(3)为了与土体孔隙体积联系起来,非饱和土的含水量采用体积含水量。体积含水量是土体中孔隙水的体积与土体体积的比值,土体中孔隙水
6、的体积可以用土体孔隙体积表示。根据式(2),介于孔径r※r+dr中孔隙水体积可以表示为24πrdΛ=dN(4)3式中:Λ=θ-θr,是相对体积含水量;θ和θr分别是体积含水量和残余体积含水量。残余体积含水量是土体中被土粒吸附的水分占孔隙水的百分含量。由式(4)可以得到非饱和土的相对含水量和饱和土的相对含水量,分别表示为3-D3-DΛ=Ar;Λs=AR(5)式中:R是土体中孔隙的最大孔径。非饱和土的水分分布往往是先充填在孔径小的孔隙中,随着饱和度增大,水分渐渐地向大孔隙中移动,当最大孔径的孔隙中充满水时,土体基本饱和,此时的吸力为非饱和土的进气值,非饱和土的进气
7、值是土体孔隙中开始出现气泡时的吸力,与土体孔隙的最大半径有关。根据Young_Laplace公式得到非饱和土的吸力和进气值与孔径的关系2σcosα2σcosαψ=;ψe=(6)rR式中:ψ和ψe分别是非饱和土的吸力和进气值;σ是表面张力;α是接触角。将式(6)代入式(5),得到非饱和土的水分特征曲线为δψSe=;δ=D-3(7)ψe式中:Se=(θ-θr)/(θs-θr),是非饱和土的有效饱和度,有效饱和度表示非饱和土孔隙中能流动的水分与饱和土孔隙中能流动的水分的比值;θs为饱和土的体积含水量,其值等于土体的孔隙率。式(7)是由土体孔隙分布的分形模型导出的水分
8、特征曲线的表达式,参数δ可以由孔隙分布
