1楼:匿名用户
我用过 equi-pf自动土壤水分特征曲线测量仪equi-pf自动土壤水分特征曲线测量仪是用于获得土壤水分含量与土壤水势之间的关系(0-100cm水势)pf曲线, 样品一放入仪器,系统自动进行所有的测量过程
2楼:匿名用户
如果用负压计(俗称陶土头),测定范围比较窄。
3楼:匿名用户
equi-pf自动土壤水分特征曲线测量仪,好贵的。测量范围也只能到100cmhuamanlou31,您好,您在南京吗。equi-pf自动土壤水分特征曲线测量仪你会用,到时候请教您,刚买回来
4楼:匿名用户
我们现在用张力计测,测量范围是比较窄,而且好麻烦
土壤水分特征曲线的测定
5楼:中地数媒
土壤水分特征曲线的测定在室内采用张力计称重法,用张力计(负压计)测定土壤负压h,用称重法测定相应的含水率θ,试验装置如图2.3.1所示。
通过试验获得了主脱湿过程的实验数据,采用van genuchten(van genuchten,1980)模型来描述主脱湿曲线(mdc)(沈荣开1993),模型如下:
图2.3.1 试验装置示意图
土壤水盐运移数值模拟
式中:s为饱和度(表示孔隙被水充满的程度,等于水的体积与孔隙体积之比,cm3/cm3);θ为含水率(cm3/cm3);θr为残留含水率(cm3/cm3);θs为饱和含水率(cm3/cm3);h(hh2o)为负压(cm);α,n,m表示土壤水分特征曲线形状的参数。
van genuchten模型含有四个参数(即α,n,m(含 n),θr,θs),所以,计算较为复杂,但一般情况下,θr和θs可由室内外试验给出,这样,模型中只剩α和 n 两个参数(其中 m 可利用 m=1-1/n求得),为求出这两个参数,一般根据最小二乘原理,用实验数据拟合的方法确定(沈荣开,1987)。
1.线性迭代法(一参数迭代)
公式(2.3.1)可变换为:
土壤水盐运移数值模拟
由于负压h的绝对值为一正值,所以可用吸力代替(这里仍取为h,推导时去掉了绝对值符号)。将上式两边取对数得:
土壤水盐运移数值模拟
令 ,b0=nlnα,b1=n,x=lnh,则上式变为一元回归模型:
y=b0+ b1x (2.3.4)
因此,可用求解一元回归方程的方法确定b0、b1,进而求出α、n:
,n=b1,从而可得
。具体计算时需要用迭代的方法来求解。首先给出初值m(0),并将实测数据含水率θ、负压h代入求解回归方程的系数b0、b1,从而可求得第一次迭代值m(1),再将m(1)代入,得出第二次迭代值m(2),……,依次迭代,直到第p+1次迭代与第p次迭代值之差的绝对值小于预先给定的常数e(e为一充分小的正数)为止。
收敛标准(迭代控制)用公式表示如下:
土壤水盐运移数值模拟
满足收敛标准时,由第p+1次迭代求出的回归系数b0、b1,即可确定出参数α、n。
2.非线性迭代法(二参数迭代)
由van genuchten模型变形为
θ-θr=(θs-θr)[1 +(αh)n]-m (2.3.6)
即(θr-θ)+(θs-θr)[1 +(αh)n]-m=0 (2.3.7)
将实测数据:含水率θi、负压hi(i=1,2,…,n,n为观测点个数)代入上式得:
(θr-θi)+(θs-θr)[1 +(αhi)n]-m=εi(i=1,2,…,n)(2.3.8)
由最小二乘原理
土壤水盐运移数值模拟
令 z=e(α,m,n),z为α、n的二元函数,m为中间变量。求多元函数z的极值:
土壤水盐运移数值模拟
令土壤水盐运移数值模拟
其中w、x为α,n的二元函数。
为了简化推导,求出w、x,令
θs-θr=θ1
θi-θr=θ2
(αhi)n=xi
推导时略去角标i,则
土壤水盐运移数值模拟
先求w:
土壤水盐运移数值模拟
因x=(αh)n,
,所以上式为
土壤水盐运移数值模拟
再求x:
土壤水盐运移数值模拟
所以土壤水盐运移数值模拟
式(2.3.11)为一组非线性方程,非线性方程的求根可运用牛顿迭代法。
牛顿迭代法其基本思想是:将非线性方程逐步归结为某种线性方程来求解。其几何解释为,方程的根,用其切线方程的根来逼近,由于这种几何背景,牛顿法亦称切线法。
由牛顿法迭代公式(李庆扬、王能超、易大义,1991.9):
土壤水盐运移数值模拟
改写为如下形式
δxf′(x)=-f(x) (2.3.20)
对于多元函数,牛顿迭代公式可扩展为:
f′x(x,y)δx+f′y(x,y)δy=-f(x,y) (2.3.21)
则式(2.3.11)的牛顿迭代公式表示为:
土壤水盐运移数值模拟
也即土壤水盐运移数值模拟
下面的任务就是求方程式(2.3.23)左端各偏导数项,根据多元复合函数的求导法则,首先求第一个方程各偏导数项。
土壤水盐运移数值模拟
考虑到,则土壤水盐运移数值模拟
所以土壤水盐运移数值模拟
土壤水盐运移数值模拟
因为nln(αh)=lnx,所以上式
土壤水盐运移数值模拟
土壤水盐运移数值模拟
因为 ,
,所以土壤水盐运移数值模拟
则土壤水盐运移数值模拟
然后求式(2.3.23)第二个方程各偏导数项。根据二阶混合偏导数在连续的条件下与求导的次序无关,则
土壤水盐运移数值模拟
而土壤水盐运移数值模拟
式(2.3.32)右端共有4项,需求4项偏导数。
第1项偏导数
土壤水盐运移数值模拟
因 ,所以上式
土壤水盐运移数值模拟
第2项偏导数
土壤水盐运移数值模拟
第3项偏导数
土壤水盐运移数值模拟
将 代入得
土壤水盐运移数值模拟
所以土壤水盐运移数值模拟
第4项偏导数
土壤水盐运移数值模拟
综合上述4项偏导数得
土壤水盐运移数值模拟
令土壤水盐运移数值模拟
则式(2.3.40)简记为:
土壤水盐运移数值模拟
所有偏导数求出后,解方程组
土壤水盐运移数值模拟
得土壤水盐运移数值模拟
迭代公式为
土壤水盐运移数值模拟
式中:p为迭代次数。
具体计算时,首先给出参数的迭代初值α(0),n(0),并由实测数据计算出相应的偏导数值,然后按照式(2.3.45)依次进行迭代。迭代控制标准为:
土壤水盐运移数值模拟
式中:e为给定的充分小的正数。满足上述迭代标准的α(p+1)和n(p+1)就是所求的参数α和n的值。
根据牛顿迭代法的局部收敛性,一般的说,牛顿法的收敛性依赖于初值的选择,如果初值偏离所求的根比较远,则牛顿法可能发散。为保证牛顿法的收敛性,首先用线性迭代法对参数进行估计,然后将线性迭代法估计的参数作为非线性牛顿迭代法的初值进行迭代。
3.土壤水分特征曲线测定结果
1998年10月于长江河口地区采集了 3 个土样,寅阳 1#(lxg-1)粉砂壤土(1998.10.23),大兴2#(stg-2)粉砂壤土(1998.
10.24),兴隆沙1#(xls-1)粉质粘壤土(1998.10.
25),取样时均去除表土20cm。所取土样的机械组成见表2.3.
1。土壤水分特征曲线的测定在室内采用张力计称重法,为了保证装土的初始含水率均匀、密度一致,土样经过粉碎和过筛(20目)处理,然后按一定的干容重装填土样,装好的土样经过充分饱和后开始脱湿试验,试验于1999年4月16日开始,5月15日结束,脱湿过程的实验数据见表2.3.
2。表2.3.1 土样颗粒分析结果(美国制)
表2.3.2 水分特征曲线实测数据
续表根据实测数据,通过vb编程计算,运用线性迭代法得到的上述三个土样,寅阳1#(lxg-1)粉砂壤土,大兴2#(stg-2)粉砂壤土,兴隆沙1#(xls-1)粉质粘壤土的van genuchten模型参数见表2.3.3。
其拟合曲线见图2.3.2。
表2.3.3 水分特征曲线参数
图2.3.2 水分特征曲线
将这些参数代入 van genuchten模型,即可根据不同的负压值来计算其相应的含水率
土壤水盐运移数值模拟
对上式求导可得到容水度c(h)
土壤水盐运移数值模拟
如果已知饱和水力传导度ks,还可得到mualem模型(mualem,1984)的非饱和水力传导度
土壤水盐运移数值模拟
根据长江河口地区土壤水分特征曲线的实测数据,选择vg(van genuchten)模型,用一个函数较好地描述了脱湿过程,比起用分段函数来描述,具有明显的优越性。通过计算表明,线性迭代法简单实用,同时也具有相当的精度,基本可以满足实际需要。本次试验由于没有电子天平,而使用普通天平(感量2g),因而给含水量的观测带来一定的误差,但通过实测值与计算值的比较含水率最大绝对误差小于2%。
6楼:况媚裴桐华
实验原理
张力计插入土样后,张力计中的纯自由水经过陶土壁与土壤水建立了水力联系。在非饱和土壤中,仪器中的自由水的势值总是高于土壤水的势值,因此,仪器中的自由水就会透过陶土管进入土壤,但因陶土材料孔隙细小,孔隙中形成的水膜不能使空气通过,而只能让水或溶质液通过(但如果压力过高水膜破裂,空气就会透过,这时的压力称为透气值),因而在仪器内形成一定的真空度,由仪器上的负压表读出。最后当仪器内外的势值趋于平衡时,仪器中水的总水势φwd与土壤中土水势φws应该相等,即:
φwd=φws
土水势的完整表述为:
φ=φm+φp+φs+φg+φt
因为陶土管为多孔透水材料,并非半透膜,故溶质也能通过,最后达到内外溶液浓度相等,内外溶质势φs相等。仪器内外温度相等,温度势φt相等。坐标0点选在陶土头中心,则陶土头中心的内外重力势φg相等。
这样仪器中和土壤中的总势平衡可表述为:
φmd+φpd=φms+φps
式中,φps为土壤水的压力势,φms为土壤水的基质势,φpd为仪器内自由水的压力势,φmd为仪器内自由水的基质势。
在非饱和土壤中,土壤水所受的压力为大气压(基准状态),故φps应为零,又仪器中自由水无基质势存在,故φmd亦为零,所以:
φms=φpd=δpd+z
式中,δpd为负压表显示的负压值(小于0),z为埋藏在土中的陶土管中心与土面以上负压表之间的静水压力即水柱高,(向上为正,大于0)。即可得到土壤水的基质势。按定义土壤水吸力为基质势的负值,因而即可测得吸力值。
s=-φms=-δpd-z
如果负压表读数记为p(大于0,即p=-δpd),则s=p-z
另外,在计算土样中水分的变化时,还应考虑集气管中水分的变化量。