快速确定非饱和土水力参数的方法

专利名称：快速确定非饱和土水力参数的方法
技术领域：
本发明涉及一种非饱和土水力参数的测试方法，尤其涉及一种快速确定非饱 和土水力参数的方法。
背景技术：
非饱和土水力参数包括土水特征曲线和渗透函数。所述的土水特征曲线描述 的是毛细压力与含水量之间的关系；所述的渗透函数描述的是渗透系数与含水量 或毛细压力之间的关系。通过它们可以计算非饱和土的强度、体应变、污染物在 土壤中的运移规律等一系列问题。因此，非饱和土水力参数在解决非饱和土及渗 流有关的问题中扮演着重要的角色。
传统确定非饱和土水力参数的试验方法主要有三大类
第一大类，利用现有的试验装置(如Temper cell、压力板等)测定非饱 和土的土水特征曲线，再根据经典模型(Mualem、 Burdine等模型)来预测非饱 和土的渗透函数；
第二大类，采用两套试验装置分别测定非饱和土的土水特征曲线和渗透函
数；
第三大类，在一套试验装置(如Constant Flow Method)上同时测定非饱 和土的土水特征曲线和渗透函数。
虽然这些传统的试验方法各具特色，但是它们存在一个共同的缺陷一一要待 达到平衡后才能进行测量，从而导致测试时间很长。
传统确定非饱和土水力参数的试验方法的原理均是建立在试样处于平衡状 态下，测定非饱和土毛细压力与含水量之间的关系。然而，试样达到平衡状态需 要一个过程。这就意味获得非饱和土水力参数需要很长时间，这主要由试样的土 性及所要测的数据的个数决定。另外，由于传统确定非饱和土水力参数的试验方
4法需要很长时间，从而又会引发另一个不可避免的问题高压气体通过水的流动 在压力较低的高进气值陶土板的背面析出，影响测试精度。 综上所述，传统的试验方法有以下缺点
① 传统确定非饱和土的水力参数的试验方法是建立在试样处于平衡状态下 测定的，从而导致其测试时间很长；
② 高压气体通过水扩散，在压力较低的高进气值陶土板背面析出，影响测试 精度。
因此，如何縮短测试时间，从而达到快速确定非饱和土水力参数的目的就成 为需要解决的关键问题。

发明内容
本发明的目的就在于克服现有技术存在的上述缺点和不足，扩大现有技术的 使用范围，提供一种快速确定非饱和土水力参数的方法。 本发明的目的是这样实现的 本发明包括下列步骤
① 开展在一定条件下的多步流动试验，由试验得到饱和度随时间变化的关系 曲线；同时，建立一个能够准确描述该条件下试样饱和度随时间变化的线性粘弹 性模型；
所述的一定条件是指在[/,,、)时间内毛细压力R,保持不变；
② 从[H,)时间段内的饱和度与时间关系曲线中任取某两时刻^、 /;2及 其对应的饱和度&(C &"2)，通过建立的线性粘弹性模型计算毛细压力;^作 用下的单位压力容水率C,及特征时间^,;
③ 利用单位压力容水率^ 、特征时间Z"d与平衡状态下的饱和度和渗 透系数&间的关系，计算毛细压力A,作用下，试样处于平衡状态时的饱和度 (S,)—及渗透系数^;
④ 重复以上步骤②③，直到试样处于平衡状态时的饱和度等于残余饱和度为止。
5本发明的工作原理
开展一定条件下的多步流动试验，由试验得到试样饱和度随时间变化的关系 曲线。从该曲线的[， )时间内任取某两时刻/n 、 ,,2及其对应的饱和度&)、
，通过建立一定条件下的线粘弹性模型计算毛细压力A,作用下的单位压
力容水率C,及特征时间^;随后，利用单位压力容水率C。特征时间^与平衡
状态下的饱和度OS,、,和渗透系数A:,间的关系，计算毛细压力;^作用下，试样
处于平衡状态时的饱和度C)^及渗透系数^ 。
本发明具有下列优点和积极效果
1、 本发明是一种用一定条件下的多步流动试验结果来确定非饱和土水力参 数的方法。该方法无需等待试样达到平衡状态后再进行测定，从而大大縮短了测 试时间，有效地避免了高压气体通过水的流动在压力较低的高进气值陶土板背面 析出的问题。
2、 本发明是一种简单、有效，而且能在较短的时间内确定非饱和土水力参 数的方法，因而突破传统确定非饱和土水力参数方法存在的耗时费力等缺陷，提 高了确定非饱和土水力参数的速度和准确性。
3、 本发明适用于对地下污染物的环境评价、降雨入渗等与渗流和非饱和土 有关的相关领域。


图1是多步流动试验装置示意其中
l一试样；
2— 装样容器；
3— 高精度伺服气泵；
4— 高精度伺服流泵；
5— 高进气值陶土板；
6— Teflon渗透膜。
图2是毛细压力和饱和度随时间变化关系示意图;其中
r,—施加第i级毛细压力增量所对应的时刻； ~—在[H》时间内任意取的时刻；
A.,—在f,时刻的毛细压力，且A,=(凡-/0,在[/,，,,+1)保持不变；
A;^—在^时刻的毛细压力增量。
^(0—试样在时刻^的饱和度；
在[(，时间内^时刻所对应的饱和度； & "2) —在[Z,, (+1 )时间内f。时刻所对应的饱和度。 图3是C, A,关系曲线示意图； 其中-
C,一毛细压力为A,时对应的单位压力容水率。
图4由于压力变化引起水的流动示意图。 其中-
/7。一施加的气压力； 施加的水压力； A/^ —降低的水压力值； L一试样的长度。
具体实施例方式
下面结合附图和实施例对本发明进一步说明 1、关于步骤①
A、开展一定条件下的多步流动试验
采用如图1所示的多步流动试验装置，即将高精度伺服气泵3、装样容器2、 试样l、 Teflon渗透膜6、高进气值陶土板5和高精度伺服流泵4依次连接。随
后，施加如图2所示毛细压力随时间变化关系(在时刻A，将毛细压力A。突然
7增加一个小增量A&至毛细压力&，并保持其不变至^时刻；随后，在^时刻， 将毛细压力Pd突然增加一个小增量Ap^至毛细压力/^，并保持其不变至 3时 刻；依此类推，在^吋刻，将毛细压力A("一D突然增加一个小增量Ap。"至毛细压
力/^，并保持其不变至/ +1时刻)，得到试样饱和度随时间变化的关系曲线，如
图2所示。
B、 一定条件下线性粘弹性模型的建立
当处于平衡状态的非饱和土受到扰动时，瞬态的毛细压力增量可以表示为
& = - ;O = 。 - " + 警+*1 (1)
结合图2中毛细压力随时间变化关系及式(1)的线性特点，可以推得在
时间段内的任意时刻^，非饱和土的饱和度为
& (fy) = ) +1 Ape4 Ct [1 - exp(- ^^)](2 )
2、 关于步骤②
从测得的饱和度随时间变化关系的[K+i)时间段内任取两时刻^和^及其 对应的饱和度& (/n)和& ((2)。根据式(2 ) ， Sr (,n)和Sr (。2)与rc,和q值之间存 在以下对应关系
WU = "O +1- eXp(—(3)
& (,, 2) = Sr (0 + ^ [1 - exp(- (4 )
从^时刻开始，利用式(3)、式(4)可以求得在毛细压力/^作用下的、,和C,。
3、 关于步骤③
A、计算在平衡状态下毛细压力为A,时所对应的饱和度(S,)^ 根据计算得到毛细压力a,下的单位压力容水率c,.，建立c, a,关系曲线，如图3所示。根据单位压力容水率c,的定义及施加较小的毛细压力增量，可以 得到试样处于平衡态下毛细压力为a,时所对应的饱和度(&)^，它可以用图3 中的带斜杠的矩形面积加上初始时刻的饱和度来表示。即
(SW" +红A;^ (5)
由于&("、 c,、 Aa,.均己知，因此，在平衡状态下毛细压力为a,时所对应
的饱和度(&)^可以通过式(5)求出。
B、计算在平衡状态下毛细压力为p 时所对应的渗透系数^ 如图4所示，在初始时刻f。试样两端施加的压力使得孔隙中的水保持平衡。 随后，在^时刻，试样右边的水压力突然降低A^，此时孔隙中的水发生流动，
并从试样的右端流出。在此过程中，孔隙中水的流动可以用达西定律来描述；此 外，孔隙中的水还必须满足质量守恒定律。将流动方程和质量守恒方程相结合， 可以得到特征时间与渗透系数之间的关系式，如下所示-
由于试样孔隙率"、试样长度丄、水的重度L、单位压力容水率G及特征 时间t,均已知，因此，在平衡状态下毛细压力为A,时所对应的渗透系数&可以 通过式(6)求出。
9
权利要求
1、一种快速确定非饱和土水力参数的方法，其特征在于包括下列步骤①开展在一定条件下的多步流动试验，由试验得到饱和度随时间变化的关系曲线；同时，建立一个能够准确描述该条件下试样饱和度随时间变化的线性粘弹性模型；所述的一定条件是指在[ti，ti+1)时间内毛细压力pci保持不变；②从[ti，ti+1)时间段内的饱和度与时间关系曲线中任取某两时刻ti1、ti2及其对应的饱和度Sr(ti1)、Sr(ti2)，通过建立的线性粘弹性模型计算毛细压力pci作用下的单位压力容水率Ci及特征时间τci；③利用单位压力容水率Ci、特征时间τci与平衡状态下的饱和度(Sr)eq，i和渗透系数ki间的关系，计算毛细压力pci作用下，试样处于平衡状态时的饱和度(Sr)eq，i及渗透系数ki；④重复以上步骤②③，直到试样处于平衡状态时的饱和度等于残余饱和度为止。
2、 按权利要求1所述的一种快速确定非饱和土水力参数的方法，其特征在 于步骤①A、 开展一定条件下的多步流动试验采用多步流动试验装置，即将高精度伺服气泵(3)、装样容器(2)、试样(1)、 Teflon渗透膜(6)、高进气值陶土板(5)和高精度伺服流泵(4)依次连接； 随后，施加一定条件的毛细压力，得到饱和度随时间变化的关系曲线；b、 一定条件下线性粘弹性模型的建立当处于平衡状态的非饱和土受到扰动时，瞬态的毛细压力增量可以表示为r * 必r *电=Wpa - a ) = W a —化、+甘s—^+甘s, (1)在[H。时间内的任意时刻^.，非饱和土的饱和度为& (~) = ^) +1 、& [1 - exp(- (2 )"1 、
3、 按权利要求1所述的一种快速确定非饱和土水力参数的方法，其特征在于步骤②在[H。时间内任取两时刻^、 ^及其对应的饱和度^(q)、 S,(/,2)，根据式(2)可以计算得到毛细压力为A,时所对应的^,和C,值。
4、 按权利要求1所述的一种快速确定非饱和土水力参数的方法，其特征在于步骤③计算在平衡状态下毛细压力为A,时所对应的饱和度(&)^(&H(0 +红A;^ (3) 计算在平衡状态下毛细压力为/^时所对应的渗透系数A:,(4)
全文摘要
本发明公开了一种快速确定非饱和土水力参数的方法，涉及一种非饱和土水力参数的测试方法。本发明主要包括下列步骤①开展在一定条件下的多步流动试验，由试验得到饱和度随时间变化的关系曲线；同时，建立一个能够准确描述该条件下试样饱和度随时间变化的线性粘弹性模型；②计算毛细压力p_ci作用下的单位压力容水率C_i及特征时间τ_ci；③计算毛细压力p_ci作用下，试样处于平衡状态时的饱和度(S_r)_eq，i及渗透系数k_i。本发明是一种用一定条件下的多步流动试验结果来确定非饱和土水力参数的方法，该方法简单、有效，而且能在较短的时间内确定非饱和土水力参数。本发明也适用于对地下污染物的环境评价、降雨入渗等与渗流和非饱和土有关的相关领域。
文档编号G01N15/08GK101487787SQ20091006073
公开日2009年7月22日 申请日期2009年2月12日 优先权日2009年2月12日
发明者辉 陈, 韦昌富, 魏厚振 申请人:中国科学院武汉岩土力学研究所