正冻土冻结过程水-热耦合模型的建立与应用.docxVIP

正冻土冻结过程水-热耦合模型的建立与应用 土壤中的水分运动和含水量分布与热流和温度分布有关。一般来说，土壤含水量的分布和变化会导致不同位置土壤段各点的热特性参数（比热容积热导率等），并随着时间的推移而变化。相反，土壤中的热流和温度分布也会影响水流的运动。温度变化会导致水的物理性质发生变化，土壤中的水分运动参数会发生变化。 早在20世纪50年代末到60年代初, 原苏联的冻土学研究者就将一维的热传导方程和一维的渗流方程进行联列, 提出了研究土壤冻结过程的水热输运模型. 该模型假定冻土和未冻土中的土热水分扩散系数、 土壤内的初始温度和含水量均为常数, 而实际上, 这些参数都是变化的, 因此, 该模型计算结果难以反映实际情况. 70年代初, Harlan根据当时一些新的试验结果和观测事实提出了Harlan模型, 之后又产生了Taylor与Luthin模型, Harlan认为水分迁移的动力是土水势梯度, 并考虑了水分在迁移过程所携带的热量, 但并没给出确定土水势的方法. Taylor与Luthin模型是Harlan模型在特定条件下的简化. 本文通过实验, 观测水分场和温度场耦合变化过程, 并采用相应的数学模型进行模拟计算. 1 试验研究、观测和分析 1.1 试验设备及试样方法 在以往的试验中, 无法测试土样在冻结过程中其含水率及温度随时间的变化, 只能通过破损性试验测试终态(经过冻结过后)的土样含水率, 不能给出温度梯度和含水率随时间的变化规律, 因此不能很好的分析水分迁移过程中水分场和温度场的耦合过程. 本试验采用由长安大学研究开发的水分迁移测试系统, 在土样无破损的条件下, 对土柱中的点位进行温度和含水率的动态观测. 试验测试系统共由3大部分组成, 即含水率与温度的测试装置(温度传感器及水分传感器等)、 试件的温度控制装置(冷浴、 顶板及底板等)和试件的绝热装置(隔温壁). 土柱试件高为22 cm, 直径为16 cm(图1). 试样用土取自青藏公路雁石坪段. 首先将烘干碾碎的土配制成水分均匀的土样, 测定其初始含水率, 本次试验的初始含水率为15.7%(重量含水率), 然后将其分层装入试件模中, 使其达到密实状态, 不要有过多的气泡存在. 然后将温度传感器和水分传感器插入试件中, 试件制好后, 静置8 h, 待试件将多余的气体排出, 达到均匀状态后, 通过制冷装置使-10 ℃的冷冻液通过试件的顶板, 使2 ℃的冷冻液通过试件的底板, 对土柱施加温度梯度. 试验持续72 h, 每隔3 h记录一组温度传感器和水分传感器的读数. 1.2 土的深度和厚度随试件深度的变化特征 通过试验观测得到, 不同时刻温度随试件深度的变化曲线(图2)及不同时刻水分随试件深度的变化曲线(图3). 土柱冻结后从表层向下土体含水率将会提高. 当土壤发生单向冻结时, 土柱从上到下产生了较大的温度梯度, 从而破坏了土体中的水量平衡, 使其水分场发生了重新分布, 水分从土样的暖端向冷端迁移, 进而土上层的含水率较冻结前有所提高. 由图2不同时刻温度随试件深度的变化曲线可以看出, 试件在24 h之内温度变化较大, 随着时间的增加后期温度降低幅度不大;试件上部的温度梯度大于下部的温度梯度. 从图3可以看出, 含水量冻结前后的变化量随土层的深度加大而有减小趋势, 其变化与土的冻结过程有关. 土柱在初始的温度与室温比较接近, 当把土柱上下两个端面降温后, 尤其在土柱的顶面0~4 cm的范围内, 降温的速度很大, 导致土样迅速冻结, 水分迁移量不是很大;但土样在4 cm左右含水率达到最大, 是因为该处产生冻结锋面(从试件的表面看到该处产生较密集的裂纹), 致使试件下部水分向冻结锋面迁移; 在8 cm处含水量最小, 主要由于在该处的水分向上迁移, 而其下部的水分由于温度梯度小, 致使其不能及时的补给土柱的中层. 因此, 可以看出温度梯度是导致水量梯度产生的一个直接原因, 在温度梯度作用下, 土中未冻水沿着温度降低的方向迁移, 迁移量随温度梯度的增大而增加. 2 数学模型 在引用Flerchingeret al.建立的垂直一维冻土系统水、 热流耦合模型基础上, 通过限制一定的边界条件进行计算. 2.1 冻土水、热迁移 (1) 假定土壤介质不可压缩、 均质连续且各向同性. (2) 无溶质迁移. (3) 水分迁移符合达西定律. (4) 冻土中未冻水含量与土壤负温处于动态平衡中. (5) 水、 热的运动主要发生在垂直方向, 近似为一维问题. 2.2 冻土系统内的水通量及冰、化反应 冻土系统内的热通量方程为: ??z[ks?Τ?z]+ρlcl?(qlΤ)?z+S=Cs?Τ?t-ρiLf?θi?t+Lν[?ρν?t+?qν?z](1)??z[ks?T?z]+ρlcl?(qlT)?z+S=C