冻结条件下冻土物理力学特性的数值模拟.docxVIP

冻结条件下冻土物理力学特性的数值模拟 中国是世界上最大的冷冻土壤分布区域。多年冻土和季节冻土的面积分别占全国面积的21.5%和53.5%。53.5%.中国位于中国东部30n的以北，川陕地区33n的以北。土壤的冻结是一个十分复杂的过程, 在冻结过程中, 土壤中水分状况、热状况和溶质状况的变化并不是独立的, 而是相互制约、相互影响的, 是一个耦合变化过程。 综合性的描述和研究土壤中水、热、溶质耦合运移的文献较为少见。20世纪90年代有些学者从不同侧面对土壤中水、热、盐运移问题进行了理论及试验研究, Nassar利用土壤中水、热、盐运动方程和连续方程, 提出了描述土壤中水、热 盐运动的耦合方程。 本文通过对棕壤土在冻结条件下水、热、溶质耦合运移的模拟, 得出水、热、溶质迁移规律, 为后续研究含多种溶质的水、热、溶质耦合运移以及田间水、热、溶质耦合运移奠定基础。 1 冻结时间和温度 室内试验在中科院寒区旱区环境与工程研究所进行, 所采用的土壤是棕壤土, 其主要物理指标见表1, 加入5%氯化钾溶液, 土样高度和直径均为10 cm, 从上端降温进行冻结, 降温速率为-0.1 ℃/h, 冻结时间为120 h, 底板控制温度在1℃, 土样上、下两端密封。试验结束后分层取样测定含水率的分布状况。 实验土柱初始含水率为40%, 初始溶液浓度为5%, 初始温度分布如表2, 土壤冻结特征曲线如图1和图2。 2 数学模型 2.1 一维和一维热流运移基本方程 假定土壤介质不可压缩、均质且各向同性;冻结过程中, 水分以液态形式迁移;忽略水汽迁移和热量对流作用的影响及溶质势梯度对水分迁移的影响;冻土中未冻水含量与土壤负温处于动态平衡中;水、热、溶质的运动主要发生在垂向, 可近似为一维问题。 一维情况下, 根据Darcy定律和质量守恒原理, 得非饱和土壤水分运动基本方程: ?θ?t=??z[D(θ)?θ?z]-?Κ(θ)?z-ρΙρw?θΙ?t(1) 式中θ为液态水含水率;θI为含冰率;ρI为冰密度;ρw为水密度;D为土壤水分扩散系数;K为土壤导水系数。 一维情况下, 根据Fourier定律和质量守恒原理, 得热流基本方程: cv?Τ?t=??z[Κh?Τ?z]+LΙρΙ?θΙ?t(2) 式中cv为土壤的热容量;T为温度;Kh为土壤热传导系数;LI为冰的溶解潜热。 一维情况下, 根据Fick定律和质量守恒原理, 得溶质运移基本方程: ?(θc)?t=??z[Dsh(ν,θ)?c?z]-?(qc)?z(3) 式中c为溶质浓度;Dsh为水动力弥散系数。 2.2 土体含水率的确定 冻土中, 未冻水含量与土体负温保持动态平衡, 即土体中未冻水含水率和温度的关系: θ≤θmax(Τ)(4) 式中θmax为相应土体负温下可能的最大未冻水含量 (土壤冻结特征曲线) 。 3 热流差分方程 冻结土柱的水、热、溶质运动状况是由土柱两端的土壤水、热、溶质条件决定的。由实验可知, 热流差分方程上、下边界条件均为已知, 属于第1类边界条件。水流差分方程和溶质差分方程上下边界为零通量面, 属于第2类边界条件。 3.1 土柱冷端差分方程 水流差分方程计算时, 采用半控制容积法处理, 即原差分网格不变, 对边界处的半个单元, 直接根据水量平衡列出差分方程。如图3, 对土柱冷端i为0～0.5处列平衡方程: (R-q1+1/2)△t=(ρΙρw△θΙ+△θu)△z12(5) 式中表面通量R=0。 q1+1/2=-D(θu)?θu?z+Κ(θu)(6) 对土柱暖端采取同样的方法处理。 3.2 b土壤蒸发通量 溶质差分方程边界采用与水流差分方程边界相似的处理, 如图4, 在i为0～0.5处列平衡方程: (Es-J1+1/2)△t=△c△z12(7) 式中表面蒸发通量Es=0。 J1+1/2=-Dsh?c?z+qc(8) 下边界采取同样处理方法。 4 土壤特征参数 4.1 抗冻土法 未冻土的导水率和扩散率可由实验测得, 而冻土的导水率和扩散率测定非常困难, 目前国内尚无测定先例。确定冻土导水率和扩散率的方法有很多种, 本文采用阻抗系数法, 即用冻土区的导水率和扩散率除以阻抗系数I来表示,I的大小取决于土壤中含冰量的多少。Taycor和Luthin提出I由下式确定: Ι=1010θΙ(9) 未冻区的导水率和扩散率采用实验值: Κ(θ)=1.5089×10-10e38.666θ(10)D(θ)=0.0017e17.596θ(11) 冻结区: Κ(θ)=1.5089×10-10e38.666θ/Ι(12)D(θ)=0.0017e17.596θ/Ι(13) 4.2 土壤中各组成物的热容量的计算 估算土壤热导率是十分困难的, 它不但取决于土壤中各组成物的比例, 还和各组成物的形状等因素有关。De Vries提