本文聚焦于干旱浅层地下水区域（AASGT）冻融循环中土壤水-热-盐耦合动态的模拟改进。研究团队通过引入毛细管束模型（Capillary Bundle Model, CBM）优化了传统SHAW-FH模型，开发出SHAW-FH-CBM新模型，其核心创新在于突破传统冰阻抗模型（Ice Impedance Model, IIM）的局限性，通过实验数据驱动的参数化方法解决冻土水迁移与盐分积累的耦合模拟难题。以下从研究背景、方法改进、实验验证和实际应用四个维度进行解读：
### 一、研究背景与问题提出
全球变暖加速了极地和高海拔地区冻融循环的时空异质性，导致浅层地下水与土壤水交换机制紊乱。在AASGT区域，冻融期土壤体积膨胀（冻胀）引发水分迁移路径改变，盐分在冻土中的迁移呈现非线性特征，但现有模型未能有效耦合这些过程。以SHAW模型为代表的传统冻融模型存在两大缺陷：
1. **早期冻融期表层土壤水分估算偏差**：未考虑冻胀引起的土壤孔隙重构对渗透率的影响，导致冰层形成初期水分迁移被低估。
2. **稳定冻结期水盐迁移停滞假设**：传统模型采用固定临界孔隙率阈值（0.13）控制冻土渗透率，忽略盐分-冻胀协同作用下的动态孔隙结构变化。
### 二、模型改进与技术突破
#### （一）土壤变形耦合机制
通过引入土壤体积重构参数（式9-11），首次将冻胀导致的孔隙率变化纳入模型：
- **孔隙率动态计算**：根据冻胀程度（冰体积占比超过临界值0.99时）实时调整孔隙率，孔隙率增量与初始干密度成反比（式9），体现冻胀引起的土壤结构重组。
- **热传导参数修正**：基于重构孔隙率调整热扩散系数，当冻胀导致有效孔隙率提升20%时，导热系数下降15%-25%，更符合实际冻土热惯性特征。
#### （二）毛细管束模型替代冰阻抗模型
CBM模型通过三个创新解决IIM的不足：
1. **非零低渗特性**：在临界温度以下仍保留约1.96×10⁻⁶ cm/s的渗透率（式19），避免完全停滞的物理假设。
2. **盐分-温度耦合效应**：通过拟合参数β（-0.417至-0.716）量化盐分对冻融曲线斜率的影响（式23），盐分浓度每增加1%，临界温度下降0.3-0.5℃。
3. **动态孔隙网络重构**：采用有效孔隙率θ*替代固定阈值，θ* = θ_s + Δθ（Δθ为冻胀导致的孔隙增量），使渗透率计算更符合实际冻融过程。
#### （三）多参数耦合优化策略
针对AASGT区域特殊的水文地质条件，开发了三阶段参数优化机制：
1. **初始水盐状态匹配**：通过核磁共振（NMR）建立冻融曲线（SFCC）与土壤质地、盐分浓度的映射关系（R²>0.83）。
2. **冻融期动态校准**：采用分层次逆向参数优化法，先冻结层后冻结层调整参数，确保模型能捕捉分层冻融特征。
3. **季节循环记忆模块**：引入冻融循环周期记忆参数（CP），使模型具备跨年度参数迁移能力。
### 三、实验验证与性能对比
#### （一）实验室验证体系
1. **多盐分梯度实验**：设置0.1%-1.0% NaCl溶液梯度（5组），通过等温法获取15种不同初始水盐组合的SFCC曲线。
2. **温度-时间响应测试**：采用-8℃至0℃的变温实验，每0.2℃间隔记录NMR信号，构建冻融动力学数据库。
3. **临界孔隙率验证**：当θ*达到0.13时，模型预测渗透率突降至0（IIM）或维持1.0×10⁻⁶ cm/s（CBM），与显微镜观测的孔隙连通性变化吻合度达92%。
#### （二）田间模型验证结果
基于2019-2021年贺兰灌区观测数据（图3-9），SHAW-FH-CBM模型在三个维度实现突破：
1. **冻融期水分动态**：
- 表层（0-5cm）冻结期水分增量达60%饱和含水率（θ_s=0.472 m³/m³）
- 深层（>20cm）稳定冻结期水分波动幅度提高3倍（RMSE从0.04降至0.01 m³/m³）
2. **盐分迁移模拟精度**：
- 表层盐分浓度预测误差<8%（RMSE=0.09 g/kg）
- 中层盐分峰值捕捉误差<15%（NSE=0.43）
3. **冻胀-热传导耦合**：
- 冻胀导致导热系数下降幅度与观测值匹配度达89%
- 冻融循环周期预测误差<12%（NSE=0.72）
### 四、模型应用与工程指导
#### （一）农业水资源管理
1. **灌溉制度优化**：通过模拟显示，当表层冻胀导致孔隙率提升15%时，灌溉水利用效率提高22%
2. **排水系统设计**：预测深层冻土（>60cm）盐分累积量比表层高40%，指导排水系统深度优化
3. **盐碱地改良**：模拟表明，冻融期盐分表聚效应使表层0-20cm区域盐分浓度年上升速率达0.8 g/kg
#### （二）生态工程评估
1. **冻融生态位重建**：模型显示，冻胀引起的土壤位移可使植被根系分布调整幅度达30-50cm
2. **冻融通道形成**：预测冻融期形成的水通道密度达4.2条/m²，影响微生物群落分布
3. **冻融灾害预警**：结合土壤水分-温度耦合模型，可提前72小时预警冻胀裂缝（准确率89%）
#### （三）气候变化应对
1. **冻融期延长模拟**：预测升温1℃可使冻融期缩短15-20天
2. **碳封存潜力评估**：冻融循环增强的有机质分解速率提升18-25%
3. **冻土碳汇计算**：修正模型使碳封存量预测误差<10%
### 五、理论创新与工程实践启示
1. **冻融动力学理论**：提出"盐分-冻胀-渗透率"三元耦合机制，解释了AASGT区域"越冻越盐"现象（式23）
2. **孔隙重构理论**：建立冻胀导致的孔隙率动态变化方程（式9），为冻土工程寿命预测提供新参数
3. **多尺度验证体系**：构建实验室（0.5m³）-田间（573,300ha）-区域（中国西北）三级验证框架
### 六、未来研究方向
1. **多相流耦合**：扩展至气-液-固三相流耦合模拟
2. **冻融疲劳累积**：建立冻融循环次数与渗透率衰减的量化关系
3. **机器学习融合**：将深度学习算法用于冻融参数实时反演
该研究为寒区农业可持续发展提供了新的技术工具，其模型框架已应用于中国西部5个农业大县的盐碱地治理工程，使灌溉水利用率提升18.7%，土壤有机质含量年增加0.32%。模型突破传统冻土模拟中"水-热-盐"分立处理的局限，建立了多物理场耦合的动态仿真体系，为全球变暖背景下寒区农业风险管理提供了科学基础。
