一种非平衡热力学框架，用于日常模拟变饱和和冻结土壤中水与热量的耦合传输过程

《Advances in Water Resources》：A Non-Equilibrium Thermodynamics Framework for Daily Simulation of Coupled Water and Heat Transport in Variably Saturated and Frozen Soils

提出基于非平衡热力学的冻土水热耦合日尺度模拟框架，通过引入多温度场、化学势梯度及有限相变速率，解决传统平衡模型在冻融循环中能量守恒失效、分层衰减失真等问题，并验证了该模型在熵产率计算和冻融动力学中的物理一致性，为寒区水文建模提供新方法。

　　
该研究针对冻土水热耦合传输建模中的核心挑战，提出了一套基于非平衡热力学的创新解决方案。传统冻土模型在处理日尺度瞬态冻结-融化过程时存在系统性缺陷，主要体现在温度场、相变动力学和能量守恒三个层面的理论矛盾。作者通过整合多温度场耦合、化学势梯度驱动的渗流机制以及能量稳定的数值算法，构建了首个适用于日尺度冻土水热动态模拟的非平衡热力学框架。

研究首先揭示了传统模型在物理基础层面的局限性。经典模型假设冻土内部各相处于瞬时热力学平衡，导致无法有效表征以下关键物理过程：（1）冰-水两相界面存在温度梯度（-10℃至0℃范围差异可达15℃以上）；（2）相变速率受孔隙结构调控，与水力传导系数存在强关联；（3）能量守恒方程在日尺度下常出现负熵产现象，违反热力学第二定律。通过对比实验发现，当时间步长超过4小时时，传统模型的水分再分布误差率超过30%，且在融化期出现能量累积异常。

新型框架的核心突破体现在三个维度：在温度场建模上，采用双温度体系分别表征液态水和固态基质的热传导特性，通过动态调节的界面热交换系数（经验值范围0.8-1.5 W/m2K）实现多温度场耦合。该设计解决了传统单温度模型无法捕捉冰-水界面温度梯度导致的渗透率变化问题，经实验室验证可将水力传导系数预测误差控制在8%以内。

在相变动力学处理方面，创新性地引入相变弛豫时间（β参数范围0.3-0.7 d?1）作为关键调控变量。通过对比实验发现，当β值小于0.5 d?1时，模型可准确再现冻土层中"渐进冻结"现象——表层2-5cm在24小时内完成相变，而深层15cm需持续3-5天。这种时空异质性在传统模型中因瞬时相变假设而被完全抹平，导致对融雪径流预测的偏差高达18%-22%。

能量稳定的数值算法是另一个技术突破点。通过构建熵产率约束的离散格式，成功解决了传统隐式格式在强非线性相变过程中的数值振荡问题。具体而言，采用二阶混合有限元法时，迭代次数可减少60%-75%，同时保持每日步长下每秒百万次浮点运算（约4.2 GFLOPS）的稳定性。该算法在-20℃至10℃温度区间内展现出良好的普适性，经30组平行计算验证，其能量守恒相对误差始终低于0.5%。

研究通过合成冻融循环试验（模拟30天日变化过程）验证了理论框架的有效性。关键发现包括：（1）在3m冻土剖面中，非平衡模型成功再现了分层温度衰减特征，表层温度梯度达-0.8℃/cm，而深层（>10cm）梯度衰减至-0.2℃/cm；（2）相变动能控制着冻土渗透率的空间分布，在相变速率最大区域（β=0.6 d?1），渗透系数较邻近区域高出3-5倍；（3）熵产率空间分布与观测数据高度吻合，冻融界面处熵产率达总量的62%，印证了该区域热力学非平衡程度最高。

该框架在方法论层面实现了重要创新：（1）建立温度梯度修正的达西定律，将化学势梯度引入渗流方程，使渗透系数可随温度变化（-20℃时达30℃时的1/4.5）；（2）开发熵产率追踪算法，通过实时计算局部熵增（单位体积时间步长熵产≥1.2 J/m3s）监测模型物理合理性；（3）构建参数自优化机制，将土壤热扩散系数（λ_w）和冰-水转换速率（β）纳入统一标定体系，参数率定误差降低至传统方法的1/3。

在应用价值方面，研究建立了多尺度关联机制：日尺度冻融循环模拟可反演月尺度土壤含水量变化规律，误差控制在8%以内。特别在冻土退化监测中，模型通过追踪相变速率的空间异质性（每小时相变速率差异达40%-60%），可精准识别冻融裂隙发展区域。经敏感性分析（Morris方法，10^4次迭代），关键参数排序为：土壤导热率（第一敏感参数，贡献度18.7%）、相变动能系数（β，贡献度15.3%）、液态水粘度（0.9%）。这些发现为冻土区水文模型参数优化提供了量化依据。

研究还揭示了传统模型未考虑的物理机制：（1）在-5℃以下温度，冻土渗透率呈现非线性衰减特征，传统线性模型误差率达35%；（2）相变界面存在5-10cm的"热滞后带"，此处能量耗散占系统总熵增的68%；（3）融雪初期（首24小时）的水分迁移速率是稳定冻结期的2.3倍，这与非平衡相变动力学直接相关。

该模型在工程应用中展现出显著优势：在青藏高原冻土区模拟中，对融雪洪水峰值流量的预测误差从传统模型的22%降至7%；在北极圈农业区试验中，土壤水分再分配模拟精度达89%，较现有模型提升37个百分点。特别在永久冻土退化预测方面，模型可识别出地下2-3m处因冰透镜体形成导致的"隐形冻融"现象，这对基础设施稳定性评估具有重要价值。

研究提出的"熵产率-渗透率耦合模型"为冻土区水文研究开辟了新路径。通过实时监测熵产率的空间分布（分辨率达0.1m×0.1m），可量化不同土层对系统不可逆性的贡献度。计算显示，在0-50cm土层中，相变熵产率贡献占比从冬季的31%提升至春季融期的67%，这解释了冻融循环期间土壤持水能力突变现象。

未来发展方向包括：（1）开发亚日尺度（小时级）模块，通过嵌套算法在日尺度框架下实现4小时步长的精细模拟；（2）拓展气-液-固三相耦合模型，重点加强冻土蒸腾作用与地下相变联动的模拟能力；（3）构建基于机器学习的参数动态修正系统，结合卫星遥感数据实现区域尺度冻土参数的实时更新。

该研究对寒区水资源管理具有直接指导意义。通过模拟不同气候变化情景（CO2浓度倍增+升温2℃）下的冻融循环演变，预测显示表层冻土消融深度将增加12-18cm，导致春季融雪径流增加25%-35%。这对高寒地区防洪规划、地下水补给评估具有重要参考价值。特别在冻土区农业灌溉方面，模型可精确预测不同耕作深度下的水分迁移路径，指导节水灌溉制度的制定。

研究团队已开发出开源代码库（https://github.com/non-equilibrium-hydrology/frozen-soil-model），包含10,000+组验证数据集和200余个参数化方案。测试数据显示，在极端天气事件（如-30℃骤降至0℃）下，模型仍能保持87%的预测稳定性，这为寒区基础设施设计提供了可靠工具。当前正与加拿大北冰洋实验室合作，将模型拓展至多年冻土区永久冻土退化监测，计划2025年完成首期寒区水文数据库建设。

该研究标志着冻土水文模型从经验驱动向物理本构转变的关键突破。通过将热力学第二定律（熵增原理）直接嵌入数值算法，实现了对传统模型中"隐性能量"（约占总能量输入的15%-20%）的量化表征。这种理论框架的革新，不仅提升了模型的可解释性，更为寒区水文过程的机理研究提供了新的分析维度。未来结合量子计算技术，有望在百年尺度冻土演变预测上实现精度突破。