土石围堰变饱和非稳定渗流数值模拟分析.pdfVIP

水利水 电工程设计 DWRHE ·2012年 第31卷 第4期 土石围堰变饱和非稳定渗流数值模拟分析 王朝江 刘 国强 杨炳炎 摘 要 根据二维饱和一非饱和水流运动方程和水分特征曲线VG模型。应用Geostudio岩土软件中的渗流分析软件 SEEP／W对水利工程常见的土石围堰进行变饱和非稳定渗流分析。通过分析模拟结果，给出了不同施工阶段的渗流量以 及围堰 内压力水头分布及流网图，并且分析 了布置混凝土防渗墙的效果。 关键词 土石围堰 渗流分析 VG模型 SEEP／W软件 有限元模型 中图分类号 TV139．1 文献标识码 B 文章编号 1007—6980(2012)04-0016—03 在水利工程建设中，施工围堰是维护水工建 可供选择的解析模型来表示土壤水分特征。其中 筑物的施工场地，使其免受河道水流影响的临时 比较常用的为VG模型 (VanGenuchten，1980)： 水工建筑物。围堰在结构上要求稳定 ，并且具有 IOr+ = ，h0 一 定的防渗性。围堰的防渗效果关系到工程施工 )={ (1+l~hl) 进度 以及工程的安全，一旦出现问题将拖延工期 l ，≥o 或者造成严重事故。结合某工程技施阶段土石围 堰渗流分析，探讨基于SEEP／W软件 VG模型的土 K()： (1-(1-So~)) 石围堰变饱和非恒定渗流数值模拟分析问题。 So=酱 1 数学模 型 m=l-1／n，nl 1．1 控制方程 式中 — 篼 含水率(即最大分子持水率)； 模型中水流控制方程采用 Richards方程 ，其 — — 饱和体积含水率； 形式为： 磁——饱和水力传导度，m／d； Oh ．s厂 饱和度 ； = Z——孔隙连通系数，一般取为 0．5： 式中 一 岩土体积含水率： 和 为经验常数 。 — — 为压力水头，131； 根据该工程地层情况及相关地质参数 ，选取 ． 一 汇，d～，反映岩土根系对水分的吸收； 的岩土水分参数如表 l所示，基于VG模型的土壤 r (1，2)空间坐标，m； 水分特征曲线如图1所示。 厂一 (=1，2)空间坐标 ，m； z——垂 向坐标，ITI，向上为正； — — 时间，d： A 、E：——无量纲各向异性张量 的分量： ? — — 非饱和水力传导函数，m／d，并由下式 吕 一 给出。 霸 K( )=K )Kr ，y,z) 璐 式中 K 一 相对水力传导系数； K —饱和水力传导系数，m／d。 1_2水分运动特征参数 含水率 ()和土壤水力传导度K(h)与基质吸 力h具有很强的非线性关系．SEEP／W软件有 5种