冻土中水热耦合运移的理论与实验研究.docxVIP

冻土中水热耦合运移的理论与实验研究 土壤冻结是一种常见的自然现象。根据冻土的生长期，将冻土分为三类：多年冻土；冬小麦甜土；立即冷冻的土壤溶液。我国是世界第三大冻土国,多年冻土分布面积约占世界多年冻土分布的10%,占我国国土面积的21.5%,季节性冻土分布面积约占国土面积的53.5%。由于冻土区有丰富的土地、森林和矿物资源,它的存在及其演变对人类的生存环境、生产活动方式和可持续发展具有重要影响。 土壤的冻结是一个非常复杂的过程,它伴随着许多物理、化学的现象和过程,其中最主要的是三种现象:一是水分的迁移及水的相变(冰晶的形成);二是热量的传输;三是溶质含量的变化。这三者的变化并不是独立的,而是相互制约、相互影响的,是耦合变化的过程。用于描述冻结过程中土壤水热耦合运移的模型主要有两类:基于多孔介质中液态水分的粘性流动、热平衡原理的机理模型和应用不可逆热力学原理描述土壤水热通量的热力学模型。 国外对于冻土的研究起步较早,Philip DeVries和DeVries根据土壤水热耦合运动方程给出非等温条件下土壤水热运动的理论,Taylor Cary最先运用不可逆热力学定律建立了温度梯度下土壤水分运动的理论,Shao等利用变量代换和Fourier变换得到了特殊的初始和边界条件下水热耦合运移方程的解析解,求解步骤非常直接,因而可用于验证水热耦合模型的数值求解过程;国内土壤水热耦合数值模型的研究起步于20世纪80年代,隋红建、康绍忠、孙景生、郭庆荣分别对不同覆盖条件下的田间水热运移进行了数值模拟研究,实现不同覆盖层下非均质土壤水热分布的定量分析,并辅以田间试验验证;胡和平、黄兴法以我国广泛分布的冻土为背景,分析了冻结条件下非饱和土壤水、热、盐耦合的运动,建立了相应的数学模型,讨论其分布的一般规律,杨诗秀、雷志栋等模拟了水平、垂直土柱冻结过程;叶佰生、陈肖柏、胡和平等引入了Clapeyron方程来研究冻土中水热迁移问题。 实验所采用的土壤取自沈阳农业大学所属实验场,土壤为棕壤土,质地为壤土,其主要物理指标:容重1.4 g/cm3;有效粒径0.01 mm;均匀系数3.00;曲率系数0.96。 采用NMR法测定土壤冻结特征曲线和土壤融化特征曲线。采用水平土柱实验法测定土壤的导水系数(K)和扩散系数(D)。 式中,θ为液态水含水率(m3/m3)。 假定土壤介质不可压缩、均质各向同性,忽略水气迁移、热量对流作用的情况下,只考虑一维情况,根据达西定律和质量守恒原理,得非饱和土壤水分运动基本方程: 式中,θI为含冰率(m3/m3),ρI为冰密度(103kg/m3),ρw为水密度(103kg/m3),D(θ)为土壤水分扩散系数(m2/s),K(θ)为土壤导水系数(m/s)。 只考虑一维情况,根据Fourier定律和质量守恒原理得热流基本方程,得非饱和土壤热流运动基本方程: 式中,cv为土壤的热容量(J/m3/K),T为地温(K),Kh为土壤热传导系数(W/m/K),LI为冰的溶解潜热(J/m3)。 补充联系方程,即土壤中未冻水含水率θ和温度的关系: 4 有限差分法的推求 由于土壤冻结过程中水热耦合运动方程的非线性及初始条件、边界条件的复杂性,用解析法或半解析法推求是非常困难的,目前常用的方法是数值计算方法,笔者采用有限差分法。利用一元函数Taylor展开公式,采用Grank-Nicolson格式,在点列方程,并代入上式整理得: 式6、7即为求解差分方程,同时和式5联合组成方程组,在MATLAB 7.01中进行求解。 5 土壤冻结前后温度变化的规律 初始含水率分别为15%、25%、35%、47%,土柱下表面温度维持在1℃左右,总共冻结时间分别为9 862、9 862、6 661、8 408 min。 图1为土壤冻结后土壤水分模拟值的分布状况,随着初始含水率的增加,土壤的冻结锋面在逐渐增加,当含水率接近饱和含水率的时候,土壤锋面出现在8~10 cm区间之间。图2、3、4和图5为土壤冻结后,温度模拟值的分布状况。由图可知,冻结后的温度随土柱高度呈现线性变化规律,特别是1 600~2 000 min之后,线性变化规律明显。图2中,冻后高度为9.9、0~3 cm未冻,3.0~9.9 cm已冻,土体表面无明显冰脉,土体中无明显冰晶;图3中冻后高度为9.9、0~3 cm未冻,3.0~9.9 cm已冻,土体表面无明显冰脉,土体中无明显冰晶;图4中,0~1、1~2 cm未冻,3 cm以上已冻,3~9 cm有大量冰晶,2.5~9.4 cm有大量冰脉和裂纹;图5中,0~1 cm未冻,1~2 cm微冻,3.0~11.1 cm已冻,3.0~8.8 cm有大量裂纹、3~9 cm有大量冰晶,9~10 cm有稍少量的冰晶。 3 数学模型 2 实验设计 6 结论