热脉冲法测定土壤热性质的研究进展

广州市深华生物技术有限公司热脉冲法测定土壤热性质的研究进展本文综述了土壤热性质的计算模型及研究现状，重点针对近年来国内外研究土壤热性质的新方法热脉冲法的理论和技术发展，及其在土壤水和其他物理性质应用方面的进展。土壤热性质是决定土壤热状况的内在因素，研究土壤热和温度的变化规律以及调节土壤热状况时必须首先了解土壤的热性质。已有研究表明，导热过程受土壤水和其他物理化学特性的影响[1，2]，土壤热性质与土壤水分状况之间存在明确的定量关系。因此研究土壤热性质还有助于从土壤热运动规律方面获取土壤水分信息。目前已有不少研究者对土壤水和土壤热性质的关系进行了探讨，并找出了一些规律，这对于土壤热性质的深入研究及其在土壤水分管理的应用方面无疑具有重要作用。本文就近年来国内外土壤热性质研究的新方法热脉冲方法在理论和实验上的发展及应用进展进行介绍，并与常规方法进行对比，以期对其推广和进一步研究起到参考作用。1土壤热性质参数及其模型11反映土壤热性质的相关参数不同土壤吸收一定热量后，其温度增减的幅度不同，即各种土壤贮热和导热能力不同，这是因为土壤的热性质不同所致。土壤热性质指标主要有土壤热容量、土壤导热系数、土壤热扩散率等。单位体积土壤的热容量Cv可用下式计算:Cv=XsCsXwCwXaCa(1)式中，Xs、Xw和Xa分别是土壤中固体物质、水和气体的体积;Cs、Cw和Ca分别是它们的比热容。土壤导热率是在标准条件下通过土壤传导热量的量度，为各组分导热率的加权平均。导热率与体积热容量之比即为土壤热扩散率:Dq=/Cv(2)式中，Cv为土壤体积热容，Jm-3K-1;为土壤导热率，Jm-1K-1s-1;Dq为土壤热扩散率，m2s-1。1.2土壤导热率的计算模型土壤热性质与土壤中各相的组成及比例有关。在上述三个土壤热性质中，土壤热容量的计算较简单，可根据土壤总孔隙度及土壤中不同组成成分的热容量求得[1];而土壤热扩散率的计算基于式(2)求出，因此获取土壤热性质的方法便集中在了导热率的计算模式上。导热率计算模型可分为物理模型和经广州市深华生物技术有限公司验模型两类。其中物理模型以deVries(1963)的最具代表性，该模型以临界含水量为界，不同的土壤含水量范围有不同的表达形式。在此，临界含水量是当土壤中液态水失去连续性时的含水量，也有研究者把压力势为-55kPa时的含水量作为临界含水量计算值[1]。若以k表示临界含水量，当k时，导热率由下式计算[1]:=!ni=1(kixii)/!ni=1kixi(3)式中，xi为组成成分i的体积比例;i为组分i的导热率;ki表示颗粒组分i的权重系数，它与颗粒形状和接触角以及各成分导热率有关;模型中考虑了5种组分(n=5)，包括:液态土壤水，湿润土壤空气，石英，其他土壤矿物和土壤有机质。1.3土壤热性质常规研究方法的研究进展国内外对热性质的研究已表明，土壤热性质受土壤水分状况的影响。Parikh(1969)用非稳定方法测定了250m的玻璃珠和粉壤土的导热率和扩散率随含水量变化特征[3]。Wierenga等(1969)分析了Yolo粉壤土的热性质，发现表观导热率为土壤含水量的函数，且用deVries模型进行计算的结果与测定值较吻合[4]。Ghuman等(1985)研究发现，土壤导热率和热扩散率都随质地和初始含水量的改变而变化，含水量增大则导热率也增大，而且在不同含水量范围内，粘性土壤的导热率比沙土的低[5]。Persaud和Chang(1985)根据两个深度的温度值计算了表观导热率，并采用四种不同方法进行了对比，结果表明不同方法计算的导热率是有差异的[6]。Kaune等(1993)测定了扰动的结构性黄土的温度，研究了团聚体对土壤热性质的影响[7]。除含水量之外，土壤中含盐水平、容重及有机质含量也影响导热性质。早在vanRooyen和Winterkorn(1959)的研究中，浓度为0.18molkg-1的CaCl2溶液，或是浓度达到0.34molkg-1的NaCl溶液对石英的导热率并无明显影响[8]。但Globus和Rozenshtok(1989)对0.25molkg-1的KOH湿润过的石英进行导热率测定，结果表明其导热率比水湿润过的石英砂的低[9]。Noborio和McInnes(1993)发现从0.1molkg-1到溶解度范围内，土壤表观导热率随CaCl2、MgCl2、NaCl、Na2SO4等盐分浓度的增加而降低[10]。在运用不同模型及其与实测值的对比方面，Bachmann等(2001)将斥水土壤与吸水土壤做对比，利用模型计算出的结果和实测值进行比较发现，吸水土壤的导热率比斥水土壤的大得多，实测法与计算模型等不同方法得出的导热率差异很大，其中deVries模型对吸水土壤导水率的计算值比实测值低0.5Wm-1K-1以上;对斥水土壤的计算值也比实测值小，而Campbell模型计算的导热率在低饱和度时偏小。此外，两模型计算干土或者饱和含水量下的导热率值均很精确[11，12]。广州市深华生物技术有限公司贺康宁等(2000)采用自记土壤温度计连续测定土温，并根据温度观测值计算土壤热性质，给出了光滑坡面、光滑地衣坡面及自然坡面等不同坡面处理的土壤热性质平均值[13]。李毅等(2003)在研制了热性质实验系统的基础上，根据实测土壤温度资料，用非稳态法测定土壤温度，并采用有限差分法离散热传导方程来计算热扩散率[14]:Dq=(Tk1，j-2Tk，jTk-1，j)!t(!x)2(Tk，j1-Tk，j)(6)式中，T为温度;x为距离坐标;t为时间坐标;下标k和j分别代表差分离散后的位置点和时间点。李毅等根据土壤孔隙分布确定土壤热容量，得到了不同质地土壤的导热率。同时还将导热率与土壤水吸力和土壤盐分浓度建立联系，研究了土壤水、盐、热性质的内在关系。2热脉冲法测定土壤热性质的理论基础及其应用2.1测定原理土壤热性质的测量有不同的方法，传统方法是在观测位置上设置输入热源，根据温度的升降直接测定。近年来研究者提出了土壤热性质测定的新方法热脉冲方法。该方法自Byrne等(1967)提出后，最初多应用于矿业，引入土壤研究中仅10多年，但在国外已引起重视并进行了一系列相关研究[15~24]。该方法成本低，对土壤扰动小，测试时间短，不易引起非饱和土壤中水分的重新分布，而且所需样品体积小，能够在田间进行自动连续定位测量，并且不断有研究证实其测量精度较高[20~24]，因而在国外得到了大量应用。自热脉冲方法提出以来，其技术方面在不断改进和更新，理论研究也逐渐深入，并且在土壤水分测定方面得到了较多应用。由于经济力量和技术等各方面原因，热脉冲技术在我国的应用还很有限。根据热传导定律，在一个无限大的均匀等温介质中，线性热源发射的热脉冲呈放射状向周围传导。使用热脉冲方法测量时，最简单的是用单探头方法进行测量。对于初始温度场均匀的介质，基于热传导方程，可将温度表示为时间的函数。假设温度与时间的对数值呈线性关系，则对两个时刻温度可得出下式[15]:=(q/4)ln[(t2/t1)(t1-th)/(t2-th)]/[T(t2)-T(t1)](7)式中，T为温度;q为热源的强度，定义为q=q/Cv，其中q是单位长度加热丝在单位时间内释放的热量;th为传感器冷却时间;t1、t2分别为测定的两个时刻。实际测定和应用中采用更多的是双探头热脉冲法。其测量设施上装有两个距离为r的平行不锈钢探针，其中一个探针含有线性加热源，另一个装有温度测量元件(如传感器或热电偶)。将双探头设备插入土壤时，通电后加热探头产生热脉冲，而另一探头可记录温度的时间变化。这些观测资料可直接用于确定包括导热率在内的热性质参数。对于土壤中的某一点，其温度随时间的变化表示为[11，广州市深华生物技术有限公司16]:!T(t，r)=q4CvDqEi-r24Dq(t-t0)-Ei-r24Dqttt0(8)式中，!T(t，r)为温度变化量，t为加热设备开启后的时间，t0为开始发射热脉冲的时刻，r为线性热源的径向距;Cv和Dq分别是介质的体积热容和热扩散率，-Ei(-x)为指数积分。用热脉冲装置实际测出的是时间增加量，热容量和热扩散率可根据上式用非线性回归方法求解，热容量和热扩散率的乘积即为导热率。热时域反射仪(ThermoTDR)，简称热TDR，是采用热脉冲方法进行土壤热和其他物理性质测量的一种新仪器，该装置结合TDR和热脉冲方法为一体，采用三个探针等间距平行排列，通常中间的探针发射热脉冲，而两边的探针(热源上、下游)监测温度随时间变化过程[17~21]。2.2理论研究进展由于热脉冲方法提供了测定温度、热性质、土壤含水量等参数的改进方法，因此自该法提出以来，在土壤物理研究中很快得到了认可和应用，尤其在近几年其理论和技术上的发展非常迅速。Byrne等(1967)首先提出用热作为示踪测定水分通量，他们使用的装置中装有温度传感器和点(或线)热源，水分通量根据仪器附近热场的变化确定。但该方法测定水分通量存在几个局限性，因而限制了仪器的应用推广。局限性之一是它在达到热平衡之前需要较长时间的稳定热输入(在平均流速下需30min)，而非饱和土壤中存在热梯度将导致土壤水分再分布;其二是仪器的校正需要把水分通量和仪器的响应方式关联起来;此外，仪器尺寸偏大改变了仪器周围土壤水流区域。而且其实验结果也表明，点热源方法测定的结果与理论值相差较大[22]。为此，不少学者致力于改进Byrne等(1967)研究中存在的问题。Ren等(2000)采用热TDR技术提出了测定土壤水通量密度的新方法，基于热迁移理论推导出了最大无量纲温度差(MDTD)的表达式:MDTD=#tmtm-t0S-1exp-(xd-Vs)24Cvs-exp-(xuVs)24Cvsds(11)式中，tm是无量纲温度差达最大值的时间，V为热脉冲速度或热锋对流速度，x为位置，下标u代表热源探头的上游，下标d指下游，s为时间差。Ren等采用热TDR方法测定MDTD，利用上述关系确定V，根据(12)式确定土壤水分通量[23]:V=VwCvwCv=JCvwCv(12)式中，J为土壤水分通量，Cvw为土壤溶液的体积热容。若测定了Cvw和Cv，则可进一步计算水分通量，Ren等的研究中探头距离xd和xu经人工调整，可使Cvw实测值为4.18MJm-3-1。Ren等对沙土、砂壤和黏壤等不同质地土壤测定MDTD，结果表明对上述土壤MDTD与通量之间的关系均为近似线性关系，砂土实测值和计算值较接近，而砂壤和粘壤则误差较大。此外，由于热TDR可同时测定土壤含水量，因此在获得土壤水分通量J后，可进一步得到孔隙水流速值(J/)。Ren等的方法克服了前人测定土壤水分通量时不能同时测定热性质的缺陷，但关于该方法在不同土壤上的适用性方面，今后还需更多的工作去检验;此外，式(11)的积分形式使得MDTD的计算偏于复杂而不便于应用。广州市深华生物技术有限公司为了解决Ren等(2000)的方法中表达式过于复杂的问题，Kluitenberg(2001)引入渗漏含水层井函数的概念，并将该函数中的指数部分进行麦克劳林级数展开，得到了形式上较为简单的近似解，从而仅用简单函数就表达了Ren等的解析解[24]:W(u，#)=!nm=0(-1)mm!#24umEm1(u)(13)Em1(u)=1m[exp(-u)-uEm(u)]m=1，2，%(14)式中，Em1(u)为井函数。但Kluitenberg的方法在实际计算中也是不易应用的。为寻求更简便的方法，Wang等(2002)对热脉冲测定土壤水分通量进行了理论分析，找出了线性热源上、下游温度增加比的自然对数与土壤水分通量之间的简单线性关系，由于其表达式在形式上较简单，因此在土壤水分通量的测定方面，更适于应用热脉冲传感器方法。热流方程写为[25]:Tt=Dq2Tx22Ty2-VTx(15)式中V可根据式(12)确定。Wang等基于热流方程得出的上、下游温度增加比与时间的关系为:上游:V2=4Dqt0lntu-t0tux2u(tu-t0)tu(16)下游:V2=4Dqt0lntd